function yp=f(t,y)
yp=-y; %yp=y'