function y=energie(u);

global Asat vitesse_rotation

% Variable d'entrée u, le vecteur rotation inertiel
%	u(1)=p, u(2)=q, u(3)=r

% Fonction donnant y l'énergie de rotation autour du centre d'inertie
% mersurée dans le repère inertiel


energie=0.5*u'*Asat*u;  % Energie de rotation

y=energie;




end