function y=mat_card(u);

% En entrée, les angles de Cardan u = angles (roulis tangage lacet) en rd

angle_roulis=u(1);
angle_tangage=u(2);
angle_lacet=u(3);

% En sortie la matrice de transformation ( pas de passage !!! ) permettant le calcul de la rotation
% du satellite par rapport au repere orbital XYZ
% Cette matrice servira aussi à la transformation des couples 

 
P(1,1)=1;
P(1,2)=0;
P(1,3)=-sin(angle_tangage);

P(2,1)=0;
P(2,2)=cos(angle_roulis);;
P(2,3)=cos(angle_tangage)*sin(angle_roulis);


P(3,1)=0;
P(3,2)=-sin(angle_roulis);;
P(3,3)=cos(angle_tangage)*cos(angle_roulis);

y=P;

end