x0=0;
xm=0.5;
f=@(x,y) -y.*sin(x)-y.^2;
y0=1;
