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

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


end