function yp = f1(t,y)
yp = -y;