function [xdot] =etat(t, x, A, u)
//  Calcul de la solution du systeme d'etat   
 
 xdot = A*x+B*u(t)

endfunction

function ut=u(t)
// Calcul de la commande
    
//    Réponse impulsionnelle
//    ut=0;
//    if (t == 0) then ut = 1;  end;

//    Réponse indicielle
      ut = 1;
endfunction
