function dxy=esesepm(t,x)
gam=0.1;
dxy=[-x(2) ; +gam*x(2)-1+x(1)^2];
