四阶Runge-Kutta方法求解初值问题.doc

《四阶Runge-Kutta方法求解初值问题.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、四阶Runge-Kutta方法求解初值问题方法四阶Runge-Kutta计算过程PROGRANMFOUR_RUNGE_KUTTAIMPLICITDOUBLEPRECISION(A-H，O-Z)PARAMETER(N=100，NT=5)DIMENSIONX（N），Y（NT，N），RK（4）EXTERNAL::RGFUNH=0.1DOI=1，NY(1，I)=1.0X（I）=0.01*iENDDODOIT=1，NT-1DOI=1，NRK(1)=RGFUN(X(I),Y(IT,I))RK(2)=RGFUN(X(I)+H/2.0，Y(I

2、T,I)+H/2.0*RK(1))RK(3)=RGFUN(X(I)+H/2.0，Y(IT,I)+H/2.0*RK(2))RK(4)=RGFUN(X(I),H,Y(IT,I)+H*RK(3))Y(IT+1，I)=Y(IT,I),H/6.0*(RK(1)+2.0*RK(2)&+2.0*RK(3)+RK(4))ENDDOENDDOOPEN(11，FILE=’Y.DAT’)DOIT=1，NTDOI=1，NWRIT(11，*)X(I),Y(IT,I)ENDDOENDDOCLOSE(11)ENDFUNCTIONRGFUN(X,Y)IMPL

3、ICITDOUBLEPRECISION(A-H,O-Z)RGFUN=Y-2.0*X/YRETURNEND