% 
% Funçao rotação
% T=Rot(a,b)
% em que T será a matriz transformação rotação em 3D,
% a será caracter correspondente ao eixo de rotação ( x,y ou z )
% b será o angulo em graus.

function T=Rot(c,a);
b=0;
a=deg2rad(a);
while(b==0);
    if (c=='x');
    T=[1  0        0    0
       0 cos(a) -sin(a) 0
       0 sin(a) cos(a)  0
       0  0        0    1];
    b=1;
    elseif (c=='y');
     T=[cos(a) 0 sin(a) 0
         0     1   0    0
       -sin(a) 0 cos(a) 0
         0     0   0    1];
     b=1;
    elseif (c=='z');
    T=[cos(a) -sin(a) 0 0
       sin(a)  cos(a) 0 0
        0       0     1 0
        0       0     0 1];
    b=1;
    else
    disp('Selecione eixo correcto, (x,y ou z)');
    c=input('Novo valor do eixo ');
    end
end
end