%Variatia momentului pt. k=0.5 si G=10
[x,yn]=ode23('g1a',0,.0101-.0001,1,.00001,1)
grid;plot(x,yn,'-r');title('Variatia lui M in functie de s');xlabel('s');ylabel('M');

hold on
[x,yn]=ode23('g1a',0,-0.1,1,.000001,1)
grid;plot(x,yn,'*b');
