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