optimazion of R=((θ_1*sin(20*θ_2)+ θ_2*sin(20*θ_1))^2)*cosh(θ_1*sin(10*θ_1))+(( θ_1*cos(10*θ_2)- θ_2*sin(10*θ_1))^2)*cosh(θ_2*cos(20*θ_2)))) with powel

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