function f=ode1(x,y)
f=x.*y;
