function yp=f1(x,y)
yp=-2*y;
%yp=(2-y)*y;