syms y(t) r 2 V odeToVectorField(diff(y, 2) + 0.5*diff(y, 1) - y^3+y^5 r*cos(t)) M matlabFunction(V, vars , { t , Y }) -syms y(t) r 2 V odeToVectorField(diff(y, 2)+ 0.5*diff(y, 1)- y^3+y^5 r*cos(t)) M matlabFunction(V, vars