che01.m

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

function dx=ode45(t,x)
global ee n alpha dx0 h gamma


dx=[x(2);
    -x(1)+ee*(-2*n*x(2)-alpha*sign(x(2)-dx0)+h*cos(gamma*t))];