function f=fsin(t,y)
   f=sin(t*y);