
R=1;
L=0.001;
J=0.03;
b=0.001;
ke=0.5;
t0=1;
V0=100;
figure(1)
clf
for Ce=[0:20:100]
   
   Tfin=5;
   sim('motccmdl',Tfin)
   
   subplot(211)
   plot(t,Ia)
   hold on
   ylabel('Corrente Ia [A]')
   
   subplot(212)
   plot(t,wm)
   hold on
   ylabel('Velocita Wm [rad/s]')
   xlabel('Tempo t  [s]')
   
end

