vdp_3.m

来自「提供积分算法变步长Gill法代码」· M 代码 · 共 7 行

function dx = vdp_3(t,x,para)
% 为题3 Val-der-Pol的微分方程
A  = para.A;
B  = para.B;
dx = [x(2);
    A*(1-x(1)^2)*x(2)-B*x(1)];