function [M]=rotz(a)
%Matriz transformação 3D
%   Matriz de transformação homogénea no espaço 

M=[cos(a)    -sin(a)       0       0
   sin(a)     cos(a)       0       0
   0          0            1       0
   0          0            0       1];


end