function z = f(x,y); % % I want to compute the ODE dy/dx = f(x,y) = cos(x)*y z = cos(x)*y;