[讨论]怎样在matlab中用四阶Runge-kutta法求解四阶常微分方程

[smiles860710]

以二阶微分方程为例
y1''-(1-y1^2)*y1'+y1=0
降阶令y2=y1'
y2'=(1-y1^2)*y2-y1
M文件
function dydx=vdp1(x,y)

dydx=[y(2);(1-y(1)^2)*y(2)-y(1)];

end在命令窗口输入[x,y]=ode45(@vdp1,[0 20],[20,5])就可得到结果