48 HPF³ÌÐòÕªÒª - II

---

• subroutine runga-kutte(h, ns)
• use nbodyvars
• real h; integer ns
• INTERFACE
• subroutine Grav(X, M, A)
• use nbodyvars
• real, dimension ( :, : ), intent(in) :: X, M
• real, dimension ( :, : ), intent(out) :: A
• end interface
• . . .
• do k = 1, ns
• Xdelta1 = V
• call Grav(X, M, A)
• Vdelta1 = A
• . . . .
• end subroutine

---

Copyright: NPACT