特殊微分方程的数值解——刚性微分方程（ode15s）
我们在微分方程的标准型转换中已经提到的Van der Pol方程就是一种刚性微分方程。刚性微分方程可能有几个状态变量，其中某个状态变量变化的特别快，远远的快于其他的状态变量，这样就会使得常规算法中，选择步长出现问题。再谈Van der Pol方程：计算mu=1000时的Van der Pol方程的数值解。h_opt=odeset;h_opt.RelTol=1e-6;x0=[2;0];t_final=3000;f=@(t,x,mu)[x(2);-mu*(x(1)^2-1)*x(2)-x(1)
