duffing01.m

来自「非线性动力学中lyapunov、duffing等方程的求解」· M 代码 · 共 7 行

function dx=ode45(t,x)
global m c k w f a


dx=[x(2);
    -c/m*x(2)+a/m*x(1)-k/m*x(1)^3+f/m*cos(w*t)];