vdp1000.m

MATLAB科学计算与工程分析源代码源程序6

function dy = vdp1000(t,y)
dy = zeros(2,1);    
%dy(1) = -10*y(1)+9*y(2);
%dy(2) = 10*y(1)-11*y(2);
dy(1) = y(2);
dy(2) = 1000*(1 - y(1)^2)*y(2) - y(1);