
, 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
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 |