function y=B_SIM(Asat,K_SAT,H_SAT,L_SAT,Q_SAT);

 % Calcul de la matrice B_SIM de l''équation dX/dt=A_SIM X+B_SIM U

 M(1:3,1:3)=inv(Asat);
 
 M(4:6,1:3)=zeros(3,3);

 y=M;  

 End;