function [M] = rotz(a)
%Rotaçao em z ; a em radianos!
%   Detailed explanation goes here
M=[cos(a) -sin(a)   0      0
   sin(a)  cos(a)   0      0
    0       0       1      0
    0       0       0      1];
end
