29 Fortran90的N体Runge Kutta子程序-I
子程序runge-kutta(X,V,M,h,ns)
C 使用Runge-Kutta方法，解决牛顿的运动方程
C 是全局4阶的，X和V是初始的位置，并且
C 速度,系统对时间间隔h*ns逐渐逼近
C X和V包含每次修改的状态
real,dimension(1:3,1:N) :: X,V,Xdelta1,Vdelta1,Xdelta2,Vdelta2,Xdelta3,Vdelta3,Xdelta4,Vdelta4
real h,h2
integer ns,k
INTERFACE
function Grav(X,M)
real array(1:3,1:N) :: X,M,Grav
end function
end INTERFACE
C Grav是难解决的并行算法，后面将给出 ！
Copyright: NPACT
