fortran下龙格-库塔法求解常微分方程(实例)

1.二阶龙格—库塔法实例：，程序代码：DIMENSION x(0:20),y(0:20),z(0:20),p(0:21)DOUBLE PRECISION x,y,z,p,h,x0,y0,z0,p0,n

1.二阶龙格—库塔法 实例： ， 程序代码： DIMENSION x(0:20),y(0:20),z(0:20),p(0:21) DOUBLEPRECISION x,y,z,p,h,x0,y0,z0,p0,n,k1,k2 f(x,y)=x-y h=1.0/20 n=20 x(0)=0.0 y(0)=0.0 DO i=0,n-1 x(i+1)=x(i)+h k1=f(x(i),y(i)) k2=f(x(i+1),y(i)+h*k1) y(i+1)=y(i)+h*(k1+k2)/2 ENDDO p(0)=0 DO i=0,n z(i)=p(i)+(-p(i))-1 exp p(i+1)=p(i)+h ENDDO open (10,file=) '1.txt' WRITE (10,10)(x(i),y(i),z(i),i=0,20) WRITE (*,10)(x(i),y(i),z(i),i=0,20) FORMAT (1x,f10.7,2x,f10.8,2x,f10.8/) 10 END 输出结果： 0.0.0. 0.0.0. 0.0.0. 0.0.0. 0.0.0. 0.0.0. 0.0.0. 0.0.0. 0.0.0. 0.0.0. 0.0.0. 0.0.0. 0.0.0. 0.0.0.