function yp=myfun(t,y)
yp(1)=y(2);
yp(2)=sin(y(1)*y(2));
yp=yp(:);