function f=yy1(t,x)