[求助]请教奇异摄动法解微分方程用Matlab如何实现

[gujian1211]

本人因毕业设计需要用奇异摄动法解微分方程，请教高手如何用Matlab实现

例如：Van der pol 方程 ：y''+(y*2-1)y'+y=0 应该如何实现呢，高手请指教

奇异摄动法的基本思想是:对于含有未知u(x,t)的方程，欲求的u(x,t)，则可以将函数u(x,t)按小参数t展开
u(x,t)=u0(x)+u1(x) t+u2(x) t*2+u3(x) t*3+……
其中：ui(x)是与t无关的未知函数，u0(x)是当t=t0时方程的初值，可以认为已知的。