ajout derniere version

BE_Beam.m

@@ -203,7 +203,6 @@ for Ts = Ts2:0.001:Ts3
     if all(abs(eig(feedback(c2d(sys,Ts,'zoh'),k))) < 1)
         ;
     else
-        Ts;
         break;
     end
     
@@ -355,6 +354,8 @@ if SAVE
 end
 
 %%
+
+for Ts=
 sys_integral_d = c2d(sys_integral, 0.5, 'zoh');
 t_d = 0:0.5:10;
 cmd_d = double (t_d >= 1) ;