資源簡介
龍格-庫塔(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
- 上一篇:最大lyapunov指數matlab程序
- 下一篇:A*算法的matlab程序
評論
共有 條評論