05-24.txt

来自「《MATLAB R2007基础教程》 清华大学出版社,包含电子课件,实例原文件,」· 文本 代码 · 共 10 行

function dydt = vdp1000(t,y)
dydt = [y(2); 1000*(1-y(1)^2)*y(2)-y(1)];


>> [t,y] = ode15s(@vdp1000,[0 3000],[2; 0]);
>> plot(t,y(:,1),'-');
>> title('Solution of van der Pol Equation, \mu = 1000');
>> xlabel('time t');
>> ylabel('solution y_1');