龙格-库塔（Runge-Kutta）法的matlab程序

龙格-库塔(Runge-Kutta)方法是一种在工程上应用广泛的高精度单步算法。本程序为4阶龙格-库塔法的matlab文件，用于求解微分方程。

fun=‘y-2*x./y‘
a=0;b=1;
y0=1;
h=0.2;
%四阶R-K方法
n=（b-a）/h;X=a:h:b;Y=zeros（1n+1）;
X（1）=a;Y（1）=y0;
for i=1:n
    x=X（i）;y=Y（i）;
    k1=eval（fun）;
    x=X（i）+h/2;y=Y（i）+h*k1/2;
    k2=eval（fun）;
    y=Y（i）+h*k2/2;
    k3=eval（fun）;
    x=X（i）+h;y=Y（i）+h*k3;
    k4=eval（fun）;
    Y（i+1）=Y（i）+h*（k1+2*k2+2*k3+k4）/6;
end
[X‘Y‘]

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件        346  2013-02-23 16:49  rk.m

     文件      31744  2014-03-22 18:58  4Runge-Kutta.doc

----------- ---------  ---------- -----  ----

                32090                    2