9 ANNANANNMINANANAANAANAAANN 
70 
4@0 
39 
A(K) =BLANK 
IF (MOD (K910) EGol) A(K) 200T 
IF (MOD (J9S )cEQel) A(K)=D0T 
CONTINUE 
DO 80 I=1,NFUN 
LOC=((F (1, J) -VFIRST(I)) STEP (1) 1.5) 
C=A(LOC) 
A(LOC) =CH(T) 
IF (C.NEeBLANKeAND.C.NE.DOT) A(LOC)=ASTER 
CONTINUE 
IF (MOD (J910).EQ.1)GO TO 95 
WRITE (6985) A 
FORMAT (25X9101A1) 
GO TO 100 
WRITE (6915)BoA 
FORMAT (12X1PEL2.451Xo101A1) 
CONTINUE 
RETURN 
END 
SUBROUTINE KUTMER(ND9T oHo V0 ,EPSE sAgHCX oF IRST) 
DIMENSION Y0(6) sV1 (6) 9¥2(6) 9F0(6) 9F 1 (6) oF 2 (6) g9EPSE (6) 9A (6) 
COMMON /O0UT /NPRINT oNPLOT END 
COMMON /ACCEL / XACCLoBWACL »CGACL oBL 
DATA NAM1L»NAM2 /2HY1ls2HY2 / 
ND = NUMBER OF EQUATIONS, NO, UF COMPONENTS OF YO 
T = INDEPENDENT VARIABLE 
H = INCREMENT FOR WHICH SULUTION IS TO BE RETURNED + OR = 
yO = THE VECTOR OF DEPENDENT VARIABLES. ENTER WITH INITIAL 
VALUES AT T AND RETURN WITH VALUES AT TeH 
EPSE = RELATIVE ERROR CRITERION FUR COMPONENTS OF YO .GT ABS(A) 
A = ABSOLUTE ERROR SRITERION FUR COMPONENTS OF YO .LT. ABS(A) 
NUTE-- EPSE AND A MUST BE SPECIFIED FOR EACH COMPONENT OF THE SYSTEM 
HCX = THE SMALLEST STEP SIZE USED IN THE INTEGRATION 
IF Iv IS ENTERED WITH A CHANGED H 
IF (FIRST) 20910520 
eee e ese oe FIRST ENTRY 
ees 227 S20 OTHER ENTRY. 
r 
o 
(o) 
(—) 
HCX HC 
IF (HC.NE.9.) GU TO 30 
WRITE (69800) 
FIRST = 2. 
RETURN 
e eee e ce = S CALLS TO DAUX 
0 CALL DAUX(T»Y09F0) 
IF (NPRINT.EQo5) WRITE (69400) YOoT»FO 
FORMAT (6 (2X 9F 10.4%) o4HTIME 92X oF 10.4%) 
IF (NPRINT.F£Q.5) WRITE (69400)HC 
DO 40 I=1sND 
D5) 
FIRST SHOULD BE 0 WHEN KUTMER 1S ENTERED FOR THE FIRST TIME 
AFTER THAT FIRST IS 1 IF KUTMER IS ENTERED WITH THE SAME H OR 
0 FORMAT (SXo4S5HKUTMER ENTERED WITH ZERO INTEGRATION INTERVAL ) 
PLOT2 88 
PLOT2 89 
PLoT2 90 
PLOTe 91 
PLOT2 92 
PLOT2 93 
PLOT2 94 
PLOT2 95 
PLOT2 96 
PLOT2 97 
PLOT2 98 
PLOT2 99 
PLOT2100 
PLOT2101 
PLOT2102 
PLOT2103 
PLOT2104 
PLOT2105 
PLOT2106 
KUTMER 
KUTMER 
KUTMER 
KUTMER 
KUTMER 
KUTMER 
KUTMER 
KUTMER 
KUTMERI1O 
KUTMERI1 
KUTMERI2 
KUTMERI13 
KUTMERIS 
KUTMERIS 
KUTMERI6 
KUTMERI7 
KUTMERI8 
KUTMERIO 
ODBNANFSF WN 
IF FIRST IS 2 THE ERROR SRITERIA CANNOT BE MEET AND THE STEP SIZE IKUTMER20 
REDUCED TO 47128. 
KUTMER21 
KUTMER22 
KUTMER23 
KUTMER24 
KUTMER2S 
KUTMER26 
KUTMER27 
KUTMER28 
KUTMER29 
KUTMER30 
KUTMER31 
KUTMER32 
KUTMER33 
KUTMER34 
KUTMER3S 
KUTMER36 
KUTMER37 
KUTMER38 
KUTMER39 
KUTMER4S0 
KUTMER41 
