function f = myfunc(x,y) % This is the function f(x,y) in the ODE % dy/dt = y*(2-y) f = y*(2-y);