CC PAGE NUMERICAL MATHEMATICS
Found at: ftp.icm.edu.pl:70/packages/netlib/cheney-kincaid/taylor.f
C
C PAGE 307: NUMERICAL MATHEMATICS AND COMPUTING, CHENEY/KINCAID, 1985
C
C FILE: TAYLOR.FOR
C
C TAYLOR SERIES METHOD (ORDER 4) FOR SOLVING AN ORDINARY DIFFERENTIAL EQUATION
C
DATA T,X/2*0.0/, H/7.8125E-3/
PRINT 3,T,X
DO 2 K=1,128
X1 = 1.0 + X*X + T**3
X2 = 2.0*X*X1 + 3.0*T*T
X3 = 2.0*X*X2 + 2.0*X1*X1 + 6.0*T
X4 = 2.0*X*X3 + 6.0*X1*X2 + 6.0
X = X + H*(X1 + H*(X2/2.0 + H*(X3/6.0 + H*X4/24.0)))
T = REAL(K)*H
PRINT 3,T,X
2 CONTINUE
3 FORMAT(5X,F10.5,5X,E22.14)
STOP
END
.