vdp_3.m

提供积分算法变步长Gill法代码

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)];