C
C PAGE 306: NUMERICAL MATHEMATICS AND COMPUTING, CHENEY/KINCAID, 1985
C
C FILE: EULER.FOR
C
C EULER'S METHOD FOR SOLVING AN ORDINARY DIFFERENTIAL EQUATION
C
      DATA  A/0.0/, B/1.0/, S/0.0/, N/128/
      F(T,X) = 1.0 + X*X + T**3
      H = (B - A)/REAL(N)
      X = S
      T = A
      PRINT *,T,X
      DO 2 K = 1,N
        X = X + H*F(T,H)
        T = T + H
        PRINT *,T,X
   2  CONTINUE
      STOP
      END