GENLS.FORވMS Fk LLIBFORE0LO CV7DGROUP GENLS_TEXTCODE_DATADATACONST_BSSBSSFAR_DATA GENLS5_DATA$SECDERSECDER$A$EXTRAEXTRA$A$DATADATA$A$ININ$A$WTIWTI$A$DONDON$A$EXTRYEXTRY$A(f,HHHN h x xØx՘x٘xp Fxx V @ESFJSRQQFISRQQFIERQQFIDRQQFIWRQQ __acrtused __FF_LDIN __FFpPAUSE __FFsSTOP__FFrdPWGMMCALCPRINTOUT_mainOBSWTSFUNC__FIlog __FCclenv__FFopn __FF_NRECTMI __FF_NRECTMO __FF_NRECFSI __FF_NRECFSO __FF_FMTOUTDERXE__FFwr__FHfexp __FF_DECFLTPLTR__FIsqrt __FF_ENCFLTPEARSN__FFcls__FIsign__ftol QGENOUTC$1BAD Y DATA FOR LN TRANSFORMATIONw WEIGHT CHOICE IS NOT VALID!! GENOUTR! GENOUTT& GENOUTF;X"  ) N T T n)T8]"!@@ MENU #1: DEFAULT CHOICES @@]">1: ACCEPT ALL DEFAULTS AND CONTINUE; MIN OUTPUT; ST.-LINE FIT]"2: CHANGE CONVERGENCE CHOICES]"3: CHANGE PRINTING CHOICES]"*4: CONTINUE TO MENU #2 - SET INPUT VALUESo v ' ~ )N[T7%]"MENU 1: DEFAULT CHOICE =H" )T84]"$MENU #1A: CHANGE CONVERGENCE CHOICES]"'1: RESET MAX # DEMING ITERATIONS, MAXD]"(2: RESET MAX # GENERAL ITERATIONS, MAXG]",3: RESET DEMING PARAMETERS CONV. LIMIT, EPS]"+4: RESET GEN. PARAMETERS CONV. LIMIT, EPS1]"65: TURN ON/OFF LEVENBERG-MARQUARDT CONVERGENCE METHOD]"6: RESET CONVERGENCE FACTOR, G]"17: RESET RELATIVE INCREMENT FOR DERIVATIVES, FAC]"&8: RESET DERIVATIVE CONV. LIMIT, EPS2]".9: RESET CONV. LIMIT FOR SUM-OF-SQUARES, EPS3]"010: RESET ALL DEFAULTS IN 1A TO ORIGINAL VALUES]"11: RETURN TO DEFAULT MENU 1{  ڜ  )8NT7%8]"MENU 1: DEFAULT CHOICE =H͠ )ZT84Z]")RESET MAX # OF DEMING ITERATIONS (DEF=10)  T 4 )ۜT85]"*RESET MAX # OF GENERAL ITERATIONS (DEF=36)(  0^ )T8@]"5RESET DEMING PARAMETERS CONVERGENCE LIMIT (DEF=1.D-6) 1@T 4 )MT8>]"3RESET GENERAL PARAMETERS CONVERG. LIMIT (DEF=1.D-8) 1D )2T8,2]"!ENTER 1 TO TURN ON; 0 TO TURN OFFΠ  +7 )ZʜT86Z]"+RESET CONVERGENCE FACTOR G (DEF=1; 05 FOR THIS CHOICE AND PICK]"&PROPER NOUT VALUES TO DEFINE FUNCTION:]"3[1 + A1X +...+ A5(X**5)]/[A6 + A7X +...+ A10(X**5)]  )T T7&T ]"LN TRANSFORMATION OF Y DATA )v T8gv ]"'NOW SELECT NEW FIT FUNCTION, e.g. NF =1]"/USUALLY UNITY WEIGHTING CHOICE APPROPRIATE HERE ) T8] ]"4MAKE SURE THAT THERE ARE NP INITIAL VALUES SPECIFIED]"ENTER NEW NP VALUE (< 11)  T 4  )4 NΜT7(4 ]"NEW TOTAL # OF PARAMETERS =H"Ҡ )X T8"X ]"ENTER NEW VALUE OF NOUTޠ  pT 4  )v NrT7(v ]"NEW # OF FIXED PARAMETERS =H" ) T81 ]"&ENTER INDICES OF NOUT FIXED PARAMETERSb  ]h !t ) t) #TT TT ? ]")SET/RESET INITIAL VALUES OF NP PARAMETERS]#  Π,  ) ) NNNNNT TT TT 8 ]"GENLS PROGRAM V.1.0]FN]"NUMBER OF OBSERVATIONS =H]"NUMBER OF PARAMETERS =H]"NUMBER OF FIXED PARAMETERS =H]"FUNCTION CHOICE NUMBER =H)") 111111 TT T T ]"*PARAMETER CONVERGENCE TOLERANCE (DEMING) =# ]"+PARAMETER CONVERGENCE TOLERANCE (GENERAL) = ]"G = CONVERGENCE FACTOR = ]"-H = RELATIVE INCREMENT FOR NUM. DERIVATIVES = ]""DERIVATIVE CONVERGENCE TOLERANCE = ]"S CONVERGENCE TOLERANCE = H )) TT]"NOTE:]X )N3)N3ɜ TTl ))[ TT!4]" PARAMETERH" IS FIXED AT THE VALUE #.Ѡ |)T8)]"NAT. LOG. TRANSFORM OF Y DATA֠ )gT8>]"3INITIAL OR CURRENT FITTING VALUES OF THE PARAMETERS)X G T TX]#  x )|T8 )dT8jd]"MENU #3: SELECT NEW WEIGHTING]":0: RESET X-Y DATA TO INPUT VALUES AND RETURN TO THIS MENU]"1: RETURN TO DEFAULT MENU #1]"#2: USE LAST CURRENT VALUES FOR FIT]"73: ALL SIG(X) CONSTANT; ALL SIG(Y) CONSTANT: SIGX,SIGY]"44: RATIO OF VARIABLE X AND Y SIGMAS CONSTANT: IDWMC]")5: GENERAL WEIGHTING CHOICES FOR X AND Y]"6: EXIT PROGRAM&  I  )N1T7(]"MENU 3: WEIGHTING CHOICE =H"b  x )NT8]"=DETERMINE FORM OF GENERAL ITERATIVE WEIGHTING; LET Q = X OR Y]"9THEN FOR EACH POINT: WGTQ = 1/[UQ*UQ + PQ*{|Q|**(2*XIQ)}] )rT8/r]"$ENTER ADDITIVE ERROR TERMS UX AND UY"1015 )1015}T=!]"UX =# ]"UY =#  )QT8/]"$ENTER CONSTANT PX; ENTER CONSTANT PY1+1+) TGG]" Tà1]"STANDARD DEVIATION OF FIT, SF =#+ R)՜T8]"MENU #5: AFTER-FIT CHOICES]"&1: RETURN TO MENU #1: DEFAULT CHOICES]"$2: RETURN TO MENU #2: INPUT CHOICES]"%3: RETURN TO MENU #3: WEIGHT CHOICES]"94: PLOT RESIDUALS ON SCREEN AND IN GENOUTR; RTRN TO MENU]"45: SHOW RESIDUAL-RELATED QUANTITIES; RETURN TO MENU]":6: MAKE OR ADD TO RESIDUALS FILE(GENOUTT); RETURN TO MENU]"87: MAKE OR ADD TO OUTPUT FILE (GENOUTF); RETURN TO MENU]"8: EXIT FROM PROGRAM Z ܜ b)T8']"$PRODUCE VARIOUS RESIDUAL ROUGH PLOTS]",USE POSITIVE VALUES FOR SCREEN-ONLY PLOTTING]"9* USE NEG. VALUES FOR SCREEN AND FILE (GENOUTR) OUTPUT * ]":DEPENDING ON WEIGHTING CHOICES, NOT ALL PLOT CHOICES BELOW] "ARE ALWAYS OPERATIVE]"1: X-RESIDUALS (RX) VS X]"2: Y-RESIDUALS (RY) VS YCALC]"3: Y-RESIDUALS VS X-RESIDUALS]"#4: WEIGHTED X-RESIDUALS (FNX) VS X]"'5: WEIGHTED Y-RESIDUALS (FNY) VS YCALC]"06: WEIGHTED Y-RESIDUALS VS WEIGHTED X-RESIDUALS]"%7: WEIGHTED COMPOSITE RESIDUALS VS X]"8: RETURN TO MENU 5 j @r" ) N T T) NԜ T T ]FN&" ) N T T )N+T7 )DNT7a)1I1I)1015)01+1I1+1I)))N3333O1STKTCT;61,'TT T))N3333J TT4))N3333 TTG]".TABLES OF X-WEIGHTS:SX; X-DATA:x; X-ADJ:X; AND] " COMPOSITE WEIGHTED RESIDUALS:DI]"6TABLES OF Y-WEIGHTS:SY; Y-DATA:y; Y-ADJ:Y; AND SQUARED] " WEIGHTED RESIDUALS: SI = DI*DI]"2TABLES OF X-RESIDUALS:RX; Y-RESIDUALS:RY; WEIGHTED] ". X-RESIDUALS:FNX; AND WEIGHTED Y-RESIDUALS:FNY]"I]"SX] "x]"X] "DIH]"I]"SY] "y]"Y] "SI]"I]"RX] "RY]"FNX] "FNY-]H# J"  ) NNNNN T T3d ) 111111 )N˜)T!TT T T  )DNΜT7L )1I1I )1015 )01+1I1+1I )B)CT;61,'TT TB]" *-*-*-*-*-*-*-*-*   )^T71^] "%** ITERATIVE RE-WEIGHTING SOLUTION **  )N$T7[]" ~T T )h1&|T )1&VTlh]"*CORR. COEFF. OF X-Y UNWEIGHTED RESIDUALS =2]"*CORR. COEFF. OF X-Y *WEIGHTED* RESIDUALS =2{ )B )B  TTX N@:0yE>& .>|=-C6??ư>?}:z8@#B ;MbP?eAKH9QrT hTTT TT rUjWVFNRPHRPFQPXP^P+&&&&,+&&4,2,NN&&++&&FVF H F H 9V}9Fs6F VIw^6F Gp Vڥ^6F V&w^6F p Vw^6F (V~w^F H p. &(&0&0$1& &I$II I[NT T ڟԟVT VȜVT ſVŽźVŸŴVŲTŮVŬũVŧţVšŞVŜřVŗŔVŒŏVōŊVňŅV}xfaTYGBT:(#TT TߟT֜ΟĸĴįĩğT ĚT ĖĒďČĉĆăĀ}zwtqnkheb_\WSPMJ@V<4V0+&V#VVT T T )>Nu >Nu&&&& +&r&p82(:b d b d >d :|>b w1b ݟ*6b  V*^뵎 2(:nPvP66~P u=t$ u=u u=u u=uS&& +&r&p&&,+&&+00&6&6&6r&6p&6 &6&6&6P"P&p& ru>0}  >0wH+&&PP66P u=tx u=u u=u u=u u=u u=u u=u u=u u= u u= u u= u4PPUΥV ơΙVƕ VV Vޟ۟ӝʟžŵT ŰT ŬͣVş͗VœŐT ŋT ņT ŁT |T wT rT mT hda\T VT QNKFBT >T 6T /T *VV ՝VɝŸĿT ĸĵT IJħğT ěĕđVďČVĊćVāzqlga\XVVSVQNVLIVGDVB=9T 5T -T &T T T  A]1PPҸPP밸PPPI&6& &I&601 P?P"-*P2NN+&&:PBP66JP u=t0 u=t u=t, u=t: u=t@ u=u&&+o+,`+&&&&,+:&>t&>t&&FVF H F H 9V|s9FrlF ݇I@؛~𐛊fw1P P+PPP6F ݄IĜPWTVNEV 8V/&V"VVVVT T T T ǜŽźŷŴűŮũťŢŘŒŋňŁ{sm0,$V V VڝVΝȝVļĶ̮VĪĤĠ̛V#ĘT ĔVĒďVčT ĈT Ą~V#{T wVuqnVgVc[V WOVKF>V:5-V)$V Ve&_F H qPPPP&6 &6PTP\P66dP u=t- u=t u=u u=u u=uznP&>u&>t+vP&6&6~P&>t-&>t"P PP&6r&6pP&p& ruP&p&r&(F H -"$>(|>&wF H +$"6$6"P&p&r F H -> |>wF H +O_\UOH?;7.*'#T T V ٟ՟џȟğŽŴT ůT ū͢VŞŖT őT ō̈́Vŀ}T xT sjV f^VZOVKBT 7T 2)V%"T T V T T VĪĦ̞VĚėē̊V Ć~VzrVnkT fT aXV TLVH=V94.%V!V V TVVg 66P&& b d ->|>wb d +66P&& b d -  >|>wb d + 6 6 P&& FVb d b d qb d 9V}9Fsb  H&G&_6b 捜 V*w^똸PP&p&rFVF H 4b ㍇ -RP6d 6b bP F H pcF H 9V|W9FrPF t&G&Wb d >0|>0vb ㍇ -RP6d 6b X눃>0|>0v lPt"CPHVD>5/,(T T  T ՟џ˟ƟͽVŹŶŲɩT ŦT ŞŖŐņT ŁT }tVp]XT RNVIVDVB?V=T 91 T T V۝؟ԟП͟ƟĹİĬĨğěĘĔċT ĆT ĂyV urnjg`ZSJFB952.%T T V \ P>0~| >0r&&&+& &p&rVFF&>(*,>u>u+&&+&&>u9u |PP&& b d ->|>wb d +66PPPP6I6IP&&>IuG>Iu@PPPP+PPP>Iu >IufI It&&FVF H F H qF H 9V}y9Fso6F VI^6F p GVڥ^6F 捜V& JVQTK84T+T '!T ܟԟΟT ſT ŻŲŮťŞ͒V ͅV!ŁyV!umV!iaV!]WPJD;V740'V #VV VݟԟП̟ßĿļĸįT ĪT Ħ̝Vę̑VčĈā}zwtold]YVWTVROVMJVFV>T 9T 4T /T *T %T ! V!w^6F p V^6F 捜(V~^+7>Iu >Iu>It>It+$@"@PPP00*55:PPPP++P++++I++IPP"PIX+IX",*P&&F܉VF H F I*ݟ6I^ZŪŧVťŢVŠŝVśŘVŒʼnŃyT tT pgVc`V^[VYVVTQVOLVJGVEBV;V7/V +#VVVV V VVVVVVVVVV ՝VɝVĽĺVĸĵVijİVĮīVĩĦVĤġVğĜV̕V đ̉Vą}Vyvsh^TMFC91T+"T/@؛~䐛ft~F GŽhڥ&PFPF ㍿&6F ݄&",~қnқ:`Iݜ:F I:ݟ:F H prF H 9V|f9Fr_I@؛~䐛fu6F &G+~қnқ*`Iݜ6++؛~䐛fu'II؛~䐛fuf h +h f >Iu>It>Iu>Iu +55m>Iu>IuFZŻŶŰŨšŚœŐňŁzsolc]SVOVLVJGVEBV:V6V3V1.V,)V$V"VVVV VVVVVVTVVVҜV͟VļĸIJĭĩVħĤVĢğVĝĚVĔĎċVĉĆVĄāV|VxVtVpVmjVgVebV^RJDV;.T+"T VVVV_zPNP6&6&VP&>u[&>uS&>uK&> uC&p& ru8>Iu>It>Iu#>Iu>&u>&t55+55>Itr>Ith&>u.&>u&&> |&>v>&t>&tx`PhPPIIPIItPI6IɛPP&&FƉVF H F H qF H 9V}9FsF p IF ㍟(I6F p GVڥ^6F V&I^g>It*>It &&F‰VF H F H qF ̜QYPKC=3T .T * TT ٜӟ˟TſŷűTŭťŏŋł}uoeT `T \SVOGVC?V=:V7V4V2/V-*V#VVVV V VVVV VޟԟT T ĻT ijT ĮĤĚĖēĊĄ}vohaZTT OT HT @T 8T 0T +"VV VH 9V}9Fs6F p V^6F 捜(V~^6F 捜(&Gp &wݜ>H }{ >F wo6F ݄ܤj؛~䐛fw>PP?PP&&FVF H 6F p &Gɛݜ66F 捜(&Gɛݜ:6F r(ܴ6&_6F 捜(ܴ:&_F H plF H 9V|`9FrY>It[>ItQ6F ݄6p &_6F ݄:bVVVTҜVɟVŶŬŜŘŒōʼnVńVłV}zVxT rnViVgdVb_V]T YQLVJGVDVAVt>tPP&&FVF H F H 9V|L9FrE6F 捄(-RP&Pp -RP-RPPF H q>5ub>5u[P*,P&6P:PRPRPp RPIP RPRPP2,P>&u >&u+,@*@+PP&P&P+P&P"@PIP&6P:PIP+PIP",P*P:PIPIPPPP*,P*@P0P@P@PPRPj!P?PP5P+P>*@u >,@uFF@P@PP RPRPj!P?PP5PFPRP+0[E2+V &T  T T ޟןV ̟ǟŸŽŸųŮũŤşŚŕŐŋņŁ|wrmhc^YTOJE@;61,)&V T T T T T T T T ؟ӟΟɟßļĵİ̧VģěTĒTČȂTTwgc^XNT IT E;V96V2.V,)V'$V"VVV V TVu6F ݄!܌0++6F ݄܌.,++6F ݄.,܌0**F ݇&@&&F ݇D&&F ݇b JJF H qF H 9V}9FsF GŽhڥRP&PF P&&6F ݄6&&܄:"&PFPF ㍿&6F &Gܤ&J ZoVmjVheV`V^TZTH@:V1.V,)V'$V"VVV V VVVVTTT؜T˜ßĭĩĠěėVĕĒVĐčVċĈVĄĀV~{VyvVtqVmiVgdVb_V]ZVVRVPMVKHVFCVA>V:6V41V/,V*'V%"VVVV V V6F r&G&wt؛~䐛fr4F J 3&&(&6?b!&?R ?2Z ?B?I>t?>t5PPPP>t>t&&FV::9V}9Fs::9u9u P>u>u+PP>u>u+PP>u>u+PƟſŹŲͩVťŠřœŌ̓Vzsmf]VYTNJF1-*&# VVVڝVΝşĻķVĵIJVİĭVīĨVĦģVġĞVĜęVėĔVĒďVčĊVĈąVăĀV~{VyvVtqVolVjT fVd`]V[XVVSVQNVE?T ;5(V$V!VVVT V T P PpG*,&6*,&:*,+6z&:&+> (P4PFPRP+&&$@"@6ZP u=u& u=u u=u u=uO u=u u=u' u=uv t=u &~қ(nқ&&|؛~䐛fs::&?&|؛~䐛fs::&::&&&@؛~o\VV؜VӟV̟ƟźŲŬŦŜVŘVŕVœŐVŎŋVŅyoVkVhVfcVa^V\YVVVRVPMVIVEVBT >V<ҟΟV ĿVĽĺVķVĵIJįĬ̢VĞ̖VĒ̊VĆ~VzwVurVnVkhVfcV_V\YVWTVPVMJVHEVC@V>;V96V41V/,V*'V%"V VVMwft(6&+&@؛~䐛ft(6&+(6I&&FV::>u>uF ݇D+>u>uF ݇&I&&@؛~䐛fu +&~қnқ&++F H pCF H 9V|79Fr0>t=>t3F ݇&@+&_&(ݟ+YqnVliVgdVb_VYRVPMVKHVFCV=3) VVVVVVVVVޟVV՟VVVVVľĻVĹĶVĴıVįĬVĪħVġĚēČĉVćĄV~wpieb^UQNVLIVGDVB?V=:V85V30V(V$V!VVVV V VV$6!q9V}9Fs&&FVF H K>G|>w>p &G㍟(&G먛؛~䐛fs @@؛~䐛fs &@(@>Iu>It/؛~䐛fuf h +h f >@}PrPxPP@@ }؃ u=tT u=u u=u u=u u=u u=u u=u u=uu@ @u"RPRPRP@PbPjP+@@(@&@CX0-V+(V%V# V VVT T T T TTٟ՟jf^VZRVNFVB=730'!VVV V VVVVۜV֟VϟɟĿVĻVĸVĶijVıĮVĪVĨĥVģĠVěVęTēďVčĊVĈąVĀV~Tzrkb]XOI?T :T 62V0-V*V VʠAL# &@ (@u&P땡@ @u&@ (@uRPx@ @uf h twRP6,U&@ (@t\f h tP&P1ڡ@ @t:&@ (@t.f h t"6,иRPj 뤸PPP6I6IP>Iu>It>Iu>Iu6&6&P>Iu>Iu P>Iu>IuPPPP&&FVF H F H 9V|Z9FrS6F b P-RPp -RPp -RP6H 6F PF H qP&&FVF H F H 9V|P9FrI6F 5H;+'"T T VVݟTTT T ɼTŹTųŪŚŖőŋŁT |T xoVkcV_WVSKVGA:1V-' V VܟVϝVÝľĹȴTıTĬĤĠĘĔČĈă~vrjf_ZTWTOKFBIu>It>Iu>Iu6&6&P>Iu>Iu P>Iu>IuPPP>5ue>5u^2,4,P*,P&6P:PRPRPp RPIP RPRPP2,P+4,ĜNKDV ?:5T 0T -T (#T T T T T  VٝV͝VŻŴͫVŧšŚ͑VōŊņxqjaV]ZVMGA;2V.)%VT T T T T T T T VܝVΝ˟ǟľT ȵT IJT ĬĥĜĆĂytlf\T WT SJVF@;2V.+'TT?'2,>Ium>IufP019$@|9"@v6$@6"@P019+|9+v6+6+6+6+&&FVF H F H 9V|y9Frr6F (&Gp &wݜ>H |7>F v.6F ݄ܤj؛~䐛fw,F H p{ PP019+|9+v 6+6+" 6+6+,P@P@PP RPRPj!P?PP5P0PRP6P>PPPXPlKVޝVҝVƝ;Vź͵V ŰT ūŦšŜŗŒɍT ŊT Ņŀ{sVolhb_[UMIE=V91V-*V(%V# VVVVVVVVVVٟПVʟVVTĹVķTijīěėĒČĂT }T yrokda]WOKG?V;84.&"V )&&FV::9V}9Fs::9u9u pP>u>u+PxP>u>u+PP>u>u+PPpG?JJ&JJ*,D&`؛~䐛fv'&`؛~䐛fv00+000 1t)r&@&@r&&r&&&&PPLĝͼVŸͰVŬťŢřŖőVŏŌVŊŇVŅłV|ypmhda^UOEVAV>V<9V74V,V(V%V# VVVV V VVVVVT VVVߟVڟVϟʟVĽĸıīĤ̛VėĒċą~uVqle_XOVKF@<8#+P>0un>0ugRPRPRP&P @PPP6,P1PRP&P @PPP+00P@PPPP+PPPFF^_]eV V!Ľ̵V!ı̩V!ĥ̝V!ęē̋VćĄāwVsnV id_ZT WT RMEVA