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 , { t , Y })

sol　 ode45(M,[0 20],[2 0])

x　 linspace(0,20,100)

y1　 d (sol,x,1)

y2　 d (sol,x,2)

plot(x,y1, r- ,x,y2, b: )

legend( y1 , y2 )

xlabel( 时间 )

ylabel( y )

grid on