Numerical Methods

Approximating the Second Derivative

Verlet's Method

Substituting the approximation above into , we get

In Fortran, this would be implemented as

DO I = 1, NSTEP
   X(I+1) = 2*X(I) - X(I-1) + H**2 * F(X(I))
ENDDO

Summed Form

Substitute DX(I) = (X(I+1) - X(I))/H into Verlet's method to get:
DO I = 1, NSTEP
   DX(I) = DX(I-1) + H*F(X(I))
   X(I+1) = X(I) + H*DX(I)
ENDDO


Applications | Mathematical Models | Numerical Methods | High Performance Computing | Animation


Maintained by H. Edward Donley <hedonley@grove.iup.edu>
Last Modified on Saturday, 24-May-1997 12:55:12 EDT