请问用四阶龙格库塔法解二阶微分方程的思想是什么？最近我遇到了一个难题。就是求一个二阶微分方程，形式如：y''=f(y)，初始条件是y(0)=0，y'(0)=0，y是t的函数.要求用四阶龙格库塔法求解。请问思路是什么？由于f(y)函数复杂，所以解析解是求不出来的，必须用四阶龙格库塔法解之，如果您有实例，希望能够和我分享一下。
请问用四阶龙格库塔法解二阶微分方程的思想是什么？ 例 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 vectordy(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按二阶龙格库塔法求解微分方程组，急用于毕业设计！ 先写状态方程function output=fangcheng(t，state)output=…再用龙格库塔求解state=[0 0 0]；t=0；for t=0：step：1k1=fangcheng(t，state)。
MATLAB中已知系统微分方程及初始值用欧拉法和龙格库塔法解一阶微分方程 function Euler欧拉法和龙格库塔算法解一阶常微分方程源代码例子dy/dx=-y+x+1f=inline('-y+x+1'，'x'，'y')；微分方程的右边项dx=0.5；x方向步长xleft=0；区域的左边界xright=10；区域的右边界xx=xleft：dx：xright；一系列离散的点n=length(xx)；点的个数y0=1；(1)欧拉法Euler=y0；for i=2：nEuler(i)=Euler(i-1)+dx*f(xx(i-1)，Euler(i-1))；end(2)龙格库塔法RK=y0；for i=2：nk1=f(xx(i-1)，RK(i-1))；k2=f(xx(i-1)+dx/2，RK(i-1)+k1*dx/2)；k3=f(xx(i-1)+dx/2，RK(i-1)+k2*dx/2)；k4=f(xx(i-1)+dx，RK(i-1)+k3*dx)；RK(i)=RK(i-1)+dx*(k1+2*k2+2*k3+k4)/6；endEuler和Rk法结果比较plot(xx，Euler，xx，RK)hold on精确解用作图syms xrightsolve=dsolve('Dy=-y+x+1'，'y(0)=1'，'x')；求出解析解rightdata=subs(rightsolve，xx)；将xx代入解析解，得到解析解对应的数值plot(xx，rightdata，'r*')legend('Euler'，'Runge-Kutta'，'analytic')
用MATLAB按二阶龙格库塔法求解微分方程组，大神速来，急急急 用MATLAB按二阶龙格库塔法求解微分方程组，在0时刻，x（0），y(0)，Z(0)，Vx（0），Vy（0），Vz（0）是各个方向的速度且已知，μ为常数，r等于x，y，z的平方和开根，求xyz。
