在matlab中用四阶龙格-库塔法解二阶微分方程怎么做？？最好有代码！！

Question

如题所述

Answer 1

例 y’= - y+x+1，y(0) = 1
首先建立M-文件 （weif.m）
function f = weif(x,y)
f=-y+x+1;
求解：[x，y]=ode45('weif',[0,1],1)
再如：
建立文件：
function dy = rigid(t,y)
dy = zeros(3,1); % a column vector
dy(1) = y(2) * y(3);
dy(2) = -y(1) * y(3);
dy(3) = -0.51 * y(1) * y(2);

命令窗口：
options = odeset('RelTol',1e-4,'AbsTol',[1e-4 1e-4 1e-5]);
[T,Y] = ode45(@rigid,[0 12],[0 1 1],options);追问

这个是一阶的，我要的二阶微分方程，而且是要自己编的程序 ，不是调用MATLAB自带函数。