40 
Cc 
50 
Cc 
60 
Cc 
70 
Cc 
80 
Gien= 
Cee 
85 
Giles 
87 
Creu 
Cee 
90 
Cumin 
91 
92 
Cee 
95 
710 
720 
VICI) = YO (1) *(HC/3.) #FO(1) KUTMERS2 
IF (NPRINT .€Q.5) WRITE (65400) YI oT KUTMER43 
KUTMER44 
CALL DAUX(T*HC/3¢ V1 oF 1) KUTMER4S 
IF (NPRINT o£ Qo5) WRITE (69400)F 1 9T KUTMER4G6 
Q0 50 I=l»ND . KUTMER47 
YI(I) = YO(1) *(HC/6.) #FO(1) + (HC/60) #F 1 (1) KUTMER48 
IF (NPRINT 0 Q05) WRITE (65400) V1 oT KUTMERS9 
KUTMERSO 
CALL DAUX(T*¢HC/3.9YV19F 1) KUTMERS1 
IF (NPRINT o=Q05) WRITE (69400)F19T KUTMERS2 
DO 60 I=l»ND KUTMERS3 
YI(I) = YO(1)+(HC/8.) #F0(1) +. 375#HC#F] (I) KUTMERSG 
IF (NPRINT £EQ05) WRITE (65400) VI oT KUTMERSS 
KUTMERS6 
CALL DAUX(T*HC/269VlsF2) KUTMERS7 
IF (NPRINT o=Q05) WRITE (69400)F20T KUTMERS8 
00 70 I=l»ND KUTMERSS 
YL(T) = vO(1)*(HC/2.) #F0(1) =) SeHCeFl (1) +2, ehnC#Fe (1) KUTMER60 
IF (NPRINT 5=Q05) WRITE (65400) V1 5T KUTMER61 
KUTMER62 
CALL DAUX(T+HCoV1 oF 1) KUTMER63 
IF (NPRINT .=Q05) WRITE (69400)F 1 9T KUTMER64 
DO 80 I=1l>5ND KUTMER65 
Y2(1) = YO(I) *HC/6.#F0(1) + (20/36) *HC#F2 (1) + (HC/6.) #F 1 (T) KUTMER66 
IF (NPRINT .£Q.5) WRITE (69400) Y20T KUTMER67 
INC = 0 KUTMER68 
ese eer eee CHECK ERROR CRITERIA KUTMER69 
DO 110 1=1.NO KUTMER7O 
ZZZ = ABS(YIL(I))=-A(I) KUTMERT1 
IF (222) 85587,87 KUTMERT2 
soso = 2 = = ABSOLUTE ERRUR KUTMER73 
ERROR = ABS(.2#(Y1(I)-Y2(1))) KUTMERT4 
IF (ERROR=A(1)) 100,100,590 KUTMERTS 
cee cee eee RELATIVE ERKOR KUTMER76 
ERROR = ABS(e2-e2"Y2(1)/VI(1)) KUTMERT7 
IF (ERROR=EPSE(I)) 1009100990 KUTMER78 
Fe a ie Sule SINCE ERROR eGT. ERROR CRITERIA CHECK IF HC.GT H/KUTMERT9 
ee eee ee IF YES THEN HALVE INTERVAL. OTHERWISE STOP, KUTMERBO 
X = 128.#AS (HC) =ABS (H) KUTMERB1 
IF (X) 91,95,95 KUTMERB2 
-2f2ee-+-ee ERROR TOU LARGE KUTMER83 
WRITE (6992) I9T sERROR®HC KUTMER84 
FORMAT (/18H FOR EQUATION NO. I2927Hs THE RELATIVE ERROR AT T = KUTMERBS 
e E15e8>5 4H IS 9£15.8913H STEP SIZE = 5615-8) KUTMER86 
FIRST = 2. KUTMERB7 
RETURN KUTMER88 
=<s-e--- = HALVE INTERVAL KUTMERB9 
HC = HC/2, KUTMER9O 
IPLOC = 2#IPLUC KUTMERS1 
LOC = 2*LUC KUTMER92 
HCX = HC KUTMER93 
WRITE (29719) Ts I s9ERROR9HC KUTMER94 
FORMAT (/8H TIME = 9Ff10¢395Xs26HHALVE INTERVAL. EQUATION 9135 KUTMERSS 
6ol3H HAS ERDOR = 5£16.896X917H STEP SIZE NOW = 5615.8) KUTMER96 
WRITE (29729) NAM2 9 (Y¥2(J) 9 J=1 ND) KUTMEROT 
WRITE (29720) NAM» (Y1 (JU) 9 J=19ND) KUTMER98 
0 FORMAT( 2XeA2 / 3(10E1305/)) KUTMER99 
Go TO 390 KUTME100 
56 
