Google
 

Trailing-Edge - PDP-10 Archives - decuslib20-04 - decus/20-0101/nonlin.exe
There is 1 other file named nonlin.exe in the archive. Click here to see a list.
?p00 `0`?hX=X?x " 
 
<++2	`}P0s}o =Xz2o((1H1,38H***CONTROL CARD ERROR*** MODEL NUMBER ,I2,10H REQUESTED)(//1H0,' *** PRINT OUTPUT FILE FOR COMPLETE INFORMATION ***  '////)(1H1,18HEND OF PROBLEM    ,5A4)+ -<oe $
 D=3 D=D &U F=E & F=0 ( H=/$(	 H=- ( H=/$( H=2( D=A F=1 $ D=E@d=C@d=B <>,n} <>c,Y? $2D=C+> =03$0&+=j+=m <>o, <>r,
+=e =00D+=o,<>U,K~,=00D+=r,<>H,K~,=00D+=t <>:,K~ =00D+=w <>-,K~ =00D+=y <>,K~6D=B+=d <>,dl =E $ \ D=W =E $&
+=e <>,o~ <>,
 0ZL=1=2;(=E=3=D=E~`;*=2 \=*=B=.=W=W=E=3=E=-=A=0=1=2=D 0Z=*2 \=B=.=E=0=1=E~00Z;*==*L=E`t2 \T8=B=2=3=3=E=-=A=0=1=2=D=D~00Z;*==*L=E`t2 \_C=B=2=3=3=E=-=A=0=1=2=D=D~00Z;*==*L=E`t2 \_=B=2=3=3=E=-=A=0=1=2=D=D~00Z;*==*L=E`t2 \^[=B=2=3=3=E=-=A=0=1=2=D=D~00Z;*==*L=E`t2 \Ss=B=2=3=3=E=-=A=0=1=2=D=D~@=E0Z;*==*L`t=0=C=1=3=2=3=E=-=A=0=1=2=D=DX=E=F$=0X=LX=E=S@=*(//30H0NONLINWOOD - DEC VERSION 1979/17H0NONLINEAR LEAST-,29HSQUARES CURVE-FITTING PROGRAM//)(31H NONLINEAR ESTIMATION,         ,5A4/23H NUMBER OF OBSERVATIONS,I5/23H NUMBER OF COEFFICIENTS,I5/16H STARTING LAMBDA,F16.3/15H STARTING NU   ,F17.3/22H MAX NO. OF ITERATIONS,I6//)(31H1NONLINEAR ESTIMATION,         ,5A4/23H NUMBER OF OBSERVATIONS,I5/23H NUMBER OF COEFFICIENTS,I5/16H STARTING LAMBDA,F16.3/15H STARTING NU   ,F17.3/22H MAX NO. OF ITERATIONS,I6//)(35H0INITIAL VALUES OF THE COEFFICIENTS)(8(4X,I2,6X))(10(4X,I2,6X))(1H ,1P6E12.5)(1H0,1P10E12.5)(53H0PROPORTIONS USED IN CALCULATING DIFFERENCE QUOTIENTS)(24H0INITIAL SUM OF SQUARES:,1PE12.4)(///35X,13HITERATION NO.,I4)(52H0EIGENVALUES OF MOMENT MATRIX - PRELIMINARY ANALYSIS)(13H0DETERMINANT:,1PE12.4)(25H0ANGLE IN SCALED. COORD,=,F5.2,8H DEGREES)(23H0VALUES OF COEFFICIENTS)(1H ,1P10E12.5)(16H0SUM OF SQUARES:,1PE12.4)(8H0LAMBDA:,1PE10.3,20X,27HSUM OF SQUARES AFTER LEAST-,12HSQUARES FIT:,1PE12.4)(54H0ITERATION STOPS - RELATIVE CHANGE IN EACH COEFFICIENT,11H LESS THAN:,1PE12.4)(63H0ITERATION STOPS ---------- CHANGE IN SUM OF SQUARES LESS THAN:,1PE12.4)(40H0NUMBER OF ITERATIONS EXCEEDS MAXIMUM OF,I4)(19H0CORRELATION MATRIX)(3X,7(4X,I2,5X))(3X,10(4X,I2,5X))(I3,7(2X,F7.4,2X))(I3,10(2X,F7.4,2X))(1H1,28X,45HNONLINEAR LEAST-SQUARES CURVE-FITTING PROGRAM//1X,5A4,2X,9HDEP.VAR.:,1X,7HMIN Y =,1PE10.3,2X,8HMAX Y = ,1PE10.3,2X,10HRANGE Y = ,1PE10.3/)(1H ,18X,18A4)(19X,22HOBSERVATIONS DELETED: ,7(A6,2X),7(/41X,7(A6,2X)))(19X,24HOBSERVATIONS NOT FOUND: ,7(A6,2X),7(/43X,7(A6,2X)))(1H0,62X,21H95% CONFIDENCE LIMITS/1X,10HIND.VAR(I),3X,4HNAME,5X,9HCOEF.B(I),6X,10HS.E. COEF ,4X,7HT-VALUE,4X,22HLOWER            UPPER)(1H ,4X,I2,6X,A4,A2,2X,1PE12.5,5X,1PE9.2,4X,0PF6.1,5X,1PE9.2,5X,1PE9.2)(20H0NO. OF OBSERVATIONS,11X,I5)(20H NO. OF COEFFICIENTS,11X,I5)(28H RESIDUAL DEGREES OF FREEDOM,4X,I4)(26H RESIDUAL ROOT MEAN SQUARE,3X,F16.8)(21H RESIDUAL MEAN SQUARE,8X,F16.8)(24H RESIDUAL SUM OF SQUARES,5X,F16.8)(1H1,18HEND OF PROBLEM    ,5A4)(1H0,11H***ERROR***/32H0EXCESSIVE NUMBER OF VARIABLES. ,36H PROGRAM DIMENSIONED FOR MAXIMUM OF ,I3)(1H0,11H***ERROR***/35H0EXCESSIVE NUMBER OF COEFFICIENTS. ,36H PROGRAM DIMENSIONED FOR MAXIMUM OF ,I3)(1H0,11H***ERROR***/26H0LESS THAN 2 COEFFICIENTS.)(1H0,11H***ERROR***/35H0EXCESSIVE NUMBER OF OBSERVATIONS. ,36H PROGRAM DIMENSIONED FOR MAXIMUM OF ,I4)(1H0,11H***ERROR***/55H0NUMBER OF COEFFICIENTS EXCEEDS NUMBER OF OBSREVATIONS.)(1H0,11H***ERROR***/33H0EXCESSIVE NUMBER OF ITERATIONS. ,44H PROGRAM LIMIT SWITCH SET FOR MAXIMUM OF 100)(1H0,11H***ERROR***/36H0STOP CRITERIA IS NEGATIVE. STOPSS =,F12.7,10H, STOPRC =,F12.7)(1H0,11H***ERROR***/12H0COEFFICIENT,I4,18H HAS NEGATIVE SIGN)_P 0 P!\Rw$='E=WN(!	F]^BT=8):s\-6+8-3U<c\O.: \Kq !p @>| !p @?W !p @F  !p @?P !p @IZ !p @IX !p @>~ p @I\ !p @? !p @?U !p @IR !p @I5 !p @?[ !p @?8 p @F p @?p !p @F p @?r p	 @IR p	 @? p
 @F p
 @?O p @?O p @?Y p
 @IN <QH,sM <QC,sM <Q@,sM <Q<,sM <Q9,sM <Q5,sM <Q2,sM <Q.,sM <Q+,sM <Q',sM <Q#,sM <Q,sM <Q,sM ;( DF <QK,, <QN, <QQ,
 <QV, <QY,
 <Q^,, <Qa,, <Qd,"F $ FF DIQ <Qg,
:IQ; F+L1, <Qi,"F $ FF DIQ <Ql,
:IQ; F+L7, <Qm,"F $ FF DIQ IQ(.?N <Qp,
:IQ; F+L=, <Qr,"F $ FC DIQ IQ(.?N <Qu,
:IQ; C+LE, <Qv,, <Qy,"F $ FC DIQ <Q|,
:IQ; C+LN, <Q~,"F $ F?t DIQ IQ(.?7 <R,
:IQ; ?t+LT, ?Y2d?r+P^ F2dIN+Pa $2dF+Pc;&F $ D?T "O B?6 "O% BIV "M\ BIV F$F7 .$(dK} &?S ?O2dIR+Pe ?O3$F+Ph $23$;(+Pj7$;%+Pl5dLq7$;$+Pl5dLo "O, B?6+Ls "O% B?6+Ls7$;$+Pl5dLs "O, BIV@d?@f? <Q,?8"?O $ ( 
.
?Z .?1 .IJP4@0 (.?Z
 L 
?
?.$8&Lv <R, <R,
 <R, <R
,
 <R, <R,
 <R, <R,
Fd;'F"F $ F?r D?q(.?NI4 ?q( .?N .?7 H ?q(@h?@j? ?q$?V.?W $ DKr <Q,?8"?O $ ( 
?V$
?q.
.
?W 
.
?1,0
(
( ?q( 
(
 JKr 
Kr.
?Z ?V$?q.Kr.?W@P0 ? ?.$8&M ?q( .?N 
.
?7
0L ?h? ?q(.?NI4:?q; ?r+M
"F $ FIK DIQ"IQ $ FI6 D?q@dC@fC"?O $ ( 
?V$
?q.
.
?W ?V$IQ..?W@P0CC.$8&M7 IQ( ?q( .?N 
.
?7
0L  .?7L  .?NL ChI4 &$IQ ?q(.?_ &$?q IQ(.?_:?q; I6+M5 IQ( ($IQ. ( ?_ HKr <Q,s?7:IQ; IK+M3ZIV1BO+O1BO%+O%1BM\+M\1BO,+O,1BN+N"F $ FI2 DIQ"IQ $ *$
 IQ(. ?_C .$IQ 
(
.
0Bs .$ 
IQ(
.
0Bs.$8&M:IQ;I2+M^$RDIL@dF!F<Q FF  DIQ &$IQ IQ(. <R,
:IQ; F +Ms, <R,"F $ FF DIQ &$IQ IQ(. <R",
:IQ; F+M|, "N BIV"F $ F>} DIQ"IQ $ IQ( *$
.
 
(
 ?7P0?7
(?_l
(Bs *$
IQ (. .$ 
IQ(
.

0Bs
 Bs.$8&N:IQ; >}+N"F $ (
 ?l ?7
 ?[ (
 ?[
 I5 *$
 (.
F
 Bs
 Bs.$8&N <Q	,bZ <R#, <R&,
Ku?M"F $ (
 ?[L?Ml ?7 .?N
@
 ?.$8&N#"F $ DIQ IQ(.?57, +Pn4LN1 IQ( 
.
?N
0L ?4lN?:IQ8&N*@d?X@f?X <Q,?8"?O $ ( 
.
?Z .?1 .IJP4@0 (.?Z
 L 
?X
?X.$8&N5@dIZ@fI[INF"F $ (
 ?[L I5
IZ
IZ (
 ?[L ?[
IN
IN (
 I5L I5
F
F.$8&NBINDFKs <Q,sIZdKr <Q,o\@KyC <R(, <R+,
"F $ DIQ IQ(.?57, +Pn4LN\ IQ( 
.
?N
0L ?4lN_:IQ8&NU?Xd?Kt4dNdK{C5dNb?MdKx?M+N";'DFF+N <R,,, <R/,,"F $ ( 
.
?N ?0.$8&Nh <R2,"F $ F?5 DIQ <R5,
:IQ; ?5+No, <R7,"F $ F?3 DIQ <R:,
:IQ; ?3+Nu, <R;,"F $ F?0 DIQ IQ(.?N <R>,
:IQ; ?0+N{, <R@,"F $ F?. DIQ IQ(.?N <RC,
:IQ; ?.+O, <RD, <RG,
 <RI, <RL,
Z?61BO+O1BO%+O%1BM\+M\1BO,+O,1BN+N"F $ DIQ IQ(
 ?[L?Ml ?77 , .?N@7 4Kzl,;%5lO :IQ8&O <RN, <RQ,
 <RR, <RU,
+O0ZIV1BO+O1BO%+O%1BM\+M\1BO,+O,1BN+N$;$?$?X5dO, <RW, <RZ,
 <R[, <R^,
+O0?X?:?T3d;(+M <R`, <Rc,
?X? ?O/F D? <Rd,, <Rg,, <P~,bZ <Rj, <Rm,
 <Ro, <Rr,
"F $ ( ,$
.
 ,0?_ LKr <P},s ?7.$8&O< <Rs,"F $ F? DIQ <Rv,
:IQ; ?+OE, <Rx,"F $ F?V DIQ <R{,
:IQ; ?V+OK,"F $ F?U DIQ"&/F.IQ IQ ( *$
IQ.
 
IQ(
0?7P ?7
(?_l
(Bs ( *$
IQ.
 
IQ(
 LKr.
?N Kr.?7@P0 
.
?7P0 
.
?NP0
(?_l
(?_ *$
 IQ(. .$IQ 
(
.

0?_
 ?_ *$
 IQ(. .$IQ 
(
.

0Bs
 Bs.$8&OR <R|, <R,
"IQ $ F?R D?q &$IQ ?q(. <S,
:?q; ?R+Ot, <S, <S,
"IQ $ F?R D?q &$IQ ?q(. <S,
:?q; ?R+O},:IQ; ?U+OP7$?+Ph4DP\h? *?d>| <P|,s?/K|IP(IP?2"?O $ H?Q DIQ(.IJ $ DKr <Pz,e_:IQ; ?Q+P$IP?2F <S, <S,
7D;)+P <S,";) $ F?Q DIQ IQ(.> <S,
:IQ; ?Q+P, $3D;)+P43D;*+P+ <S,";* $ F?P D?M ?M(.IO <S,
:?M; ?P+P$, ;)3D;*+P4 $.;* D?q <S,"&/;).?q ?q FC D?M ?M(.IO <S,
:?M; C+P/, <S,, <Py,d?"F $ ( ?7H?/IT 
(
.
?N
0lIT
?ZH?? 
(
.
?N0?0F 
(
.
?N0?1 ( H?4/( H?q <S ,(.I\ 
?4(
.
I\ (.?N <S#,
.$8&P8 <S(, <S+,
 <S,, <S/,
d?K|D DF! <S1, <S4,
 <S5, <S8,
 <S:, <S=,
 <S>, <SA,
+Pr <SC, <SF,
+Pq <SH, <SK,
+P\ <SM, <SP,
+P\ <SQ,,+P\ <ST, <SW,
+P\ <SY,,+P\ <S\,,+P\ <S_, <Sb,
+P\ <Sd, <Sg,
+P\ $ DF Kq I\ Ap F Ap ?r Ap IR Ap	 ?O Ap IN Ap
,~x?hIP?2Krx>|xKrP?tF?\F!?pINxKsxKr0IX?IR?OFF ?rIRINI\PCF?\IL?pINXCF"?SFINxKr0IX>~Kr?OFF ?rIRINI\0IX>~IR?OFF ?rIRINI\PKwIM>|IMKwIRPKwIK?WIJKwIR@KwI5F I3Kw?r?IRPKw>?P>Kw?PKw>~IZI\KwFPKwIYIXIXKw?OPKwIW>~?NKwINPKw?L??7KwINPKw?7?U?5KwINPKw?4IR?1KwIRPKw?/I5?ZKwIR@Kw?Y?[?WKwIR?VINPKwIQFIOKw?pXI]	XIfDIYIX$?O$F0;&0;'$;(XI\IyDIYIX$?O$F0;&0;'$;(XJXI\JXJ$IQXI\J$IQXJ0XI\J0XI\JXI\J$IQXI\J0XJ0?XI\J0?XI\J$?TXJ$?TXJ"XI\J"XJ0BsXI\J0BsXI\J(0?pXI\J+0CXI\J0XJ0XJ$IQXI\J$IQXJ0XI\J30XJ50?XXI\J80F0?XXJ@0;%XI\J@0;%XI\JH0;$XJH0;$XJP$;(XJUXI\JUXJ(0?pXI\J(0?pXJX$IQXI\JZ$IQXJ\$IQ0BsXI\J^$IQ0BsXI\J`DIYIX0IP0?20FXI\Jo0XI\Jp0XI\Jv0XI\J|
$0 000@0IT0?Z0F0?1XI\K$?OXI\K$FXI\K$F!XI\K0?/XI\K 0>|XI\K$0?XI\K(DIYIXXI\K+
$?rXI\K5
$INXI\K@XI\KE
$IRXI\KOXI\KW
XI\Kb	0;$0;%XI\Kk$IQ[>$[ \Sp !p @Sh !p @Sl p @Si !p @Sk p @Sj p @Sp p @Sn <T,sM <T,sM <T
,sM"Si $ So 
(
.
Sm Sk$.Sl 0@ PSq 0  PSq <T	,qE@ 0.$8&S Sp,~pSqSqPSrSoShSoSrSnPSrSnSlSmSrSp@SrSmSkSlSrSjSkSp 0 }\=8(zpR^<f3Lf3^K
f3Lf3_'Lf3Lf3_/a#kGW`kBG.
[>$[
,\T. !p,@T !p@T%p@T"!p @T# p$@T.p@T+p@T)<U,sM <U
T%$
.
T&
0
T6lT2
T T%$.T&@T7pT2T T%$.T&PT3tT2T T%$.T&`T5xT2THTT"T!HTTT!HTTLT
TTHTT'PTT- T),TL T!P 
TL 
T X T"P 	,TL 

0TP ,T'L 

TP 

T TP T!T TP TT TP TP 0 T,lT,
T. <U,qET* 
(
.
T' 
T%$
.
T&T*0T. ,  LT/ <U,qET p
T,T,

(.$8&TF T.,~pT.T/pT.T$PT1T*TT)T1T)PT1T(T%T'T1T+@T1T&T#T&T1TT%T+END                 (5A4,I1,1X,3I2,2X,I2,F4.2,F4.0,1X,I3,F4.3,2F8.7,2I2,2I1)(1H1,42X,45HNONLINEAR LEAST-SQUARES CURVE-FITTING PROGRAM//31H 1979 VERSION OF THE NONLINWOOD,I3,10H VARIABLE,,I3,13H COEFFICIENT,,I4,20H OBSERVATION PROGRAM//55H REFER TO FITTING EQUATIONS TO DATA BY DANIEL AND WOOD,,33H SECOND EDITION, WILEY PUBLISHER,/49H FOR GLOSSARY OF TERMS, USER'S MANUAL, DETAILS OF,44H CALCULATIONS AND INTERPRETATION OF RESULTS.//1H0,5A4/)(1H0,25X,24HCONTROL CARD INFORMATION,39X,14HORDER OF CARDS,//1X,49HCOL.  INPUT  MAX ITEM   (NOTE: BLANK ON CARD = 0),20X,13HFIRST PROBLEM/2X,4H1-20,12X,26HIDENTIFICATION OF PROBLEM.,26X,16H1  CONTROL CARD.)(1H ,3X,2H21,I3,9X,40H0  OBSERVATIONS READ FROM CARDS OR FILE.,12X,53H2  FORMAT CARD (72 COLUMNS, ORDER: IDENT, Y, X(I)'S),/18X,36H1  REUSE DATA FROM PREVIOUS PROBLEM.,20X,31HE.G. (A6, F6.0, (NOIND*F6.0) ).)(1H ,5H23-24,I3,5X,I3,1X,39HNUMBER OF COEFFICIENTS TO BE ESTIMATED.,13X,39H3  DELETE-OBSERVATIONS CARD(S), IF ANY.)(1H ,5H25-26,I3,9X,41HFILE NUMBER IF DATA ARE TO BE READ FROM A,11X,45H4  STARTING VALUES (GUESSES) OF COEFFICIENTS /19X,42HSEPARATE FILE.  NO END CARD(S) IF ONLY ONE,13X,50H(10 COLUMNS / COEFFICIENT, 7 COEFFICIENTS / CARD)./19X,28HSET OF DATA IS ON EACH FILE.,23X,42H5  INFORMATION CARDS FOR PRINTOUT, IF ANY.)(1H ,5H27-28,I3,5X,I3,1X,46HNUMBER OF INDEPENDENT VARIABLES TO BE READ IN.,6X,54H6  NAMES OF COEFFICIENTS CARD(S) FOR PRINTOUT, IF ANY.)(1H ,5H31-32,I3,9X,30HNUMBER OF EQUATION TO BE USED.,22X,40H7  DATA CARDS (IF NOT READ FROM A FILE)./70X,51H8  END CARD (END IN COLUMNS 1 - 3 OF IDENTIFICATION)(1H ,5H33-36,F7.3,5X,44HSTARTING VALUE FOR LAMBDA, F4.2, (E.G. 0.1),,12X,47HFIELD).  THE NUMBER OF END CARDS MUST EQUAL THE/19X,42HUSED AS A MULTIPLIER TO SCALE THE SPACE OR13X,44HNUMBER OF CARDS PER OBSERVATION.  (END CARDS/19X,20HSIZE OF STEPS TAKEN.,35X,45HARE NOT NEEDED IF DATA ARE READ FROM A FILE).)(1H ,5H37-40,F7.3,5X,30HVALUE OF NU, F4.0, (E.G. 10.)./19X,40HDIVISOR AND MULTIPLIER TO CHANGE SIZE OF11X,14HSECOND PROBLEM/19X,45HLAMBDA DEPENDING ON WHETHER SUM OF SQUARES OF,9X,38HIF DATA ARE REUSED FROM FIRST PROBLEM,/19X,38HITERATION IS NEAR OR FAR FROM MINIMUM.,19X,38HDELETE THE FORMAT, DATA AND END CARDS.)(1H ,5H43-44,I3,9X,44HMAXIMUM NUMBER OF ITERATIONS, I2, (E.G. 20).,11X,18HIF DIFFERENT DATA,)(1H ,5H45-48,F7.3,5X,37HMULTIPLIER USED TO INCREMENT VALUE OF,21X,19HREPEAT 1 - 8 ABOVE./19X,31HCOEFFICIENTS, F4.3 (E.G. 0.01).//18X,48HNOTE: IF VALUES IN COLUMNS 33-48 ARE NOT DEFINED/24X,41HON THE CONTROL CARD.  THEIR LEVEL WILL BE/24X,43HSET AUTOMATICALLY TO THE ABOVE E.G. VALUES.//18X,42HCRITERIA FOR ENDING CONVERGING ITERATIONS.)(1H ,5H49-56,1PE9.1,3X31HSUM OF SQUARES CRITERION, F8.7,/19X,45H(E.G. 0.0001, A CHANGE OF LESS THAN 0.0001 IN/19X,29HTHE RESIDUAL SUM OF SQUARES).)(1H ,5H57-64,1PE9.1,3X,38HRATIO OF COEFFICIENTS CRITERION, F8.7,/19X,43H(E.G. 0.001, A CHANGE OF LESS THAN 0.001 IN/19X,43HTHE RATIOS OF ALL COMPARABLE COEFFICIENTS).//18X,47HNOTE: VALUES IN COLUMNS 49-64 CAN BE SET AT 0.0/24X,44HIF CONTROL OF EITHER OR BOTH IS NOT DESIRED./)(1H ,5H65-66,I3,5X,I3,1X,50HNUMBER OF INFORMATION CARDS TO BE READ FOR DISPLAY/19X,40HON PRINTOUT IF DESIRED, 72 COLUMNS EACH.)(1H ,3X,2H68,I3,9X44H1  READ NAMES OF COEFFICIENTS FROM CARDS FOR/21X,42HDISPLAY ON PRINTOUT, 1ST 6 OF 10 COLUMNS //21X,23HCOEFFICIENT,  7 / CARD.)(1H ,3X,2H69,I3,9X48H3  PLOT RESIDUALS VS. EACH INDEPENDENT VARIABLE.)(1H ,3X,2H70,I3,5X,I3,1X,48HNUMBER OF DELETE-OBSERVATIONS CARDS, OBSERVATION/19X,48HIDENTIFICATION IN 1ST 6 OF 10 COLUMNS, 7 / CARD.)(1H0,16X,35HSUBROUTINE ARGUMENTS AND DIMENSIONS/1H ,16HSUBROUTINE MODEL,I1,35H (NPROB, B, FY, NOB, NC, X, NVARX, ,20HNOBMAX, NCMAX, KTOU)/1X,33HDIMENSION  B(NCMAX), FY(NOBMAX), ,15HX(NVARX,NOBMAX)//)(18A4)(7(A6,4X))(1H0//1X,6H******,I2,32H DELETE CARD(S) WITH (7(A6,4X)) ,55HFORMAT REQUESTED, NO OBSERVATION IDENTIFICATION FOUND. ,6H******//)(A6)(1H0//1X,6H******,35HDATA SAME AS PREVIOUS PROBLEM, NEW ,31HDELETE CARDS READ BUT NOT USED.,6H******//)(7F10.0)(18A4)(7(A4,A2,4X))(1H ,18A4)(13H0DATA FORMAT ,18A4//24H OBSV. NO.  IDENT.  SEQ.,2X,107H 1-11-21    2-12-22    3-13-23    4-14-24    5-15-25    6-16-26    7-17-27    8-18-28    9-19-29   10-20-30)(1H0//////1X,39H******WARNING--THIS VERSION OF PROGRAM ,15HDIMENSIONED FOR,I3,11H VARIABLES,,I3,16H REQUESTED******//////)(1H ,2X,I4,5X,A6,1X,I4,10F11.3)(19X,I4,10F11.3)(1H1//24H0PROGRAM DIMENSIONED FOR,I4,16H OBSERVATIONS.  ,I4,8H READ IN/17H0PROBLEM SKIPPED.)(1H0,42HOBSERVATIONS THE SAME AS PREVIOUS PROBLEM.)(1H0//////1X,45H******THIS VERSION OF PROGRAM DIMENSIONED FOR,I3,14H COEFFICIENTS,,I3,11H REQUESTED./17H0PROBLEM SKIPPED.)(1H0//////1X,45H******THIS VERSION OF PROGRAM DIMENSIONED FOR,I2,26H DELETE OBSERVATION CARDS,,I2,42H REQUESTED.  NOT ALL OBSERVATIONS DELETED.,6H******//////)(21H0OBSERVATIONS DELETED)(10(6X,A6))(23H0OBSERVATIONS NOT FOUND) 0"_T=8(:pQua!9Lf3Lf3e
J7 \Y7 p @UJ !p @UF !p @U. !p @U; !p @U" !p @UI !p @UH !p @U  !p @U !p @U# p @U p @UB p @U. !p @U7 p @U0 p @U4 p @UF p @U p	 @U8 p	 @U! p
 @U! p
 @U/ p @UE p @UC <\2,sM <\/,sM <\+,sM <\(,sM <\$,sM <\!,sM <\,sM <\,sM <\,sM <\,sM@dU U0&& FU U&&	 FU <\6, <\9,
 <\C, <\F,
 <\J,, <\M, <\P,
 <\Q, <\T,
 <\V, <\Y,
 <\[, <\^,
 <\`, <\c,
 <\d, <\g,
 <\i, <\l,
 <\m, <\p,
 <\r, <\u,
 <\v, <\y,
 <\{, <\~,
 <\, <],
 <], <],
 <]	, <]
,
,<] <],
 <], <],
@h3H;&2J;&+Z
Y=
;&3H;'2J;'+Z
Y<
;'2D;(+Z &
 F;(3HU2JU+Z
Y;
U ;(2fUC+Z"
;( $ (.UY7  (.UU .$8,Z U@1D+ZB <], <],
 $2DU+Z# U. FU D;)3DU9+ZA U FU"3&U9+Z' U9 HU""U" $ &.f;)/& FU1 <],"
;)/*.
U1 U1 JU6 FUD 
UD(
.
UD <],
:UD; U6+Z.,.$8(Z( ;)(.UD3DU02FU1+Z:>;)5DZ4 <]!, <]$,
 U93dU+ZA:U""U9/&.U" U" <]%, <](,
.$8&Z>@d;*+ZI7DU9+ZI"U9 $ <]*, <]-,
.$8&ZD <].,, ;(3dUC+ZQ <]1,";( $ FU5 DUD <]4,
:UD; U5+ZM,+Z`@dU# $.dU#/$ DU1 U#3f;(+ZV ;( FU# <]6, U1/&/U# U1 FU4 DU U(.U <]9,
:U; U4+ZZ, U#3$;(+ZQ7D;)+Zs U FU" DU#3dU" DU" U"$$	 D;) <]:,";) $ HU3 DUD UD(.U <]=,
:UD; U3+Zh, U#3dU+Zs"U.U# DU""U" $ <]?, <]B,
.$8&Zp ;(( DU3 U 0F+[ ;(3fUC+Z~ <]C,"U3 $ FU@ DUD <]F,
:UD; U@+Z{,+[ <]H,"U3 $ FU? DUD UD(.UJ <]K,
:UD; U?+[,+["U3 $ (.UJ
U0
 .$8&[7D;)+[ <]L,";) $ FU= DUD UD(.U <]O,
:UD; U=+[, U@1D+[h <]Q, <]T,
 U42$U/+[ <]V, <]Y,
 $3dU/+[ U/ DU/+[ $ DU/:U <][,( .UB 
.
U? <]^,
"U/ $ FU< DUD U$U UD(..U= <]`,
:UD; U<+[%, U(.UB3DU$2FU$7+[c U$U.U=3DY:2FY;7+[c ;*3D;)+[=.$ DU2 U(.UB U2(.UD 3H2J+[=:;*"&.FU+[ $ DU9 <]b, U( .UB 
.
U? <]e,
"U/ $ FU< DUD U$U UD(..U= <]h,
:UD; U<+[E, $ DU2 U/3dU2+[ $.U2 DU $.U22$U/+[S $.U2 DU+[T U/ DU:U9 <]i, <]l,
"&.U/U U FU; DUD U$U UD(..U= <]n,
:UD; U;+[Y, $.DU2+[L U3DU.+[c <]o,"$.U DU!3dUF+[j <]r, <]u,
+Ya <]w,,+\ ;(3dUC+[n <]z, <]},
+Ya7D;)+\ U93dU+[s <], <^,
7D;*+[~ <^,, <^,";* $ FU DU U(.UD <^
,
:U; U+[x, ;*3D;)+\ <^,, $.;* DU2 <^,"&/;).U2 U2 FUI DU U(.UD <^,
:U; UI+\,+\ $ DUB Y7 U Ap UB Ap U! Ap
 U/ Ap
 UE Ap,~PY9UCUFUBY:UFPY9UAU.U?Y:UF@Y:U>U;U=Y:U4UUFPY9UU"UY:UPY9UUIUJY:U8PY9UHUHUGY9U!PY9UGU UY:UCPY9UUUY:UCPY9UU#UY:UCPY9UEU7UDY:U0XU.\
UKDUHUH$U@$;($U$U/$U0;&0;'$;(0U0;$0;%$;)$U$UE$U9XUJUQ%$U4$UC$UFDUHUHXUJUvXUJV$U@XUJV $;($UCXUJV+$UXUJVK$U/$U4XUJVY$UXUJVi0;&XUJW0;'XUJW'	$;(XUJW1"0UXUJWS0;$XUJWb0;%XUJW~
XU.X\0U:XU.X\0XU.X]0XU.XQ$U#XU.X^0U:XU.X^0XUJX`0XUJXa@	U%XUJXr
$U4$U/XU[aU%	0 000XUJX~$U0 $U9000XUJY$U90`UXUJY	$UF$UXUJY
$UC$;(XUJY$U$U9XUJY/XUJY20XUJY4XUJY20(///52H0***** MODEL 1 SPECIFIED BUT FORTRAN SOURCE WAS NOT ,52HMADE AVAILABLE FOR THIS RUN.  EXECUTION STOPS. *****)[>$[ \^% !p @^ !p @^ p @^ !p @^ p @^ p @^ p @^ p @^ <^?,sM <^<,sM <^8,sM7`^+^6 <^D,, <^7,o~ ^%,~P^%^^^^&^P^%^^^^&^@^&^^^^&^^^X^^
(///52H0***** MODEL 2 SPECIFIED BUT FORTRAN SOURCE WAS NOT ,52HMADE AVAILABLE FOR THIS RUN.  EXECUTION STOPS. *****)[>$[	 \^Y !p @^F !p @^F p @^L !p @^F p @^G p @^K p @^J p @^L <^s,sM <^p,sM <^l,sM7`^L+^j <^x,, <^k,o~ ^Y,~P^Y^K^F^J^Z^JP^Y^I^F^I^Z^K@^Z^H^F^H^Z^G^G^KX^L^M
(///52H0***** MODEL'3SPECIFIEDBUT FORTRANsSOURCE+WASNOT ,52HMADEAVAILABLEFORTHIS RUN.  EXECUTION STOPS. *****)[>$[	 \
(///52H0***** MODEL 4 SPECIFIED BUT FORTRAN SOURCE WAS NOT ,52HMADE AVAILABLE FOR THIS RUN.  EXECUTION STOPS. *****)[>$[
 \_A !p @_. !p @_. p @_4 !p @_. p @_/ p @_3 p @_2 p @_4 <_[,sM <_X,sM <_T,sM7`_4+_R <_`,, <_S,o~ _A,~P_A_3_._2_B_2P_A_1_._1_B_3@_B_0_._0_B_/_/_3X_4_5
(///52H0***** MODEL 5 SPECIFIED BUT FORTRAN SOURCE WAS NOT ,52HMADE AVAILABLE FOR THIS RUN.  EXECUTION STOPS. *****)[>$[
 \_u !p @_b !p @_b p @_h !p @_b p @_c p @_g p @_f p @_h <`,sM <`
,sM <`,sM7`_h+` <`,, <`,o~ _u,~P_u_g_b_f_v_fP_u_e_b_e_v_g@_v_d_b_d_v_c_c_gX_h_i
 0 P p K&<[U \`8 !p @` !p @`3p`, p @`1 p B` `1$` D`& `.$ D`+ `1/$ D`% `1$`1,s2`/"`1 $ F`1 D`7 `7/`&&`/& `7 (.`3@l'@n'.`8&`P:`7; `1+`M &/`&&`+/& $ (.`3
`9
'.`+8&`W`,d`/`  "`o B`@d`$@f`% $ D`"@f`( &/`1 $ F`. D`7 `.D`".$.D`( `"/&/`( `" (.`
'L'
`$
`$.$8&`f:`7; `.+`bZ`1B`o+`o1BaZ+aZ"` D`("`% $ F`+ D`) `.f`(.`( D` `).$ D`7 F`"&/`1.`7 `7 F`* D` `.f`.` D`*.`) F`7 `.`( D`2 `7(.`H(` 4haV `(.`` `*(.`'` `7(.`'`6,`
`
`&H`=H`6L`&

`8 <a^,s`46 `&+a$`4`4`&d`4``<H`6h`4`#`9D`>`8 <a],s`H`>h``!``D`>`-`(`
`L`>L(`#H`6`" `(.``"' `*(.``-(`" `7(.`@h@j `2(.` `/` D` `)/` D`5 ` D` `( D`5"`1 $:`:
`5 
`.L`.L`5(.`30'P`!(
.
`3
/L`

`2 `(.`3'T` `5(.`3GP`!' `5(.`3
' 
`/
 `)/$4HaU `(.`,'L`! `5(.`'H`` 
`(
.
`7T` `5(.`GP`!7 
`5(
.
`7 
`(
.
` `(.`G7 
`5(
.
`7.$8&a.:`; `*+`y:`); `++`r "aZ B`+`]`$$`,5d`o `8,~x`8x`8_P 0y`[%[ \bT !p @bK p @ab !p @ab p @aapbS p @bM <ct,sM <cp,sMbWbS bY@fac*dbd@dbR DaabVac"ab $ FbL DbH6@ac+bw"ab $ Da` a`6@ac+bv 
bQ$
a` bH(..bR
 7 ,
a_3Lac2nac3,ac+bv bH DbN a` DbJ
ac:a`8&bk:bH; bL+bi7`ac+cj bJ bN Fac3FbJ+c(bSbS"ab $ bQ$ 
bN(
.
.
bR
0
a_ bQ$.bR bN(. bJ(.P@ bQ$ bJ(..bR
@.$8(b~"aa $ 
bQ$
 bN(..bO a_ 
bQ$
.
bO bN(. bJ(.P@ bQ$ 
bJ(
.
.
bO0.$8&cbQ$bJ bJ(..bRaDbSbS:aa bQ$bJ bJ(..bRbW bQ$bJ bJ(..bRdaa"ab $ 
bQ$
 bJ(..bR
a_L 
 .$8&c,"aa $ 
bQ$
 bJ(..bO
a_L 
 .$8&c2"ab $ FbH Da`3DbJ+cT bQ$bJ(..bRa_ bQ$bJ a`(..bR@h@j"ab $ bQ$ 
a`(
.
.
bR 
bJ(
.
.
bR,0La_
(
(.$8&cC"aa $ bQ$ 
a`(
.
.
bO 
bJ(
.
.
bO,0La_
(
(.$8&cL:a`; bH+c:+bf bR b DbJ ac FbN"ab $ Da` 
bQ$
bN a`(..bR a_ a`( 
bQ$
bN.
.
bR bQ$bJ..bR@0 bQ$bJ 
a`(
.
.
bR0:a`8&cZ"$.DbR6@bR+cVaab/aa,Daa bT aa ApbSp,~@bXbSbKbRbXbMbQbM@bXbPabbObXbMbQaa!4^oNlE q4V Y:/
$nKc R
$nc)}z Qi<5_Y QMi<5_Y Q7%q~|vd Q#kBW. Q=8(zpQv Q7%qT~} PzpQ5a#kB PpQu!#kBG Pg6"(9,
 P_Y!e0 PX	]FT PQua#kBGW PKc)=ymHZ PFS{3[4 PADM)<5_ P=8(:pQua P9,
$nKcPP5Y!e P1:/
=8) P+	Rxj ?.
=8) P@A P+aoNlEa
?>;iR.=\dDdp@cy0`+d"cy(cx+d3$<2$cy+d$d+d3$3$cy+d*".cy`d	d"+d3$3$cy+d/".cy`dd"+d3".cy`d(1H0,44H---------ORDERED BY COMPUTER INPUT----------,10X,53H----------------ORDERED BY RESIDUALS-----------------/9H OBS. NO.,4X,6HOBS. Y,5X,8HFITTED Y,5X,8HRESIDUAL,10X,8HOBS. NO.,4X,6HOBS. Y,5X,8HFITTED Y,2X,20HORDERED RESID.  SEQ.)(1H ,1X,A6,F11.3,2F13.3,11X,A6,F11.3,2F13.3,3X,I4):s\-6+8-3g?@ \de !p @dG !p @d; !p @d4 !p @dF !p @dB !p @d;pd?pd7 p @dG p @d8 p @d9 p @dH <eM,sM <eI,sM <eF,sM <eB,sM <e?,sMdjd?(d?d7"d9 $ Hd> DdD(.d@ $ Dde <e=,e:dD; d>+e d; & F/d9 $ Fd= Dd5(.d6dA d5/$ Dd<"d< $ Dd8 d8.d;  Dd9 d9(.d6( dA4he:d88&e d5.d; d5 F+e d5 Dd:>d:.d;  F d:2dd8+e d8.d; d5 F:d5; d=+e	 <eP,, $.d9 Dd?"d9 $ d?/ Hd5.d; 
' Jd: <eS, ( Hde(
 Jdf de.dE (  Hdf de.dD (  Hdg de.d@ (  Hdg de.d6 (  Hdh df.dE (  Hdh df.dD (  Hdi df.d@ 
df.
d6 <eV,
.$8&e# ded?pd7p,~hd?d7dePdidFdGdEdjdGPdidEd;dDdjdGPdidCd4dCdid8PdidBdFd@djdGPdid6dBd6djdGXdHXdHd`1df1dg1dg1dh1dh1di0 00$[&uX<pe]pe^pe\e\$e]5$eee\e]e\$e^4dege\e^e]pe^p,~EX3e.9GG(P5,*UF Z^FUpzFYu2pZzq6JgVcw&xkE02*3S8Bh= ~OQ]U1j>YXL,^\yr%<uWkpzo3~
ajM,A'f?nPh2-4V47>a:J;G=zNPhLRfRTw>DW8YTA\K^|+$puCur5!s|9{uJXuw
&xx^{zVZ|7yc~5J P    O    S    I    T    I    V    E    R    E    S    I    D    U    A    L    S    N    E    G    A    T    I    V    E    +++++(1H1,1X,5A4,33X,36HCUMULATIVE DISTRIBUTION OF RESIDUALS)(1H ,17X,101H.0002  .001  .005  .01 .02   .05  .1    .2    3  .4  .5  .6  .7   .8    .9   .95  .98 .99 .995   .999)(1H ,19X,25A4,A2)(1H ,17X,A1,1X,25A4,A2)(1H ,17X,1H0,1X,25A4,A2)(1H ,5X,A1,13X,25A4,A2)(1H ,5X,A1,11X,1H0,1X,25A4,A2)(1H ,54X,33HCUMULATIVE FREQUENCY, NORMAL GRID)(1H1,1X,5A4,39X,21HRESIDUAL VS. FITTED Y)(1H ,13X,F10.3,40X,F10.3,40X,F10.3)(1H ,66X,8HFITTED Y)(1H1,1X,5A4,30X,35HRESIDUALS VS. INDEPENDENT VARIABLE ,I2)(1H ,13X,F10.3,2(14X,F10.3,16X,F10.3))(1H0,58X,20HINDEPENDENT VARIABLE,I2)8 P~0( 0	Rxj?h|h:s\-6+8-3@@!Y@( h 0
\=8({!y a'$: \g6 !p @f` !p @ej !p @ej !p @f^ !p @f !p @fH !p @f[pfdpf !p @fF !p @f p @fz p @fT p @fq p @f p @fZ p @f p @f% p	 @fH p	 @fp <kt,sM <kq,sM <km,sM <kj,sM <kf,sM <kb,sM f% Dfb & Ffz & Ff	d Dft fb.fH  Dfl(.elfv fH  Df<(.el(HgBhg?fIf$fdDgBdgCfh"fb $ Dek ek.fH  DfG fG(.el 4(gy fG(.el( fvhfI,s<.( Hfl+gz:ek8&go <k_,m^"fb $ FfY Dek.fH  FfG(.el$fvdfI,s6.$ Dfc0$0d+h2 fG(.f(fdhfh,s<.( HfQ0(0h2+h0H+h
 & FfQ fQ.& Ff;&&  Ff
7(f;.(( Hfa5Fh $ Df
"$.Dfa $$fa.fc.fF $e Dg6 <k],n`dek$ft Df2dgA+h  g=ZDfm fm f2&f+h*dh $ Dfm fm Df_+h& g=ZDfm"fm.& Ff[!(f3(f*dh! $.fm Df_ $.f_ DfS&$  Df
7&fS.&( Ffa5Dh. $ Df
"$.Dfa $$fa.fc.f $e Dg6 <k[,n`:ek; fY+g} <kx, <k{,
 <k},, <l, gA_DfWZDek g:_DfVZDf
 $$f
.ek.f <l,
:f
; fV+h;:ek; fW+h9, <l, g?_DfUZDek ek <l,
 g:_DfSZDf
 $$f
.ek.f <l	,
:f
; fS+hG:ek; fU+hD, g; &3Ffl+hW <l, g:_FfjZFf
 &$f
.&.f <l,
:f
; fj+hR,+h_ <l, g:_FfiZFf
 &$f
.&.f <l,
:f
; fi+hZ,*dhN g9 &3Ffl+hj <l, <l,
 g:_FfeZFf
 &$f
.&.f <l,
:f
; fe+he,+hs <l, <l,
 g:_FfcZFf
 &$f
.&.f <l,
:f
; fc+hn,*dh` g@ &3Ffl+h} <l , g:_Ff_ZFf
 &$f
.&.f <l#,
:f
; f_+hx,+i <l$, g:_FejZFf
 &$f
.&.f <l',
:f
; ej+i,*dht <l), g>_DfxZDek ek <l,,
 g:_DfwZDf
 $$f
.ek.f <l-,
:f
; fw+i:ek; fx+i, <l/, g<_DfuZDek g:_DfuZDf
 $$f
.ek.f <l2,
:f
; fu+i:ek; fu+i, <l3,, <l6,,f3Dfd2Ffd7+j <l9, <l<,
fddg9ff <l?, <lB,
 <lD, gA_DfsZDek g:_DfrZDf
 $$f
.ek.fF <lG,
:f
; fr+i+:ek; fs+i), <lI, g?_DfqZDek ek <lL,
 g:_Df	ZDf
 $$f
.ek.fF <lM,
:f
; f	+i7:ek; fq+i4, g; &3Ffl+iG <lO, g:_FfZFf
 &$f
.&.fF <lR,
:f
; f+iB,+iO <lS, g:_FfZFf
 &$f
.&.fF <lV,
:f
; f+iJ,*di> g9 &3Ffl+iZ <lX, <l[,
 g:_FfZFf
 &$f
.&.fF <l\,
:f
; f+iU,+ic <l^, <la,
 g:_FfZFf
 &$f
.&.fF <lb,
:f
; f+i^,*diP g@ &3Ffl+im <ld, g:_FfGZFf
 &$f
.&.fF <lg,
:f
; fG+ih,+iu <lh, g:_FfFZFf
 &$f
.&.fF <lk,
:f
; fF+ip,*did <lm, g>_DfAZDek ek <lp,
 g:_Df@ZDf
 $$f
.ek.fF <lq,
:f
; f@+i{:ek; fA+ix, <ls, g<_Df;ZDek g:_Df:ZDf
 $$f
.ek.fF <lv,
:f
; f:+j:ek; f;+j, <lw, <lz,
 <l},, $2Dfp+kQ Df	"fH $ Ff: Dfog=f(ffJ"fb $ fn$ 
fo(
.
.
fn ,0 Lg6 <kZ,q7f fn$ 
fo(
.
.
fn ,0 Lg6 <kX,q*fJ.$8(j fo(fel fo(fJf%:fo; f:+j@dfO"fH $ FfV Dfof3HfJ2JfJ7+kP:fO f<(.el(fvHg;hg?fI(elf fo(f%fJ5dj8,f
fJLg;lgC
fh+j;$ffJDg;dgCfh <kV,m^"fb $ FfU Dek.fH  FfG(.el$fvdfI,s6.$ Dfc0$0d+jZ fG$fn fO(..fn(fhfh,s<.( HfQ0(0h2+jZ0H+jN & FfQ fQ.& Ff;&&  Ff
7(f;.(( Hfa5FjV $ Df
"$.Dfa $$fa.fc.fF $e Dg6 <kT,n`:ek; fU+j> fo(/$ DfQ <m, <m,
ffJdg9fyfJdg9fwfyfdg9fW <m, <m	,
 <m
, gA_DfTZDek g:_DfRZDf
 $$f
.ek.fF <m,
:f
; fR+jk:ek; fT+ji, <m, g?_DfRZDek ek <m,
 g:_DfPZDf
 $$f
.ek.fF <m,
:f
; fP+jw:ek; fR+jt, g; &3Ffl+k <m, g:_FfOZFf
 &$f
.&.fF <m,
:f
; fO+k,+k <m, g:_FfeZFf
 &$f
.&.fF <m,
:f
; fe+k
,*dj~ g9 &3Ffl+k <m , <m#,
 g:_Ff`ZFf
 &$f
.&.fF <m$,
:f
; f`+k,+k# <m&, <m),
 g:_Ff^ZFf
 &$f
.&.fF <m*,
:f
; f^+k,*dk g@ &3Ffl+k- <m,, g:_Ff]ZFf
 &$f
.&.fF <m/,
:f
; f]+k(,+k5 <m0, g:_Ff\ZFf
 &$f
.&.fF <m3,
:f
; f\+k0,*dk$ <m5, g>_Df\ZDek ek <m8,
 g:_DftZDf
 $$f
.ek.fF <m9,
:f
; ft+k;:ek; f\+k8, <m;, g<_DfsZDek g:_DfrZDf
 $$f
.ek.fF <m>,
:f
; fr+kE:ek; fs+kD, <m?, <mB,
 <mF, <mI,
:fo; fV+j* g6 fz Ap,~h	g6f
fp`	fF	fflf	pfJg6pfg6h	g6f
fph	g6f
fp`	fF	fflf	@g8fof`fng8fTfnfqPg8fmejfkg8fPg8fjejfig8fZPg8fgf^fg8fqPg8elfelg8fqPg8ekf[fyg8fXfzf{Dekf[XfzgXfzg
(eXfzg(f:(eXfzg
(eXfzg(eXfzg(f6(eXfzg(f6(eXfzg
(eXfzg(eXfzg$f5(eXfzg
(eXfzgXfzgXfzgDekf[Xfzg"0fd0ff0fXfzg
(eXfzg(f:(eXfzg
(eXfzg(eXfzg(f6(eXfzg(f6(eXfzg
(eXfzg(eXfzg$f5(eXfzg
(eXfzg"0fd0ff0fXfzg&Xfzg(Dekf[$foXfzg.0f0fW0fy0fw0fJXfzg
(eXfzg(f:(eXfzg
(eXfzg(eXfzg(f6(eXfzg(f6(eXfzg
(eXfzg(eXfzg$f5(eXfzg
(eXfzg.0f0fW0fy0fw0fJXfzg2$fo     I     I   I--- -I   *--- -*-- --*- ---* ---- *-   --      ' ~0~@|hXOJL@ !p @mL !p @mK p @mZ p @mK6bmK0b+me+mc4me4n34nA m]DmRZDmJ m]"(.mJ (..& FmS$&.mL mS He mS$&.mL mS H*dmg:mJ; mR+mg mL mV F F0 mK mS H mT H mU H5 mT HO mT Hj H mU H mU H9 mS HT mT Hn mU H	 H# mU H> HX mU Hs mU H
 m\_DmQZDmJ m\ &$&.mJ.mL mL He &$&.mJ.mK He*dn mJ.mK mW F=:mJ; mQ+n m\ &.mL mM H &.mL mN H &.mK mM H &.mK mN H*dn m\ &$&.mZ.mL mU He &$&.mZ.mK mU He*dn$ mZ.mL mN F mZ.mL mO F mZ.mK mN F mZ.mK mO F+n\ m]_DmMZDmJ m]"(.mJ (..& FmS$&.mL mS He mS$&.mL mS H*dn5:mJ; mM+n5 mL mV F F0 m\_DmYZDmJ m\ &$&.mJ.mL mL He*dnC:mJ; mY+nC mZ0D1+nO m\ &.mL mM H &.mL mN H*dnJ+n\ m\ &.mL mM H &.mL mN H*dnP m\ &$&.mZ.mL mU He*dnU mZ.mL mN F mZ.mL mO F,~a0 \n_ p @n^ p Bn] p Bn_ @n^ <nk,oB n^ Dn^ n_ n^ Ap n] Ap n_ Ap,~hn_n]n^(A10)(' INPUT THE NAME OF THE OUTPUT FILE')(' INPUT THE NAME OF THE INPUT FILE')DSK: SEQINDSK: SEQINOU   _B-` \nw <o,, <o, <o,
 <o,, <o, <o,
:<o,<o,nw,~XnsX~no0nnXnoX~no0nmP~nx^nyPnnP~nz^n{PnmW&,~xoX	`Z(X^. >X.~(3~+oi	~+ouZ. \x <x*<x .ez+op!"	b+  "or	b+Z`.Z. h+ggM?}
&ow	b+oe,
?NO CORE AVAILABLE - CAN NOT PROCEED
a.3@1,>sZO@p9+pgR~@1,>sZ@@p9,>,>@6 w+p
"p:0B+p!"A`X#p5@p/#p!"X#p+p
\p+p
a.9'1,>sZ,>,>O@p9+p
,>
6+ _4`p"`Q@,> <p;,,^,> <p>,
,^*bp+p,> <p?,,^7@p9+p,> <pB,
,^+p,> <pD,
,^Z5"p,> <pE,
,^,7@p9,
,> <pG,,^,,> <pJ,,^,> <pM,
,^(@qd1@,,
1@#+p.0@*+p,>q<pN,,^+p,^,^,^,~('TypeGtoContinue,XtoExit,TToqtrace.'/2H*,$)(A1)(1Hq14A5/)(1HA5,1X,O6)PAUSESTOP.pXp6pXp7pbp8p"p9p"Xp0X~p5p xC~01,>sZ#p2`p+p]4@p\+ "`p @p`,> <pa,q,^p+  @p`,> <pa,p|,^7!pp`,~ p,~+Pov0pb@,~ 0 9~j Y~jxp`ACOS OF ARG > 1.0 IN MAGNITUDECNM`1,>sZ#p2bpv+ps"pp+ pv4@pq @pu,> <pw,q,^ p Bpu,> <pw,p|,^,~ pv7ap"`,~+Pov0px@,~ 9~j 0xpuASIN OF ARG > 1.0 IN MAGNITUDECR
q%#p+q  ,~>attempt to take SQRT of negative argI6 1,>sZ,>,>e|Pwp+q1 9p2``+q12@`+q13"``*|q.a|+q5 ph8 ax~+q.,^,^,~I6Md 1,>sZ,>,>e|Pwp+q> 9p3 `+q>2@`+q>2b``*|q;a|+qB ph8 ax~+q;,^,^,~IFgp =p+qGIFI,>,>,>7@p+qc4@qkqCp# $)D0dM+qV/$@#5HqV(FZ $) + qV5FqV p7 "`,>p DpqC,qw,^p+qaqC5 qY +Pov0qlqC,> <qt,ro,^qCp Dqj@qC+ qdqC,> <qt,r+,^ qj,^,^,^,~!  "+qa qj4@qc4 qh+Dov0
qgUU"+qa+Dov0qi@`+qa6 p+qa+qfAttempt to raise negative double precision number to non-integer powerxqCIFGp =p+qwIFI7@p+r
r	*,r6,p+r#`!& (f+"r!  "+rF+0r(LeL+r@+0r5Lq!
r	*,,~7 &+r5&raL&+rr	!  ",~6 p,~+Dov0
qyU U",~+q~r&!vrX 7
Ql\|((080619.jEu] 3?A#
HoZ!um."{EtH^AZIs7% 0f3
f3M"bb,7,IF1,>sZ r* r(p4@rH#2$r+rI@r_({j(7 m ~d ~))'.(F)D'Qd"X@rF ,2B0r+r:(dr9+r90&+rL3"0r=
r7/0r(L)|l"rMr"fr(r)r%fr)@r$r#r(r!`r4LrF@0r@!
r*,,~! +rFU  r&7!p   +rF(L0r+rF Ix51(3 TT.C}DyyH}|ant0bx!b{x\G}x9ThC(UL
<O'@6x4Q.+-xA&a'AK`;50,(&#! _/Qg^5<8qc8r`0R$I$I$I`*U*UjU*U+YLf3Lf3`u*U**U*U_p 0rfI2|t(1,>sZ #p7B+s	,ro@rm,~{e^b72I2|p1,>sZ rf* rhp5`rt4@s	+Pov0s!p3@rc4Bs&s/&@FM (FrZred ~)@f,s%0rV %0rVl".eb. (
0rNd`+s)+ry)|h rd&rerer[ ,r\F0.,1lrc+s@re!
rf*,,~! ,~Attempt to take DLOG of Negative Arg@`+s6 >
 WsdA psL!	Ge!8ld Zm[7 zm[74rxtrx[6m[{6m[sIO-1,>sZ s* sp5`s4@s-+Pov0s-!p 
(Jr/*@(*XJs+(Jnd ~@(? (s"(s(s"(s(s (s`Fs`@!
s**,~Attempt to take DSQRT of Negative ArgIad@ &)|h$Fs4,~S1d@,>_({5$s9$hd$~)w27 "`,^,~S1e,>_({5(s?(hd(~)w27 "`,^,~g:9d@4$sG`fmD+sK AsL,~$`fmD+sK"KsL,~",~`CR@1,>sZ,>,> p  p  Ep<bsV $Ep$ E. .&+sR".!p Ap $Ep,^,^,~sZY\s\Q<X\s\*+6s]6`s]S	YVBG9s\G9s\
GJsZGK/KsZUexhUc*4ZYYXsK	0sC
sP)G
VL`
BS	YVB-sguaesguaZs-guaZs,guaYs+
guajssPU7RMSQsqswkgs]sguaj|0sG
s	CGrfW%
rerc(Q8rN(P6[rm+w3rVP)yP)G
W@
BS	YVB	
jrdq{rtcN6MrZqzwsqzvsqzvsqzuryqzuruq{ro
c>friPU7RMSQwkgrfar[`r\`r]r^r^r`^ra]rb]rcq{z0rNCGr=l?rsJrF4|kr74|er:4|er9
D,rrfr#P)yP)Gjr%
W@
VL`
Bklr)S	YVBcN4crd^Zr+=GrF`2rr$PU7Klr(rr"wkgr&^r!d^Z{ r;grCGr	;-9r
X!.r24:rP)yP)GS	YVBSo&B
_=Gqu
_=L%qwrr
r
q
<<Vqw	q~PU7+T2rRMSQwkgr	=L%{`quXqYd^OqV
_=G"qE
=L&qG0kgqjPU7RMSQTEt
d:q7dfqBdfq>deq>deq;PU7d:~q6
dzq*d&q5d&q1d%q1d%q.PU7dz~q)G
q$aeqTvFq!aZqaYqajqP3?K~qBqqqDAp}DRp|I;qI:q
Y]@p9
Y?p
\!p{-^p
{-Xp{-Op{-Np.{-Np{-Mp{-Mpp0S	YVBSo&BTE)|o{\6.AwBP@VKV4PbXlVkBI
YPbKYYYh[|ggMA
JTUR/}F
/|N	@\:T(q5vX_<F

zR(TsT	`u6
(J"ov
&*
 R	WJ
>
_Bb
_2v	^,6RW	`
K^
2b
2RYmtL(bIzoe)	Gob[o_ood	qGoZFtfoY0KoZ05o[06o[\x]"oXoX
\x]"o\\x]"~0o\GoT[oQooW	qGoM\y@(oLFtfoK0KoM05oM06oNoJ
\y@(oO\y@(~@oN'&oF
_.
+oBu=oI_-$K@oBGo2[o.oo?o;o7+o3	qGo$Y1YMo$w04o#w05o#mCMo"w05o"w06o!w06o!w07o w07o w08o0Ko$	!ow0Ho
u^c;o&u^c;} o&zAnozUnoGozins[n~oooo+o?oSo
_o
so	qGno
Hnn
Cnm0Knw2xn}2x}Pn}Gnf[ndonk	qGn_n_n^0Kn_n^
n]3*nn`3*n~Pn`nJnP
an#'nU)n;=nAQnCynLGn\!nO[mb5nRomc!nUmfmg+mg?n
Sn
_n
sn	qGm[m[9mZPHmZPImYPImYPJmXPJmXPKmWU6JmW"KmSmSP4mRP5mRP5mQP6mQP6mPP7mPP7mO"KmMP8mMUJmL
g`}mL
ZomK
g`~mKmJ"Lm^"Lxm]h04g6i]3kj>om&Ag
ijT%skf9KlOl3ti^3jEm)IgyihW
?Cheiou7m
m/
ah'ia Ak'h
?WlipuKjim0ag;ik+k]?kl%ahjci	uk#,'h@hmlI&[iwzjkKm3
tAj$,Oh!zSl?i4'lp{msk7
k5,ck[zghtSlL'iz{jsk8&1hs3jo,wkxz{l?i6'+i{{/jtsm8hk;h$Gjz-3l{7l#{i>,?j{kjw/m9'i&AjV[k2h9{Kl$lO,Sj{mCm;)h%oklisiI-7lz;mm>Gh>,Qg"kSk?3hCiiJ-Kl}Omm?!h.[hK?k>3+hD/l,slV-jck	'mBohJ,yiGSkJ3?lCi
iP-sjwk
;mF!h2hU	hAkgkH3ShFWilX2_jmOmIh]-!iO8?hGkl-/l[2skZk+hh kQ8Sl	l/CiT3kX3m ?hq2!iZ
Ak- #gY8ghNiWiU3j*m#qf{Sh{ 7kt8{likl\3/j*3k
_i2Iicj8 Kkq9hQil^3CkVGkg
si
Ag2^km9#hR'ila3Wj=[m$	qGf{*ACfzB}hfyw0JfyP]fxY9fwP]fw
":fpw04fo	fow05fnw05fnfmw06fmfl9flPHfkw06fkPIfjw07fjPIfiw07fiWLfhPJfgw08fgH+hffPJfeQ,feHthfd+vGfcPKfcQ-fbYJITfb
jf^Q.f]Q/f]Q/f\Q0f\f[mf[k2LfZP4fZP5fYP5fYP6fXXfWP6fWQfVP7fVQfUP7fUQfT0Kg6	!fT
fQQfPQfPfOQfO4fKCEMfJ]lfI
g`}fF*#fAQfAQf@&Bf<ybf<Qf;Xzf
"f

Zof	Ppf	CFhfPqf
g`~fPqfPrfPrfPsfPsfPtfdfHsMf':fSBLfw0HfPhemw0Ielw0Ielw0JekekP\ej
=?1Me=?1M~@eh08dih07dhh07dhh06dgh06dgh05dfh05dfh04de
ae
uee8Ge9[dxoeMeIeF+eB?e?Se
_e
se=	qGdH*ACdHh=OdGY1YMdGdFw04dFw05dEw05dEdDw06dDw06dCw07dCdBw07dB;-9dAw08d@Hthd?0d?d>P4d>P5d=P5d=P6d<0Kde&E
ad 
ud"Gd3[d	qGdc3d0Kd05d]cyycyc3dc3}`d!c2oAck
acCcZcZtactAcccjac[b`octcp{%bwbd
acT+bh?biSbk
_bkabz
sb~	qGbT_5bSw04bS	bRw05bRw05bQbQw06bP9a`PHbPw06bOPIbOPIbN&AbNQbMP4bMP5bLP5bLbK:"	a_P6bKP6bJXH7bJP7bIP7bI0KbTP8bH
bH{y}adOMacabza_09jab
a`faaa_
PJaa
8YvPbZ8YvPv0bYh05`9h04`8
aaU`oa)aGa\[`Eo`M`M`P+`W?`bS`b
_`f
s`q	qG`8`7	`7.`6`59`5`4`3m`2.`2m`1P4`1P5`0P5`0X!`/P6`.`-P6`-`,`+P7`+P7`*.`*0K`8P8`)
`)})`(`(`'`&)`&E`%E`$9LA`#`"})`"9LN`!` -`	`i`
`%`&&`
`a`n``Re`-`Rr`
z	+J`?z	+Ju@`?5`I_i@`3`G`[_o``
`+`	qG_i*AC_hf_hY1YM_gw04_gw05_fmCM_fw05_ew06_ew06_dw07_dw07_c0K_u	!_cA|p_v5_OI_5@_R3_SG_R[_Ko_[_X_T+_`	qG_5*AC_4f_4Y1YM_3w04_3w05_2mCM_2w05_1w06_1w06_0w07_0w07_/0K_A	!_/
ACA|`B5I@3G[o'$ +,	qG*ACfY1YM^w04^w05^~mCM^~w05^}w06^}w06^|w07^|w07^{0K
A_A|`_5^gI^M@^j3^kG^j[^co^s^p^l+^x	qG^M*AC^Lf^LY1YM^Kw04^Kw05^JmCM^Jw05^Iw06^Iw06^Hw07^Hw07^G0K^Y	!^G
A^[A|`^Z5^3I^@^63^7G^6[^/o^?^<^8+^D	qG^*AC^f^Y1YM^w04^w05^mCM^w05^w06^w06^w07^w07^0K^%	!^A|p^&2qXra[M)Y
X!	X\3[%\y9KZMO]`IZ3[ZZ
?CZc
a['\']+]?kZgQX~c\,']@]=[c
Z*u	[j\C,cZzg]BS]u
XC e[Q3\F,wZz{]C?]w
)Z4 y[SbEYu1Y\J-Z{ZzS]zbAX.
=XD!\M-]{#Z{g]}
QZ:G\P-3]{7]F{][\Q2Z#{K]H^
yZ>&![Lo\T23Z'{_[#^
XQ\V2GZ({s[7^
!ZA{[n\Y2[]]KK[x+\[2oZ.s[[xzAUK
!ZB3Z1{)Y?\^3Z.]Ls^
zUUQ
5ZDGZN{=[sS\`3][^
IXR[Z\?\c3+]!/[s^oZjz}Zp,y[/#ry\S\d3?]$C]O\
j[]K\r9]-]^
s[{ _\u9#].'[%	qGUK*ACUJw0JUJP]UIyUIVZq;$w0KUHzUHw0KUGw0LUGf;(h=OUFY1YMUFqCUEw04UE	UDw05UDmCMUCw05UC'3UBw06UBPHUAw06UA
SZ3U@PIU@w07U?PIU?w07U>PJU>w08U=PJU=PKU<PKU<U;PLU;k;(\7U:YG>U9&EYU9k2LU8P4U8
L<kU7P5U7P5U6P6U6P6U5P7U5P7U40KY7	!U4P8U3yU3NU2
U2U1go;'UJU0@G|kU0 +';&LU/dGU/U.*@HU.g[U%GU$LU#U#U"bU"zU!XzU!QYU U |k;)w0\U"Uw0]U63wUw0]Uw0^U/R;)w0^UT
/G;*@G~>U
U*B`UFRrUsU%(U%2USBLUw0HUw0IU*"Uw0IUVZnR;%w0JUP\U
aUGU[T@oUUT.UT-T,Y1YMT+T*w04T*w05T)mCMT)w05T(T~T'w06T'w06T&w07T&QT%w07T%6T$T#TiT")T!P4T!T T0KT.	!T)T
T*TzTTjTT*TT~T+TTjT
AT8AzT8h05Sqh04Sq
aTGT[S{oTTT
+S?T		qGSpY1YMSpw04Sow05SomCMSnw05Snw06Smw06Smw07SlQSlw07SkSkP4Sj0KSp	!Sj
SizSiSh
ASsA|Ssh06Ksh05Ksh05Krh04KryNSdDAIyO'0AJaMPuMK5Nh&WLD?R
aLk!uJ@'L'Qv?WMs1uJ^;QH+Qy?kR
OQC,LM?RR>&GR{usS(%aLoQuJ\cQ@,'LN@M{R@&[OOzRS+R{/S4Lh{RL,?S{kS:&AO0biPKuP+[Q#2R{KR#O,SS{S=&!N"PaJvAKoQ23R{R&#O,gS|S>,LnQ2GR{sN#7RN,{O}SAQK2[R
ON*KRQ-O}SCz3$Pr+QN2oRsN*sRR-#S'SF3L2zAMF)K?QQ3RQsRU-7P|;SHGL8SQV3RN5SRW-KP
WQR^2_SSQLVPc	hUN1{Qd8?QkR(/R`2sSST,AM\1Pe+Mw	hiN-
aPL L18SMR+CRc3S3SWEPh2!Ls?M	h}N?tUN\)Nd	O< #L18gM2NUWRd3PSYYPjtiNX=N_SNp3]K( 7Qg8{M3NUkRg3/P3S\mPlN#QNb
_Nv KQi9M4R,P~3CSGS_PneN`
sN~JL79#M5'R/Rj3WS[Sb	qGI]*ACI\w0JI\P]I[:-IZyIZVZq;$w0KIYP]IYzIXw0KIXP^IWw0LIWw|A(IVqkIVP^IUf;(PIUM/@TITPISP`ISY1YMIRIRw04IQ	IQsphIPw05IO:-INmCMINw05IMw06IMPHILILw06IKPIIKw07IJ%1I6PII6I5w07I5;-9I4PJI3w08I3PJI2F"PKF"
L<kFP5FP5FCP6C?tQ?sP7?s0KKq	!?rP8?r
?qgo;'Y{?p@G|k?pdG?Yw0q?Y:2?Xw0q?W?Ww0r?VQ?VQ?UL?U?T0]?L?8w0]?7w0^?7&>h?6Pp?6/R;)w0^?5Pq?5T
/G;*w0_?4
?4Pq?3soM?2[?1w02?1Pr?09+|?/w0`?/Pr?.Ps?.R?w0I>w0I>a>~VZnR;%w0J>~P\>}2Z>|&.K~&.A0K~(R=Xh05=Wh04=W
_=t
s>:	qG=F*AC=EVZq;$y=z=*f;(h=O0Y1YM=EqC=DmCM=D'3=C2Hth=Bk;(Y:4=Bk2L=A
L<k=3	!=3go;'@G|k=2Ltb;*z=1Xz=0L|k;)63w=0/R;)T
/G;*@G~>=/FRr=/%2`HsM=.SBL=-VZnR;%
\|@v;*J\|@v (;*
\|@u;$J\|@uh;$
\|@u0J\|@u.@0(Ru`=X\B:/\>X\>5'\7N+\2g'r'I(\)8*\)7\	\	)\+,\6.\.,[rO.O[mNO[mO[m	O[l|[l{&O[l+O[kO[kO[j]O[jRO[i%r%H([Z[X%[Qgr%0[Ms/[M"-[M	+[L,[FG[@v*r#eOD^w'dw]!+{0