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

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


end