函数Dydt = vdp1000(t,y)求mu = 1000时的van der Pol ode。%%参见ODE15S, ODE23S, ODE23T, ODE23TB。Jacek Kierzenka和Lawrence F. Shampine。版权所有The MathWorks, Inc.Dydt = [y(2);1000 * (1 y (1) ^ 2) * y (2) - y (1)];
用数值默认的相对和绝对容错(1 e - 3而且1 e-6,分别)非常慢,需要几分钟来求解和绘制解。数值需要数百万的时间步骤来完成集成,由于它的刚度领域,它努力满足公差。