x0=[0;0];
h=0.01;
ti=0;
tf=300;
%[t,x]=feuler(@masares,x0,ti,tf,h);
[t,x]=rk23(@masares,x0,ti,tf,1e-3,1e-9);
plot(t,x)