home *** CD-ROM | disk | FTP | other *** search
open in: 
MacOS 8.1
|
Win98
|
DOS
browse contents |
view JSON data
|
view as text
This file was processed as: LBR Archive
(archive/lbr).
| Confidence | Program | Detection | Match Type | Support
|
|---|
| 100%
| dexvert
| LBR Archive (archive/lbr)
| magic
| Supported |
| 1%
| dexvert
| Micro Focus COBOL library (other/microFocusCobolLibrary)
| ext
| Unsupported |
| 1%
| dexvert
| Micro Focus COBOL Library (other/microFocusCOBOLLibrary)
| ext
| Unsupported |
| 100%
| file
| LBR archive data
| default
| |
| 99%
| file
| data
| default
| |
| 87%
| TrID
| LU library
| default
| |
| 12%
| TrID
| C64 Hires bitmap
| default (weak)
| |
| 100%
| lsar
| LBR
| default
|
|
hex view+--------+-------------------------+-------------------------+--------+--------+
|00000000| 00 20 20 20 20 20 20 20 | 20 20 20 20 00 00 03 00 |. | ....|
|00000010| f9 ed 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00000020| 00 41 4e 59 46 49 54 20 | 20 42 41 53 03 00 4a 00 |.ANYFIT | BAS..J.|
|00000030| f3 fb 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00000040| 00 4c 49 4e 46 49 54 20 | 20 42 41 53 4d 00 1d 00 |.LINFIT | BASM...|
|00000050| 96 d1 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00000060| 00 4d 4f 53 54 46 49 54 | 20 42 41 53 6a 00 44 00 |.MOSTFIT| BASj.D.|
|00000070| 0a 3e 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |.>......|........|
|00000080| 00 50 4f 4c 59 46 49 54 | 20 42 41 53 ae 00 2c 00 |.POLYFIT| BAS..,.|
|00000090| 11 2c 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |.,......|........|
|000000a0| 00 43 55 52 56 46 49 54 | 20 44 4f 43 da 00 4a 00 |.CURVFIT| DOC..J.|
|000000b0| b6 0d 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|000000c0| 00 4d 54 45 53 54 31 20 | 20 44 41 54 24 01 01 00 |.MTEST1 | DAT$...|
|000000d0| a2 18 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|000000e0| 00 43 55 52 54 45 53 54 | 20 44 41 54 25 01 01 00 |.CURTEST| DAT%...|
|000000f0| cd 22 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |."......|........|
|00000100| 00 4c 49 4e 54 45 53 54 | 32 44 41 54 26 01 01 00 |.LINTEST|2DAT&...|
|00000110| 64 a4 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |d.......|........|
|00000120| 00 4d 54 45 53 54 32 20 | 20 44 41 54 27 01 01 00 |.MTEST2 | DAT'...|
|00000130| 1f 26 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |.&......|........|
|00000140| 00 50 4f 4c 59 54 45 53 | 54 44 41 54 28 01 01 00 |.POLYTES|TDAT(...|
|00000150| ab e1 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00000160| 00 4c 49 4e 54 45 53 54 | 31 44 41 54 29 01 01 00 |.LINTEST|1DAT)...|
|00000170| d8 55 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |.U......|........|
|00000180| 31 30 30 30 20 27 20 20 | 20 20 20 20 20 20 20 20 |1000 ' | |
|00000190| 20 20 20 20 20 20 2a 2a | 2a 20 20 20 50 52 4f 47 | **|* PROG|
|000001a0| 52 41 4d 20 41 4e 59 46 | 49 54 20 20 20 2a 2a 2a |RAM ANYF|IT ***|
|000001b0| 0d 0a 31 30 31 30 20 27 | 0d 0a 31 30 32 30 20 27 |..1010 '|..1020 '|
|000001c0| 20 20 41 20 44 52 49 56 | 45 52 20 50 52 4f 47 52 | A DRIV|ER PROGR|
|000001d0| 41 4d 20 46 4f 52 20 54 | 48 45 20 43 55 52 46 49 |AM FOR T|HE CURFI|
|000001e0| 54 20 53 55 42 52 4f 55 | 54 49 4e 45 53 20 54 48 |T SUBROU|TINES TH|
|000001f0| 41 54 20 57 49 4c 4c 20 | 46 49 4e 44 0d 0a 31 30 |AT WILL |FIND..10|
|00000200| 33 30 20 27 20 54 48 45 | 20 50 41 52 41 4d 45 54 |30 ' THE| PARAMET|
|00000210| 45 52 53 20 49 4e 20 41 | 20 47 49 56 45 4e 20 46 |ERS IN A| GIVEN F|
|00000220| 55 4e 43 54 49 4f 4e 20 | 54 48 41 54 20 47 49 56 |UNCTION |THAT GIV|
|00000230| 45 20 54 48 45 20 42 45 | 53 54 20 46 49 54 0d 0a |E THE BE|ST FIT..|
|00000240| 31 30 34 30 20 27 20 28 | 49 4e 20 41 20 4c 45 41 |1040 ' (|IN A LEA|
|00000250| 53 54 20 53 51 55 41 52 | 45 53 20 53 45 4e 53 45 |ST SQUAR|ES SENSE|
|00000260| 29 20 54 4f 20 41 20 53 | 45 54 20 4f 46 20 49 4e |) TO A S|ET OF IN|
|00000270| 50 55 54 20 44 41 54 41 | 2e 20 20 58 20 26 20 59 |PUT DATA|. X & Y|
|00000280| 20 41 52 45 0d 0a 31 30 | 35 30 20 27 20 54 48 45 | ARE..10|50 ' THE|
|00000290| 20 49 4e 44 45 50 45 4e | 44 45 4e 54 20 41 4e 44 | INDEPEN|DENT AND|
|000002a0| 20 44 45 50 45 4e 44 45 | 4e 54 20 56 41 52 49 41 | DEPENDE|NT VARIA|
|000002b0| 42 4c 45 53 20 52 45 53 | 50 45 43 54 49 56 45 4c |BLES RES|PECTIVEL|
|000002c0| 59 2e 20 20 54 48 45 0d | 0a 31 30 36 30 20 27 20 |Y. THE.|.1060 ' |
|000002d0| 46 55 4e 43 54 49 4f 4e | 20 54 48 41 54 20 44 45 |FUNCTION| THAT DE|
|000002e0| 54 45 52 4d 49 4e 45 53 | 20 59 20 47 49 56 45 4e |TERMINES| Y GIVEN|
|000002f0| 20 41 20 56 41 4c 55 45 | 20 4f 46 20 58 20 4d 55 | A VALUE| OF X MU|
|00000300| 53 54 20 42 45 20 53 50 | 45 43 49 46 49 45 44 0d |ST BE SP|ECIFIED.|
|00000310| 0a 31 30 37 30 20 27 20 | 49 4e 20 54 48 45 20 53 |.1070 ' |IN THE S|
|00000320| 55 42 52 4f 55 54 49 4e | 45 20 46 55 4e 43 54 4e |UBROUTIN|E FUNCTN|
|00000330| 20 28 53 45 45 20 45 4e | 44 20 4f 46 20 53 4f 55 | (SEE EN|D OF SOU|
|00000340| 52 43 45 20 43 4f 44 45 | 29 2e 20 20 53 45 45 20 |RCE CODE|). SEE |
|00000350| 43 4f 4d 4d 45 4e 54 53 | 0d 0a 31 30 38 30 20 27 |COMMENTS|..1080 '|
|00000360| 20 49 4e 20 43 55 52 46 | 49 54 20 53 55 42 52 4f | IN CURF|IT SUBRO|
|00000370| 55 54 49 4e 45 20 46 4f | 52 20 45 58 50 4c 41 4e |UTINE FO|R EXPLAN|
|00000380| 41 54 49 4f 4e 20 4f 46 | 20 56 41 52 49 41 42 4c |ATION OF| VARIABL|
|00000390| 45 53 2e 0d 0a 31 30 39 | 30 20 27 0d 0a 31 31 30 |ES...109|0 '..110|
|000003a0| 30 20 27 20 20 20 50 72 | 6f 67 72 61 6d 6d 65 72 |0 ' Pr|ogrammer|
|000003b0| 3a 20 20 47 2e 20 4d 2e | 20 52 65 73 63 68 2c 20 |: G. M.| Resch, |
|000003c0| 20 46 65 62 2e 20 37 39 | 0d 0a 31 31 31 30 20 27 | Feb. 79|..1110 '|
|000003d0| 0d 0a 31 31 32 30 20 27 | 09 44 49 4d 45 4e 53 49 |..1120 '|.DIMENSI|
|000003e0| 4f 4e 20 41 4c 4c 20 41 | 52 52 41 59 53 20 48 45 |ON ALL A|RRAYS HE|
|000003f0| 52 45 0d 0a 31 31 33 30 | 20 27 20 54 48 49 53 20 |RE..1130| ' THIS |
|00000400| 56 45 52 53 49 4f 4e 20 | 49 53 20 53 45 54 20 55 |VERSION |IS SET U|
|00000410| 50 20 54 4f 20 48 41 4e | 44 4c 45 20 55 50 20 54 |P TO HAN|DLE UP T|
|00000420| 4f 20 31 30 30 20 44 41 | 54 41 20 50 4f 49 4e 54 |O 100 DA|TA POINT|
|00000430| 53 20 41 4e 44 0d 0a 31 | 31 34 30 20 27 20 41 53 |S AND..1|140 ' AS|
|00000440| 20 4d 41 4e 59 20 41 53 | 20 31 30 20 50 41 52 41 | MANY AS| 10 PARA|
|00000450| 4d 45 54 45 52 53 2e 0d | 0a 31 31 35 30 20 44 49 |METERS..|.1150 DI|
|00000460| 4d 20 58 28 31 30 30 29 | 2c 20 59 28 31 30 30 29 |M X(100)|, Y(100)|
|00000470| 2c 20 53 49 47 59 28 31 | 30 30 29 2c 20 59 46 49 |, SIGY(1|00), YFI|
|00000480| 54 28 31 30 30 29 2c 20 | 57 47 48 54 28 31 30 30 |T(100), |WGHT(100|
|00000490| 29 0d 0a 31 31 36 30 20 | 44 49 4d 20 41 28 31 30 |)..1160 |DIM A(10|
|000004a0| 29 2c 20 44 45 4c 41 28 | 31 30 29 2c 20 49 4b 28 |), DELA(|10), IK(|
|000004b0| 31 30 29 2c 20 4a 4b 28 | 31 30 29 2c 20 42 45 54 |10), JK(|10), BET|
|000004c0| 41 28 31 30 29 2c 20 42 | 28 31 30 29 2c 20 44 45 |A(10), B|(10), DE|
|000004d0| 52 49 56 28 31 30 29 0d | 0a 31 31 37 30 20 44 49 |RIV(10).|.1170 DI|
|000004e0| 4d 20 41 4c 50 41 23 28 | 31 30 2c 31 30 29 2c 20 |M ALPA#(|10,10), |
|000004f0| 41 52 59 23 28 31 30 2c | 31 30 29 0d 0a 31 31 38 |ARY#(10,|10)..118|
|00000500| 30 20 27 0d 0a 31 31 39 | 30 20 4d 41 58 4e 3d 20 |0 '..119|0 MAXN= |
|00000510| 31 30 30 20 20 3a 20 4d | 41 58 50 3d 31 30 09 27 |100 : M|AXP=10.'|
|00000520| 53 45 54 20 4d 41 58 49 | 4d 55 4d 20 23 20 4f 46 |SET MAXI|MUM # OF|
|00000530| 20 50 4f 49 4e 54 20 41 | 4e 44 20 50 41 52 4d 53 | POINT A|ND PARMS|
|00000540| 0d 0a 31 32 30 30 20 27 | 0d 0a 31 32 31 30 20 27 |..1200 '|..1210 '|
|00000550| 09 20 4f 50 45 4e 20 26 | 20 52 45 41 44 20 54 48 |. OPEN &| READ TH|
|00000560| 45 20 49 4e 50 55 54 20 | 44 41 54 41 20 53 45 54 |E INPUT |DATA SET|
|00000570| 0d 0a 31 32 31 35 20 50 | 52 49 4e 54 22 20 50 72 |..1215 P|RINT" Pr|
|00000580| 6f 67 72 61 6d 20 41 4e | 59 46 49 54 2c 20 20 47 |ogram AN|YFIT, G|
|00000590| 2e 20 4d 2e 20 52 65 73 | 63 68 2c 20 46 65 62 2e |. M. Res|ch, Feb.|
|000005a0| 20 37 39 22 0d 0a 31 32 | 32 30 20 49 4e 50 55 54 | 79"..12|20 INPUT|
|000005b0| 20 22 4e 61 6d 65 20 6f | 66 20 49 6e 70 75 74 20 | "Name o|f Input |
|000005c0| 64 61 74 61 20 73 65 74 | 3d 22 3b 4e 24 0d 0a 31 |data set|=";N$..1|
|000005d0| 32 33 30 20 4f 50 45 4e | 22 49 22 2c 32 2c 4e 24 |230 OPEN|"I",2,N$|
|000005e0| 0d 0a 31 32 34 30 20 46 | 4f 52 20 49 3d 31 20 54 |..1240 F|OR I=1 T|
|000005f0| 4f 20 4d 41 58 4e 0d 0a | 31 32 35 30 20 20 20 49 |O MAXN..|1250 I|
|00000600| 46 20 45 4f 46 28 32 29 | 20 20 47 4f 54 4f 20 31 |F EOF(2)| GOTO 1|
|00000610| 32 38 30 0d 0a 31 32 36 | 30 20 20 20 49 4e 50 55 |280..126|0 INPU|
|00000620| 54 23 32 2c 58 28 49 29 | 2c 59 28 49 29 2c 53 49 |T#2,X(I)|,Y(I),SI|
|00000630| 47 59 28 49 29 0d 0a 31 | 32 37 30 20 4e 45 58 54 |GY(I)..1|270 NEXT|
|00000640| 20 49 0d 0a 31 32 38 30 | 20 4e 50 54 53 3d 49 2d | I..1280| NPTS=I-|
|00000650| 31 20 20 3a 20 43 4c 4f | 53 45 23 32 0d 0a 31 32 |1 : CLO|SE#2..12|
|00000660| 39 30 20 50 52 49 4e 54 | 22 54 48 45 52 45 20 41 |90 PRINT|"THERE A|
|00000670| 52 45 20 22 3b 4e 50 54 | 53 3b 22 20 44 41 54 41 |RE ";NPT|S;" DATA|
|00000680| 20 50 4f 49 4e 54 53 22 | 0d 0a 31 33 30 30 20 27 | POINTS"|..1300 '|
|00000690| 0d 0a 31 33 31 30 20 27 | 09 20 49 4e 50 55 54 20 |..1310 '|. INPUT |
|000006a0| 54 48 45 20 41 50 52 49 | 4f 52 49 20 56 41 4c 55 |THE APRI|ORI VALU|
|000006b0| 45 53 20 46 4f 52 20 54 | 48 45 20 50 41 52 41 4d |ES FOR T|HE PARAM|
|000006c0| 45 54 45 52 53 0d 0a 31 | 33 32 30 20 50 52 49 4e |ETERS..1|320 PRIN|
|000006d0| 54 22 20 20 4d 4f 44 45 | 3d 2b 31 2c 20 57 47 48 |T" MODE|=+1, WGH|
|000006e0| 54 28 49 29 3d 31 2f 53 | 49 47 59 28 49 29 5e 32 |T(I)=1/S|IGY(I)^2|
|000006f0| 22 0d 0a 31 33 33 30 20 | 50 52 49 4e 54 22 20 20 |"..1330 |PRINT" |
|00000700| 20 20 20 20 3d 20 30 2c | 20 57 47 48 54 28 49 29 | = 0,| WGHT(I)|
|00000710| 3d 31 22 0d 0a 31 33 34 | 30 20 50 52 49 4e 54 22 |=1"..134|0 PRINT"|
|00000720| 20 20 20 20 20 20 3d 2d | 31 2c 20 57 47 48 54 28 | =-|1, WGHT(|
|00000730| 49 29 3d 31 2f 59 28 49 | 29 22 0d 0a 31 33 35 30 |I)=1/Y(I|)"..1350|
|00000740| 20 50 52 49 4e 54 22 20 | 20 20 20 20 20 3d 2d 32 | PRINT" | =-2|
|00000750| 2c 20 57 47 48 54 28 49 | 29 3d 53 49 47 59 28 49 |, WGHT(I|)=SIGY(I|
|00000760| 29 22 0d 0a 31 33 36 30 | 20 49 4e 50 55 54 22 53 |)"..1360| INPUT"S|
|00000770| 50 45 43 49 46 59 20 54 | 48 45 20 57 45 49 47 48 |PECIFY T|HE WEIGH|
|00000780| 54 49 4e 47 20 4d 4f 44 | 45 20 28 2b 31 2c 30 2c |TING MOD|E (+1,0,|
|00000790| 2d 31 2c 20 4f 52 20 2d | 32 29 20 2d 22 3b 4d 4f |-1, OR -|2) -";MO|
|000007a0| 44 45 0d 0a 31 33 37 30 | 20 49 46 20 28 4d 4f 44 |DE..1370| IF (MOD|
|000007b0| 45 3c 3d 2b 31 20 41 4e | 44 20 4d 4f 44 45 3e 3d |E<=+1 AN|D MODE>=|
|000007c0| 2d 32 29 20 20 47 4f 54 | 4f 20 31 33 39 30 0d 0a |-2) GOT|O 1390..|
|000007d0| 31 33 38 30 20 50 52 49 | 4e 54 22 45 52 52 4f 52 |1380 PRI|NT"ERROR|
|000007e0| 20 49 4e 20 4d 4f 44 45 | 20 53 50 45 43 49 46 49 | IN MODE| SPECIFI|
|000007f0| 43 41 54 49 4f 4e 22 20 | 20 3a 20 20 47 4f 54 4f |CATION" | : GOTO|
|00000800| 20 31 33 32 30 0d 0a 31 | 33 39 30 20 50 52 49 4e | 1320..1|390 PRIN|
|00000810| 54 0d 0a 31 34 30 30 20 | 49 4e 50 55 54 22 48 4f |T..1400 |INPUT"HO|
|00000820| 57 20 4d 41 4e 59 20 50 | 41 52 41 4d 45 54 45 52 |W MANY P|ARAMETER|
|00000830| 53 3f 20 4e 54 52 4d 3d | 22 3b 4e 54 52 4d 0d 0a |S? NTRM=|";NTRM..|
|00000840| 31 34 31 30 20 49 46 20 | 28 4e 54 45 52 4d 3c 3d |1410 IF |(NTERM<=|
|00000850| 31 30 20 41 4e 44 20 4e | 54 52 4d 3e 30 29 20 20 |10 AND N|TRM>0) |
|00000860| 47 4f 54 4f 20 31 34 33 | 30 0d 0a 31 34 32 30 20 |GOTO 143|0..1420 |
|00000870| 50 52 49 4e 54 22 49 4c | 4c 45 47 41 4c 20 4e 55 |PRINT"IL|LEGAL NU|
|00000880| 4d 42 45 52 20 4f 46 20 | 50 41 52 41 4d 45 54 45 |MBER OF |PARAMETE|
|00000890| 52 53 22 20 20 3a 20 20 | 47 4f 54 4f 20 31 33 39 |RS" : |GOTO 139|
|000008a0| 30 0d 0a 31 34 33 30 20 | 50 52 49 4e 54 0d 0a 31 |0..1430 |PRINT..1|
|000008b0| 34 34 30 20 46 4f 52 20 | 49 3d 31 20 54 4f 20 4e |440 FOR |I=1 TO N|
|000008c0| 54 52 4d 0d 0a 31 34 35 | 30 20 20 20 50 52 49 4e |TRM..145|0 PRIN|
|000008d0| 54 22 53 50 45 43 49 46 | 59 20 54 48 45 20 53 54 |T"SPECIF|Y THE ST|
|000008e0| 41 52 54 49 4e 47 20 56 | 41 4c 55 45 20 41 4e 44 |ARTING V|ALUE AND|
|000008f0| 20 44 45 4c 54 41 20 46 | 4f 52 20 4e 54 52 4d 3d | DELTA F|OR NTRM=|
|00000900| 22 3b 49 0d 0a 31 34 36 | 30 20 20 20 49 4e 50 55 |";I..146|0 INPU|
|00000910| 54 22 20 20 20 20 4e 54 | 52 4d 3d 22 3b 41 28 49 |T" NT|RM=";A(I|
|00000920| 29 0d 0a 31 34 37 30 20 | 20 20 49 4e 50 55 54 22 |)..1470 | INPUT"|
|00000930| 20 20 20 20 44 45 4c 54 | 41 3d 22 3b 44 45 4c 41 | DELT|A=";DELA|
|00000940| 28 49 29 0d 0a 31 34 38 | 30 20 4e 45 58 54 20 49 |(I)..148|0 NEXT I|
|00000950| 0d 0a 31 34 39 30 20 27 | 0d 0a 31 35 30 30 20 27 |..1490 '|..1500 '|
|00000960| 09 20 53 45 54 55 50 20 | 54 48 45 20 43 41 4c 4c |. SETUP |THE CALL|
|00000970| 20 54 4f 20 43 55 52 46 | 49 54 0d 0a 31 35 31 30 | TO CURF|IT..1510|
|00000980| 20 46 49 52 53 54 24 3d | 22 54 52 55 45 22 09 27 | FIRST$=|"TRUE".'|
|00000990| 54 48 49 53 20 49 53 20 | 54 48 45 20 46 49 52 53 |THIS IS |THE FIRS|
|000009a0| 54 20 43 41 4c 4c 0d 0a | 31 35 32 30 20 4c 53 54 |T CALL..|1520 LST|
|000009b0| 43 48 49 3d 31 45 2b 30 | 38 09 27 53 45 54 20 54 |CHI=1E+0|8.'SET T|
|000009c0| 4f 20 52 45 44 49 43 55 | 4c 4f 55 53 4c 59 20 48 |O REDICU|LOUSLY H|
|000009d0| 49 47 48 20 56 41 4c 55 | 45 20 46 4f 52 20 53 54 |IGH VALU|E FOR ST|
|000009e0| 41 52 54 0d 0a 31 35 33 | 30 20 46 4c 4d 44 3d 2e |ART..153|0 FLMD=.|
|000009f0| 30 31 09 09 27 47 4f 4f | 44 20 56 41 4c 55 45 20 |01..'GOO|D VALUE |
|00000a00| 46 4f 52 20 53 54 41 52 | 54 0d 0a 31 35 34 30 20 |FOR STAR|T..1540 |
|00000a10| 4c 4f 4f 50 3d 30 09 09 | 27 49 4e 49 54 49 41 4c |LOOP=0..|'INITIAL|
|00000a20| 49 5a 45 20 4c 4f 4f 50 | 20 43 4f 55 4e 54 0d 0a |IZE LOOP| COUNT..|
|00000a30| 31 35 35 30 20 4c 4f 4f | 50 3d 4c 4f 4f 50 2b 31 |1550 LOO|P=LOOP+1|
|00000a40| 0d 0a 31 35 36 30 20 20 | 47 4f 53 55 42 20 32 30 |..1560 |GOSUB 20|
|00000a50| 34 30 20 20 09 27 43 41 | 4c 4c 20 43 55 52 46 49 |40 .'CA|LL CURFI|
|00000a60| 54 0d 0a 31 35 37 30 20 | 20 20 50 52 49 4e 54 22 |T..1570 | PRINT"|
|00000a70| 4c 4f 4f 50 20 23 22 3b | 4c 4f 4f 50 2c 22 43 48 |LOOP #";|LOOP,"CH|
|00000a80| 49 53 51 3d 22 3b 43 48 | 49 53 51 2c 22 46 4c 41 |ISQ=";CH|ISQ,"FLA|
|00000a90| 4d 44 41 3d 22 2c 46 4c | 4d 44 0d 0a 31 35 38 30 |MDA=",FL|MD..1580|
|00000aa0| 20 20 20 49 46 20 28 4c | 53 54 43 48 49 2d 43 48 | IF (L|STCHI-CH|
|00000ab0| 49 53 51 29 3c 2e 30 30 | 31 20 20 47 4f 54 4f 20 |ISQ)<.00|1 GOTO |
|00000ac0| 31 36 34 30 20 20 27 45 | 58 49 54 20 49 46 20 4e |1640 'E|XIT IF N|
|00000ad0| 4f 20 43 48 41 4e 47 45 | 20 49 4e 20 43 48 49 53 |O CHANGE| IN CHIS|
|00000ae0| 51 0d 0a 31 35 39 30 20 | 20 20 4c 53 54 43 48 49 |Q..1590 | LSTCHI|
|00000af0| 3d 43 48 49 53 51 0d 0a | 31 36 30 30 20 20 20 49 |=CHISQ..|1600 I|
|00000b00| 46 20 46 4c 4d 44 3c 2e | 30 30 30 30 30 31 20 20 |F FLMD<.|000001 |
|00000b10| 54 48 45 4e 20 46 4c 4d | 44 3d 2e 30 30 30 30 31 |THEN FLM|D=.00001|
|00000b20| 0d 0a 31 36 31 30 20 20 | 49 46 20 46 4c 4d 44 3c |..1610 |IF FLMD<|
|00000b30| 31 45 2b 30 36 20 20 47 | 4f 54 4f 20 31 35 35 30 |1E+06 G|OTO 1550|
|00000b40| 09 27 4b 45 45 50 20 4c | 4f 4f 50 49 4e 47 20 49 |.'KEEP L|OOPING I|
|00000b50| 46 20 46 4c 4d 44 20 49 | 53 20 4e 4f 54 20 54 4f |F FLMD I|S NOT TO|
|00000b60| 4f 20 4c 41 52 47 45 0d | 0a 31 36 32 30 20 27 0d |O LARGE.|.1620 '.|
|00000b70| 0a 31 36 33 30 20 27 09 | 57 45 20 46 4f 55 4e 44 |.1630 '.|WE FOUND|
|00000b80| 20 41 4e 20 41 43 43 45 | 50 54 41 42 4c 45 20 46 | AN ACCE|PTABLE F|
|00000b90| 49 54 2c 20 53 55 4d 4d | 41 52 49 5a 45 20 54 48 |IT, SUMM|ARIZE TH|
|00000ba0| 45 20 52 45 53 55 4c 54 | 53 0d 0a 31 36 34 30 20 |E RESULT|S..1640 |
|00000bb0| 50 52 49 4e 54 0d 0a 31 | 36 35 30 20 46 4f 52 20 |PRINT..1|650 FOR |
|00000bc0| 4a 3d 31 20 54 4f 20 4e | 54 52 4d 0d 0a 31 36 36 |J=1 TO N|TRM..166|
|00000bd0| 30 20 20 50 52 49 4e 54 | 22 41 28 22 3b 4a 3b 22 |0 PRINT|"A(";J;"|
|00000be0| 29 3d 22 3b 41 28 4a 29 | 2c 22 53 49 47 4d 41 41 |)=";A(J)|,"SIGMAA|
|00000bf0| 28 22 3b 4a 3b 22 29 3d | 22 3b 53 49 47 41 28 4a |(";J;")=|";SIGA(J|
|00000c00| 29 0d 0a 31 36 37 30 20 | 4e 45 58 54 20 4a 0d 0a |)..1670 |NEXT J..|
|00000c10| 31 36 38 30 20 50 52 49 | 4e 54 22 43 48 49 20 53 |1680 PRI|NT"CHI S|
|00000c20| 51 55 41 52 45 3d 22 3b | 43 48 49 53 51 0d 0a 31 |QUARE=";|CHISQ..1|
|00000c30| 36 39 30 20 50 52 49 4e | 54 0d 0a 31 37 30 30 20 |690 PRIN|T..1700 |
|00000c40| 50 52 49 4e 54 22 58 28 | 49 29 22 2c 22 59 28 49 |PRINT"X(|I)","Y(I|
|00000c50| 29 22 2c 22 59 46 49 54 | 22 2c 22 59 2d 59 46 49 |)","YFIT|","Y-YFI|
|00000c60| 54 22 0d 0a 31 37 31 30 | 20 46 4f 52 20 49 3d 31 |T"..1710| FOR I=1|
|00000c70| 20 54 4f 20 4e 50 54 53 | 0d 0a 31 37 32 30 20 20 | TO NPTS|..1720 |
|00000c80| 50 52 49 4e 54 20 58 28 | 49 29 2c 59 28 49 29 2c |PRINT X(|I),Y(I),|
|00000c90| 59 46 49 54 28 49 29 2c | 28 59 28 49 29 2d 59 46 |YFIT(I),|(Y(I)-YF|
|00000ca0| 49 54 28 49 29 29 0d 0a | 31 37 33 30 20 4e 45 58 |IT(I))..|1730 NEX|
|00000cb0| 54 20 49 0d 0a 31 37 34 | 30 20 45 4e 44 0d 0a 31 |T I..174|0 END..1|
|00000cc0| 37 35 30 20 52 45 4d 2a | 2a 2a 2a 2a 2a 2a 2a 2a |750 REM*|********|
|00000cd0| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00000ce0| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00000cf0| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00000d00| 2a 2a 2a 2a 2a 2a 0d 0a | 31 37 36 30 20 27 20 20 |******..|1760 ' |
|00000d10| 53 55 42 52 4f 55 54 49 | 4e 45 20 43 55 52 46 49 |SUBROUTI|NE CURFI|
|00000d20| 54 28 58 2c 20 59 2c 20 | 53 49 47 59 2c 20 4e 50 |T(X, Y, |SIGY, NP|
|00000d30| 54 53 2c 20 4e 54 52 4d | 2c 20 4d 4f 44 45 2c 20 |TS, NTRM|, MODE, |
|00000d40| 41 2c 20 44 45 4c 41 2c | 0d 0a 31 37 37 30 20 27 |A, DELA,|..1770 '|
|00000d50| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00000d60| 20 20 20 20 53 49 47 41 | 2c 20 46 4c 4d 44 2c 20 | SIGA|, FLMD, |
|00000d70| 59 46 49 54 2c 20 43 48 | 49 53 51 2c 20 24 46 49 |YFIT, CH|ISQ, $FI|
|00000d80| 52 53 54 29 0d 0a 31 37 | 38 30 20 27 0d 0a 31 37 |RST)..17|80 '..17|
|00000d90| 39 30 20 27 09 41 20 52 | 4f 55 54 49 4e 45 20 54 |90 '.A R|OUTINE T|
|00000da0| 48 41 54 20 57 49 4c 4c | 20 46 49 4e 44 20 54 48 |HAT WILL| FIND TH|
|00000db0| 45 20 4c 45 41 53 54 2d | 53 51 55 41 52 45 53 20 |E LEAST-|SQUARES |
|00000dc0| 46 49 54 20 54 4f 20 41 | 20 4e 4f 4e 2d 4c 49 4e |FIT TO A| NON-LIN|
|00000dd0| 45 41 52 0d 0a 31 38 30 | 30 20 27 20 46 55 4e 43 |EAR..180|0 ' FUNC|
|00000de0| 54 49 4f 4e 20 57 49 54 | 48 20 41 20 4c 49 4e 45 |TION WIT|H A LINE|
|00000df0| 41 52 49 5a 41 54 49 4f | 4e 20 4f 46 20 54 48 45 |ARIZATIO|N OF THE|
|00000e00| 20 46 49 54 54 49 4e 47 | 20 46 55 4e 43 54 49 4f | FITTING| FUNCTIO|
|00000e10| 4e 2e 0d 0a 31 38 33 30 | 20 27 0d 0a 31 38 34 30 |N...1830| '..1840|
|00000e20| 20 27 20 50 41 52 41 4d | 45 54 45 52 53 3a 0d 0a | ' PARAM|ETERS:..|
|00000e30| 31 38 35 30 20 27 20 20 | 58 3d 41 52 52 41 59 20 |1850 ' |X=ARRAY |
|00000e40| 4f 46 20 49 4e 44 45 50 | 45 4e 44 2e 20 56 41 52 |OF INDEP|END. VAR|
|00000e50| 49 41 42 4c 45 2c 20 20 | 59 3d 41 52 52 41 59 20 |IABLE, |Y=ARRAY |
|00000e60| 4f 46 20 44 45 50 45 4e | 44 45 4e 54 20 56 41 52 |OF DEPEN|DENT VAR|
|00000e70| 49 41 42 4c 45 0d 0a 31 | 38 36 30 20 27 20 20 53 |IABLE..1|860 ' S|
|00000e80| 49 47 59 3d 53 54 4e 44 | 20 44 45 56 2e 20 4f 46 |IGY=STND| DEV. OF|
|00000e90| 20 59 2c 20 20 4e 50 54 | 53 3d 23 20 4f 46 20 50 | Y, NPT|S=# OF P|
|00000ea0| 41 49 52 53 20 4f 46 20 | 44 41 54 41 20 50 4f 49 |AIRS OF |DATA POI|
|00000eb0| 4e 54 53 0d 0a 31 38 37 | 30 20 27 20 20 4e 54 52 |NTS..187|0 ' NTR|
|00000ec0| 4d 3d 23 20 4f 46 20 50 | 41 52 4d 53 2c 20 20 41 |M=# OF P|ARMS, A|
|00000ed0| 3d 41 52 52 41 59 20 4f | 46 20 50 41 52 4d 53 0d |=ARRAY O|F PARMS.|
|00000ee0| 0a 31 38 38 30 20 27 20 | 20 4d 4f 44 45 3d 2b 31 |.1880 ' | MODE=+1|
|00000ef0| 2c 20 57 47 48 54 28 49 | 29 3d 31 2f 53 49 47 59 |, WGHT(I|)=1/SIGY|
|00000f00| 28 49 29 5e 32 0d 0a 31 | 38 39 30 20 27 09 20 20 |(I)^2..1|890 '. |
|00000f10| 20 20 3d 20 30 2c 20 57 | 47 48 54 28 49 29 3d 31 | = 0, W|GHT(I)=1|
|00000f20| 0d 0a 31 39 30 30 20 27 | 20 20 20 20 20 20 3d 2d |..1900 '| =-|
|00000f30| 31 2c 20 57 47 48 54 28 | 49 29 3d 31 2f 59 28 49 |1, WGHT(|I)=1/Y(I|
|00000f40| 29 0d 0a 31 39 31 30 20 | 27 09 20 20 20 20 3d 2d |)..1910 |'. =-|
|00000f50| 32 2c 20 57 47 48 54 28 | 49 29 3d 53 49 47 59 28 |2, WGHT(|I)=SIGY(|
|00000f60| 49 29 0d 0a 31 39 32 30 | 20 27 20 20 44 45 4c 41 |I)..1920| ' DELA|
|00000f70| 3d 41 52 52 41 59 20 4f | 46 20 49 4e 43 52 45 4d |=ARRAY O|F INCREM|
|00000f80| 45 4e 54 53 20 4f 46 20 | 41 2c 20 20 53 49 47 41 |ENTS OF |A, SIGA|
|00000f90| 41 3d 53 54 4e 44 20 44 | 45 56 20 4f 46 20 41 0d |A=STND D|EV OF A.|
|00000fa0| 0a 31 39 33 30 20 27 20 | 20 46 4c 4d 44 3d 50 52 |.1930 ' | FLMD=PR|
|00000fb0| 4f 50 4f 52 54 49 4f 4e | 20 4f 46 20 47 52 41 44 |OPORTION| OF GRAD|
|00000fc0| 49 45 4e 54 20 53 45 41 | 52 43 48 20 49 4e 43 4c |IENT SEA|RCH INCL|
|00000fd0| 55 44 45 44 0d 0a 31 39 | 34 30 20 27 20 20 43 48 |UDED..19|40 ' CH|
|00000fe0| 49 53 51 3d 52 45 44 55 | 43 45 44 20 43 48 49 20 |ISQ=REDU|CED CHI |
|00000ff0| 53 51 55 41 52 45 0d 0a | 31 39 35 30 20 27 20 20 |SQUARE..|1950 ' |
|00001000| 46 49 52 53 54 24 3d 20 | 41 20 46 4c 41 47 20 54 |FIRST$= |A FLAG T|
|00001010| 48 41 54 20 49 53 20 54 | 52 55 45 20 57 48 45 4e |HAT IS T|RUE WHEN|
|00001020| 20 43 55 52 46 49 54 20 | 49 53 20 46 49 52 53 54 | CURFIT |IS FIRST|
|00001030| 20 43 41 4c 4c 45 44 0d | 0a 31 39 36 30 20 27 0d | CALLED.|.1960 '.|
|00001040| 0a 31 39 37 30 20 27 20 | 53 55 42 52 4f 55 54 49 |.1970 ' |SUBROUTI|
|00001050| 4e 45 53 20 43 41 4c 4c | 45 44 3b 0d 0a 31 39 38 |NES CALL|ED;..198|
|00001060| 30 20 27 20 20 20 46 55 | 4e 43 54 4e 28 58 2c 49 |0 ' FU|NCTN(X,I|
|00001070| 2c 41 29 2c 20 20 46 43 | 48 49 53 51 28 59 2c 53 |,A), FC|HISQ(Y,S|
|00001080| 49 47 59 2c 4e 50 54 53 | 2c 4e 46 52 45 2c 4d 4f |IGY,NPTS|,NFRE,MO|
|00001090| 44 45 2c 59 46 49 54 29 | 0d 0a 31 39 39 30 20 27 |DE,YFIT)|..1990 '|
|000010a0| 20 20 20 46 44 45 52 49 | 56 28 58 2c 49 2c 41 2c | FDERI|V(X,I,A,|
|000010b0| 44 45 4c 41 2c 4e 54 52 | 4d 2c 44 45 52 49 56 29 |DELA,NTR|M,DERIV)|
|000010c0| 2c 20 20 4d 41 54 49 4e | 56 28 41 52 59 23 2c 4e |, MATIN|V(ARY#,N|
|000010d0| 54 52 4d 2c 44 45 54 29 | 0d 0a 32 30 30 30 20 27 |TRM,DET)|..2000 '|
|000010e0| 0d 0a 32 30 31 30 20 27 | 2d 2d 2d 20 53 54 41 52 |..2010 '|--- STAR|
|000010f0| 54 20 4f 46 20 50 52 4f | 47 52 41 4d 20 43 4f 44 |T OF PRO|GRAM COD|
|00001100| 45 20 2d 2d 2d 0d 0a 32 | 30 32 30 20 27 0d 0a 32 |E ---..2|020 '..2|
|00001110| 30 33 30 20 27 20 20 52 | 4f 55 54 49 4e 45 20 49 |030 ' R|OUTINE I|
|00001120| 53 20 46 49 52 53 54 20 | 43 41 4c 4c 45 44 0d 0a |S FIRST |CALLED..|
|00001130| 32 30 34 30 20 49 46 20 | 46 49 52 53 54 24 3d 22 |2040 IF |FIRST$="|
|00001140| 46 41 4c 53 45 22 20 20 | 47 4f 54 4f 20 32 31 35 |FALSE" |GOTO 215|
|00001150| 30 20 20 27 53 4b 49 50 | 20 49 46 20 4e 4f 54 20 |0 'SKIP| IF NOT |
|00001160| 46 49 52 53 54 20 43 41 | 4c 4c 0d 0a 32 30 35 30 |FIRST CA|LL..2050|
|00001170| 20 46 49 52 53 54 24 3d | 22 46 41 4c 53 45 22 0d | FIRST$=|"FALSE".|
|00001180| 0a 32 30 36 30 20 4e 46 | 52 45 3d 4e 50 54 53 2d |.2060 NF|RE=NPTS-|
|00001190| 4e 54 52 4d 20 20 3a 20 | 43 48 49 53 51 3d 30 20 |NTRM : |CHISQ=0 |
|000011a0| 20 3a 20 49 46 20 4e 46 | 52 45 3c 31 20 54 48 45 | : IF NF|RE<1 THE|
|000011b0| 4e 20 52 45 54 55 52 4e | 0d 0a 32 30 37 30 20 27 |N RETURN|..2070 '|
|000011c0| 09 20 45 56 41 4c 55 41 | 54 45 20 57 45 49 47 48 |. EVALUA|TE WEIGH|
|000011d0| 54 53 0d 0a 32 30 38 30 | 20 46 4f 52 20 49 3d 31 |TS..2080| FOR I=1|
|000011e0| 20 54 4f 20 4e 50 54 53 | 0d 0a 32 30 39 30 20 20 | TO NPTS|..2090 |
|000011f0| 49 46 20 4d 4f 44 45 3d | 31 20 20 54 48 45 4e 20 |IF MODE=|1 THEN |
|00001200| 57 47 48 54 28 49 29 3d | 31 2f 28 53 49 47 59 28 |WGHT(I)=|1/(SIGY(|
|00001210| 49 29 2a 53 49 47 59 28 | 49 29 29 20 20 3a 20 47 |I)*SIGY(|I)) : G|
|00001220| 4f 54 4f 20 32 31 33 30 | 0d 0a 32 31 30 30 20 20 |OTO 2130|..2100 |
|00001230| 49 46 20 4d 4f 44 45 3d | 30 20 20 54 48 45 4e 20 |IF MODE=|0 THEN |
|00001240| 57 47 48 54 28 49 29 3d | 31 20 20 3a 20 47 4f 54 |WGHT(I)=|1 : GOT|
|00001250| 4f 20 32 31 33 30 0d 0a | 32 31 31 30 20 20 49 46 |O 2130..|2110 IF|
|00001260| 20 4d 4f 44 45 3d 2d 31 | 20 54 48 45 4e 20 57 47 | MODE=-1| THEN WG|
|00001270| 48 54 28 49 29 3d 31 2f | 41 42 53 28 59 28 49 29 |HT(I)=1/|ABS(Y(I)|
|00001280| 29 20 20 3a 20 47 4f 54 | 4f 20 32 31 33 30 0d 0a |) : GOT|O 2130..|
|00001290| 32 31 32 30 20 20 49 46 | 20 4d 4f 44 45 3d 2d 32 |2120 IF| MODE=-2|
|000012a0| 20 54 48 45 4e 20 57 47 | 48 54 28 49 29 3d 53 49 | THEN WG|HT(I)=SI|
|000012b0| 47 59 28 49 29 0d 0a 32 | 31 33 30 20 4e 45 58 54 |GY(I)..2|130 NEXT|
|000012c0| 20 49 0d 0a 32 31 34 30 | 20 27 09 20 45 56 41 4c | I..2140| '. EVAL|
|000012d0| 20 41 4c 50 48 41 20 26 | 20 42 45 54 41 20 4d 41 | ALPHA &| BETA MA|
|000012e0| 54 52 49 43 45 53 0d 0a | 32 31 35 30 20 46 4f 52 |TRICES..|2150 FOR|
|000012f0| 20 4a 3d 31 20 54 4f 20 | 4e 54 52 4d 0d 0a 32 31 | J=1 TO |NTRM..21|
|00001300| 36 30 20 20 20 42 45 54 | 41 28 4a 29 3d 30 0d 0a |60 BET|A(J)=0..|
|00001310| 32 31 37 30 20 46 4f 52 | 20 4b 3d 31 20 54 4f 20 |2170 FOR| K=1 TO |
|00001320| 4a 0d 0a 32 31 38 30 20 | 20 20 41 4c 50 41 23 28 |J..2180 | ALPA#(|
|00001330| 4a 2c 4b 29 3d 30 0d 0a | 32 31 39 30 20 4e 45 58 |J,K)=0..|2190 NEX|
|00001340| 54 20 4b 2c 4a 0d 0a 32 | 32 30 30 20 27 0d 0a 32 |T K,J..2|200 '..2|
|00001350| 32 31 30 20 46 4f 52 20 | 49 3d 31 20 54 4f 20 4e |210 FOR |I=1 TO N|
|00001360| 50 54 53 0d 0a 32 32 32 | 30 20 20 47 4f 53 55 42 |PTS..222|0 GOSUB|
|00001370| 20 32 38 35 30 20 20 3a | 52 45 4d 2a 20 43 41 4c | 2850 :|REM* CAL|
|00001380| 4c 20 46 44 45 52 49 56 | 28 58 2c 49 2c 41 2c 44 |L FDERIV|(X,I,A,D|
|00001390| 45 4c 41 2c 4e 54 52 4d | 2c 44 45 52 49 56 29 0d |ELA,NTRM|,DERIV).|
|000013a0| 0a 32 32 33 30 20 20 47 | 4f 53 55 42 20 33 37 38 |.2230 G|OSUB 378|
|000013b0| 30 20 20 3a 52 45 4d 2a | 20 43 41 4c 4c 20 46 55 |0 :REM*| CALL FU|
|000013c0| 4e 43 54 4e 28 58 2c 49 | 2c 41 29 0d 0a 32 32 34 |NCTN(X,I|,A)..224|
|000013d0| 30 20 46 4f 52 20 4a 3d | 31 20 54 4f 20 4e 54 52 |0 FOR J=|1 TO NTR|
|000013e0| 4d 0d 0a 32 32 35 30 20 | 20 20 42 45 54 41 28 4a |M..2250 | BETA(J|
|000013f0| 29 3d 42 45 54 41 28 4a | 29 20 2b 57 47 48 54 28 |)=BETA(J|) +WGHT(|
|00001400| 49 29 2a 28 59 28 49 29 | 2d 46 55 4e 43 54 4e 29 |I)*(Y(I)|-FUNCTN)|
|00001410| 2a 44 45 52 49 56 28 4a | 29 0d 0a 32 32 36 30 20 |*DERIV(J|)..2260 |
|00001420| 46 4f 52 20 4b 3d 31 20 | 54 4f 20 4a 0d 0a 32 32 |FOR K=1 |TO J..22|
|00001430| 37 30 20 20 20 41 4c 50 | 41 23 28 4a 2c 4b 29 3d |70 ALP|A#(J,K)=|
|00001440| 41 4c 50 41 23 28 4a 2c | 4b 29 20 2b 57 47 48 54 |ALPA#(J,|K) +WGHT|
|00001450| 28 49 29 2a 44 45 52 49 | 56 28 4a 29 2a 44 45 52 |(I)*DERI|V(J)*DER|
|00001460| 49 56 28 4b 29 0d 0a 32 | 32 38 30 20 4e 45 58 54 |IV(K)..2|280 NEXT|
|00001470| 20 4b 2c 4a 0d 0a 32 32 | 39 30 20 4e 45 58 54 20 | K,J..22|90 NEXT |
|00001480| 49 0d 0a 32 33 30 30 20 | 27 0d 0a 32 33 31 30 20 |I..2300 |'..2310 |
|00001490| 46 4f 52 20 4a 3d 31 20 | 54 4f 20 4e 54 52 4d 0d |FOR J=1 |TO NTRM.|
|000014a0| 0a 32 33 32 30 20 46 4f | 52 20 4b 3d 31 20 54 4f |.2320 FO|R K=1 TO|
|000014b0| 20 4a 0d 0a 32 33 33 30 | 20 20 20 41 4c 50 41 23 | J..2330| ALPA#|
|000014c0| 28 4b 2c 4a 29 3d 41 4c | 50 41 23 28 4a 2c 4b 29 |(K,J)=AL|PA#(J,K)|
|000014d0| 0d 0a 32 33 34 30 20 4e | 45 58 54 20 4b 2c 4a 0d |..2340 N|EXT K,J.|
|000014e0| 0a 32 33 35 30 20 27 09 | 20 45 56 41 4c 20 43 48 |.2350 '.| EVAL CH|
|000014f0| 49 20 53 51 55 41 52 45 | 20 41 54 20 53 54 41 52 |I SQUARE| AT STAR|
|00001500| 54 20 50 4f 49 4e 54 0d | 0a 32 33 36 30 20 46 4f |T POINT.|.2360 FO|
|00001510| 52 20 49 3d 31 20 54 4f | 20 4e 50 54 53 0d 0a 32 |R I=1 TO| NPTS..2|
|00001520| 33 37 30 20 20 20 47 4f | 53 55 42 20 33 37 38 30 |370 GO|SUB 3780|
|00001530| 20 20 3a 52 45 4d 2a 20 | 43 41 4c 4c 20 46 55 4e | :REM* |CALL FUN|
|00001540| 43 54 4e 28 58 2c 49 2c | 41 29 0d 0a 32 33 38 30 |CTN(X,I,|A)..2380|
|00001550| 20 20 20 20 59 46 49 54 | 28 49 29 3d 46 55 4e 43 | YFIT|(I)=FUNC|
|00001560| 54 4e 0d 0a 32 33 39 30 | 20 20 20 20 47 4f 53 55 |TN..2390| GOSU|
|00001570| 42 20 32 39 39 30 20 20 | 3a 52 45 4d 2a 20 43 41 |B 2990 |:REM* CA|
|00001580| 4c 4c 20 46 43 48 49 53 | 51 28 59 2c 53 49 47 59 |LL FCHIS|Q(Y,SIGY|
|00001590| 2c 4e 50 54 53 2c 4e 46 | 52 45 2c 4d 4f 44 45 2c |,NPTS,NF|RE,MODE,|
|000015a0| 59 46 49 54 29 0d 0a 32 | 34 30 30 20 20 20 43 48 |YFIT)..2|400 CH|
|000015b0| 49 3d 46 43 48 49 53 51 | 0d 0a 32 34 31 30 20 4e |I=FCHISQ|..2410 N|
|000015c0| 45 58 54 20 49 0d 0a 32 | 34 32 30 20 27 09 20 49 |EXT I..2|420 '. I|
|000015d0| 4e 56 45 52 54 20 4d 4f | 44 49 46 49 45 44 20 43 |NVERT MO|DIFIED C|
|000015e0| 55 52 56 41 54 55 52 45 | 20 4d 41 54 52 49 58 20 |URVATURE| MATRIX |
|000015f0| 54 4f 20 46 49 4e 44 20 | 4e 45 57 20 50 41 52 4d |TO FIND |NEW PARM|
|00001600| 53 0d 0a 32 34 33 30 20 | 46 4f 52 20 4a 3d 31 20 |S..2430 |FOR J=1 |
|00001610| 54 4f 20 4e 54 52 4d 0d | 0a 32 34 34 30 20 46 4f |TO NTRM.|.2440 FO|
|00001620| 52 20 4b 3d 31 20 54 4f | 20 4e 54 52 4d 0d 0a 32 |R K=1 TO| NTRM..2|
|00001630| 34 35 30 20 20 20 41 52 | 59 23 28 4a 2c 4b 29 3d |450 AR|Y#(J,K)=|
|00001640| 41 4c 50 41 23 28 4a 2c | 4b 29 2f 53 51 52 28 41 |ALPA#(J,|K)/SQR(A|
|00001650| 4c 50 41 23 28 4a 2c 4a | 29 2a 41 4c 50 41 23 28 |LPA#(J,J|)*ALPA#(|
|00001660| 4b 2c 4b 29 29 0d 0a 32 | 34 36 30 20 4e 45 58 54 |K,K))..2|460 NEXT|
|00001670| 20 4b 0d 0a 32 34 37 30 | 20 20 20 41 52 59 23 28 | K..2470| ARY#(|
|00001680| 4a 2c 4a 29 3d 31 2b 46 | 4c 4d 44 0d 0a 32 34 38 |J,J)=1+F|LMD..248|
|00001690| 30 20 4e 45 58 54 20 4a | 0d 0a 32 34 39 30 20 27 |0 NEXT J|..2490 '|
|000016a0| 0d 0a 32 35 30 30 20 4e | 4f 52 3d 4e 54 52 4d 20 |..2500 N|OR=NTRM |
|000016b0| 20 3a 20 47 4f 53 55 42 | 20 33 31 34 30 20 20 3a | : GOSUB| 3140 :|
|000016c0| 52 45 4d 2a 20 43 41 4c | 4c 20 4d 41 54 49 4e 56 |REM* CAL|L MATINV|
|000016d0| 28 41 52 59 23 2c 4e 54 | 52 4d 2c 44 45 54 29 0d |(ARY#,NT|RM,DET).|
|000016e0| 0a 32 35 31 30 20 46 4f | 52 20 4a 3d 31 20 54 4f |.2510 FO|R J=1 TO|
|000016f0| 20 4e 54 52 4d 0d 0a 32 | 35 32 30 20 20 20 42 28 | NTRM..2|520 B(|
|00001700| 4a 29 3d 41 28 4a 29 0d | 0a 32 35 33 30 20 46 4f |J)=A(J).|.2530 FO|
|00001710| 52 20 4b 3d 31 20 54 4f | 20 4e 54 52 4d 0d 0a 32 |R K=1 TO| NTRM..2|
|00001720| 35 34 30 20 20 20 42 28 | 4a 29 3d 42 28 4a 29 20 |540 B(|J)=B(J) |
|00001730| 2b 42 45 54 41 28 4b 29 | 2a 41 52 59 23 28 4a 2c |+BETA(K)|*ARY#(J,|
|00001740| 4b 29 2f 53 51 52 28 41 | 4c 50 41 23 28 4a 2c 4a |K)/SQR(A|LPA#(J,J|
|00001750| 29 2a 41 4c 50 41 23 28 | 4b 2c 4b 29 29 0d 0a 32 |)*ALPA#(|K,K))..2|
|00001760| 35 35 30 20 4e 45 58 54 | 20 4b 2c 4a 0d 0a 32 35 |550 NEXT| K,J..25|
|00001770| 36 30 20 27 09 20 49 46 | 20 43 48 49 20 53 51 55 |60 '. IF| CHI SQU|
|00001780| 41 52 45 20 49 4e 43 52 | 45 41 53 45 44 2c 20 49 |ARE INCR|EASED, I|
|00001790| 4e 43 52 45 41 53 45 20 | 46 4c 4d 44 20 41 4e 44 |NCREASE |FLMD AND|
|000017a0| 20 54 52 59 20 41 47 41 | 49 4e 0d 0a 32 35 37 30 | TRY AGA|IN..2570|
|000017b0| 20 46 4f 52 20 4a 3d 31 | 20 54 4f 20 4e 54 52 4d | FOR J=1| TO NTRM|
|000017c0| 0d 0a 32 35 38 30 20 20 | 20 54 45 4d 50 3d 41 28 |..2580 | TEMP=A(|
|000017d0| 4a 29 20 20 3a 20 41 28 | 4a 29 3d 42 28 4a 29 20 |J) : A(|J)=B(J) |
|000017e0| 20 3a 20 42 28 4a 29 3d | 54 45 4d 50 0d 0a 32 35 | : B(J)=|TEMP..25|
|000017f0| 39 30 20 4e 45 58 54 20 | 4a 0d 0a 32 36 30 30 20 |90 NEXT |J..2600 |
|00001800| 27 0d 0a 32 36 31 30 20 | 46 4f 52 20 49 3d 31 20 |'..2610 |FOR I=1 |
|00001810| 54 4f 20 4e 50 54 53 0d | 0a 32 36 32 30 20 20 20 |TO NPTS.|.2620 |
|00001820| 47 4f 53 55 42 20 33 37 | 38 30 20 20 3a 52 45 4d |GOSUB 37|80 :REM|
|00001830| 2a 20 43 41 4c 4c 20 46 | 55 4e 43 54 4e 28 58 2c |* CALL F|UNCTN(X,|
|00001840| 49 2c 41 29 0d 0a 32 36 | 33 30 20 20 20 59 46 49 |I,A)..26|30 YFI|
|00001850| 54 28 49 29 3d 46 55 4e | 43 54 4e 0d 0a 32 36 34 |T(I)=FUN|CTN..264|
|00001860| 30 20 4e 45 58 54 20 49 | 0d 0a 32 36 35 30 20 46 |0 NEXT I|..2650 F|
|00001870| 4f 52 20 4a 3d 31 20 54 | 4f 20 4e 54 52 4d 0d 0a |OR J=1 T|O NTRM..|
|00001880| 32 36 36 30 20 54 45 4d | 50 3d 42 28 4a 29 20 20 |2660 TEM|P=B(J) |
|00001890| 3a 20 42 28 4a 29 3d 41 | 28 4a 29 20 20 3a 20 41 |: B(J)=A|(J) : A|
|000018a0| 28 4a 29 3d 54 45 4d 50 | 0d 0a 32 36 37 30 20 4e |(J)=TEMP|..2670 N|
|000018b0| 45 58 54 20 4a 0d 0a 32 | 36 38 30 20 27 0d 0a 32 |EXT J..2|680 '..2|
|000018c0| 36 39 30 20 47 4f 53 55 | 42 20 32 39 39 30 20 20 |690 GOSU|B 2990 |
|000018d0| 3a 52 45 4d 2a 20 43 41 | 4c 4c 20 46 43 48 49 53 |:REM* CA|LL FCHIS|
|000018e0| 51 28 59 2c 53 49 47 59 | 2c 4e 50 54 53 2c 4e 46 |Q(Y,SIGY|,NPTS,NF|
|000018f0| 52 45 2c 4d 4f 44 45 2c | 59 46 49 54 29 0d 0a 32 |RE,MODE,|YFIT)..2|
|00001900| 37 30 30 20 43 48 49 53 | 51 3d 46 43 48 49 53 51 |700 CHIS|Q=FCHISQ|
|00001910| 0d 0a 32 37 31 30 20 49 | 46 20 28 43 48 49 2d 43 |..2710 I|F (CHI-C|
|00001920| 48 49 53 51 29 3e 30 20 | 47 4f 54 4f 20 32 37 34 |HISQ)>0 |GOTO 274|
|00001930| 30 0d 0a 32 37 32 30 20 | 46 4c 4d 44 3d 31 30 2a |0..2720 |FLMD=10*|
|00001940| 46 4c 4d 44 20 20 3a 20 | 47 4f 54 4f 20 32 34 33 |FLMD : |GOTO 243|
|00001950| 30 0d 0a 32 37 33 30 20 | 27 09 20 45 56 41 4c 20 |0..2730 |'. EVAL |
|00001960| 50 41 52 4d 53 20 26 20 | 55 4e 43 45 52 54 41 49 |PARMS & |UNCERTAI|
|00001970| 4e 54 49 45 53 0d 0a 32 | 37 34 30 20 46 4f 52 20 |NTIES..2|740 FOR |
|00001980| 4a 3d 31 20 54 4f 20 4e | 54 52 4d 0d 0a 32 37 35 |J=1 TO N|TRM..275|
|00001990| 30 20 20 20 41 28 4a 29 | 3d 42 28 4a 29 20 20 3a |0 A(J)|=B(J) :|
|000019a0| 20 53 49 47 41 28 4a 29 | 3d 20 53 51 52 28 41 52 | SIGA(J)|= SQR(AR|
|000019b0| 59 23 28 4a 2c 4a 29 2f | 41 4c 50 41 23 28 4a 2c |Y#(J,J)/|ALPA#(J,|
|000019c0| 4a 29 29 0d 0a 32 37 36 | 30 20 4e 45 58 54 20 4a |J))..276|0 NEXT J|
|000019d0| 0d 0a 32 37 37 30 20 46 | 4c 4d 44 3d 46 4c 4d 44 |..2770 F|LMD=FLMD|
|000019e0| 2f 31 30 0d 0a 32 37 38 | 30 20 52 45 54 55 52 4e |/10..278|0 RETURN|
|000019f0| 0d 0a 32 37 39 30 20 52 | 45 4d 2a 2a 2a 2a 2a 2a |..2790 R|EM******|
|00001a00| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00001a10| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00001a20| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00001a30| 2a 2a 2a 2a 2a 2a 2a 2a | 0d 0a 32 38 30 30 20 27 |********|..2800 '|
|00001a40| 09 20 20 53 55 42 52 4f | 55 54 49 4e 45 20 46 44 |. SUBRO|UTINE FD|
|00001a50| 45 52 49 56 28 58 2c 20 | 49 2c 20 41 2c 20 44 45 |ERIV(X, |I, A, DE|
|00001a60| 4c 41 20 4e 54 52 4d 2c | 20 44 45 52 49 56 29 0d |LA NTRM,| DERIV).|
|00001a70| 0a 32 38 31 30 20 27 0d | 0a 32 38 32 30 20 27 20 |.2810 '.|.2820 ' |
|00001a80| 41 20 52 4f 55 54 49 4e | 45 20 54 4f 20 45 56 41 |A ROUTIN|E TO EVA|
|00001a90| 4c 55 41 54 45 20 54 48 | 45 20 44 45 52 49 56 41 |LUATE TH|E DERIVA|
|00001aa0| 54 49 56 45 20 4f 46 20 | 41 20 46 55 4e 43 54 49 |TIVE OF |A FUNCTI|
|00001ab0| 4f 4e 20 57 49 54 48 0d | 0a 32 38 33 30 20 27 20 |ON WITH.|.2830 ' |
|00001ac0| 52 45 53 50 45 43 54 20 | 54 4f 20 54 48 45 20 50 |RESPECT |TO THE P|
|00001ad0| 41 52 41 4d 45 54 45 52 | 53 20 41 28 4a 29 0d 0a |ARAMETER|S A(J)..|
|00001ae0| 32 38 34 30 20 27 0d 0a | 32 38 35 30 20 46 4f 52 |2840 '..|2850 FOR|
|00001af0| 20 4a 4a 3d 31 20 54 4f | 20 4e 54 52 4d 0d 0a 32 | JJ=1 TO| NTRM..2|
|00001b00| 38 36 30 20 20 41 4a 3d | 41 28 4a 4a 29 20 20 3a |860 AJ=|A(JJ) :|
|00001b10| 20 44 45 4c 3d 44 45 4c | 41 28 4a 4a 29 20 3a 20 | DEL=DEL|A(JJ) : |
|00001b20| 41 28 4a 4a 29 3d 41 4a | 2b 44 45 4c 0d 0a 32 38 |A(JJ)=AJ|+DEL..28|
|00001b30| 37 30 20 20 20 47 4f 53 | 55 42 20 33 37 38 30 20 |70 GOS|UB 3780 |
|00001b40| 20 3a 52 45 4d 2a 20 43 | 41 4c 4c 20 46 55 4e 43 | :REM* C|ALL FUNC|
|00001b50| 54 4e 28 58 2c 49 2c 41 | 29 0d 0a 32 38 38 30 20 |TN(X,I,A|)..2880 |
|00001b60| 20 20 20 59 46 3d 46 55 | 4e 43 54 4e 20 20 3a 20 | YF=FU|NCTN : |
|00001b70| 41 28 4a 4a 29 3d 41 4a | 2d 44 45 4c 0d 0a 32 38 |A(JJ)=AJ|-DEL..28|
|00001b80| 39 30 20 20 20 20 47 4f | 53 55 42 20 33 37 38 30 |90 GO|SUB 3780|
|00001b90| 0d 0a 32 39 30 30 20 20 | 20 44 45 52 49 56 28 4a |..2900 | DERIV(J|
|00001ba0| 4a 29 3d 28 59 46 2d 46 | 55 4e 43 54 4e 29 2f 28 |J)=(YF-F|UNCTN)/(|
|00001bb0| 32 2a 44 45 4c 29 0d 0a | 32 39 31 30 20 20 41 28 |2*DEL)..|2910 A(|
|00001bc0| 4a 4a 29 3d 41 4a 0d 0a | 32 39 32 30 20 4e 45 58 |JJ)=AJ..|2920 NEX|
|00001bd0| 54 20 4a 4a 0d 0a 32 39 | 33 30 20 52 45 54 55 52 |T JJ..29|30 RETUR|
|00001be0| 4e 0d 0a 32 39 34 30 20 | 52 45 4d 2a 2a 2a 2a 2a |N..2940 |REM*****|
|00001bf0| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00001c00| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00001c10| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00001c20| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 0d 0a 32 39 35 30 20 |********|*..2950 |
|00001c30| 27 09 20 20 53 55 42 52 | 4f 55 54 49 4e 45 20 46 |'. SUBR|OUTINE F|
|00001c40| 43 48 49 53 51 28 59 2c | 20 57 47 48 54 2c 20 4e |CHISQ(Y,| WGHT, N|
|00001c50| 50 54 53 2c 20 4e 46 52 | 45 2c 20 59 46 49 54 29 |PTS, NFR|E, YFIT)|
|00001c60| 0d 0a 32 39 36 30 20 27 | 0d 0a 32 39 37 30 20 27 |..2960 '|..2970 '|
|00001c70| 20 20 41 20 52 4f 55 54 | 49 4e 45 20 54 4f 20 43 | A ROUT|INE TO C|
|00001c80| 41 4c 43 55 4c 41 54 45 | 20 54 48 45 20 52 45 44 |ALCULATE| THE RED|
|00001c90| 55 43 45 44 20 43 48 49 | 20 53 51 55 41 52 45 20 |UCED CHI| SQUARE |
|00001ca0| 46 49 54 20 54 4f 20 44 | 41 54 41 0d 0a 32 39 38 |FIT TO D|ATA..298|
|00001cb0| 30 20 27 0d 0a 32 39 39 | 30 20 43 48 53 3d 30 20 |0 '..299|0 CHS=0 |
|00001cc0| 20 3a 20 49 46 20 4e 46 | 52 45 3e 30 20 47 4f 54 | : IF NF|RE>0 GOT|
|00001cd0| 4f 20 33 30 32 30 0d 0a | 33 30 30 30 20 46 43 48 |O 3020..|3000 FCH|
|00001ce0| 49 53 51 3d 30 20 20 3a | 20 52 45 54 55 52 4e 0d |ISQ=0 :| RETURN.|
|00001cf0| 0a 33 30 31 30 20 27 20 | 20 41 43 43 55 4d 55 4c |.3010 ' | ACCUMUL|
|00001d00| 41 54 45 20 43 48 49 20 | 53 51 55 41 52 45 0d 0a |ATE CHI |SQUARE..|
|00001d10| 33 30 32 30 20 46 4f 52 | 20 49 49 3d 31 20 54 4f |3020 FOR| II=1 TO|
|00001d20| 20 4e 50 54 53 0d 0a 33 | 30 33 30 20 20 20 43 48 | NPTS..3|030 CH|
|00001d30| 53 3d 43 48 53 20 2b 57 | 47 48 54 28 49 49 29 2a |S=CHS +W|GHT(II)*|
|00001d40| 28 59 28 49 49 29 2d 59 | 46 49 54 28 49 49 29 29 |(Y(II)-Y|FIT(II))|
|00001d50| 2a 28 59 28 49 49 29 2d | 59 46 49 54 28 49 49 29 |*(Y(II)-|YFIT(II)|
|00001d60| 29 0d 0a 33 30 34 30 20 | 4e 45 58 54 20 49 49 0d |)..3040 |NEXT II.|
|00001d70| 0a 33 30 35 30 20 46 43 | 48 49 53 51 3d 43 48 53 |.3050 FC|HISQ=CHS|
|00001d80| 2f 4e 46 52 45 0d 0a 33 | 30 36 30 20 52 45 54 55 |/NFRE..3|060 RETU|
|00001d90| 52 4e 0d 0a 33 30 37 30 | 20 52 45 4d 2a 2a 2a 2a |RN..3070| REM****|
|00001da0| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00001db0| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00001dc0| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00001dd0| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 0d 0a 33 30 38 30 |********|**..3080|
|00001de0| 20 27 20 20 20 20 20 20 | 20 20 20 20 53 55 42 52 | ' | SUBR|
|00001df0| 4f 55 54 49 4e 45 20 4d | 41 54 49 4e 56 28 41 52 |OUTINE M|ATINV(AR|
|00001e00| 59 2c 4e 4f 52 2c 44 45 | 54 29 0d 0a 33 30 39 30 |Y,NOR,DE|T)..3090|
|00001e10| 20 27 0d 0a 33 31 30 30 | 20 27 20 20 41 20 52 4f | '..3100| ' A RO|
|00001e20| 55 54 49 4e 45 20 54 4f | 20 49 4e 56 45 52 54 20 |UTINE TO| INVERT |
|00001e30| 41 20 53 59 4d 4d 45 54 | 52 49 43 20 4d 41 54 52 |A SYMMET|RIC MATR|
|00001e40| 49 58 20 41 4e 44 20 43 | 41 4c 43 55 4c 41 54 45 |IX AND C|ALCULATE|
|00001e50| 20 49 54 53 0d 0a 33 31 | 31 30 20 27 20 20 44 45 | ITS..31|10 ' DE|
|00001e60| 54 45 52 4d 49 4e 41 4e | 54 2e 20 20 41 52 59 28 |TERMINAN|T. ARY(|
|00001e70| 49 2c 4a 29 20 49 53 20 | 54 48 45 20 49 4e 50 55 |I,J) IS |THE INPU|
|00001e80| 54 20 4d 41 54 52 49 58 | 20 4f 46 20 4f 52 44 45 |T MATRIX| OF ORDE|
|00001e90| 52 3d 4e 4f 52 0d 0a 33 | 31 32 30 20 27 20 20 41 |R=NOR..3|120 ' A|
|00001ea0| 4e 44 20 49 53 20 52 45 | 50 4c 41 43 45 44 20 42 |ND IS RE|PLACED B|
|00001eb0| 59 20 49 54 53 20 49 4e | 56 45 52 53 45 0d 0a 33 |Y ITS IN|VERSE..3|
|00001ec0| 31 33 30 20 27 0d 0a 33 | 31 34 30 20 44 45 54 3d |130 '..3|140 DET=|
|00001ed0| 31 21 0d 0a 33 31 35 30 | 20 27 09 09 53 54 41 52 |1!..3150| '..STAR|
|00001ee0| 54 20 54 48 45 20 4b 25 | 20 4c 4f 4f 50 0d 0a 33 |T THE K%| LOOP..3|
|00001ef0| 31 36 30 20 46 4f 52 20 | 4b 25 3d 31 20 54 4f 20 |160 FOR |K%=1 TO |
|00001f00| 4e 4f 52 0d 0a 33 31 37 | 30 20 20 41 4d 58 23 3d |NOR..317|0 AMX#=|
|00001f10| 30 21 0d 0a 33 31 38 30 | 09 46 4f 52 20 49 25 3d |0!..3180|.FOR I%=|
|00001f20| 4b 25 20 54 4f 20 4e 4f | 52 0d 0a 33 31 39 30 09 |K% TO NO|R..3190.|
|00001f30| 46 4f 52 20 4a 25 3d 4b | 25 20 54 4f 20 4e 4f 52 |FOR J%=K|% TO NOR|
|00001f40| 0d 0a 33 32 30 30 20 20 | 09 20 49 46 20 41 42 53 |..3200 |. IF ABS|
|00001f50| 28 41 4d 58 23 29 3e 41 | 42 53 28 41 52 59 23 28 |(AMX#)>A|BS(ARY#(|
|00001f60| 49 25 2c 4a 25 29 29 20 | 47 4f 54 4f 20 33 32 32 |I%,J%)) |GOTO 322|
|00001f70| 30 0d 0a 33 32 31 30 20 | 20 09 20 41 4d 58 23 3d |0..3210 | . AMX#=|
|00001f80| 41 52 59 23 28 49 25 2c | 4a 25 29 20 20 3a 20 49 |ARY#(I%,|J%) : I|
|00001f90| 4b 28 4b 25 29 3d 49 25 | 20 20 3a 20 4a 4b 28 4b |K(K%)=I%| : JK(K|
|00001fa0| 25 29 3d 4a 25 0d 0a 33 | 32 32 30 09 4e 45 58 54 |%)=J%..3|220.NEXT|
|00001fb0| 20 4a 25 2c 49 25 0d 0a | 33 32 33 30 20 20 49 46 | J%,I%..|3230 IF|
|00001fc0| 20 41 4d 58 23 3c 3e 30 | 21 20 47 4f 54 4f 20 33 | AMX#<>0|! GOTO 3|
|00001fd0| 32 35 30 0d 0a 33 32 34 | 30 20 20 44 45 54 3d 30 |250..324|0 DET=0|
|00001fe0| 21 20 20 3a 20 52 45 54 | 55 52 4e 09 27 52 45 54 |! : RET|URN.'RET|
|00001ff0| 55 52 4e 20 49 46 20 44 | 45 54 45 52 4d 49 4e 41 |URN IF D|ETERMINA|
|00002000| 4e 54 20 49 53 20 5a 45 | 52 4f 0d 0a 33 32 35 30 |NT IS ZE|RO..3250|
|00002010| 20 20 49 25 3d 49 4b 28 | 4b 25 29 0d 0a 33 32 36 | I%=IK(|K%)..326|
|00002020| 30 20 20 49 46 20 49 25 | 3c 4b 25 20 47 4f 54 4f |0 IF I%|<K% GOTO|
|00002030| 20 33 31 38 30 0d 0a 33 | 32 37 30 20 20 49 46 20 | 3180..3|270 IF |
|00002040| 49 25 3d 4b 25 20 47 4f | 54 4f 20 33 33 31 30 0d |I%=K% GO|TO 3310.|
|00002050| 0a 33 32 38 30 09 46 4f | 52 20 4a 25 3d 31 20 54 |.3280.FO|R J%=1 T|
|00002060| 4f 20 4e 4f 52 0d 0a 33 | 32 39 30 20 20 20 20 20 |O NOR..3|290 |
|00002070| 53 41 56 23 3d 41 52 59 | 23 28 4b 25 2c 4a 25 29 |SAV#=ARY|#(K%,J%)|
|00002080| 20 3a 20 41 52 59 23 28 | 4b 25 2c 4a 25 29 3d 41 | : ARY#(|K%,J%)=A|
|00002090| 52 59 23 28 49 25 2c 4a | 25 29 20 3a 20 41 52 59 |RY#(I%,J|%) : ARY|
|000020a0| 23 28 49 25 2c 4a 25 29 | 3d 2d 53 41 56 23 0d 0a |#(I%,J%)|=-SAV#..|
|000020b0| 33 33 30 30 09 4e 45 58 | 54 20 4a 25 0d 0a 33 33 |3300.NEX|T J%..33|
|000020c0| 31 30 20 20 4a 25 3d 4a | 4b 28 4b 25 29 0d 0a 33 |10 J%=J|K(K%)..3|
|000020d0| 33 32 30 20 20 49 46 20 | 4a 25 3c 4b 25 20 47 4f |320 IF |J%<K% GO|
|000020e0| 54 4f 20 33 31 38 30 0d | 0a 33 33 33 30 20 20 49 |TO 3180.|.3330 I|
|000020f0| 46 20 4a 25 3d 4b 25 20 | 47 4f 54 4f 20 33 33 37 |F J%=K% |GOTO 337|
|00002100| 30 0d 0a 33 33 34 30 09 | 46 4f 52 20 49 25 3d 31 |0..3340.|FOR I%=1|
|00002110| 20 54 4f 20 4e 4f 52 0d | 0a 33 33 35 30 09 20 53 | TO NOR.|.3350. S|
|00002120| 41 56 23 3d 41 52 59 23 | 28 49 25 2c 4b 25 29 20 |AV#=ARY#|(I%,K%) |
|00002130| 3a 20 41 52 59 23 28 49 | 25 2c 4b 25 29 3d 41 52 |: ARY#(I|%,K%)=AR|
|00002140| 59 23 28 49 25 2c 4a 25 | 29 20 3a 20 41 52 59 23 |Y#(I%,J%|) : ARY#|
|00002150| 28 49 25 2c 4a 25 29 3d | 2d 53 41 56 23 0d 0a 33 |(I%,J%)=|-SAV#..3|
|00002160| 33 36 30 09 4e 45 58 54 | 20 49 25 0d 0a 33 33 37 |360.NEXT| I%..337|
|00002170| 30 09 46 4f 52 20 49 25 | 3d 31 20 54 4f 20 4e 4f |0.FOR I%|=1 TO NO|
|00002180| 52 0d 0a 33 33 38 30 09 | 20 49 46 20 49 25 3d 4b |R..3380.| IF I%=K|
|00002190| 25 20 47 4f 54 4f 20 33 | 34 30 30 0d 0a 33 33 39 |% GOTO 3|400..339|
|000021a0| 30 09 20 41 52 59 23 28 | 49 25 2c 4b 25 29 3d 2d |0. ARY#(|I%,K%)=-|
|000021b0| 41 52 59 23 28 49 25 2c | 4b 25 29 2f 41 4d 58 23 |ARY#(I%,|K%)/AMX#|
|000021c0| 0d 0a 33 34 30 30 09 4e | 45 58 54 20 49 25 0d 0a |..3400.N|EXT I%..|
|000021d0| 33 34 31 30 09 46 4f 52 | 20 49 25 3d 31 20 54 4f |3410.FOR| I%=1 TO|
|000021e0| 20 4e 4f 52 0d 0a 33 34 | 32 30 09 46 4f 52 20 4a | NOR..34|20.FOR J|
|000021f0| 25 3d 31 20 54 4f 20 4e | 4f 52 0d 0a 33 34 33 30 |%=1 TO N|OR..3430|
|00002200| 09 20 20 49 46 20 49 25 | 3d 4b 25 20 47 4f 54 4f |. IF I%|=K% GOTO|
|00002210| 20 33 34 36 30 0d 0a 33 | 34 34 30 09 20 20 49 46 | 3460..3|440. IF|
|00002220| 20 4a 25 3d 4b 25 20 47 | 4f 54 4f 20 33 34 36 30 | J%=K% G|OTO 3460|
|00002230| 0d 0a 33 34 35 30 09 20 | 20 41 52 59 23 28 49 25 |..3450. | ARY#(I%|
|00002240| 2c 4a 25 29 3d 41 52 59 | 23 28 49 25 2c 4a 25 29 |,J%)=ARY|#(I%,J%)|
|00002250| 2b 41 52 59 23 28 49 25 | 2c 4b 25 29 2a 41 52 59 |+ARY#(I%|,K%)*ARY|
|00002260| 23 28 4b 25 2c 4a 25 29 | 0d 0a 33 34 36 30 09 4e |#(K%,J%)|..3460.N|
|00002270| 45 58 54 20 4a 25 2c 49 | 25 0d 0a 33 34 37 30 09 |EXT J%,I|%..3470.|
|00002280| 46 4f 52 20 4a 25 3d 31 | 20 54 4f 20 4e 4f 52 0d |FOR J%=1| TO NOR.|
|00002290| 0a 33 34 38 30 09 20 20 | 49 46 20 4a 25 3d 4b 25 |.3480. |IF J%=K%|
|000022a0| 20 47 4f 54 4f 20 33 35 | 30 30 0d 0a 33 34 39 30 | GOTO 35|00..3490|
|000022b0| 09 20 20 41 52 59 23 28 | 4b 25 2c 4a 25 29 3d 41 |. ARY#(|K%,J%)=A|
|000022c0| 52 59 23 28 4b 25 2c 4a | 25 29 2f 41 4d 58 23 0d |RY#(K%,J|%)/AMX#.|
|000022d0| 0a 33 35 30 30 09 4e 45 | 58 54 20 4a 25 0d 0a 33 |.3500.NE|XT J%..3|
|000022e0| 35 31 30 20 20 41 52 59 | 23 28 4b 25 2c 4b 25 29 |510 ARY|#(K%,K%)|
|000022f0| 3d 31 21 2f 41 4d 58 23 | 0d 0a 33 35 32 30 20 20 |=1!/AMX#|..3520 |
|00002300| 44 45 54 3d 44 45 54 2a | 41 4d 58 23 0d 0a 33 35 |DET=DET*|AMX#..35|
|00002310| 33 30 20 4e 45 58 54 20 | 4b 25 0d 0a 33 35 34 30 |30 NEXT |K%..3540|
|00002320| 20 27 09 09 54 45 52 4d | 49 4e 41 54 45 20 54 48 | '..TERM|INATE TH|
|00002330| 45 20 4b 25 20 4c 4f 4f | 50 0d 0a 33 35 35 30 20 |E K% LOO|P..3550 |
|00002340| 46 4f 52 20 4c 25 3d 31 | 20 54 4f 20 4e 4f 52 0d |FOR L%=1| TO NOR.|
|00002350| 0a 33 35 36 30 20 20 4b | 25 3d 4e 4f 52 2d 4c 25 |.3560 K|%=NOR-L%|
|00002360| 20 2b 31 20 20 3a 20 4a | 25 3d 49 4b 28 4b 25 29 | +1 : J|%=IK(K%)|
|00002370| 0d 0a 33 35 37 30 20 20 | 49 46 20 4a 25 3c 3d 4b |..3570 |IF J%<=K|
|00002380| 25 20 47 4f 54 4f 20 33 | 36 31 30 0d 0a 33 35 38 |% GOTO 3|610..358|
|00002390| 30 20 20 20 20 46 4f 52 | 20 49 25 3d 31 20 54 4f |0 FOR| I%=1 TO|
|000023a0| 20 4e 4f 52 0d 0a 33 35 | 39 30 20 20 20 20 20 53 | NOR..35|90 S|
|000023b0| 41 56 23 3d 41 52 59 23 | 28 49 25 2c 4b 25 29 20 |AV#=ARY#|(I%,K%) |
|000023c0| 20 3a 20 41 52 59 23 28 | 49 25 2c 4b 25 29 3d 2d | : ARY#(|I%,K%)=-|
|000023d0| 41 52 59 23 28 49 25 2c | 4a 25 29 20 20 3a 20 41 |ARY#(I%,|J%) : A|
|000023e0| 52 59 23 28 49 25 2c 4a | 25 29 3d 53 41 56 23 0d |RY#(I%,J|%)=SAV#.|
|000023f0| 0a 33 36 30 30 20 20 20 | 20 4e 45 58 54 20 49 25 |.3600 | NEXT I%|
|00002400| 0d 0a 33 36 31 30 20 20 | 49 25 3d 4a 4b 28 4b 25 |..3610 |I%=JK(K%|
|00002410| 29 0d 0a 33 36 32 30 20 | 20 49 46 20 49 25 3c 3d |)..3620 | IF I%<=|
|00002420| 4b 25 20 47 4f 54 4f 20 | 33 36 36 30 0d 0a 33 36 |K% GOTO |3660..36|
|00002430| 33 30 20 20 20 20 46 4f | 52 20 4a 25 3d 31 20 54 |30 FO|R J%=1 T|
|00002440| 4f 20 4e 4f 52 0d 0a 33 | 36 34 30 20 20 20 20 20 |O NOR..3|640 |
|00002450| 53 41 56 23 3d 41 52 59 | 23 28 4b 25 2c 4a 25 29 |SAV#=ARY|#(K%,J%)|
|00002460| 20 20 3a 20 41 52 59 23 | 28 4b 25 2c 4a 25 29 3d | : ARY#|(K%,J%)=|
|00002470| 2d 41 52 59 23 28 49 25 | 2c 4a 25 29 20 20 3a 20 |-ARY#(I%|,J%) : |
|00002480| 41 52 59 23 28 49 25 2c | 4a 25 29 3d 53 41 56 23 |ARY#(I%,|J%)=SAV#|
|00002490| 0d 0a 33 36 35 30 20 20 | 20 20 4e 45 58 54 20 4a |..3650 | NEXT J|
|000024a0| 25 0d 0a 33 36 36 30 20 | 4e 45 58 54 20 4c 25 0d |%..3660 |NEXT L%.|
|000024b0| 0a 33 36 37 30 20 52 45 | 54 55 52 4e 0d 0a 33 36 |.3670 RE|TURN..36|
|000024c0| 38 30 20 27 0d 0a 33 36 | 39 30 20 52 45 4d 2a 2a |80 '..36|90 REM**|
|000024d0| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|000024e0| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|000024f0| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00002500| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 0d 0a 33 37 |********|****..37|
|00002510| 30 30 20 27 09 20 20 53 | 55 42 52 4f 55 54 49 4e |00 '. S|UBROUTIN|
|00002520| 45 20 46 55 4e 43 54 4e | 28 58 2c 20 49 2c 20 41 |E FUNCTN|(X, I, A|
|00002530| 29 0d 0a 33 37 31 30 20 | 27 0d 0a 33 37 32 30 20 |)..3710 |'..3720 |
|00002540| 27 20 45 56 41 4c 55 41 | 54 45 20 41 4e 20 41 52 |' EVALUA|TE AN AR|
|00002550| 42 49 54 52 41 52 59 20 | 46 55 4e 43 54 49 4f 4e |BITRARY |FUNCTION|
|00002560| 20 46 4f 52 20 43 55 52 | 46 49 54 0d 0a 33 37 33 | FOR CUR|FIT..373|
|00002570| 30 20 27 20 20 49 4e 20 | 54 48 49 53 20 43 41 53 |0 ' IN |THIS CAS|
|00002580| 45 20 54 48 45 20 46 55 | 4e 43 54 49 4f 4e 20 46 |E THE FU|NCTION F|
|00002590| 52 4f 4d 20 42 59 54 45 | 20 4d 41 47 41 5a 49 4e |ROM BYTE| MAGAZIN|
|000025a0| 45 2c 20 4d 41 59 20 38 | 34 2c 20 50 33 34 31 0d |E, MAY 8|4, P341.|
|000025b0| 0a 33 37 34 30 20 27 0d | 0a 33 37 35 30 20 27 09 |.3740 '.|.3750 '.|
|000025c0| 20 20 58 3d 41 52 52 41 | 59 20 4f 46 20 44 41 54 | X=ARRA|Y OF DAT|
|000025d0| 41 20 50 4f 49 4e 54 53 | 2c 20 20 49 3d 49 4e 44 |A POINTS|, I=IND|
|000025e0| 45 58 20 4f 46 20 44 41 | 54 41 20 50 4f 49 4e 54 |EX OF DA|TA POINT|
|000025f0| 53 0d 0a 33 37 36 30 20 | 27 09 20 20 41 3d 41 52 |S..3760 |'. A=AR|
|00002600| 52 41 59 20 4f 46 20 43 | 4f 45 46 46 49 43 49 45 |RAY OF C|OEFFICIE|
|00002610| 4e 54 53 0d 0a 33 37 37 | 30 20 27 0d 0a 33 37 38 |NTS..377|0 '..378|
|00002620| 30 20 46 55 4e 43 54 4e | 3d 20 41 28 31 29 2a 58 |0 FUNCTN|= A(1)*X|
|00002630| 28 49 29 2f 28 58 28 49 | 29 2b 41 28 32 29 29 0d |(I)/(X(I|)+A(2)).|
|00002640| 0a 33 37 39 30 20 52 45 | 54 55 52 4e 0d 0a 1a 1a |.3790 RE|TURN....|
|00002650| 1a 1a 1a 1a 1a 1a 1a 1a | 1a 1a 1a 1a 1a 1a 1a 1a |........|........|
|00002660| 1a 1a 1a 1a 1a 1a 1a 1a | 1a 1a 1a 1a 1a 1a 1a 1a |........|........|
|00002670| 1a 1a 1a 1a 1a 1a 1a 1a | 1a 1a 1a 1a 1a 1a 1a 1a |........|........|
|00002680| 31 30 30 30 20 27 20 20 | 20 20 20 20 20 20 20 20 |1000 ' | |
|00002690| 20 20 20 20 20 20 2a 2a | 2a 20 20 20 50 52 4f 47 | **|* PROG|
|000026a0| 52 41 4d 20 4c 49 4e 46 | 49 54 20 20 20 2a 2a 2a |RAM LINF|IT ***|
|000026b0| 0d 0a 31 30 31 30 20 27 | 0d 0a 31 30 32 30 20 27 |..1010 '|..1020 '|
|000026c0| 20 20 41 20 44 52 49 56 | 45 52 20 50 52 4f 47 52 | A DRIV|ER PROGR|
|000026d0| 41 4d 20 46 4f 52 20 54 | 48 45 20 4c 49 4e 46 49 |AM FOR T|HE LINFI|
|000026e0| 54 20 53 55 42 52 4f 55 | 54 49 4e 45 20 54 48 41 |T SUBROU|TINE THA|
|000026f0| 54 20 57 49 4c 4c 20 46 | 49 4e 44 0d 0a 31 30 33 |T WILL F|IND..103|
|00002700| 30 20 27 20 54 48 45 20 | 50 41 52 41 4d 45 54 45 |0 ' THE |PARAMETE|
|00002710| 52 53 20 49 4e 20 41 20 | 4c 49 4e 45 41 52 20 46 |RS IN A |LINEAR F|
|00002720| 55 4e 43 54 49 4f 4e 20 | 54 48 41 54 20 47 49 56 |UNCTION |THAT GIV|
|00002730| 45 20 54 48 45 20 42 45 | 53 54 20 46 49 54 0d 0a |E THE BE|ST FIT..|
|00002740| 31 30 34 30 20 27 20 28 | 49 4e 20 41 20 4c 45 41 |1040 ' (|IN A LEA|
|00002750| 53 54 20 53 51 55 41 52 | 45 53 20 53 45 4e 53 45 |ST SQUAR|ES SENSE|
|00002760| 29 20 54 4f 20 41 20 53 | 45 54 20 4f 46 20 49 4e |) TO A S|ET OF IN|
|00002770| 50 55 54 20 44 41 54 41 | 2e 20 20 58 20 26 20 59 |PUT DATA|. X & Y|
|00002780| 20 41 52 45 0d 0a 31 30 | 35 30 20 27 20 54 48 45 | ARE..10|50 ' THE|
|00002790| 20 49 4e 44 45 50 45 4e | 44 45 4e 54 20 41 4e 44 | INDEPEN|DENT AND|
|000027a0| 20 44 45 50 45 4e 44 45 | 4e 54 20 56 41 52 49 41 | DEPENDE|NT VARIA|
|000027b0| 42 4c 45 53 20 52 45 53 | 50 45 43 54 49 56 45 4c |BLES RES|PECTIVEL|
|000027c0| 59 2e 20 20 0d 0a 31 30 | 36 30 20 27 20 53 45 45 |Y. ..10|60 ' SEE|
|000027d0| 20 43 4f 4d 4d 45 4e 54 | 53 20 49 4e 20 4c 49 4e | COMMENT|S IN LIN|
|000027e0| 46 49 54 20 53 55 42 52 | 4f 55 54 49 4e 45 20 46 |FIT SUBR|OUTINE F|
|000027f0| 4f 52 20 45 58 50 4c 41 | 4e 41 54 49 4f 4e 20 4f |OR EXPLA|NATION O|
|00002800| 46 20 56 41 52 49 41 42 | 4c 45 53 2e 0d 0a 31 30 |F VARIAB|LES...10|
|00002810| 37 30 20 27 20 20 50 52 | 4f 47 52 41 4d 4d 45 52 |70 ' PR|OGRAMMER|
|00002820| 3a 20 20 47 65 6f 72 67 | 65 20 4d 2e 20 52 65 73 |: Georg|e M. Res|
|00002830| 63 68 2c 20 20 46 65 62 | 2e 20 31 39 37 39 0d 0a |ch, Feb|. 1979..|
|00002840| 31 30 38 30 20 27 0d 0a | 31 30 39 30 20 27 09 44 |1080 '..|1090 '.D|
|00002850| 49 4d 45 4e 53 49 4f 4e | 20 41 4c 4c 20 41 52 52 |IMENSION| ALL ARR|
|00002860| 41 59 53 20 48 45 52 45 | 0d 0a 31 31 30 30 20 27 |AYS HERE|..1100 '|
|00002870| 20 54 48 49 53 20 56 45 | 52 53 49 4f 4e 20 49 53 | THIS VE|RSION IS|
|00002880| 20 53 45 54 20 55 50 20 | 54 4f 20 48 41 4e 44 4c | SET UP |TO HANDL|
|00002890| 45 20 55 50 20 54 4f 20 | 31 30 30 20 44 41 54 41 |E UP TO |100 DATA|
|000028a0| 20 50 4f 49 4e 54 53 0d | 0a 31 31 31 30 20 44 49 | POINTS.|.1110 DI|
|000028b0| 4d 20 58 28 31 30 30 29 | 2c 20 59 28 31 30 30 29 |M X(100)|, Y(100)|
|000028c0| 2c 20 53 49 47 59 28 31 | 30 30 29 2c 20 59 46 49 |, SIGY(1|00), YFI|
|000028d0| 54 28 31 30 30 29 2c 20 | 57 47 48 54 28 31 30 30 |T(100), |WGHT(100|
|000028e0| 29 0d 0a 31 31 32 30 20 | 27 0d 0a 31 31 33 30 20 |)..1120 |'..1130 |
|000028f0| 4d 41 58 4e 3d 20 31 30 | 30 09 09 27 53 45 54 20 |MAXN= 10|0..'SET |
|00002900| 4d 41 58 49 4d 55 4d 20 | 23 20 4f 46 20 50 4f 49 |MAXIMUM |# OF POI|
|00002910| 4e 54 53 0d 0a 31 31 34 | 30 20 27 0d 0a 31 31 35 |NTS..114|0 '..115|
|00002920| 30 20 27 09 20 4f 50 45 | 4e 20 26 20 52 45 41 44 |0 '. OPE|N & READ|
|00002930| 20 54 48 45 20 49 4e 50 | 55 54 20 44 41 54 41 20 | THE INP|UT DATA |
|00002940| 53 45 54 0d 0a 31 31 35 | 35 20 50 52 49 4e 54 22 |SET..115|5 PRINT"|
|00002950| 20 50 72 6f 67 72 61 6d | 20 4c 49 4e 46 49 54 2c | Program| LINFIT,|
|00002960| 20 20 47 2e 20 4d 2e 20 | 52 65 73 63 68 2c 20 46 | G. M. |Resch, F|
|00002970| 65 62 2e 20 37 39 22 0d | 0a 31 31 36 30 20 49 4e |eb. 79".|.1160 IN|
|00002980| 50 55 54 20 22 4e 61 6d | 65 20 6f 66 20 49 6e 70 |PUT "Nam|e of Inp|
|00002990| 75 74 20 64 61 74 61 20 | 73 65 74 3d 22 3b 4e 24 |ut data |set=";N$|
|000029a0| 0d 0a 31 31 37 30 20 4f | 50 45 4e 22 49 22 2c 32 |..1170 O|PEN"I",2|
|000029b0| 2c 4e 24 0d 0a 31 31 38 | 30 20 46 4f 52 20 49 3d |,N$..118|0 FOR I=|
|000029c0| 31 20 54 4f 20 4d 41 58 | 4e 0d 0a 31 31 39 30 20 |1 TO MAX|N..1190 |
|000029d0| 20 20 49 46 20 45 4f 46 | 28 32 29 20 20 47 4f 54 | IF EOF|(2) GOT|
|000029e0| 4f 20 31 32 32 30 0d 0a | 31 32 30 30 20 20 20 49 |O 1220..|1200 I|
|000029f0| 4e 50 55 54 23 32 2c 58 | 28 49 29 2c 59 28 49 29 |NPUT#2,X|(I),Y(I)|
|00002a00| 2c 53 49 47 59 28 49 29 | 0d 0a 31 32 31 30 20 4e |,SIGY(I)|..1210 N|
|00002a10| 45 58 54 20 49 0d 0a 31 | 32 32 30 20 4e 50 54 53 |EXT I..1|220 NPTS|
|00002a20| 3d 49 2d 31 20 20 3a 20 | 43 4c 4f 53 45 23 32 0d |=I-1 : |CLOSE#2.|
|00002a30| 0a 31 32 33 30 20 50 52 | 49 4e 54 22 54 48 45 52 |.1230 PR|INT"THER|
|00002a40| 45 20 41 52 45 20 22 3b | 4e 50 54 53 3b 22 20 44 |E ARE ";|NPTS;" D|
|00002a50| 41 54 41 20 50 4f 49 4e | 54 53 22 0d 0a 31 32 34 |ATA POIN|TS"..124|
|00002a60| 30 20 27 0d 0a 31 32 35 | 30 20 27 09 20 49 4e 50 |0 '..125|0 '. INP|
|00002a70| 55 54 20 54 48 45 20 41 | 50 52 49 4f 52 49 20 56 |UT THE A|PRIORI V|
|00002a80| 41 4c 55 45 53 20 46 4f | 52 20 54 48 45 20 50 41 |ALUES FO|R THE PA|
|00002a90| 52 41 4d 45 54 45 52 53 | 0d 0a 31 32 35 35 20 50 |RAMETERS|..1255 P|
|00002aa0| 52 49 4e 54 0d 0a 31 32 | 36 30 20 50 52 49 4e 54 |RINT..12|60 PRINT|
|00002ab0| 22 20 20 4d 4f 44 45 3d | 2b 31 2c 20 57 47 48 54 |" MODE=|+1, WGHT|
|00002ac0| 28 49 29 3d 31 2f 53 49 | 47 59 28 49 29 5e 32 22 |(I)=1/SI|GY(I)^2"|
|00002ad0| 0d 0a 31 32 37 30 20 50 | 52 49 4e 54 22 20 20 20 |..1270 P|RINT" |
|00002ae0| 20 20 20 3d 20 30 2c 20 | 57 47 48 54 28 49 29 3d | = 0, |WGHT(I)=|
|00002af0| 31 22 0d 0a 31 32 38 30 | 20 50 52 49 4e 54 22 20 |1"..1280| PRINT" |
|00002b00| 20 20 20 20 20 3d 2d 31 | 2c 20 57 47 48 54 28 49 | =-1|, WGHT(I|
|00002b10| 29 3d 31 2f 59 28 49 29 | 22 0d 0a 31 32 39 30 20 |)=1/Y(I)|"..1290 |
|00002b20| 50 52 49 4e 54 22 20 20 | 20 20 20 20 3d 2d 32 2c |PRINT" | =-2,|
|00002b30| 20 57 47 48 54 28 49 29 | 3d 53 49 47 59 28 49 29 | WGHT(I)|=SIGY(I)|
|00002b40| 22 0d 0a 31 33 30 30 20 | 49 4e 50 55 54 22 53 50 |"..1300 |INPUT"SP|
|00002b50| 45 43 49 46 59 20 54 48 | 45 20 57 45 49 47 48 54 |ECIFY TH|E WEIGHT|
|00002b60| 49 4e 47 20 4d 4f 44 45 | 20 28 2b 31 2c 30 2c 2d |ING MODE| (+1,0,-|
|00002b70| 31 2c 20 4f 52 20 2d 32 | 29 20 2d 22 3b 4d 4f 44 |1, OR -2|) -";MOD|
|00002b80| 45 0d 0a 31 33 31 30 20 | 27 0d 0a 31 33 32 30 20 |E..1310 |'..1320 |
|00002b90| 20 47 4f 53 55 42 20 31 | 37 35 30 20 20 09 27 43 | GOSUB 1|750 .'C|
|00002ba0| 41 4c 4c 20 4c 49 4e 46 | 49 54 0d 0a 31 33 33 30 |ALL LINF|IT..1330|
|00002bb0| 20 27 0d 0a 31 33 34 30 | 20 27 09 53 55 4d 4d 41 | '..1340| '.SUMMA|
|00002bc0| 52 49 5a 45 20 54 48 45 | 20 52 45 53 55 4c 54 53 |RIZE THE| RESULTS|
|00002bd0| 0d 0a 31 33 35 30 20 50 | 52 49 4e 54 0d 0a 31 33 |..1350 P|RINT..13|
|00002be0| 36 30 20 50 52 49 4e 54 | 22 20 59 20 49 4e 54 45 |60 PRINT|" Y INTE|
|00002bf0| 52 43 45 50 54 3d 22 3b | 41 2c 22 20 53 49 47 4d |RCEPT=";|A," SIGM|
|00002c00| 41 3d 22 3b 53 49 47 4d | 41 41 0d 0a 31 33 37 30 |A=";SIGM|AA..1370|
|00002c10| 20 50 52 49 4e 54 22 20 | 20 20 20 20 20 20 53 4c | PRINT" | SL|
|00002c20| 4f 50 45 3d 22 3b 42 2c | 22 20 53 49 47 4d 41 3d |OPE=";B,|" SIGMA=|
|00002c30| 22 3b 53 49 47 4d 41 42 | 0d 0a 31 33 38 30 20 50 |";SIGMAB|..1380 P|
|00002c40| 52 49 4e 54 22 20 4c 49 | 4e 45 41 52 20 43 4f 52 |RINT" LI|NEAR COR|
|00002c50| 52 45 4c 41 54 49 4f 4e | 20 43 4f 45 46 49 43 49 |RELATION| COEFICI|
|00002c60| 45 4e 54 3d 22 3b 52 0d | 0a 31 33 39 30 20 50 52 |ENT=";R.|.1390 PR|
|00002c70| 49 4e 54 0d 0a 31 34 30 | 30 20 50 52 49 4e 54 22 |INT..140|0 PRINT"|
|00002c80| 58 28 49 29 22 2c 22 59 | 28 49 29 22 2c 22 59 46 |X(I)","Y|(I)","YF|
|00002c90| 49 54 22 2c 22 59 2d 59 | 46 49 54 22 0d 0a 31 34 |IT","Y-Y|FIT"..14|
|00002ca0| 31 30 20 46 4f 52 20 49 | 3d 31 20 54 4f 20 4e 50 |10 FOR I|=1 TO NP|
|00002cb0| 54 53 0d 0a 31 34 32 30 | 20 20 59 46 49 54 28 49 |TS..1420| YFIT(I|
|00002cc0| 29 3d 20 41 20 2b 42 2a | 58 28 49 29 0d 0a 31 34 |)= A +B*|X(I)..14|
|00002cd0| 33 30 20 20 50 52 49 4e | 54 20 58 28 49 29 2c 59 |30 PRIN|T X(I),Y|
|00002ce0| 28 49 29 2c 59 46 49 54 | 28 49 29 2c 28 59 28 49 |(I),YFIT|(I),(Y(I|
|00002cf0| 29 2d 59 46 49 54 28 49 | 29 29 0d 0a 31 34 34 30 |)-YFIT(I|))..1440|
|00002d00| 20 4e 45 58 54 20 49 0d | 0a 31 34 35 30 20 45 4e | NEXT I.|.1450 EN|
|00002d10| 44 0d 0a 31 34 36 30 20 | 52 45 4d 2a 2a 2a 2a 2a |D..1460 |REM*****|
|00002d20| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00002d30| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00002d40| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00002d50| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 0d 0a 31 34 37 30 |********|**..1470|
|00002d60| 20 27 20 20 53 55 42 52 | 4f 55 54 49 4e 45 20 4c | ' SUBR|OUTINE L|
|00002d70| 49 4e 46 49 54 28 58 2c | 20 59 2c 20 53 49 47 59 |INFIT(X,| Y, SIGY|
|00002d80| 2c 20 4e 50 54 53 2c 20 | 4d 4f 44 45 2c 20 41 2c |, NPTS, |MODE, A,|
|00002d90| 20 53 49 47 4d 41 41 2c | 0d 0a 31 34 38 30 20 27 | SIGMAA,|..1480 '|
|00002da0| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00002db0| 20 20 20 20 42 2c 20 53 | 49 47 4d 41 42 2c 20 52 | B, S|IGMAB, R|
|00002dc0| 29 0d 0a 31 34 39 30 20 | 27 0d 0a 31 35 30 30 20 |)..1490 |'..1500 |
|00002dd0| 27 09 41 20 52 4f 55 54 | 49 4e 45 20 54 48 41 54 |'.A ROUT|INE THAT|
|00002de0| 20 57 49 4c 4c 20 46 49 | 4e 44 20 54 48 45 20 4c | WILL FI|ND THE L|
|00002df0| 45 41 53 54 2d 53 51 55 | 41 52 45 53 20 46 49 54 |EAST-SQU|ARES FIT|
|00002e00| 20 54 4f 20 41 20 4c 49 | 4e 45 41 52 0d 0a 31 35 | TO A LI|NEAR..15|
|00002e10| 31 30 20 27 20 46 55 4e | 43 54 49 4f 4e 2e 0d 0a |10 ' FUN|CTION...|
|00002e20| 31 35 32 30 20 27 09 09 | 59 20 3d 20 41 20 2b 20 |1520 '..|Y = A + |
|00002e30| 42 2a 58 0d 0a 31 35 33 | 30 20 27 0d 0a 31 35 35 |B*X..153|0 '..155|
|00002e40| 30 20 27 20 20 50 52 4f | 47 52 41 4d 4d 45 52 3a |0 ' PRO|GRAMMER:|
|00002e50| 20 20 47 2e 20 4d 2e 20 | 52 65 73 63 68 2c 20 20 | G. M. |Resch, |
|00002e60| 46 65 62 2e 20 31 39 37 | 39 0d 0a 31 35 36 30 20 |Feb. 197|9..1560 |
|00002e70| 27 0d 0a 31 35 37 30 20 | 27 20 50 41 52 41 4d 45 |'..1570 |' PARAME|
|00002e80| 54 45 52 53 3a 0d 0a 31 | 35 38 30 20 27 20 20 58 |TERS:..1|580 ' X|
|00002e90| 3d 41 52 52 41 59 20 4f | 46 20 49 4e 44 45 50 45 |=ARRAY O|F INDEPE|
|00002ea0| 4e 44 2e 20 56 41 52 49 | 41 42 4c 45 2c 20 20 59 |ND. VARI|ABLE, Y|
|00002eb0| 3d 41 52 52 41 59 20 4f | 46 20 44 45 50 45 4e 44 |=ARRAY O|F DEPEND|
|00002ec0| 45 4e 54 20 56 41 52 49 | 41 42 4c 45 0d 0a 31 35 |ENT VARI|ABLE..15|
|00002ed0| 39 30 20 27 20 20 53 49 | 47 59 3d 53 54 41 4e 44 |90 ' SI|GY=STAND|
|00002ee0| 41 52 44 20 44 45 56 2e | 20 4f 46 20 59 2c 20 20 |ARD DEV.| OF Y, |
|00002ef0| 4e 50 54 53 3d 23 20 4f | 46 20 50 41 49 52 53 20 |NPTS=# O|F PAIRS |
|00002f00| 4f 46 20 44 41 54 41 20 | 50 4f 49 4e 54 53 0d 0a |OF DATA |POINTS..|
|00002f10| 31 36 30 30 20 27 20 20 | 4d 4f 44 45 3d 2b 31 2c |1600 ' |MODE=+1,|
|00002f20| 20 57 47 48 54 28 49 29 | 3d 31 2f 53 49 47 59 28 | WGHT(I)|=1/SIGY(|
|00002f30| 49 29 5e 32 0d 0a 31 36 | 31 30 20 27 09 20 20 20 |I)^2..16|10 '. |
|00002f40| 20 3d 20 30 2c 20 57 47 | 48 54 28 49 29 3d 31 0d | = 0, WG|HT(I)=1.|
|00002f50| 0a 31 36 32 30 20 27 20 | 20 20 20 20 20 3d 2d 31 |.1620 ' | =-1|
|00002f60| 2c 20 57 47 48 54 28 49 | 29 3d 31 2f 59 28 49 29 |, WGHT(I|)=1/Y(I)|
|00002f70| 0d 0a 31 36 33 30 20 27 | 09 20 20 20 20 3d 2d 32 |..1630 '|. =-2|
|00002f80| 2c 20 57 47 48 54 28 49 | 29 3d 53 49 47 59 28 49 |, WGHT(I|)=SIGY(I|
|00002f90| 29 0d 0a 31 36 34 30 20 | 27 20 20 41 3d 20 59 20 |)..1640 |' A= Y |
|00002fa0| 49 4e 54 45 52 43 45 50 | 54 20 4f 46 20 46 49 54 |INTERCEP|T OF FIT|
|00002fb0| 54 45 44 20 53 54 52 41 | 49 47 48 54 20 4c 49 4e |TED STRA|IGHT LIN|
|00002fc0| 45 0d 0a 31 36 35 30 20 | 27 09 53 49 47 4d 41 41 |E..1650 |'.SIGMAA|
|00002fd0| 3d 53 54 4e 44 20 44 45 | 56 20 4f 46 20 41 0d 0a |=STND DE|V OF A..|
|00002fe0| 31 36 36 30 20 27 20 20 | 42 3d 20 53 4c 4f 50 45 |1660 ' |B= SLOPE|
|00002ff0| 20 4f 46 20 46 49 54 54 | 45 44 20 53 54 52 41 49 | OF FITT|ED STRAI|
|00003000| 47 48 54 20 4c 49 4e 45 | 0d 0a 31 36 37 30 20 27 |GHT LINE|..1670 '|
|00003010| 20 20 53 49 47 4d 41 42 | 3d 20 53 54 4e 44 20 44 | SIGMAB|= STND D|
|00003020| 45 56 20 4f 46 20 42 0d | 0a 31 36 38 30 20 27 20 |EV OF B.|.1680 ' |
|00003030| 20 52 3d 20 4c 49 4e 45 | 41 52 20 43 4f 52 52 45 | R= LINE|AR CORRE|
|00003040| 4c 41 54 49 4f 4e 20 43 | 4f 45 46 49 43 49 45 4e |LATION C|OEFICIEN|
|00003050| 54 0d 0a 31 36 39 30 20 | 27 0d 0a 31 37 30 30 20 |T..1690 |'..1700 |
|00003060| 27 20 53 55 42 52 4f 55 | 54 49 4e 45 53 20 43 41 |' SUBROU|TINES CA|
|00003070| 4c 4c 45 44 3b 0d 0a 31 | 37 31 30 20 27 20 20 20 |LLED;..1|710 ' |
|00003080| 4e 4f 4e 45 20 0d 0a 31 | 37 32 30 20 27 0d 0a 31 |NONE ..1|720 '..1|
|00003090| 37 33 30 20 27 2d 2d 2d | 20 53 54 41 52 54 20 4f |730 '---| START O|
|000030a0| 46 20 50 52 4f 47 52 41 | 4d 20 43 4f 44 45 20 2d |F PROGRA|M CODE -|
|000030b0| 2d 2d 0d 0a 31 37 34 30 | 20 27 0d 0a 31 37 35 30 |--..1740| '..1750|
|000030c0| 20 53 55 4d 3d 30 20 20 | 3a 20 53 55 4d 58 3d 30 | SUM=0 |: SUMX=0|
|000030d0| 20 20 3a 20 53 55 4d 59 | 3d 30 0d 0a 31 37 36 30 | : SUMY|=0..1760|
|000030e0| 20 53 55 4d 58 32 3d 30 | 20 20 3a 20 53 55 4d 58 | SUMX2=0| : SUMX|
|000030f0| 59 3d 30 20 20 3a 20 53 | 55 4d 59 32 3d 30 0d 0a |Y=0 : S|UMY2=0..|
|00003100| 31 37 37 30 20 46 4f 52 | 20 49 3d 31 20 54 4f 20 |1770 FOR| I=1 TO |
|00003110| 4e 50 54 53 0d 0a 31 37 | 38 30 20 20 49 46 20 4d |NPTS..17|80 IF M|
|00003120| 4f 44 45 3d 31 20 20 54 | 48 45 4e 20 57 47 48 54 |ODE=1 T|HEN WGHT|
|00003130| 28 49 29 3d 31 2f 28 53 | 49 47 59 28 49 29 2a 53 |(I)=1/(S|IGY(I)*S|
|00003140| 49 47 59 28 49 29 29 20 | 20 3a 20 47 4f 54 4f 20 |IGY(I)) | : GOTO |
|00003150| 31 38 32 30 0d 0a 31 37 | 39 30 20 20 49 46 20 4d |1820..17|90 IF M|
|00003160| 4f 44 45 3d 30 20 20 54 | 48 45 4e 20 57 47 48 54 |ODE=0 T|HEN WGHT|
|00003170| 28 49 29 3d 31 20 20 3a | 20 47 4f 54 4f 20 31 38 |(I)=1 :| GOTO 18|
|00003180| 32 30 0d 0a 31 38 30 30 | 20 20 49 46 20 4d 4f 44 |20..1800| IF MOD|
|00003190| 45 3d 2d 31 20 54 48 45 | 4e 20 57 47 48 54 28 49 |E=-1 THE|N WGHT(I|
|000031a0| 29 3d 31 2f 41 42 53 28 | 59 28 49 29 29 20 20 3a |)=1/ABS(|Y(I)) :|
|000031b0| 20 47 4f 54 4f 20 31 38 | 32 30 0d 0a 31 38 31 30 | GOTO 18|20..1810|
|000031c0| 20 20 49 46 20 4d 4f 44 | 45 3d 2d 32 20 54 48 45 | IF MOD|E=-2 THE|
|000031d0| 4e 20 57 47 48 54 28 49 | 29 3d 53 49 47 59 28 49 |N WGHT(I|)=SIGY(I|
|000031e0| 29 0d 0a 31 38 32 30 20 | 20 53 55 4d 3d 20 53 55 |)..1820 | SUM= SU|
|000031f0| 4d 20 2b 20 57 47 48 54 | 28 49 29 0d 0a 31 38 33 |M + WGHT|(I)..183|
|00003200| 30 20 20 53 55 4d 58 3d | 20 53 55 4d 58 20 2b 57 |0 SUMX=| SUMX +W|
|00003210| 47 48 54 28 49 29 2a 58 | 28 49 29 0d 0a 31 38 34 |GHT(I)*X|(I)..184|
|00003220| 30 20 20 53 55 4d 59 3d | 20 53 55 4d 59 20 2b 57 |0 SUMY=| SUMY +W|
|00003230| 47 48 54 28 49 29 2a 59 | 28 49 29 0d 0a 31 38 35 |GHT(I)*Y|(I)..185|
|00003240| 30 20 20 53 55 4d 58 32 | 3d 20 53 55 4d 58 32 20 |0 SUMX2|= SUMX2 |
|00003250| 2b 57 47 48 54 28 49 29 | 2a 58 28 49 29 2a 58 28 |+WGHT(I)|*X(I)*X(|
|00003260| 49 29 0d 0a 31 38 36 30 | 20 20 53 55 4d 58 59 3d |I)..1860| SUMXY=|
|00003270| 20 53 55 4d 58 59 20 2b | 57 47 48 54 28 49 29 2a | SUMXY +|WGHT(I)*|
|00003280| 58 28 49 29 2a 59 28 49 | 29 0d 0a 31 38 37 30 20 |X(I)*Y(I|)..1870 |
|00003290| 20 53 55 4d 59 32 3d 20 | 53 55 4d 59 32 20 2b 57 | SUMY2= |SUMY2 +W|
|000032a0| 47 48 54 28 49 29 2a 59 | 28 49 29 2a 59 28 49 29 |GHT(I)*Y|(I)*Y(I)|
|000032b0| 0d 0a 31 38 38 30 20 4e | 45 58 54 20 49 0d 0a 31 |..1880 N|EXT I..1|
|000032c0| 38 39 30 20 27 0d 0a 31 | 39 30 30 20 27 09 43 41 |890 '..1|900 '.CA|
|000032d0| 4c 43 20 43 4f 45 46 46 | 49 43 49 45 4e 54 53 20 |LC COEFF|ICIENTS |
|000032e0| 26 20 53 54 4e 44 20 44 | 45 56 49 41 54 49 4f 4e |& STND D|EVIATION|
|000032f0| 53 0d 0a 31 39 31 30 20 | 44 45 4c 54 41 3d 20 53 |S..1910 |DELTA= S|
|00003300| 55 4d 2a 53 55 4d 58 32 | 20 2d 53 55 4d 58 2a 53 |UM*SUMX2| -SUMX*S|
|00003310| 55 4d 58 0d 0a 31 39 32 | 30 20 41 3d 20 28 53 55 |UMX..192|0 A= (SU|
|00003320| 4d 58 32 2a 53 55 4d 59 | 20 2d 53 55 4d 58 2a 53 |MX2*SUMY| -SUMX*S|
|00003330| 55 4d 58 59 29 2f 44 45 | 4c 54 41 0d 0a 31 39 33 |UMXY)/DE|LTA..193|
|00003340| 30 20 42 3d 20 28 53 55 | 4d 58 59 2a 53 55 4d 20 |0 B= (SU|MXY*SUM |
|00003350| 2d 53 55 4d 58 2a 53 55 | 4d 59 29 2f 44 45 4c 54 |-SUMX*SU|MY)/DELT|
|00003360| 41 0d 0a 31 39 34 30 20 | 49 46 20 4d 4f 44 45 3d |A..1940 |IF MODE=|
|00003370| 30 20 20 47 4f 54 4f 20 | 31 39 36 30 0d 0a 31 39 |0 GOTO |1960..19|
|00003380| 35 30 20 56 41 52 4e 43 | 45 3d 31 20 20 3a 20 47 |50 VARNC|E=1 : G|
|00003390| 4f 54 4f 20 31 39 38 30 | 0d 0a 31 39 36 30 20 43 |OTO 1980|..1960 C|
|000033a0| 3d 20 4e 50 54 53 2d 32 | 0d 0a 31 39 37 30 20 56 |= NPTS-2|..1970 V|
|000033b0| 41 52 4e 43 45 3d 20 28 | 53 55 4d 59 32 2b 41 2a |ARNCE= (|SUMY2+A*|
|000033c0| 41 2a 53 55 4d 2b 42 2a | 42 2a 53 55 4d 58 32 20 |A*SUM+B*|B*SUMX2 |
|000033d0| 2d 32 2a 28 41 2a 53 55 | 4d 59 2b 42 2a 53 55 4d |-2*(A*SU|MY+B*SUM|
|000033e0| 58 59 2d 41 2a 42 2a 53 | 55 4d 58 29 29 2f 43 0d |XY-A*B*S|UMX))/C.|
|000033f0| 0a 31 39 38 30 20 53 49 | 47 4d 41 41 3d 20 53 51 |.1980 SI|GMAA= SQ|
|00003400| 52 28 56 41 52 4e 43 45 | 2a 53 55 4d 58 32 2f 44 |R(VARNCE|*SUMX2/D|
|00003410| 45 4c 54 41 29 0d 0a 31 | 39 39 30 20 53 49 47 4d |ELTA)..1|990 SIGM|
|00003420| 41 42 3d 20 53 51 52 28 | 56 41 52 4e 43 45 2a 53 |AB= SQR(|VARNCE*S|
|00003430| 55 4d 2f 44 45 4c 54 41 | 29 0d 0a 32 30 30 30 20 |UM/DELTA|)..2000 |
|00003440| 52 3d 20 28 53 55 4d 2a | 53 55 4d 58 59 2d 53 55 |R= (SUM*|SUMXY-SU|
|00003450| 4d 58 2a 53 55 4d 59 29 | 2f 53 51 52 28 44 45 4c |MX*SUMY)|/SQR(DEL|
|00003460| 54 41 2a 28 53 55 4d 2a | 53 55 4d 59 32 2d 53 55 |TA*(SUM*|SUMY2-SU|
|00003470| 4d 59 2a 53 55 4d 59 29 | 29 0d 0a 32 30 31 30 20 |MY*SUMY)|)..2010 |
|00003480| 52 45 54 55 52 4e 0d 0a | 1a 1a 1a 1a 1a 1a 1a 1a |RETURN..|........|
|00003490| 1a 1a 1a 1a 1a 1a 1a 1a | 1a 1a 1a 1a 1a 1a 1a 1a |........|........|
|000034a0| 1a 1a 1a 1a 1a 1a 1a 1a | 1a 1a 1a 1a 1a 1a 1a 1a |........|........|
|000034b0| 1a 1a 1a 1a 1a 1a 1a 1a | 1a 1a 1a 1a 1a 1a 1a 1a |........|........|
|000034c0| 1a 1a 1a 1a 1a 1a 1a 1a | 1a 1a 1a 1a 1a 1a 1a 1a |........|........|
|000034d0| 1a 1a 1a 1a 1a 1a 1a 1a | 1a 1a 1a 1a 1a 1a 1a 1a |........|........|
|000034e0| 1a 1a 1a 1a 1a 1a 1a 1a | 1a 1a 1a 1a 1a 1a 1a 1a |........|........|
|000034f0| 1a 1a 1a 1a 1a 1a 1a 1a | 1a 1a 1a 1a 1a 1a 1a 1a |........|........|
|00003500| 31 30 30 30 20 27 20 20 | 20 20 20 20 20 20 20 20 |1000 ' | |
|00003510| 20 20 20 20 20 20 2a 2a | 2a 20 20 20 50 52 4f 47 | **|* PROG|
|00003520| 52 41 4d 20 4d 4f 53 54 | 46 49 54 20 20 20 2a 2a |RAM MOST|FIT **|
|00003530| 2a 0d 0a 31 30 31 30 20 | 27 0d 0a 31 30 32 30 20 |*..1010 |'..1020 |
|00003540| 27 20 20 41 20 44 52 49 | 56 45 52 20 50 52 4f 47 |' A DRI|VER PROG|
|00003550| 52 41 4d 20 46 4f 52 20 | 54 48 45 20 52 45 47 52 |RAM FOR |THE REGR|
|00003560| 45 53 53 20 53 55 42 52 | 4f 55 54 49 4e 45 20 54 |ESS SUBR|OUTINE T|
|00003570| 48 41 54 20 57 49 4c 4c | 20 46 49 4e 44 0d 0a 31 |HAT WILL| FIND..1|
|00003580| 30 33 30 20 27 20 54 48 | 45 20 50 41 52 41 4d 45 |030 ' TH|E PARAME|
|00003590| 54 45 52 53 20 49 4e 20 | 41 20 47 49 56 45 4e 20 |TERS IN |A GIVEN |
|000035a0| 46 55 4e 43 54 49 4f 4e | 20 54 48 41 54 20 47 49 |FUNCTION| THAT GI|
|000035b0| 56 45 20 54 48 45 20 42 | 45 53 54 20 46 49 54 0d |VE THE B|EST FIT.|
|000035c0| 0a 31 30 34 30 20 27 20 | 28 49 4e 20 41 20 4c 45 |.1040 ' |(IN A LE|
|000035d0| 41 53 54 20 53 51 55 41 | 52 45 53 20 53 45 4e 53 |AST SQUA|RES SENS|
|000035e0| 45 29 20 54 4f 20 41 20 | 53 45 54 20 4f 46 20 49 |E) TO A |SET OF I|
|000035f0| 4e 50 55 54 20 44 41 54 | 41 2e 20 20 58 20 26 20 |NPUT DAT|A. X & |
|00003600| 59 20 41 52 45 0d 0a 31 | 30 35 30 20 27 20 54 48 |Y ARE..1|050 ' TH|
|00003610| 45 20 49 4e 44 45 50 45 | 4e 44 45 4e 54 20 41 4e |E INDEPE|NDENT AN|
|00003620| 44 20 44 45 50 45 4e 44 | 45 4e 54 20 56 41 52 49 |D DEPEND|ENT VARI|
|00003630| 41 42 4c 45 53 20 52 45 | 53 50 45 43 54 49 56 45 |ABLES RE|SPECTIVE|
|00003640| 4c 59 2e 20 20 54 48 45 | 0d 0a 31 30 36 30 20 27 |LY. THE|..1060 '|
|00003650| 20 46 55 4e 43 54 49 4f | 4e 20 54 48 41 54 20 44 | FUNCTIO|N THAT D|
|00003660| 45 54 45 52 4d 49 4e 45 | 53 20 59 20 47 49 56 45 |ETERMINE|S Y GIVE|
|00003670| 4e 20 41 20 56 41 4c 55 | 45 20 4f 46 20 58 20 4d |N A VALU|E OF X M|
|00003680| 55 53 54 20 42 45 20 53 | 50 45 43 49 46 49 45 44 |UST BE S|PECIFIED|
|00003690| 0d 0a 31 30 37 30 20 27 | 20 49 4e 20 54 48 45 20 |..1070 '| IN THE |
|000036a0| 53 55 42 52 4f 55 54 49 | 4e 45 20 46 43 54 4e 20 |SUBROUTI|NE FCTN |
|000036b0| 28 53 45 45 20 45 4e 44 | 20 4f 46 20 53 4f 55 52 |(SEE END| OF SOUR|
|000036c0| 43 45 20 43 4f 44 45 29 | 20 41 4e 44 20 4d 55 53 |CE CODE)| AND MUS|
|000036d0| 54 20 42 45 20 4c 49 4e | 45 41 52 0d 0a 31 30 38 |T BE LIN|EAR..108|
|000036e0| 30 20 27 20 49 4e 20 54 | 48 45 20 43 4f 45 46 46 |0 ' IN T|HE COEFF|
|000036f0| 49 43 49 45 4e 54 53 2e | 0d 0a 31 30 39 30 20 27 |ICIENTS.|..1090 '|
|00003700| 0d 0a 31 31 30 30 20 27 | 09 50 52 4f 47 52 41 4d |..1100 '|.PROGRAM|
|00003710| 4d 45 52 3a 20 20 47 65 | 6f 72 67 65 20 4d 2e 20 |MER: Ge|orge M. |
|00003720| 52 65 73 63 68 2c 20 20 | 46 65 62 2e 20 31 39 37 |Resch, |Feb. 197|
|00003730| 39 0d 0a 31 31 31 30 20 | 27 0d 0a 31 31 32 30 20 |9..1110 |'..1120 |
|00003740| 27 09 44 49 4d 45 4e 53 | 49 4f 4e 20 41 4c 4c 20 |'.DIMENS|ION ALL |
|00003750| 41 52 52 41 59 53 20 48 | 45 52 45 0d 0a 31 31 33 |ARRAYS H|ERE..113|
|00003760| 30 20 27 20 54 48 49 53 | 20 56 45 52 53 49 4f 4e |0 ' THIS| VERSION|
|00003770| 20 49 53 20 53 45 54 20 | 55 50 20 54 4f 20 48 41 | IS SET |UP TO HA|
|00003780| 4e 44 4c 45 20 55 50 20 | 54 4f 20 31 30 30 20 44 |NDLE UP |TO 100 D|
|00003790| 41 54 41 20 50 4f 49 4e | 54 53 20 41 4e 44 0d 0a |ATA POIN|TS AND..|
|000037a0| 31 31 34 30 20 27 20 41 | 53 20 4d 41 4e 59 20 41 |1140 ' A|S MANY A|
|000037b0| 53 20 31 30 20 50 41 52 | 41 4d 45 54 45 52 53 2e |S 10 PAR|AMETERS.|
|000037c0| 0d 0a 31 31 35 30 20 44 | 49 4d 20 20 58 28 31 30 |..1150 D|IM X(10|
|000037d0| 30 29 2c 20 59 28 31 30 | 30 29 2c 20 53 49 47 59 |0), Y(10|0), SIGY|
|000037e0| 28 31 30 30 29 2c 20 59 | 46 49 54 28 31 30 30 29 |(100), Y|FIT(100)|
|000037f0| 2c 20 57 47 48 54 28 31 | 30 30 29 0d 0a 31 31 36 |, WGHT(1|00)..116|
|00003800| 30 20 44 49 4d 20 20 41 | 28 31 30 29 2c 20 49 4b |0 DIM A|(10), IK|
|00003810| 28 31 30 29 2c 20 4a 4b | 28 31 30 29 2c 20 58 4d |(10), JK|(10), XM|
|00003820| 4e 23 28 31 30 29 2c 20 | 53 49 47 58 23 28 31 30 |N#(10), |SIGX#(10|
|00003830| 29 2c 20 41 52 59 23 28 | 31 30 2c 31 30 29 0d 0a |), ARY#(|10,10)..|
|00003840| 31 31 37 30 20 27 0d 0a | 31 31 38 30 20 4d 41 58 |1170 '..|1180 MAX|
|00003850| 4e 3d 20 31 30 30 20 20 | 3a 20 4d 41 58 50 3d 31 |N= 100 |: MAXP=1|
|00003860| 30 09 27 53 45 54 20 4d | 41 58 49 4d 55 4d 20 23 |0.'SET M|AXIMUM #|
|00003870| 20 4f 46 20 50 4f 49 4e | 54 20 41 4e 44 20 50 41 | OF POIN|T AND PA|
|00003880| 52 4d 53 0d 0a 31 31 39 | 30 20 27 0d 0a 31 32 30 |RMS..119|0 '..120|
|00003890| 30 20 27 09 20 4f 50 45 | 4e 20 26 20 52 45 41 44 |0 '. OPE|N & READ|
|000038a0| 20 54 48 45 20 49 4e 50 | 55 54 20 44 41 54 41 20 | THE INP|UT DATA |
|000038b0| 53 45 54 0d 0a 31 32 30 | 35 20 50 52 49 4e 54 22 |SET..120|5 PRINT"|
|000038c0| 20 50 72 6f 67 72 61 6d | 20 4d 4f 53 54 46 49 54 | Program| MOSTFIT|
|000038d0| 2c 20 20 47 2e 20 4d 2e | 20 52 65 73 63 68 2c 20 |, G. M.| Resch, |
|000038e0| 46 65 62 2e 20 37 39 22 | 0d 0a 31 32 31 30 20 49 |Feb. 79"|..1210 I|
|000038f0| 4e 50 55 54 20 22 4e 61 | 6d 65 20 6f 66 20 49 6e |NPUT "Na|me of In|
|00003900| 70 75 74 20 64 61 74 61 | 20 73 65 74 3d 22 3b 4e |put data| set=";N|
|00003910| 24 0d 0a 31 32 32 30 20 | 4f 50 45 4e 22 49 22 2c |$..1220 |OPEN"I",|
|00003920| 32 2c 4e 24 0d 0a 31 32 | 33 30 20 46 4f 52 20 49 |2,N$..12|30 FOR I|
|00003930| 3d 31 20 54 4f 20 4d 41 | 58 4e 0d 0a 31 32 34 30 |=1 TO MA|XN..1240|
|00003940| 20 20 20 49 46 20 45 4f | 46 28 32 29 20 20 47 4f | IF EO|F(2) GO|
|00003950| 54 4f 20 31 32 37 30 0d | 0a 31 32 35 30 20 20 20 |TO 1270.|.1250 |
|00003960| 49 4e 50 55 54 23 32 2c | 58 28 49 29 2c 59 28 49 |INPUT#2,|X(I),Y(I|
|00003970| 29 2c 53 49 47 59 28 49 | 29 0d 0a 31 32 36 30 20 |),SIGY(I|)..1260 |
|00003980| 4e 45 58 54 20 49 0d 0a | 31 32 37 30 20 4e 50 54 |NEXT I..|1270 NPT|
|00003990| 53 3d 49 2d 31 20 20 3a | 20 43 4c 4f 53 45 23 32 |S=I-1 :| CLOSE#2|
|000039a0| 0d 0a 31 32 38 30 20 50 | 52 49 4e 54 22 54 48 45 |..1280 P|RINT"THE|
|000039b0| 52 45 20 41 52 45 20 22 | 3b 4e 50 54 53 3b 22 20 |RE ARE "|;NPTS;" |
|000039c0| 44 41 54 41 20 50 4f 49 | 4e 54 53 22 0d 0a 31 32 |DATA POI|NTS"..12|
|000039d0| 39 30 20 50 52 49 4e 54 | 0d 0a 31 33 30 30 20 27 |90 PRINT|..1300 '|
|000039e0| 0d 0a 31 33 31 30 20 27 | 09 20 49 4e 50 55 54 20 |..1310 '|. INPUT |
|000039f0| 56 41 4c 55 45 53 20 46 | 4f 52 20 54 48 45 20 50 |VALUES F|OR THE P|
|00003a00| 41 52 41 4d 45 54 45 52 | 53 0d 0a 31 33 32 30 20 |ARAMETER|S..1320 |
|00003a10| 50 52 49 4e 54 22 20 20 | 4d 4f 44 45 3d 2b 31 20 |PRINT" |MODE=+1 |
|00003a20| 49 4e 53 54 52 2e 20 57 | 47 48 54 28 49 29 3d 31 |INSTR. W|GHT(I)=1|
|00003a30| 2f 53 49 47 59 28 49 29 | 5e 32 22 0d 0a 31 33 33 |/SIGY(I)|^2"..133|
|00003a40| 30 20 50 52 49 4e 54 22 | 20 20 20 20 20 20 3d 20 |0 PRINT"| = |
|00003a50| 30 20 4e 4f 20 20 20 20 | 20 57 47 48 54 28 49 29 |0 NO | WGHT(I)|
|00003a60| 3d 31 22 0d 0a 31 33 34 | 30 20 50 52 49 4e 54 22 |=1"..134|0 PRINT"|
|00003a70| 20 20 20 20 20 20 3d 2d | 31 20 53 54 41 54 2e 20 | =-|1 STAT. |
|00003a80| 20 57 47 48 54 28 49 29 | 3d 31 2f 59 28 49 29 22 | WGHT(I)|=1/Y(I)"|
|00003a90| 0d 0a 31 33 35 30 20 50 | 52 49 4e 54 22 20 20 20 |..1350 P|RINT" |
|00003aa0| 20 20 20 3d 2d 32 20 41 | 52 42 49 54 2e 20 57 47 | =-2 A|RBIT. WG|
|00003ab0| 48 54 28 49 29 3d 53 49 | 47 59 28 49 29 22 0d 0a |HT(I)=SI|GY(I)"..|
|00003ac0| 31 33 36 30 20 49 4e 50 | 55 54 22 53 50 45 43 49 |1360 INP|UT"SPECI|
|00003ad0| 46 59 20 54 48 45 20 57 | 45 49 47 48 54 49 4e 47 |FY THE W|EIGHTING|
|00003ae0| 20 4d 4f 44 45 20 28 2b | 31 2c 30 2c 2d 31 2c 20 | MODE (+|1,0,-1, |
|00003af0| 4f 52 20 2d 32 29 20 2d | 22 3b 4d 4f 44 45 0d 0a |OR -2) -|";MODE..|
|00003b00| 31 33 37 30 20 49 4e 50 | 55 54 22 48 4f 57 20 4d |1370 INP|UT"HOW M|
|00003b10| 41 4e 59 20 54 45 52 4d | 53 3f 20 20 4e 54 52 4d |ANY TERM|S? NTRM|
|00003b20| 3d 22 3b 4e 54 52 4d 0d | 0a 31 33 37 31 20 49 46 |=";NTRM.|.1371 IF|
|00003b30| 20 28 4e 54 52 4d 3c 4e | 50 54 53 20 41 4e 44 20 | (NTRM<N|PTS AND |
|00003b40| 4e 54 52 4d 3c 3d 4d 41 | 58 50 29 20 20 47 4f 54 |NTRM<=MA|XP) GOT|
|00003b50| 4f 20 31 33 38 30 0d 0a | 31 33 37 32 20 50 52 49 |O 1380..|1372 PRI|
|00003b60| 4e 54 22 20 54 68 65 20 | 6e 75 6d 62 65 72 20 6f |NT" The |number o|
|00003b70| 66 20 74 65 72 6d 73 20 | 6d 75 73 74 20 62 65 20 |f terms |must be |
|00003b80| 6c 65 73 73 20 74 68 61 | 6e 20 31 31 22 0d 0a 31 |less tha|n 11"..1|
|00003b90| 33 37 33 20 47 4f 54 4f | 20 31 33 37 30 0d 0a 31 |373 GOTO| 1370..1|
|00003ba0| 33 38 30 20 4e 43 3d 30 | 0d 0a 31 33 39 30 20 27 |380 NC=0|..1390 '|
|00003bb0| 46 4f 52 20 4a 3d 31 20 | 54 4f 20 31 30 30 0d 0a |FOR J=1 |TO 100..|
|00003bc0| 31 34 30 30 20 27 20 50 | 52 49 4e 54 22 20 55 53 |1400 ' P|RINT" US|
|00003bd0| 45 20 58 5e 22 3b 4a 2c | 0d 0a 31 34 31 30 20 27 |E X^";J,|..1410 '|
|00003be0| 20 49 4e 50 55 54 22 20 | 59 20 4f 52 20 4e 20 22 | INPUT" |Y OR N "|
|00003bf0| 3b 41 24 0d 0a 31 34 32 | 30 20 27 20 49 46 20 41 |;A$..142|0 ' IF A|
|00003c00| 24 3d 22 4e 22 20 20 47 | 4f 54 4f 20 31 34 34 30 |$="N" G|OTO 1440|
|00003c10| 0d 0a 31 34 33 30 20 27 | 20 4e 43 3d 4e 43 2b 31 |..1430 '| NC=NC+1|
|00003c20| 20 20 3a 20 4a 54 52 4d | 28 4e 43 29 3d 4a 20 20 | : JTRM|(NC)=J |
|00003c30| 3a 20 49 46 20 4e 43 3d | 4e 54 52 4d 20 20 47 4f |: IF NC=|NTRM GO|
|00003c40| 54 4f 20 31 34 36 30 0d | 0a 31 34 34 30 20 27 4e |TO 1460.|.1440 'N|
|00003c50| 45 58 54 20 4a 0d 0a 31 | 34 35 30 20 27 0d 0a 31 |EXT J..1|450 '..1|
|00003c60| 34 36 30 20 20 47 4f 53 | 55 42 20 31 38 32 30 20 |460 GOS|UB 1820 |
|00003c70| 20 09 27 43 41 4c 4c 20 | 52 45 47 52 45 53 53 0d | .'CALL |REGRESS.|
|00003c80| 0a 31 34 37 30 20 27 0d | 0a 31 34 38 30 20 27 09 |.1470 '.|.1480 '.|
|00003c90| 53 55 4d 4d 41 52 49 5a | 45 20 54 48 45 20 52 45 |SUMMARIZ|E THE RE|
|00003ca0| 53 55 4c 54 53 0d 0a 31 | 34 39 30 20 50 52 49 4e |SULTS..1|490 PRIN|
|00003cb0| 54 0d 0a 31 35 30 30 20 | 50 52 49 4e 54 22 41 28 |T..1500 |PRINT"A(|
|00003cc0| 30 29 3d 22 3b 41 30 2c | 22 53 49 47 4d 41 41 28 |0)=";A0,|"SIGMAA(|
|00003cd0| 30 29 3d 22 3b 53 49 47 | 41 30 0d 0a 31 35 31 30 |0)=";SIG|A0..1510|
|00003ce0| 20 46 4f 52 20 4a 3d 31 | 20 54 4f 20 4e 54 52 4d | FOR J=1| TO NTRM|
|00003cf0| 0d 0a 31 35 32 30 20 20 | 50 52 49 4e 54 22 41 28 |..1520 |PRINT"A(|
|00003d00| 22 3b 4a 3b 22 29 3d 22 | 3b 41 28 4a 29 2c 22 53 |";J;")="|;A(J),"S|
|00003d10| 49 47 4d 41 41 28 22 3b | 4a 3b 22 29 3d 22 3b 53 |IGMAA(";|J;")=";S|
|00003d20| 49 47 41 28 4a 29 2c 22 | 52 28 22 3b 4a 3b 22 29 |IGA(J),"|R(";J;")|
|00003d30| 3d 22 3b 52 28 4a 29 0d | 0a 31 35 33 30 20 4e 45 |=";R(J).|.1530 NE|
|00003d40| 58 54 20 4a 0d 0a 31 35 | 34 30 20 50 52 49 4e 54 |XT J..15|40 PRINT|
|00003d50| 22 43 48 49 20 53 51 55 | 41 52 45 3d 22 3b 43 48 |"CHI SQU|ARE=";CH|
|00003d60| 53 51 0d 0a 31 35 35 30 | 20 50 52 49 4e 54 0d 0a |SQ..1550| PRINT..|
|00003d70| 31 35 36 30 20 50 52 49 | 4e 54 22 58 28 49 29 22 |1560 PRI|NT"X(I)"|
|00003d80| 2c 22 59 28 49 29 22 2c | 22 59 46 49 54 22 2c 22 |,"Y(I)",|"YFIT","|
|00003d90| 59 2d 59 46 49 54 22 0d | 0a 31 35 37 30 20 46 4f |Y-YFIT".|.1570 FO|
|00003da0| 52 20 49 3d 31 20 54 4f | 20 4e 50 54 53 0d 0a 31 |R I=1 TO| NPTS..1|
|00003db0| 35 38 30 20 20 50 52 49 | 4e 54 20 58 28 49 29 2c |580 PRI|NT X(I),|
|00003dc0| 59 28 49 29 2c 59 46 49 | 54 28 49 29 2c 28 59 28 |Y(I),YFI|T(I),(Y(|
|00003dd0| 49 29 2d 59 46 49 54 28 | 49 29 29 0d 0a 31 35 39 |I)-YFIT(|I))..159|
|00003de0| 30 20 4e 45 58 54 20 49 | 0d 0a 31 36 30 30 20 45 |0 NEXT I|..1600 E|
|00003df0| 4e 44 0d 0a 31 36 31 30 | 20 52 45 4d 2d 2d 2d 2d |ND..1610| REM----|
|00003e00| 2d 2d 2d 2d 2d 2d 2d 2d | 2d 2d 2d 2d 2d 2d 2d 2d |--------|--------|
|00003e10| 2d 2d 2d 2d 2d 2d 2d 2d | 2d 2d 2d 2d 2d 2d 2d 2d |--------|--------|
|00003e20| 2d 2d 2d 2d 2d 2d 2d 2d | 2d 2d 2d 2d 2d 2d 2d 2d |--------|--------|
|00003e30| 2d 2d 2d 2d 2d 2d 2d 2d | 2d 2d 2d 2d 2d 2d 2d 2d |--------|--------|
|00003e40| 2d 2d 0d 0a 31 36 32 30 | 20 27 2a 2a 2a 20 20 53 |--..1620| '*** S|
|00003e50| 55 42 52 4f 55 54 49 4e | 45 20 52 45 47 52 45 53 |UBROUTIN|E REGRES|
|00003e60| 28 58 2c 20 59 2c 20 53 | 49 47 59 2c 20 4e 50 54 |(X, Y, S|IGY, NPT|
|00003e70| 53 2c 20 4e 54 52 4d 2c | 20 4d 4f 44 45 2c 20 59 |S, NTRM,| MODE, Y|
|00003e80| 46 49 54 2c 0d 0a 31 36 | 33 30 20 27 20 20 20 20 |FIT,..16|30 ' |
|00003e90| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00003ea0| 20 20 20 41 30 2c 20 41 | 2c 20 53 49 47 41 30 2c | A0, A|, SIGA0,|
|00003eb0| 20 53 49 47 41 41 2c 20 | 52 2c 20 52 4d 55 4c 2c | SIGAA, |R, RMUL,|
|00003ec0| 20 43 48 53 51 2c 20 46 | 54 53 54 29 0d 0a 31 36 | CHSQ, F|TST)..16|
|00003ed0| 33 35 20 27 0d 0a 31 36 | 34 30 20 27 20 20 50 55 |35 '..16|40 ' PU|
|00003ee0| 52 50 4f 53 45 2d 20 4d | 41 4b 45 20 41 20 46 49 |RPOSE- M|AKE A FI|
|00003ef0| 54 20 54 4f 20 44 41 54 | 41 20 57 49 54 48 20 41 |T TO DAT|A WITH A|
|00003f00| 20 53 50 45 43 49 46 49 | 45 44 20 46 55 4e 43 54 | SPECIFI|ED FUNCT|
|00003f10| 49 4f 4e 20 54 48 41 54 | 0d 0a 31 36 35 30 20 27 |ION THAT|..1650 '|
|00003f20| 09 09 20 49 53 20 4c 49 | 4e 45 41 52 20 49 4e 20 |.. IS LI|NEAR IN |
|00003f30| 54 48 45 20 43 4f 45 46 | 46 49 43 49 45 4e 54 53 |THE COEF|FICIENTS|
|00003f40| 2e 0d 0a 31 36 36 30 20 | 27 20 20 58 3d 49 4e 44 |...1660 |' X=IND|
|00003f50| 45 50 45 4e 44 45 4e 54 | 20 56 41 52 49 41 42 4c |EPENDENT| VARIABL|
|00003f60| 45 2c 20 20 59 3d 44 45 | 50 45 4e 44 45 4e 54 20 |E, Y=DE|PENDENT |
|00003f70| 56 41 52 49 41 42 4c 45 | 0d 0a 31 36 37 30 20 27 |VARIABLE|..1670 '|
|00003f80| 20 20 53 49 47 59 3d 53 | 49 47 4d 41 53 20 4f 46 | SIGY=S|IGMAS OF|
|00003f90| 20 59 2c 20 20 4e 50 54 | 53 3d 23 20 4f 46 20 44 | Y, NPT|S=# OF D|
|00003fa0| 41 54 41 20 50 4f 49 4e | 54 53 0d 0a 31 36 38 30 |ATA POIN|TS..1680|
|00003fb0| 20 27 20 20 4e 54 52 4d | 3d 23 20 4f 46 20 43 4f | ' NTRM|=# OF CO|
|00003fc0| 45 46 46 2e 2c 20 20 4d | 3d 41 52 52 41 59 20 4f |EFF., M|=ARRAY O|
|00003fd0| 46 20 49 4e 43 4c 2f 52 | 45 4a 45 43 54 20 43 52 |F INCL/R|EJECT CR|
|00003fe0| 49 54 45 52 49 41 20 46 | 4f 52 20 46 43 54 4e 0d |ITERIA F|OR FCTN.|
|00003ff0| 0a 31 36 39 30 20 27 20 | 20 4d 4f 44 45 3d 2b 31 |.1690 ' | MODE=+1|
|00004000| 20 49 4e 53 54 52 2e 20 | 57 47 48 54 28 49 29 3d | INSTR. |WGHT(I)=|
|00004010| 31 2f 53 49 47 59 28 49 | 29 5e 32 0d 0a 31 37 30 |1/SIGY(I|)^2..170|
|00004020| 30 20 27 20 20 20 20 20 | 20 3d 20 30 20 4e 4f 20 |0 ' | = 0 NO |
|00004030| 20 20 20 20 57 47 48 54 | 28 49 29 3d 31 0d 0a 31 | WGHT|(I)=1..1|
|00004040| 37 31 30 20 27 20 20 20 | 20 20 20 3d 2d 31 20 53 |710 ' | =-1 S|
|00004050| 54 41 54 2e 20 20 57 47 | 48 54 28 49 29 3d 31 2f |TAT. WG|HT(I)=1/|
|00004060| 59 28 49 29 0d 0a 31 37 | 32 30 20 27 20 20 20 20 |Y(I)..17|20 ' |
|00004070| 20 20 3d 2d 32 20 41 52 | 42 49 54 2e 20 57 47 48 | =-2 AR|BIT. WGH|
|00004080| 54 28 49 29 3d 53 49 47 | 59 28 49 29 0d 0a 31 37 |T(I)=SIG|Y(I)..17|
|00004090| 33 30 20 27 20 20 59 46 | 49 54 3d 43 41 4c 43 20 |30 ' YF|IT=CALC |
|000040a0| 56 41 4c 55 45 53 20 4f | 46 20 59 2c 20 20 41 30 |VALUES O|F Y, A0|
|000040b0| 3d 43 4f 4e 53 54 41 4e | 54 20 54 45 52 4d 0d 0a |=CONSTAN|T TERM..|
|000040c0| 31 37 34 30 20 27 20 20 | 41 3d 41 52 52 41 59 20 |1740 ' |A=ARRAY |
|000040d0| 4f 46 20 43 4f 45 46 46 | 2e 2c 20 20 53 49 47 41 |OF COEFF|., SIGA|
|000040e0| 30 3d 53 54 4e 44 20 44 | 45 56 20 4f 46 20 41 30 |0=STND D|EV OF A0|
|000040f0| 0d 0a 31 37 35 30 20 27 | 20 20 53 49 47 41 41 3d |..1750 '| SIGAA=|
|00004100| 53 54 4e 44 20 44 45 56 | 20 4f 46 20 41 2c 20 20 |STND DEV| OF A, |
|00004110| 52 3d 41 52 52 41 59 20 | 4f 46 20 4c 49 4e 45 41 |R=ARRAY |OF LINEA|
|00004120| 52 20 43 4f 52 52 2e 20 | 43 4f 45 46 46 2e 0d 0a |R CORR. |COEFF...|
|00004130| 31 37 36 30 20 27 20 20 | 52 4d 55 4c 3d 4d 55 4c |1760 ' |RMUL=MUL|
|00004140| 54 49 50 4c 45 20 4c 49 | 4e 45 41 52 20 43 4f 52 |TIPLE LI|NEAR COR|
|00004150| 52 2e 20 43 4f 45 46 46 | 2e 2c 20 20 43 48 53 51 |R. COEFF|., CHSQ|
|00004160| 3d 52 45 44 55 43 45 44 | 20 43 48 49 20 53 51 55 |=REDUCED| CHI SQU|
|00004170| 41 52 45 0d 0a 31 37 37 | 30 20 27 20 20 46 54 45 |ARE..177|0 ' FTE|
|00004180| 53 54 3d 56 41 4c 55 45 | 20 4f 46 20 46 20 46 4f |ST=VALUE| OF F FO|
|00004190| 52 20 54 45 53 54 20 4f | 46 20 46 49 54 0d 0a 31 |R TEST O|F FIT..1|
|000041a0| 37 38 30 20 27 20 20 53 | 55 42 52 4f 55 54 49 4e |780 ' S|UBROUTIN|
|000041b0| 45 53 20 52 45 51 55 49 | 52 45 44 2d 2d 2d 2d 2d |ES REQUI|RED-----|
|000041c0| 0d 0a 31 37 39 30 20 27 | 20 20 20 20 20 46 43 54 |..1790 '| FCT|
|000041d0| 4e 28 58 2c 49 2c 4a 2c | 4d 29 20 26 20 4d 41 54 |N(X,I,J,|M) & MAT|
|000041e0| 49 4e 56 28 41 52 52 41 | 59 2c 20 4e 54 52 4d 2c |INV(ARRA|Y, NTRM,|
|000041f0| 20 44 45 54 29 0d 0a 31 | 38 30 30 20 27 0d 0a 31 | DET)..1|800 '..1|
|00004200| 38 31 30 20 27 09 20 49 | 4e 49 54 20 53 55 4d 53 |810 '. I|NIT SUMS|
|00004210| 20 26 20 41 52 52 41 59 | 53 0d 0a 31 38 32 30 20 | & ARRAY|S..1820 |
|00004220| 53 55 4d 23 3d 30 21 20 | 20 3a 20 59 4d 4e 23 3d |SUM#=0! | : YMN#=|
|00004230| 30 21 20 20 3a 20 53 49 | 47 23 3d 30 21 20 20 3a |0! : SI|G#=0! :|
|00004240| 20 43 48 49 23 3d 30 21 | 20 20 3a 20 52 4d 55 4c | CHI#=0!| : RMUL|
|00004250| 3d 30 21 0d 0a 31 38 33 | 30 20 46 4f 52 20 49 3d |=0!..183|0 FOR I=|
|00004260| 31 20 54 4f 20 4e 50 54 | 53 20 20 3a 20 59 46 49 |1 TO NPT|S : YFI|
|00004270| 54 28 49 29 3d 30 20 20 | 3a 20 4e 45 58 54 20 49 |T(I)=0 |: NEXT I|
|00004280| 0d 0a 31 38 34 30 20 46 | 4f 52 20 4a 3d 31 20 54 |..1840 F|OR J=1 T|
|00004290| 4f 20 4e 54 52 4d 0d 0a | 31 38 35 30 20 20 20 58 |O NTRM..|1850 X|
|000042a0| 4d 4e 23 28 4a 29 3d 30 | 21 20 20 3a 20 53 49 47 |MN#(J)=0|! : SIG|
|000042b0| 58 23 28 4a 29 3d 30 21 | 20 20 3a 20 52 28 4a 29 |X#(J)=0!| : R(J)|
|000042c0| 3d 30 21 20 20 3a 20 41 | 28 4a 29 3d 30 21 20 20 |=0! : A|(J)=0! |
|000042d0| 3a 20 53 49 47 41 28 4a | 29 3d 30 21 0d 0a 31 38 |: SIGA(J|)=0!..18|
|000042e0| 36 30 20 20 46 4f 52 20 | 4b 3d 31 20 54 4f 20 4e |60 FOR |K=1 TO N|
|000042f0| 54 52 4d 0d 0a 31 38 37 | 30 20 20 20 20 41 52 59 |TRM..187|0 ARY|
|00004300| 23 28 4a 2c 4b 29 3d 30 | 21 0d 0a 31 38 38 30 20 |#(J,K)=0|!..1880 |
|00004310| 4e 45 58 54 20 4b 2c 4a | 0d 0a 31 38 39 30 20 27 |NEXT K,J|..1890 '|
|00004320| 09 2a 2a 2a 20 41 43 43 | 55 4d 55 4c 41 54 45 20 |.*** ACC|UMULATE |
|00004330| 57 45 49 47 48 54 45 44 | 20 53 55 4d 53 20 2a 2a |WEIGHTED| SUMS **|
|00004340| 2a 0d 0a 31 39 30 30 20 | 46 4f 52 20 49 3d 31 20 |*..1900 |FOR I=1 |
|00004350| 54 4f 20 4e 50 54 53 0d | 0a 31 39 31 30 20 20 49 |TO NPTS.|.1910 I|
|00004360| 46 20 4d 4f 44 45 3c 3e | 2b 31 20 47 4f 54 4f 20 |F MODE<>|+1 GOTO |
|00004370| 31 39 33 30 0d 0a 31 39 | 32 30 20 20 20 57 47 48 |1930..19|20 WGH|
|00004380| 54 28 49 29 3d 31 2f 28 | 53 49 47 59 28 49 29 2a |T(I)=1/(|SIGY(I)*|
|00004390| 53 49 47 59 28 49 29 29 | 20 20 3a 20 47 4f 54 4f |SIGY(I))| : GOTO|
|000043a0| 20 32 30 31 30 0d 0a 31 | 39 33 30 20 20 49 46 20 | 2010..1|930 IF |
|000043b0| 4d 4f 44 45 3c 3e 30 20 | 47 4f 54 4f 20 31 39 35 |MODE<>0 |GOTO 195|
|000043c0| 30 0d 0a 31 39 34 30 20 | 20 20 57 47 48 54 28 49 |0..1940 | WGHT(I|
|000043d0| 29 3d 31 21 20 20 3a 20 | 47 4f 54 4f 20 32 30 31 |)=1! : |GOTO 201|
|000043e0| 30 0d 0a 31 39 35 30 20 | 20 49 46 20 4d 4f 44 45 |0..1950 | IF MODE|
|000043f0| 3c 3e 2d 31 20 20 47 4f | 54 4f 20 31 39 38 30 0d |<>-1 GO|TO 1980.|
|00004400| 0a 31 39 36 30 20 20 20 | 49 46 20 59 28 49 29 3c |.1960 |IF Y(I)<|
|00004410| 30 20 54 48 45 4e 20 57 | 47 48 54 28 49 29 3d 31 |0 THEN W|GHT(I)=1|
|00004420| 2f 28 2d 59 28 49 29 29 | 20 45 4c 53 45 20 20 57 |/(-Y(I))| ELSE W|
|00004430| 47 48 54 28 49 29 3d 31 | 2f 59 28 49 29 0d 0a 31 |GHT(I)=1|/Y(I)..1|
|00004440| 39 37 30 20 20 20 47 4f | 54 4f 20 32 30 31 30 0d |970 GO|TO 2010.|
|00004450| 0a 31 39 38 30 20 20 49 | 46 20 4d 4f 44 45 3d 2d |.1980 I|F MODE=-|
|00004460| 32 20 20 47 4f 54 4f 20 | 32 30 30 30 0d 0a 31 39 |2 GOTO |2000..19|
|00004470| 39 30 20 20 20 50 52 49 | 4e 54 22 4d 4f 44 45 20 |90 PRI|NT"MODE |
|00004480| 53 50 45 43 49 46 49 43 | 41 54 49 4f 4e 20 45 52 |SPECIFIC|ATION ER|
|00004490| 52 4f 52 22 20 20 3a 20 | 52 45 54 55 52 4e 0d 0a |ROR" : |RETURN..|
|000044a0| 32 30 30 30 20 20 57 47 | 48 54 28 49 29 3d 53 49 |2000 WG|HT(I)=SI|
|000044b0| 47 59 28 49 29 0d 0a 32 | 30 31 30 20 20 53 55 4d |GY(I)..2|010 SUM|
|000044c0| 23 3d 53 55 4d 23 2b 57 | 47 48 54 28 49 29 0d 0a |#=SUM#+W|GHT(I)..|
|000044d0| 32 30 32 30 20 20 59 4d | 4e 23 3d 59 4d 4e 23 2b |2020 YM|N#=YMN#+|
|000044e0| 57 47 48 54 28 49 29 2a | 59 28 49 29 0d 0a 32 30 |WGHT(I)*|Y(I)..20|
|000044f0| 33 30 20 46 4f 52 20 4a | 3d 31 20 54 4f 20 4e 54 |30 FOR J|=1 TO NT|
|00004500| 52 4d 0d 0a 32 30 34 30 | 20 20 20 20 47 4f 53 55 |RM..2040| GOSU|
|00004510| 42 20 33 34 34 30 20 20 | 3a 52 45 4d 2a 20 43 41 |B 3440 |:REM* CA|
|00004520| 4c 43 20 46 43 54 4e 0d | 0a 32 30 35 30 20 20 20 |LC FCTN.|.2050 |
|00004530| 20 58 4d 4e 23 28 4a 29 | 3d 58 4d 4e 23 28 4a 29 | XMN#(J)|=XMN#(J)|
|00004540| 2b 57 47 48 54 28 49 29 | 2a 46 43 54 4e 0d 0a 32 |+WGHT(I)|*FCTN..2|
|00004550| 30 36 30 20 4e 45 58 54 | 20 4a 2c 49 0d 0a 32 30 |060 NEXT| J,I..20|
|00004560| 37 30 20 59 4d 4e 23 3d | 59 4d 4e 23 2f 53 55 4d |70 YMN#=|YMN#/SUM|
|00004570| 23 0d 0a 32 30 38 30 20 | 27 0d 0a 32 30 39 30 20 |#..2080 |'..2090 |
|00004580| 46 4f 52 20 4a 3d 31 20 | 54 4f 20 4e 54 52 4d 20 |FOR J=1 |TO NTRM |
|00004590| 20 3a 20 58 4d 4e 23 28 | 4a 29 3d 58 4d 4e 23 28 | : XMN#(|J)=XMN#(|
|000045a0| 4a 29 2f 53 55 4d 23 20 | 20 3a 20 20 4e 45 58 54 |J)/SUM# | : NEXT|
|000045b0| 20 4a 0d 0a 32 31 30 30 | 20 57 4d 4e 3d 53 55 4d | J..2100| WMN=SUM|
|000045c0| 23 2f 4e 50 54 53 0d 0a | 32 31 31 30 20 46 4f 52 |#/NPTS..|2110 FOR|
|000045d0| 20 49 3d 31 20 54 4f 20 | 4e 50 54 53 20 20 3a 20 | I=1 TO |NPTS : |
|000045e0| 57 47 48 54 28 49 29 3d | 57 47 48 54 28 49 29 2f |WGHT(I)=|WGHT(I)/|
|000045f0| 57 4d 4e 20 20 3a 20 20 | 4e 45 58 54 20 49 0d 0a |WMN : |NEXT I..|
|00004600| 32 31 32 30 20 27 09 2a | 2a 2a 20 41 43 43 55 4d |2120 '.*|** ACCUM|
|00004610| 55 4c 41 54 45 20 4d 41 | 54 52 49 43 45 53 20 52 |ULATE MA|TRICES R|
|00004620| 20 26 20 41 52 59 23 0d | 0a 32 31 33 30 20 46 4f | & ARY#.|.2130 FO|
|00004630| 52 20 49 3d 31 20 54 4f | 20 4e 50 54 53 0d 0a 32 |R I=1 TO| NPTS..2|
|00004640| 31 34 30 20 20 53 49 47 | 23 3d 53 49 47 23 2b 57 |140 SIG|#=SIG#+W|
|00004650| 47 48 54 28 49 29 2a 28 | 59 28 49 29 2d 59 4d 4e |GHT(I)*(|Y(I)-YMN|
|00004660| 23 29 2a 28 59 28 49 29 | 2d 59 4d 4e 23 29 0d 0a |#)*(Y(I)|-YMN#)..|
|00004670| 32 31 35 30 20 46 4f 52 | 20 4a 3d 31 20 54 4f 20 |2150 FOR| J=1 TO |
|00004680| 4e 54 52 4d 0d 0a 32 31 | 36 30 20 20 20 47 4f 53 |NTRM..21|60 GOS|
|00004690| 55 42 20 33 34 34 30 20 | 20 3a 52 45 4d 2a 20 43 |UB 3440 | :REM* C|
|000046a0| 41 4c 43 20 46 43 54 4e | 28 58 2c 49 2c 4a 2c 4d |ALC FCTN|(X,I,J,M|
|000046b0| 29 0d 0a 32 31 37 30 20 | 20 20 53 49 47 58 23 28 |)..2170 | SIGX#(|
|000046c0| 4a 29 3d 53 49 47 58 23 | 28 4a 29 20 2b 57 47 48 |J)=SIGX#|(J) +WGH|
|000046d0| 54 28 49 29 2a 28 46 43 | 54 4e 20 2d 58 4d 4e 23 |T(I)*(FC|TN -XMN#|
|000046e0| 28 4a 29 29 2a 28 46 43 | 54 4e 2d 58 4d 4e 23 28 |(J))*(FC|TN-XMN#(|
|000046f0| 4a 29 29 0d 0a 32 31 38 | 30 20 20 20 52 28 4a 29 |J))..218|0 R(J)|
|00004700| 3d 52 28 4a 29 20 2b 57 | 47 48 54 28 49 29 2a 28 |=R(J) +W|GHT(I)*(|
|00004710| 46 43 54 4e 20 2d 58 4d | 4e 23 28 4a 29 29 2a 28 |FCTN -XM|N#(J))*(|
|00004720| 59 28 49 29 2d 59 4d 4e | 23 29 0d 0a 32 31 39 30 |Y(I)-YMN|#)..2190|
|00004730| 20 20 20 46 4a 3d 46 43 | 54 4e 0d 0a 32 32 30 30 | FJ=FC|TN..2200|
|00004740| 20 46 4f 52 20 4b 3d 31 | 20 54 4f 20 4a 0d 0a 32 | FOR K=1| TO J..2|
|00004750| 32 31 30 20 20 20 20 4a | 4a 3d 4a 20 20 3a 20 4a |210 J|J=J : J|
|00004760| 3d 4b 20 20 3a 20 47 4f | 53 55 42 20 33 34 34 30 |=K : GO|SUB 3440|
|00004770| 0d 0a 32 32 32 30 20 20 | 20 20 46 4b 3d 46 43 54 |..2220 | FK=FCT|
|00004780| 4e 20 20 3a 20 4a 3d 4a | 4a 0d 0a 32 32 33 30 20 |N : J=J|J..2230 |
|00004790| 20 20 20 41 52 59 23 28 | 4a 2c 4b 29 3d 41 52 59 | ARY#(|J,K)=ARY|
|000047a0| 23 28 4a 2c 4b 29 20 2b | 57 47 48 54 28 49 29 2a |#(J,K) +|WGHT(I)*|
|000047b0| 28 46 4a 2d 58 4d 4e 23 | 28 4a 29 29 2a 28 46 4b |(FJ-XMN#|(J))*(FK|
|000047c0| 2d 58 4d 4e 23 28 4b 29 | 29 0d 0a 32 32 34 30 20 |-XMN#(K)|)..2240 |
|000047d0| 4e 45 58 54 20 4b 2c 4a | 2c 49 0d 0a 32 32 35 30 |NEXT K,J|,I..2250|
|000047e0| 20 27 0d 0a 32 32 36 30 | 20 46 52 45 31 3d 4e 50 | '..2260| FRE1=NP|
|000047f0| 54 53 2d 31 20 20 3a 20 | 53 49 47 23 3d 53 51 52 |TS-1 : |SIG#=SQR|
|00004800| 28 53 49 47 23 2f 46 52 | 45 31 29 0d 0a 32 32 37 |(SIG#/FR|E1)..227|
|00004810| 30 20 46 4f 52 20 4a 3d | 31 20 54 4f 20 4e 54 52 |0 FOR J=|1 TO NTR|
|00004820| 4d 0d 0a 32 32 38 30 20 | 20 53 49 47 58 23 28 4a |M..2280 | SIGX#(J|
|00004830| 29 3d 53 51 52 28 53 49 | 47 58 23 28 4a 29 2f 46 |)=SQR(SI|GX#(J)/F|
|00004840| 52 45 31 29 0d 0a 32 32 | 39 30 20 20 52 28 4a 29 |RE1)..22|90 R(J)|
|00004850| 3d 52 28 4a 29 2f 28 46 | 52 45 31 2a 53 49 47 58 |=R(J)/(F|RE1*SIGX|
|00004860| 23 28 4a 29 2a 53 49 47 | 23 29 0d 0a 32 33 30 30 |#(J)*SIG|#)..2300|
|00004870| 20 46 4f 52 20 4b 3d 31 | 20 54 4f 20 4a 0d 0a 32 | FOR K=1| TO J..2|
|00004880| 33 31 30 20 20 20 41 52 | 59 23 28 4a 2c 4b 29 3d |310 AR|Y#(J,K)=|
|00004890| 20 41 52 59 23 28 4a 2c | 4b 29 2f 28 46 52 45 31 | ARY#(J,|K)/(FRE1|
|000048a0| 2a 53 49 47 58 23 28 4a | 29 2a 53 49 47 58 23 28 |*SIGX#(J|)*SIGX#(|
|000048b0| 4b 29 29 0d 0a 32 33 32 | 30 20 20 20 41 52 59 23 |K))..232|0 ARY#|
|000048c0| 28 4b 2c 4a 29 3d 41 52 | 59 23 28 4a 2c 4b 29 0d |(K,J)=AR|Y#(J,K).|
|000048d0| 0a 32 33 33 30 20 4e 45 | 58 54 20 4b 0d 0a 32 33 |.2330 NE|XT K..23|
|000048e0| 34 30 20 4e 45 58 54 20 | 4a 0d 0a 32 33 35 30 20 |40 NEXT |J..2350 |
|000048f0| 27 09 2a 2a 2a 20 49 4e | 56 45 52 54 20 53 59 4d |'.*** IN|VERT SYM|
|00004900| 4d 45 54 52 49 43 20 4d | 41 54 52 49 58 0d 0a 32 |METRIC M|ATRIX..2|
|00004910| 33 36 30 20 4e 4f 52 3d | 4e 54 52 4d 20 20 3a 20 |360 NOR=|NTRM : |
|00004920| 47 4f 53 55 42 20 32 38 | 33 30 20 20 3a 52 45 4d |GOSUB 28|30 :REM|
|00004930| 2a 20 43 41 4c 4c 20 4d | 41 54 49 4e 56 0d 0a 32 |* CALL M|ATINV..2|
|00004940| 33 37 30 20 49 46 20 44 | 45 54 3c 3e 30 20 47 4f |370 IF D|ET<>0 GO|
|00004950| 54 4f 20 32 34 31 30 0d | 0a 32 33 38 30 20 41 30 |TO 2410.|.2380 A0|
|00004960| 3d 30 20 20 3a 20 53 49 | 47 41 30 3d 30 20 20 3a |=0 : SI|GA0=0 :|
|00004970| 20 52 4d 55 4c 3d 30 20 | 20 3a 20 43 48 53 51 3d | RMUL=0 | : CHSQ=|
|00004980| 30 20 20 3a 20 46 54 45 | 53 54 3d 30 0d 0a 32 33 |0 : FTE|ST=0..23|
|00004990| 39 30 20 52 45 54 55 52 | 4e 0d 0a 32 34 30 30 20 |90 RETUR|N..2400 |
|000049a0| 27 09 2a 2a 2a 20 43 41 | 4c 43 20 43 4f 45 46 46 |'.*** CA|LC COEFF|
|000049b0| 2c 20 46 49 54 2c 20 26 | 20 43 48 49 20 53 51 55 |, FIT, &| CHI SQU|
|000049c0| 41 52 45 0d 0a 32 34 31 | 30 20 41 30 3d 59 4d 4e |ARE..241|0 A0=YMN|
|000049d0| 23 0d 0a 32 34 32 30 20 | 46 4f 52 20 4a 3d 31 20 |#..2420 |FOR J=1 |
|000049e0| 54 4f 20 4e 54 52 4d 0d | 0a 32 34 33 30 20 20 46 |TO NTRM.|.2430 F|
|000049f0| 4f 52 20 4b 3d 31 20 54 | 4f 20 4e 54 52 4d 0d 0a |OR K=1 T|O NTRM..|
|00004a00| 32 34 34 30 20 20 20 20 | 41 28 4a 29 3d 41 28 4a |2440 |A(J)=A(J|
|00004a10| 29 20 2b 52 28 4b 29 2a | 41 52 59 23 28 4a 2c 4b |) +R(K)*|ARY#(J,K|
|00004a20| 29 0d 0a 32 34 35 30 20 | 20 4e 45 58 54 20 4b 0d |)..2450 | NEXT K.|
|00004a30| 0a 32 34 36 30 20 20 41 | 28 4a 29 3d 41 28 4a 29 |.2460 A|(J)=A(J)|
|00004a40| 2a 53 49 47 23 2f 53 49 | 47 58 23 28 4a 29 0d 0a |*SIG#/SI|GX#(J)..|
|00004a50| 32 34 37 30 20 20 41 30 | 3d 41 30 2d 41 28 4a 29 |2470 A0|=A0-A(J)|
|00004a60| 2a 58 4d 4e 23 28 4a 29 | 0d 0a 32 34 38 30 20 20 |*XMN#(J)|..2480 |
|00004a70| 46 4f 52 20 49 3d 31 20 | 54 4f 20 4e 50 54 53 0d |FOR I=1 |TO NPTS.|
|00004a80| 0a 32 34 39 30 20 20 20 | 47 4f 53 55 42 20 33 33 |.2490 |GOSUB 33|
|00004a90| 39 30 20 20 3a 52 45 4d | 2a 20 43 41 4c 43 20 46 |90 :REM|* CALC F|
|00004aa0| 43 54 4e 28 58 2c 49 2c | 4a 2c 4d 29 0d 0a 32 35 |CTN(X,I,|J,M)..25|
|00004ab0| 30 30 20 20 20 59 46 49 | 54 28 49 29 3d 59 46 49 |00 YFI|T(I)=YFI|
|00004ac0| 54 28 49 29 20 2b 41 28 | 4a 29 2a 46 43 54 4e 0d |T(I) +A(|J)*FCTN.|
|00004ad0| 0a 32 35 31 30 20 4e 45 | 58 54 20 49 2c 4a 0d 0a |.2510 NE|XT I,J..|
|00004ae0| 32 35 32 30 20 27 0d 0a | 32 35 33 30 20 46 4f 52 |2520 '..|2530 FOR|
|00004af0| 20 49 3d 31 20 54 4f 20 | 4e 50 54 53 0d 0a 32 35 | I=1 TO |NPTS..25|
|00004b00| 34 30 20 20 59 46 49 54 | 28 49 29 3d 59 46 49 54 |40 YFIT|(I)=YFIT|
|00004b10| 28 49 29 20 2b 41 30 0d | 0a 32 35 35 30 20 20 43 |(I) +A0.|.2550 C|
|00004b20| 48 49 23 3d 43 48 49 23 | 20 2b 57 47 48 54 28 49 |HI#=CHI#| +WGHT(I|
|00004b30| 29 2a 28 59 28 49 29 2d | 59 46 49 54 28 49 29 29 |)*(Y(I)-|YFIT(I))|
|00004b40| 2a 28 59 28 49 29 2d 59 | 46 49 54 28 49 29 29 0d |*(Y(I)-Y|FIT(I)).|
|00004b50| 0a 32 35 36 30 20 4e 45 | 58 54 20 49 0d 0a 32 35 |.2560 NE|XT I..25|
|00004b60| 37 30 20 27 0d 0a 32 35 | 38 30 20 46 52 45 4e 3d |70 '..25|80 FREN=|
|00004b70| 4e 50 54 53 2d 4e 54 52 | 4d 2d 31 20 20 3a 20 43 |NPTS-NTR|M-1 : C|
|00004b80| 48 53 51 3d 43 48 49 23 | 2a 57 4d 4e 2f 46 52 45 |HSQ=CHI#|*WMN/FRE|
|00004b90| 4e 0d 0a 32 35 39 30 20 | 27 09 2a 2a 2a 20 43 41 |N..2590 |'.*** CA|
|00004ba0| 4c 43 20 55 4e 43 45 52 | 54 41 49 4e 54 49 45 53 |LC UNCER|TAINTIES|
|00004bb0| 0d 0a 32 36 30 30 20 43 | 48 53 51 3d 43 48 49 23 |..2600 C|HSQ=CHI#|
|00004bc0| 0d 0a 32 36 31 30 20 49 | 46 20 4d 4f 44 45 3d 30 |..2610 I|F MODE=0|
|00004bd0| 20 20 47 4f 54 4f 20 32 | 36 33 30 0d 0a 32 36 32 | GOTO 2|630..262|
|00004be0| 30 20 56 41 52 4e 3d 31 | 2f 57 4d 4e 20 20 3a 20 |0 VARN=1|/WMN : |
|00004bf0| 47 4f 54 4f 20 32 36 34 | 30 0d 0a 32 36 33 30 20 |GOTO 264|0..2630 |
|00004c00| 56 41 52 4e 3d 43 48 53 | 51 0d 0a 32 36 34 30 20 |VARN=CHS|Q..2640 |
|00004c10| 46 4f 52 20 4a 3d 31 20 | 54 4f 20 4e 54 52 4d 0d |FOR J=1 |TO NTRM.|
|00004c20| 0a 32 36 35 30 20 20 53 | 49 47 41 28 4a 29 3d 41 |.2650 S|IGA(J)=A|
|00004c30| 52 59 23 28 4a 2c 4a 29 | 2a 56 41 52 4e 2f 28 46 |RY#(J,J)|*VARN/(F|
|00004c40| 52 45 31 2a 53 49 47 58 | 23 28 4a 29 2a 53 49 47 |RE1*SIGX|#(J)*SIG|
|00004c50| 58 23 28 4a 29 29 0d 0a | 32 36 36 30 20 20 53 49 |X#(J))..|2660 SI|
|00004c60| 47 41 28 4a 29 3d 53 51 | 52 28 53 49 47 41 28 4a |GA(J)=SQ|R(SIGA(J|
|00004c70| 29 29 0d 0a 32 36 37 30 | 20 20 52 4d 55 4c 3d 52 |))..2670| RMUL=R|
|00004c80| 4d 55 4c 20 2b 41 28 4a | 29 2a 52 28 4a 29 2a 53 |MUL +A(J|)*R(J)*S|
|00004c90| 49 47 58 23 28 4a 29 2f | 53 49 47 23 0d 0a 32 36 |IGX#(J)/|SIG#..26|
|00004ca0| 38 30 20 4e 45 58 54 20 | 4a 0d 0a 32 36 39 30 20 |80 NEXT |J..2690 |
|00004cb0| 46 52 45 4a 3d 4e 54 52 | 4d 20 20 3a 20 46 54 53 |FREJ=NTR|M : FTS|
|00004cc0| 54 3d 28 52 4d 55 4c 2f | 46 52 45 4a 29 2f 28 28 |T=(RMUL/|FREJ)/((|
|00004cd0| 31 2d 52 4d 55 4c 29 2f | 46 52 45 4e 29 0d 0a 32 |1-RMUL)/|FREN)..2|
|00004ce0| 37 30 30 20 52 4d 55 4c | 3d 53 51 52 28 52 4d 55 |700 RMUL|=SQR(RMU|
|00004cf0| 4c 29 20 20 3a 20 53 49 | 47 41 30 3d 56 41 52 4e |L) : SI|GA0=VARN|
|00004d00| 2f 4e 50 54 53 0d 0a 32 | 37 31 30 20 46 4f 52 20 |/NPTS..2|710 FOR |
|00004d10| 4a 3d 31 20 54 4f 20 4e | 54 52 4d 0d 0a 32 37 32 |J=1 TO N|TRM..272|
|00004d20| 30 20 46 4f 52 20 4b 3d | 31 20 54 4f 20 4e 54 52 |0 FOR K=|1 TO NTR|
|00004d30| 4d 0d 0a 32 37 33 30 20 | 20 53 49 47 41 30 3d 53 |M..2730 | SIGA0=S|
|00004d40| 49 47 41 30 20 2b 56 41 | 52 4e 2a 58 4d 4e 23 28 |IGA0 +VA|RN*XMN#(|
|00004d50| 4a 29 2a 58 4d 4e 23 28 | 4b 29 2a 41 52 59 23 28 |J)*XMN#(|K)*ARY#(|
|00004d60| 4a 2c 4b 29 2f 28 46 52 | 45 31 2a 53 49 47 58 23 |J,K)/(FR|E1*SIGX#|
|00004d70| 28 4a 29 2a 53 49 47 58 | 23 28 4b 29 29 0d 0a 32 |(J)*SIGX|#(K))..2|
|00004d80| 37 34 30 20 4e 45 58 54 | 20 4b 2c 4a 0d 0a 32 37 |740 NEXT| K,J..27|
|00004d90| 35 30 20 53 49 47 41 30 | 3d 53 51 52 28 53 49 47 |50 SIGA0|=SQR(SIG|
|00004da0| 41 30 29 0d 0a 32 37 36 | 30 20 52 45 54 55 52 4e |A0)..276|0 RETURN|
|00004db0| 0d 0a 32 37 37 30 20 52 | 45 4d 2d 2d 2d 2d 2d 2d |..2770 R|EM------|
|00004dc0| 2d 2d 2d 2d 2d 2d 2d 2d | 2d 2d 2d 2d 2d 2d 2d 2d |--------|--------|
|00004dd0| 2d 2d 2d 2d 2d 2d 2d 2d | 2d 2d 2d 2d 2d 2d 2d 2d |--------|--------|
|00004de0| 2d 2d 2d 2d 2d 2d 2d 2d | 2d 2d 2d 2d 2d 2d 2d 2d |--------|--------|
|00004df0| 2d 2d 2d 2d 2d 2d 2d 2d | 2d 2d 2d 2d 2d 2d 2d 2d |--------|--------|
|00004e00| 0d 0a 32 37 38 30 20 27 | 20 20 20 20 20 20 20 20 |..2780 '| |
|00004e10| 20 20 53 55 42 52 4f 55 | 54 49 4e 45 20 4d 41 54 | SUBROU|TINE MAT|
|00004e20| 49 4e 56 28 41 52 59 23 | 2c 4e 4f 52 2c 44 45 54 |INV(ARY#|,NOR,DET|
|00004e30| 29 0d 0a 32 37 39 30 20 | 27 20 20 41 20 52 4f 55 |)..2790 |' A ROU|
|00004e40| 54 49 4e 45 20 54 4f 20 | 49 4e 56 45 52 54 20 41 |TINE TO |INVERT A|
|00004e50| 20 53 59 4d 4d 45 54 52 | 49 43 20 4d 41 54 52 49 | SYMMETR|IC MATRI|
|00004e60| 58 20 41 4e 44 20 43 41 | 4c 43 55 4c 41 54 45 20 |X AND CA|LCULATE |
|00004e70| 49 54 53 0d 0a 32 38 30 | 30 20 27 20 20 44 45 54 |ITS..280|0 ' DET|
|00004e80| 45 52 4d 49 4e 41 4e 54 | 2e 20 20 41 52 59 23 28 |ERMINANT|. ARY#(|
|00004e90| 49 2c 4a 29 20 49 53 20 | 54 48 45 20 49 4e 50 55 |I,J) IS |THE INPU|
|00004ea0| 54 20 4d 41 54 52 49 58 | 20 4f 46 20 4f 52 44 45 |T MATRIX| OF ORDE|
|00004eb0| 52 3d 4e 4f 52 0d 0a 32 | 38 31 30 20 27 20 20 41 |R=NOR..2|810 ' A|
|00004ec0| 4e 44 20 49 53 20 52 45 | 50 4c 41 43 45 44 20 42 |ND IS RE|PLACED B|
|00004ed0| 59 20 49 54 53 20 49 4e | 56 45 52 53 45 0d 0a 32 |Y ITS IN|VERSE..2|
|00004ee0| 38 32 30 20 27 0d 0a 32 | 38 33 30 20 44 45 54 3d |820 '..2|830 DET=|
|00004ef0| 31 21 0d 0a 32 38 34 30 | 20 27 09 2a 2a 2a 20 46 |1!..2840| '.*** F|
|00004f00| 49 4e 44 20 4c 41 52 47 | 45 53 54 20 45 4c 45 4d |IND LARG|EST ELEM|
|00004f10| 45 4e 54 0d 0a 32 38 35 | 30 20 46 4f 52 20 4b 4b |ENT..285|0 FOR KK|
|00004f20| 25 3d 31 20 54 4f 20 4e | 4f 52 0d 0a 32 38 36 30 |%=1 TO N|OR..2860|
|00004f30| 20 41 4d 58 23 3d 30 21 | 0d 0a 32 38 37 30 20 46 | AMX#=0!|..2870 F|
|00004f40| 4f 52 20 49 49 25 3d 4b | 4b 25 20 54 4f 20 4e 4f |OR II%=K|K% TO NO|
|00004f50| 52 0d 0a 32 38 38 30 20 | 46 4f 52 20 4a 4a 25 3d |R..2880 |FOR JJ%=|
|00004f60| 4b 4b 25 20 54 4f 20 4e | 4f 52 0d 0a 32 38 39 30 |KK% TO N|OR..2890|
|00004f70| 20 20 20 49 46 20 41 42 | 53 28 41 4d 58 23 29 3e | IF AB|S(AMX#)>|
|00004f80| 41 42 53 28 41 52 59 23 | 28 49 49 25 2c 4a 4a 25 |ABS(ARY#|(II%,JJ%|
|00004f90| 29 29 20 47 4f 54 4f 20 | 32 39 31 30 0d 0a 32 39 |)) GOTO |2910..29|
|00004fa0| 30 30 20 20 20 41 4d 58 | 23 3d 41 52 59 23 28 49 |00 AMX|#=ARY#(I|
|00004fb0| 49 25 2c 4a 4a 25 29 20 | 20 3a 20 49 4b 28 4b 4b |I%,JJ%) | : IK(KK|
|00004fc0| 25 29 3d 49 49 25 20 20 | 3a 20 4a 4b 28 4b 4b 25 |%)=II% |: JK(KK%|
|00004fd0| 29 3d 4a 4a 25 0d 0a 32 | 39 31 30 20 4e 45 58 54 |)=JJ%..2|910 NEXT|
|00004fe0| 20 4a 4a 25 2c 49 49 25 | 0d 0a 32 39 32 30 20 49 | JJ%,II%|..2920 I|
|00004ff0| 46 20 41 4d 58 23 3c 3e | 30 21 20 47 4f 54 4f 20 |F AMX#<>|0! GOTO |
|00005000| 32 39 35 30 0d 0a 32 39 | 33 30 20 27 09 2a 2a 2a |2950..29|30 '.***|
|00005010| 20 49 4e 54 45 52 43 48 | 41 4e 47 45 20 52 4f 57 | INTERCH|ANGE ROW|
|00005020| 53 20 26 20 43 4f 4c 55 | 4d 4e 53 0d 0a 32 39 34 |S & COLU|MNS..294|
|00005030| 30 20 44 45 54 3d 30 21 | 20 20 3a 20 52 45 54 55 |0 DET=0!| : RETU|
|00005040| 52 4e 0d 0a 32 39 35 30 | 20 49 49 25 3d 49 4b 28 |RN..2950| II%=IK(|
|00005050| 4b 4b 25 29 0d 0a 32 39 | 36 30 20 49 46 20 49 49 |KK%)..29|60 IF II|
|00005060| 25 3c 4b 4b 25 20 47 4f | 54 4f 20 32 38 37 30 0d |%<KK% GO|TO 2870.|
|00005070| 0a 32 39 37 30 20 49 46 | 20 49 49 25 3d 4b 4b 25 |.2970 IF| II%=KK%|
|00005080| 20 47 4f 54 4f 20 33 30 | 31 30 0d 0a 32 39 38 30 | GOTO 30|10..2980|
|00005090| 20 46 4f 52 20 4a 4a 25 | 3d 31 20 54 4f 20 4e 4f | FOR JJ%|=1 TO NO|
|000050a0| 52 0d 0a 32 39 39 30 20 | 53 41 56 23 3d 41 52 59 |R..2990 |SAV#=ARY|
|000050b0| 23 28 4b 4b 25 2c 4a 4a | 25 29 20 20 3a 20 41 52 |#(KK%,JJ|%) : AR|
|000050c0| 59 23 28 4b 4b 25 2c 4a | 4a 25 29 3d 41 52 59 23 |Y#(KK%,J|J%)=ARY#|
|000050d0| 28 49 49 25 2c 4a 4a 25 | 29 20 20 3a 20 41 52 59 |(II%,JJ%|) : ARY|
|000050e0| 23 28 49 49 25 2c 4a 4a | 25 29 3d 2d 53 41 56 23 |#(II%,JJ|%)=-SAV#|
|000050f0| 0d 0a 33 30 30 30 20 4e | 45 58 54 20 4a 4a 25 0d |..3000 N|EXT JJ%.|
|00005100| 0a 33 30 31 30 20 4a 4a | 25 3d 4a 4b 28 4b 4b 25 |.3010 JJ|%=JK(KK%|
|00005110| 29 0d 0a 33 30 32 30 20 | 49 46 20 4a 4a 25 3c 4b |)..3020 |IF JJ%<K|
|00005120| 4b 25 20 47 4f 54 4f 20 | 32 38 37 30 0d 0a 33 30 |K% GOTO |2870..30|
|00005130| 33 30 20 49 46 20 4a 4a | 25 3d 4b 4b 25 20 47 4f |30 IF JJ|%=KK% GO|
|00005140| 54 4f 20 33 30 38 30 0d | 0a 33 30 34 30 20 46 4f |TO 3080.|.3040 FO|
|00005150| 52 20 49 49 25 3d 31 20 | 54 4f 20 4e 4f 52 0d 0a |R II%=1 |TO NOR..|
|00005160| 33 30 35 30 20 53 41 56 | 23 3d 41 52 59 23 28 49 |3050 SAV|#=ARY#(I|
|00005170| 49 25 2c 4b 4b 25 29 20 | 20 3a 20 41 52 59 23 28 |I%,KK%) | : ARY#(|
|00005180| 49 49 25 2c 4b 4b 25 29 | 3d 41 52 59 23 28 49 49 |II%,KK%)|=ARY#(II|
|00005190| 25 2c 4a 4a 25 29 20 20 | 3a 20 41 52 59 23 28 49 |%,JJ%) |: ARY#(I|
|000051a0| 49 25 2c 4a 4a 25 29 3d | 2d 53 41 56 23 0d 0a 33 |I%,JJ%)=|-SAV#..3|
|000051b0| 30 36 30 20 4e 45 58 54 | 20 49 49 25 0d 0a 33 30 |060 NEXT| II%..30|
|000051c0| 37 30 20 27 09 2a 2a 2a | 20 41 43 43 55 4d 55 4c |70 '.***| ACCUMUL|
|000051d0| 41 54 45 20 45 4c 45 4d | 45 4e 54 53 20 4f 46 20 |ATE ELEM|ENTS OF |
|000051e0| 49 4e 56 45 52 53 45 20 | 4d 41 54 52 49 58 20 2a |INVERSE |MATRIX *|
|000051f0| 2a 2a 0d 0a 33 30 38 30 | 20 46 4f 52 20 49 49 25 |**..3080| FOR II%|
|00005200| 3d 31 20 54 4f 20 4e 4f | 52 0d 0a 33 30 39 30 20 |=1 TO NO|R..3090 |
|00005210| 20 49 46 20 49 49 25 3d | 4b 4b 25 20 47 4f 54 4f | IF II%=|KK% GOTO|
|00005220| 20 33 31 31 30 0d 0a 33 | 31 30 30 20 20 41 52 59 | 3110..3|100 ARY|
|00005230| 23 28 49 49 25 2c 4b 4b | 25 29 3d 2d 41 52 59 23 |#(II%,KK|%)=-ARY#|
|00005240| 28 49 49 25 2c 4b 4b 25 | 29 2f 41 4d 58 23 0d 0a |(II%,KK%|)/AMX#..|
|00005250| 33 31 31 30 20 4e 45 58 | 54 20 49 49 25 0d 0a 33 |3110 NEX|T II%..3|
|00005260| 31 32 30 20 46 4f 52 20 | 49 49 25 3d 31 20 54 4f |120 FOR |II%=1 TO|
|00005270| 20 4e 4f 52 0d 0a 33 31 | 33 30 20 46 4f 52 20 4a | NOR..31|30 FOR J|
|00005280| 4a 25 3d 31 20 54 4f 20 | 4e 4f 52 0d 0a 33 31 34 |J%=1 TO |NOR..314|
|00005290| 30 20 20 20 49 46 20 49 | 49 25 3d 4b 4b 25 20 47 |0 IF I|I%=KK% G|
|000052a0| 4f 54 4f 20 33 31 37 30 | 0d 0a 33 31 35 30 20 20 |OTO 3170|..3150 |
|000052b0| 20 49 46 20 4a 4a 25 3d | 4b 4b 25 20 47 4f 54 4f | IF JJ%=|KK% GOTO|
|000052c0| 20 33 31 37 30 0d 0a 33 | 31 36 30 20 20 20 41 52 | 3170..3|160 AR|
|000052d0| 59 23 28 49 49 25 2c 4a | 4a 25 29 3d 41 52 59 23 |Y#(II%,J|J%)=ARY#|
|000052e0| 28 49 49 25 2c 4a 4a 25 | 29 2b 41 52 59 23 28 49 |(II%,JJ%|)+ARY#(I|
|000052f0| 49 25 2c 4b 4b 25 29 2a | 41 52 59 23 28 4b 4b 25 |I%,KK%)*|ARY#(KK%|
|00005300| 2c 4a 4a 25 29 0d 0a 33 | 31 37 30 20 4e 45 58 54 |,JJ%)..3|170 NEXT|
|00005310| 20 4a 4a 25 2c 49 49 25 | 0d 0a 33 31 38 30 20 46 | JJ%,II%|..3180 F|
|00005320| 4f 52 20 4a 4a 25 3d 31 | 20 54 4f 20 4e 4f 52 0d |OR JJ%=1| TO NOR.|
|00005330| 0a 33 31 39 30 20 20 49 | 46 20 4a 4a 25 3d 4b 4b |.3190 I|F JJ%=KK|
|00005340| 25 20 47 4f 54 4f 20 33 | 32 31 30 0d 0a 33 32 30 |% GOTO 3|210..320|
|00005350| 30 20 20 41 52 59 23 28 | 4b 4b 25 2c 4a 4a 25 29 |0 ARY#(|KK%,JJ%)|
|00005360| 3d 41 52 59 23 28 4b 4b | 25 2c 4a 4a 25 29 2f 41 |=ARY#(KK|%,JJ%)/A|
|00005370| 4d 58 23 0d 0a 33 32 31 | 30 20 4e 45 58 54 20 4a |MX#..321|0 NEXT J|
|00005380| 4a 25 0d 0a 33 32 32 30 | 20 20 41 52 59 23 28 4b |J%..3220| ARY#(K|
|00005390| 4b 25 2c 4b 4b 25 29 3d | 31 21 2f 41 4d 58 23 20 |K%,KK%)=|1!/AMX# |
|000053a0| 20 3a 20 44 45 54 3d 44 | 45 54 2a 41 4d 58 23 0d | : DET=D|ET*AMX#.|
|000053b0| 0a 33 32 33 30 20 4e 45 | 58 54 20 4b 4b 25 0d 0a |.3230 NE|XT KK%..|
|000053c0| 33 32 34 30 20 27 09 2a | 2a 2a 20 52 45 53 54 4f |3240 '.*|** RESTO|
|000053d0| 52 45 20 4f 52 44 45 52 | 49 4e 47 20 4f 46 20 4d |RE ORDER|ING OF M|
|000053e0| 41 54 52 49 58 0d 0a 33 | 32 35 30 20 46 4f 52 20 |ATRIX..3|250 FOR |
|000053f0| 4c 25 3d 31 20 54 4f 20 | 4e 4f 52 0d 0a 33 32 36 |L%=1 TO |NOR..326|
|00005400| 30 20 20 20 4b 4b 25 3d | 4e 4f 52 2d 4c 25 20 2b |0 KK%=|NOR-L% +|
|00005410| 31 20 20 3a 20 4a 4a 25 | 3d 49 4b 28 4b 4b 25 29 |1 : JJ%|=IK(KK%)|
|00005420| 0d 0a 33 32 37 30 20 20 | 20 49 46 20 4a 4a 25 3c |..3270 | IF JJ%<|
|00005430| 3d 4b 4b 25 20 47 4f 54 | 4f 20 33 33 31 30 0d 0a |=KK% GOT|O 3310..|
|00005440| 33 32 38 30 20 46 4f 52 | 20 49 49 25 3d 31 20 54 |3280 FOR| II%=1 T|
|00005450| 4f 20 4e 4f 52 0d 0a 33 | 32 39 30 20 20 53 41 56 |O NOR..3|290 SAV|
|00005460| 23 3d 41 52 59 23 28 49 | 49 25 2c 4b 4b 25 29 20 |#=ARY#(I|I%,KK%) |
|00005470| 20 3a 20 41 52 59 23 28 | 49 49 25 2c 4b 4b 25 29 | : ARY#(|II%,KK%)|
|00005480| 3d 2d 41 52 59 23 28 49 | 49 25 2c 4a 4a 25 29 20 |=-ARY#(I|I%,JJ%) |
|00005490| 20 3a 20 41 52 59 23 28 | 49 49 25 2c 4a 4a 25 29 | : ARY#(|II%,JJ%)|
|000054a0| 3d 53 41 56 23 0d 0a 33 | 33 30 30 20 4e 45 58 54 |=SAV#..3|300 NEXT|
|000054b0| 20 49 49 25 0d 0a 33 33 | 31 30 20 49 49 25 3d 4a | II%..33|10 II%=J|
|000054c0| 4b 28 4b 4b 25 29 0d 0a | 33 33 32 30 20 49 46 20 |K(KK%)..|3320 IF |
|000054d0| 49 49 25 3c 3d 4b 4b 25 | 20 47 4f 54 4f 20 33 33 |II%<=KK%| GOTO 33|
|000054e0| 36 30 0d 0a 33 33 33 30 | 20 46 4f 52 20 4a 4a 25 |60..3330| FOR JJ%|
|000054f0| 3d 31 20 54 4f 20 4e 4f | 52 0d 0a 33 33 34 30 20 |=1 TO NO|R..3340 |
|00005500| 20 53 41 56 23 3d 41 52 | 59 23 28 4b 4b 25 2c 4a | SAV#=AR|Y#(KK%,J|
|00005510| 4a 25 29 20 20 3a 20 41 | 52 59 23 28 4b 4b 25 2c |J%) : A|RY#(KK%,|
|00005520| 4a 4a 25 29 3d 2d 41 52 | 59 23 28 49 49 25 2c 4a |JJ%)=-AR|Y#(II%,J|
|00005530| 4a 25 29 20 20 3a 20 41 | 52 59 23 28 49 49 25 2c |J%) : A|RY#(II%,|
|00005540| 4a 4a 25 29 3d 53 41 56 | 23 0d 0a 33 33 35 30 20 |JJ%)=SAV|#..3350 |
|00005550| 4e 45 58 54 20 4a 4a 25 | 0d 0a 33 33 36 30 20 4e |NEXT JJ%|..3360 N|
|00005560| 45 58 54 20 4c 25 0d 0a | 33 33 37 30 20 52 45 54 |EXT L%..|3370 RET|
|00005570| 55 52 4e 0d 0a 33 33 38 | 30 20 52 45 4d 2d 2d 2d |URN..338|0 REM---|
|00005580| 2d 2d 2d 2d 2d 2d 2d 2d | 2d 2d 2d 2d 2d 2d 2d 2d |--------|--------|
|00005590| 2d 2d 2d 2d 2d 2d 2d 2d | 2d 2d 2d 2d 2d 2d 2d 2d |--------|--------|
|000055a0| 2d 2d 2d 2d 2d 2d 2d 2d | 2d 2d 2d 2d 2d 2d 2d 2d |--------|--------|
|000055b0| 2d 2d 2d 2d 2d 2d 2d 2d | 2d 2d 2d 2d 2d 2d 2d 2d |--------|--------|
|000055c0| 2d 2d 2d 0d 0a 33 33 39 | 30 20 27 09 53 55 42 52 |---..339|0 '.SUBR|
|000055d0| 4f 55 54 49 4e 45 20 46 | 43 54 4e 28 58 2c 20 49 |OUTINE F|CTN(X, I|
|000055e0| 2c 20 4a 2c 20 4a 54 52 | 4d 29 0d 0a 33 34 30 30 |, J, JTR|M)..3400|
|000055f0| 20 27 20 45 56 41 4c 55 | 41 54 45 20 54 45 52 4d | ' EVALU|ATE TERM|
|00005600| 53 20 4f 46 20 41 20 50 | 4f 4c 59 4e 4f 4d 49 41 |S OF A P|OLYNOMIA|
|00005610| 4c 20 46 55 4e 43 54 49 | 4f 4e 20 46 4f 52 20 52 |L FUNCTI|ON FOR R|
|00005620| 45 47 52 45 53 0d 0a 33 | 34 31 30 20 27 20 20 58 |EGRES..3|410 ' X|
|00005630| 3d 41 52 52 41 59 20 4f | 46 20 44 41 54 41 20 50 |=ARRAY O|F DATA P|
|00005640| 4f 49 4e 54 53 2c 20 20 | 49 3d 49 4e 44 45 58 20 |OINTS, |I=INDEX |
|00005650| 4f 46 20 44 41 54 41 20 | 50 4f 49 4e 54 53 0d 0a |OF DATA |POINTS..|
|00005660| 33 34 32 30 20 27 20 20 | 4a 3d 49 4e 44 45 58 20 |3420 ' |J=INDEX |
|00005670| 4f 46 20 54 45 52 4d 20 | 49 4e 20 50 4f 4c 59 4e |OF TERM |IN POLYN|
|00005680| 4f 4d 49 41 4c 2c 20 20 | 4a 54 52 4d 3d 41 52 52 |OMIAL, |JTRM=ARR|
|00005690| 41 59 20 4f 46 20 50 4f | 57 45 52 53 0d 0a 33 34 |AY OF PO|WERS..34|
|000056a0| 33 30 20 27 0d 0a 33 34 | 34 30 20 4a 45 58 50 3d |30 '..34|40 JEXP=|
|000056b0| 4a 54 52 4d 28 4a 29 20 | 20 3a 20 46 43 54 4e 3d |JTRM(J) | : FCTN=|
|000056c0| 53 49 4e 28 58 28 49 29 | 29 20 20 3a 20 52 45 54 |SIN(X(I)|) : RET|
|000056d0| 55 52 4e 0d 0a 1a 1a 1a | 1a 1a 1a 1a 1a 1a 1a 1a |URN.....|........|
|000056e0| 1a 1a 1a 1a 1a 1a 1a 1a | 1a 1a 1a 1a 1a 1a 1a 1a |........|........|
|000056f0| 1a 1a 1a 1a 1a 1a 1a 1a | 1a 1a 1a 1a 1a 1a 1a 1a |........|........|
|00005700| 31 30 30 30 20 27 20 20 | 20 20 20 20 20 20 20 20 |1000 ' | |
|00005710| 20 20 20 20 20 20 2a 2a | 2a 20 20 20 50 52 4f 47 | **|* PROG|
|00005720| 52 41 4d 20 50 4f 4c 59 | 46 49 54 20 20 20 2a 2a |RAM POLY|FIT **|
|00005730| 2a 0d 0a 31 30 31 30 20 | 27 0d 0a 31 30 32 30 20 |*..1010 |'..1020 |
|00005740| 27 20 20 41 20 44 52 49 | 56 45 52 20 50 52 4f 47 |' A DRI|VER PROG|
|00005750| 52 41 4d 20 46 4f 52 20 | 54 48 45 20 43 55 52 46 |RAM FOR |THE CURF|
|00005760| 49 54 20 53 55 42 52 4f | 55 54 49 4e 45 53 20 54 |IT SUBRO|UTINES T|
|00005770| 48 41 54 20 57 49 4c 4c | 20 46 49 4e 44 0d 0a 31 |HAT WILL| FIND..1|
|00005780| 30 33 30 20 27 20 54 48 | 45 20 50 41 52 41 4d 45 |030 ' TH|E PARAME|
|00005790| 54 45 52 53 20 49 4e 20 | 41 20 47 49 56 45 4e 20 |TERS IN |A GIVEN |
|000057a0| 46 55 4e 43 54 49 4f 4e | 20 54 48 41 54 20 47 49 |FUNCTION| THAT GI|
|000057b0| 56 45 20 54 48 45 20 42 | 45 53 54 20 46 49 54 0d |VE THE B|EST FIT.|
|000057c0| 0a 31 30 34 30 20 27 20 | 4f 46 20 41 20 50 4f 4c |.1040 ' |OF A POL|
|000057d0| 59 4e 4f 4d 49 41 4c 20 | 54 4f 20 41 20 53 45 54 |YNOMIAL |TO A SET|
|000057e0| 20 4f 46 20 49 4e 50 55 | 54 20 44 41 54 41 2e 20 | OF INPU|T DATA. |
|000057f0| 20 58 20 26 20 59 20 41 | 52 45 0d 0a 31 30 35 30 | X & Y A|RE..1050|
|00005800| 20 27 20 54 48 45 20 49 | 4e 44 45 50 45 4e 44 45 | ' THE I|NDEPENDE|
|00005810| 4e 54 20 41 4e 44 20 44 | 45 50 45 4e 44 45 4e 54 |NT AND D|EPENDENT|
|00005820| 20 56 41 52 49 41 42 4c | 45 53 20 52 45 53 50 45 | VARIABL|ES RESPE|
|00005830| 43 54 49 56 45 4c 59 2e | 20 20 0d 0a 31 30 36 30 |CTIVELY.| ..1060|
|00005840| 20 27 0d 0a 31 30 37 30 | 20 27 20 20 50 52 4f 47 | '..1070| ' PROG|
|00005850| 52 41 4d 4d 45 52 3a 20 | 20 47 2e 20 4d 2e 20 52 |RAMMER: | G. M. R|
|00005860| 45 53 43 48 2c 20 20 46 | 45 42 2e 20 31 39 37 39 |ESCH, F|EB. 1979|
|00005870| 0d 0a 31 30 38 30 20 27 | 0d 0a 31 30 39 30 20 27 |..1080 '|..1090 '|
|00005880| 09 44 49 4d 45 4e 53 49 | 4f 4e 20 41 4c 4c 20 41 |.DIMENSI|ON ALL A|
|00005890| 52 52 41 59 53 20 48 45 | 52 45 0d 0a 31 31 30 30 |RRAYS HE|RE..1100|
|000058a0| 20 27 20 54 48 49 53 20 | 56 45 52 53 49 4f 4e 20 | ' THIS |VERSION |
|000058b0| 49 53 20 53 45 54 20 55 | 50 20 54 4f 20 48 41 4e |IS SET U|P TO HAN|
|000058c0| 44 4c 45 20 55 50 20 54 | 4f 20 31 30 30 20 44 41 |DLE UP T|O 100 DA|
|000058d0| 54 41 20 50 4f 49 4e 54 | 53 20 41 4e 44 0d 0a 31 |TA POINT|S AND..1|
|000058e0| 31 31 30 20 27 20 41 53 | 20 4d 41 4e 59 20 41 53 |110 ' AS| MANY AS|
|000058f0| 20 31 30 20 50 41 52 41 | 4d 45 54 45 52 53 2e 0d | 10 PARA|METERS..|
|00005900| 0a 31 31 32 30 20 44 49 | 4d 20 20 58 28 31 30 30 |.1120 DI|M X(100|
|00005910| 29 2c 20 59 28 31 30 30 | 29 2c 20 53 49 47 59 28 |), Y(100|), SIGY(|
|00005920| 31 30 30 29 2c 20 59 46 | 49 54 28 31 30 30 29 2c |100), YF|IT(100),|
|00005930| 20 57 47 48 54 28 31 30 | 30 29 0d 0a 31 31 33 30 | WGHT(10|0)..1130|
|00005940| 20 44 49 4d 20 20 41 28 | 31 30 29 2c 20 53 55 4d | DIM A(|10), SUM|
|00005950| 58 23 28 31 30 29 2c 20 | 53 55 4d 59 23 28 31 30 |X#(10), |SUMY#(10|
|00005960| 29 2c 20 41 52 52 41 59 | 23 28 31 30 2c 31 30 29 |), ARRAY|#(10,10)|
|00005970| 0d 0a 31 31 34 30 20 27 | 0d 0a 31 31 35 30 20 4d |..1140 '|..1150 M|
|00005980| 41 58 4e 3d 20 31 30 30 | 20 20 3a 20 4d 41 58 50 |AXN= 100| : MAXP|
|00005990| 3d 31 30 09 27 53 45 54 | 20 4d 41 58 49 4d 55 4d |=10.'SET| MAXIMUM|
|000059a0| 20 23 20 4f 46 20 50 4f | 49 4e 54 20 41 4e 44 20 | # OF PO|INT AND |
|000059b0| 50 41 52 4d 53 0d 0a 31 | 31 36 30 20 50 52 49 4e |PARMS..1|160 PRIN|
|000059c0| 54 20 22 20 20 20 2a 2a | 2a 20 50 72 6f 67 72 61 |T " **|* Progra|
|000059d0| 6d 20 50 4f 4c 59 46 49 | 54 20 2d 20 47 2e 20 4d |m POLYFI|T - G. M|
|000059e0| 2e 20 52 65 73 63 68 2c | 20 46 65 62 2e 20 31 39 |. Resch,| Feb. 19|
|000059f0| 37 39 20 2a 2a 2a 22 0d | 0a 31 31 37 30 20 50 52 |79 ***".|.1170 PR|
|00005a00| 49 4e 54 0d 0a 31 31 38 | 30 20 27 0d 0a 31 31 39 |INT..118|0 '..119|
|00005a10| 30 20 27 09 20 4f 50 45 | 4e 20 26 20 52 45 41 44 |0 '. OPE|N & READ|
|00005a20| 20 54 48 45 20 49 4e 50 | 55 54 20 44 41 54 41 20 | THE INP|UT DATA |
|00005a30| 53 45 54 0d 0a 31 32 30 | 30 20 49 4e 50 55 54 20 |SET..120|0 INPUT |
|00005a40| 22 4e 61 6d 65 20 6f 66 | 20 49 6e 70 75 74 20 64 |"Name of| Input d|
|00005a50| 61 74 61 20 73 65 74 3d | 22 3b 4e 24 0d 0a 31 32 |ata set=|";N$..12|
|00005a60| 31 30 20 4f 50 45 4e 22 | 49 22 2c 32 2c 4e 24 0d |10 OPEN"|I",2,N$.|
|00005a70| 0a 31 32 32 30 20 46 4f | 52 20 49 3d 31 20 54 4f |.1220 FO|R I=1 TO|
|00005a80| 20 4d 41 58 4e 0d 0a 31 | 32 33 30 20 20 20 49 46 | MAXN..1|230 IF|
|00005a90| 20 45 4f 46 28 32 29 20 | 20 47 4f 54 4f 20 31 32 | EOF(2) | GOTO 12|
|00005aa0| 36 30 0d 0a 31 32 34 30 | 20 20 20 49 4e 50 55 54 |60..1240| INPUT|
|00005ab0| 23 32 2c 58 28 49 29 2c | 59 28 49 29 2c 53 49 47 |#2,X(I),|Y(I),SIG|
|00005ac0| 59 28 49 29 0d 0a 31 32 | 35 30 20 4e 45 58 54 20 |Y(I)..12|50 NEXT |
|00005ad0| 49 0d 0a 31 32 36 30 20 | 4e 50 54 53 3d 49 2d 31 |I..1260 |NPTS=I-1|
|00005ae0| 20 20 3a 20 43 4c 4f 53 | 45 23 32 0d 0a 31 32 37 | : CLOS|E#2..127|
|00005af0| 30 20 50 52 49 4e 54 22 | 54 48 45 52 45 20 41 52 |0 PRINT"|THERE AR|
|00005b00| 45 20 22 3b 4e 50 54 53 | 3b 22 20 44 41 54 41 20 |E ";NPTS|;" DATA |
|00005b10| 50 4f 49 4e 54 53 22 0d | 0a 31 32 38 30 20 27 0d |POINTS".|.1280 '.|
|00005b20| 0a 31 32 39 30 20 27 09 | 20 49 4e 50 55 54 20 54 |.1290 '.| INPUT T|
|00005b30| 48 45 20 4d 4f 44 45 20 | 41 4e 44 20 4e 55 4d 42 |HE MODE |AND NUMB|
|00005b40| 45 52 20 4f 46 20 54 45 | 52 4d 53 2e 0d 0a 31 33 |ER OF TE|RMS...13|
|00005b50| 30 30 20 50 52 49 4e 54 | 22 20 20 4d 4f 44 45 3d |00 PRINT|" MODE=|
|00005b60| 2b 31 20 49 4e 53 54 52 | 2e 20 57 47 48 54 28 49 |+1 INSTR|. WGHT(I|
|00005b70| 29 3d 31 2f 53 49 47 59 | 28 49 29 5e 32 22 0d 0a |)=1/SIGY|(I)^2"..|
|00005b80| 31 33 31 30 20 50 52 49 | 4e 54 22 20 20 20 20 20 |1310 PRI|NT" |
|00005b90| 20 3d 20 30 20 4e 4f 20 | 20 20 20 20 57 47 48 54 | = 0 NO | WGHT|
|00005ba0| 28 49 29 3d 31 22 0d 0a | 31 33 32 30 20 50 52 49 |(I)=1"..|1320 PRI|
|00005bb0| 4e 54 22 20 20 20 20 20 | 20 3d 2d 31 20 53 54 41 |NT" | =-1 STA|
|00005bc0| 54 2e 20 20 57 47 48 54 | 28 49 29 3d 31 2f 59 28 |T. WGHT|(I)=1/Y(|
|00005bd0| 49 29 22 0d 0a 31 33 33 | 30 20 50 52 49 4e 54 22 |I)"..133|0 PRINT"|
|00005be0| 20 20 20 20 20 20 3d 2d | 32 20 41 52 42 49 54 2e | =-|2 ARBIT.|
|00005bf0| 20 57 47 48 54 28 49 29 | 3d 53 49 47 59 28 49 29 | WGHT(I)|=SIGY(I)|
|00005c00| 22 0d 0a 31 33 34 30 20 | 49 4e 50 55 54 22 53 50 |"..1340 |INPUT"SP|
|00005c10| 45 43 49 46 59 20 54 48 | 45 20 57 45 49 47 48 54 |ECIFY TH|E WEIGHT|
|00005c20| 49 4e 47 20 4d 4f 44 45 | 20 28 2b 31 2c 30 2c 2d |ING MODE| (+1,0,-|
|00005c30| 31 2c 20 4f 52 20 2d 32 | 29 20 2d 22 3b 4d 4f 44 |1, OR -2|) -";MOD|
|00005c40| 45 0d 0a 31 33 35 30 20 | 49 4e 50 55 54 22 48 4f |E..1350 |INPUT"HO|
|00005c50| 57 20 4d 41 4e 59 20 50 | 41 52 41 4d 45 54 45 52 |W MANY P|ARAMETER|
|00005c60| 53 3f 20 4e 54 52 4d 3d | 22 3b 4e 54 52 4d 0d 0a |S? NTRM=|";NTRM..|
|00005c70| 31 33 36 30 20 49 46 20 | 4e 54 52 4d 3c 3d 4e 50 |1360 IF |NTRM<=NP|
|00005c80| 54 53 20 20 47 4f 54 4f | 20 31 34 30 30 0d 0a 31 |TS GOTO| 1400..1|
|00005c90| 33 37 30 20 50 52 49 4e | 54 22 2a 2a 2a 20 59 4f |370 PRIN|T"*** YO|
|00005ca0| 55 20 4d 55 53 54 20 48 | 41 56 45 20 4d 4f 52 45 |U MUST H|AVE MORE|
|00005cb0| 20 44 41 54 41 20 54 48 | 41 4e 20 53 4f 4c 56 45 | DATA TH|AN SOLVE|
|00005cc0| 2d 46 4f 52 20 50 41 52 | 41 4d 45 54 45 52 53 21 |-FOR PAR|AMETERS!|
|00005cd0| 22 0d 0a 31 33 38 30 20 | 47 4f 54 4f 20 31 33 35 |"..1380 |GOTO 135|
|00005ce0| 30 0d 0a 31 33 39 30 20 | 27 0d 0a 31 34 30 30 20 |0..1390 |'..1400 |
|00005cf0| 20 47 4f 53 55 42 20 31 | 37 36 30 20 20 09 27 43 | GOSUB 1|760 .'C|
|00005d00| 41 4c 4c 20 50 4f 4c 59 | 46 49 54 0d 0a 31 34 31 |ALL POLY|FIT..141|
|00005d10| 30 20 27 0d 0a 31 34 32 | 30 20 27 09 53 55 4d 4d |0 '..142|0 '.SUMM|
|00005d20| 41 52 49 5a 45 20 54 48 | 45 20 52 45 53 55 4c 54 |ARIZE TH|E RESULT|
|00005d30| 53 0d 0a 31 34 33 30 20 | 50 52 49 4e 54 0d 0a 31 |S..1430 |PRINT..1|
|00005d40| 34 34 30 20 46 4f 52 20 | 4a 3d 31 20 54 4f 20 4e |440 FOR |J=1 TO N|
|00005d50| 54 52 4d 0d 0a 31 34 35 | 30 20 20 50 52 49 4e 54 |TRM..145|0 PRINT|
|00005d60| 22 41 28 22 3b 4a 3b 22 | 29 3d 22 3b 41 28 4a 29 |"A(";J;"|)=";A(J)|
|00005d70| 0d 0a 31 34 36 30 20 4e | 45 58 54 20 4a 0d 0a 31 |..1460 N|EXT J..1|
|00005d80| 34 37 30 20 50 52 49 4e | 54 22 43 48 49 20 53 51 |470 PRIN|T"CHI SQ|
|00005d90| 55 41 52 45 3d 22 3b 43 | 48 49 53 51 52 0d 0a 31 |UARE=";C|HISQR..1|
|00005da0| 34 38 30 20 50 52 49 4e | 54 0d 0a 31 34 39 30 20 |480 PRIN|T..1490 |
|00005db0| 50 52 49 4e 54 22 58 28 | 49 29 22 2c 22 59 28 49 |PRINT"X(|I)","Y(I|
|00005dc0| 29 22 2c 22 59 46 49 54 | 22 2c 22 59 2d 59 46 49 |)","YFIT|","Y-YFI|
|00005dd0| 54 22 0d 0a 31 35 30 30 | 20 46 4f 52 20 49 3d 31 |T"..1500| FOR I=1|
|00005de0| 20 54 4f 20 4e 50 54 53 | 0d 0a 31 35 31 30 20 20 | TO NPTS|..1510 |
|00005df0| 59 46 49 54 28 49 29 3d | 20 41 28 31 29 0d 0a 31 |YFIT(I)=| A(1)..1|
|00005e00| 35 32 30 20 20 20 46 4f | 52 20 4a 3d 32 20 54 4f |520 FO|R J=2 TO|
|00005e10| 20 4e 54 52 4d 20 20 3a | 20 59 46 49 54 28 49 29 | NTRM :| YFIT(I)|
|00005e20| 3d 59 46 49 54 28 49 29 | 20 2b 41 28 4a 29 2a 28 |=YFIT(I)| +A(J)*(|
|00005e30| 58 28 49 29 5e 28 4a 2d | 31 29 29 20 3a 20 4e 45 |X(I)^(J-|1)) : NE|
|00005e40| 58 54 20 4a 0d 0a 31 35 | 33 30 20 20 50 52 49 4e |XT J..15|30 PRIN|
|00005e50| 54 20 58 28 49 29 2c 59 | 28 49 29 2c 59 46 49 54 |T X(I),Y|(I),YFIT|
|00005e60| 28 49 29 2c 28 59 28 49 | 29 2d 59 46 49 54 28 49 |(I),(Y(I|)-YFIT(I|
|00005e70| 29 29 0d 0a 31 35 34 30 | 20 4e 45 58 54 20 49 0d |))..1540| NEXT I.|
|00005e80| 0a 31 35 35 30 20 45 4e | 44 0d 0a 31 35 36 30 20 |.1550 EN|D..1560 |
|00005e90| 52 45 4d 2d 2d 2d 2d 2d | 2d 2d 2d 2d 2d 2d 2d 2d |REM-----|--------|
|00005ea0| 2d 2d 2d 2d 2d 2d 2d 2d | 2d 2d 2d 2d 2d 2d 2d 2d |--------|--------|
|00005eb0| 2d 2d 2d 2d 2d 2d 2d 2d | 2d 2d 2d 2d 2d 2d 2d 2d |--------|--------|
|00005ec0| 2d 2d 2d 2d 2d 2d 2d 2d | 2d 2d 2d 2d 2d 2d 2d 2d |--------|--------|
|00005ed0| 2d 2d 2d 2d 2d 2d 2d 2d | 2d 0d 0a 31 35 37 30 20 |--------|-..1570 |
|00005ee0| 27 2a 2a 2a 20 20 53 55 | 42 52 4f 55 54 49 4e 45 |'*** SU|BROUTINE|
|00005ef0| 20 50 4f 4c 59 46 49 54 | 28 58 2c 20 59 2c 20 53 | POLYFIT|(X, Y, S|
|00005f00| 49 47 59 2c 20 4e 50 54 | 53 2c 20 4e 54 52 4d 2c |IGY, NPT|S, NTRM,|
|00005f10| 20 4d 4f 44 45 2c 20 41 | 2c 20 43 48 49 53 51 52 | MODE, A|, CHISQR|
|00005f20| 29 0d 0a 31 35 38 30 20 | 27 0d 0a 31 35 39 30 20 |)..1580 |'..1590 |
|00005f30| 27 20 20 50 55 52 50 4f | 53 45 2d 20 4d 41 4b 45 |' PURPO|SE- MAKE|
|00005f40| 20 41 20 4c 45 41 53 54 | 20 53 51 55 41 52 45 53 | A LEAST| SQUARES|
|00005f50| 20 46 49 54 20 54 4f 20 | 44 41 54 41 20 57 49 54 | FIT TO |DATA WIT|
|00005f60| 48 20 41 20 0d 0a 31 36 | 30 30 20 27 20 20 50 4f |H A ..16|00 ' PO|
|00005f70| 4c 59 4e 4f 4d 49 41 4c | 20 46 55 4e 43 54 49 4f |LYNOMIAL| FUNCTIO|
|00005f80| 4e 2e 0d 0a 31 36 31 30 | 20 27 0d 0a 31 36 32 30 |N...1610| '..1620|
|00005f90| 20 27 20 20 58 3d 49 4e | 44 45 50 45 4e 44 45 4e | ' X=IN|DEPENDEN|
|00005fa0| 54 20 56 41 52 49 41 42 | 4c 45 2c 20 20 59 3d 44 |T VARIAB|LE, Y=D|
|00005fb0| 45 50 45 4e 44 45 4e 54 | 20 56 41 52 49 41 42 4c |EPENDENT| VARIABL|
|00005fc0| 45 0d 0a 31 36 33 30 20 | 27 20 20 53 49 47 59 3d |E..1630 |' SIGY=|
|00005fd0| 53 49 47 4d 41 53 20 4f | 46 20 59 2c 20 20 4e 50 |SIGMAS O|F Y, NP|
|00005fe0| 54 53 3d 23 20 4f 46 20 | 44 41 54 41 20 50 4f 49 |TS=# OF |DATA POI|
|00005ff0| 4e 54 53 0d 0a 31 36 34 | 30 20 27 20 20 4e 54 52 |NTS..164|0 ' NTR|
|00006000| 4d 3d 23 20 4f 46 20 43 | 4f 45 46 46 2e 2c 20 28 |M=# OF C|OEFF., (|
|00006010| 44 45 47 52 45 45 20 4f | 46 20 50 4f 4c 59 4e 4f |DEGREE O|F POLYNO|
|00006020| 4d 49 41 4c 20 2b 20 31 | 29 0d 0a 31 36 35 30 20 |MIAL + 1|)..1650 |
|00006030| 27 20 20 4d 4f 44 45 3d | 2b 31 20 49 4e 53 54 52 |' MODE=|+1 INSTR|
|00006040| 2e 20 57 47 48 54 28 49 | 29 3d 31 2f 53 49 47 59 |. WGHT(I|)=1/SIGY|
|00006050| 28 49 29 5e 32 0d 0a 31 | 36 36 30 20 27 20 20 20 |(I)^2..1|660 ' |
|00006060| 20 20 20 3d 20 30 20 4e | 4f 20 20 20 20 20 57 47 | = 0 N|O WG|
|00006070| 48 54 28 49 29 3d 31 0d | 0a 31 36 37 30 20 27 20 |HT(I)=1.|.1670 ' |
|00006080| 20 20 20 20 20 3d 2d 31 | 20 53 54 41 54 2e 20 20 | =-1| STAT. |
|00006090| 57 47 48 54 28 49 29 3d | 31 2f 59 28 49 29 0d 0a |WGHT(I)=|1/Y(I)..|
|000060a0| 31 36 37 35 20 27 20 20 | 20 20 20 20 3d 2d 32 20 |1675 ' | =-2 |
|000060b0| 41 52 42 49 54 2e 20 57 | 47 48 54 28 49 29 3d 53 |ARBIT. W|GHT(I)=S|
|000060c0| 49 47 59 28 49 29 0d 0a | 31 36 38 30 20 27 20 20 |IGY(I)..|1680 ' |
|000060d0| 59 46 49 54 3d 43 41 4c | 43 20 56 41 4c 55 45 53 |YFIT=CAL|C VALUES|
|000060e0| 20 4f 46 20 59 2c 20 20 | 41 30 3d 43 4f 4e 53 54 | OF Y, |A0=CONST|
|000060f0| 41 4e 54 20 54 45 52 4d | 0d 0a 31 36 39 30 20 27 |ANT TERM|..1690 '|
|00006100| 20 20 41 3d 41 52 52 41 | 59 20 4f 46 20 43 4f 45 | A=ARRA|Y OF COE|
|00006110| 46 46 2e 20 4f 46 20 54 | 48 45 20 50 4f 4c 59 4e |FF. OF T|HE POLYN|
|00006120| 4f 4d 49 41 4c 0d 0a 31 | 37 30 30 20 27 20 20 43 |OMIAL..1|700 ' C|
|00006130| 48 49 53 51 52 3d 52 45 | 44 55 43 45 44 20 43 48 |HISQR=RE|DUCED CH|
|00006140| 49 20 53 51 55 41 52 45 | 20 46 4f 52 20 46 49 54 |I SQUARE| FOR FIT|
|00006150| 0d 0a 31 37 31 30 20 27 | 20 20 0d 0a 31 37 32 30 |..1710 '| ..1720|
|00006160| 20 27 20 20 53 55 42 52 | 4f 55 54 49 4e 45 53 20 | ' SUBR|OUTINES |
|00006170| 52 45 51 55 49 52 45 44 | 2d 2d 2d 2d 2d 0d 0a 31 |REQUIRED|-----..1|
|00006180| 37 33 30 20 27 20 20 20 | 20 20 44 45 54 45 52 4d |730 ' | DETERM|
|00006190| 28 41 52 52 41 59 2c 4e | 4f 52 44 45 52 29 0d 0a |(ARRAY,N|ORDER)..|
|000061a0| 31 37 34 30 20 27 0d 0a | 31 37 35 30 20 27 09 20 |1740 '..|1750 '. |
|000061b0| 49 4e 49 54 20 53 55 4d | 53 20 26 20 41 52 52 41 |INIT SUM|S & ARRA|
|000061c0| 59 53 0d 0a 31 37 36 30 | 20 4e 4d 41 58 3d 20 32 |YS..1760| NMAX= 2|
|000061d0| 2a 4e 54 52 4d 2d 31 0d | 0a 31 37 37 30 20 46 4f |*NTRM-1.|.1770 FO|
|000061e0| 52 20 4e 3d 31 20 54 4f | 20 4e 4d 41 58 20 20 3a |R N=1 TO| NMAX :|
|000061f0| 20 53 55 4d 58 23 28 4e | 29 3d 30 21 20 20 3a 20 | SUMX#(N|)=0! : |
|00006200| 4e 45 58 54 20 4e 0d 0a | 31 37 38 30 20 46 4f 52 |NEXT N..|1780 FOR|
|00006210| 20 4a 3d 31 20 54 4f 20 | 4e 54 52 4d 20 20 3a 20 | J=1 TO |NTRM : |
|00006220| 53 55 4d 59 23 28 4a 29 | 3d 30 21 20 20 3a 20 4e |SUMY#(J)|=0! : N|
|00006230| 45 58 54 20 4a 0d 0a 31 | 37 39 30 20 43 48 49 53 |EXT J..1|790 CHIS|
|00006240| 51 23 3d 20 30 21 0d 0a | 31 38 30 30 20 27 09 2a |Q#= 0!..|1800 '.*|
|00006250| 2a 2a 20 41 43 43 55 4d | 55 4c 41 54 45 20 57 45 |** ACCUM|ULATE WE|
|00006260| 49 47 48 54 45 44 20 53 | 55 4d 53 20 2a 2a 2a 0d |IGHTED S|UMS ***.|
|00006270| 0a 31 38 31 30 20 46 4f | 52 20 49 3d 31 20 54 4f |.1810 FO|R I=1 TO|
|00006280| 20 4e 50 54 53 0d 0a 31 | 38 32 30 20 20 58 31 3d | NPTS..1|820 X1=|
|00006290| 58 28 49 29 20 20 3a 20 | 59 31 3d 59 28 49 29 0d |X(I) : |Y1=Y(I).|
|000062a0| 0a 31 38 33 30 20 20 49 | 46 20 4d 4f 44 45 3c 3e |.1830 I|F MODE<>|
|000062b0| 2b 31 20 20 47 4f 54 4f | 20 31 38 35 30 0d 0a 31 |+1 GOTO| 1850..1|
|000062c0| 38 34 30 20 20 20 57 47 | 48 54 28 49 29 3d 31 2f |840 WG|HT(I)=1/|
|000062d0| 28 53 49 47 59 28 49 29 | 2a 53 49 47 59 28 49 29 |(SIGY(I)|*SIGY(I)|
|000062e0| 29 20 20 3a 20 47 4f 54 | 4f 20 31 39 33 30 0d 0a |) : GOT|O 1930..|
|000062f0| 31 38 35 30 20 20 49 46 | 20 4d 4f 44 45 3c 3e 30 |1850 IF| MODE<>0|
|00006300| 20 20 47 4f 54 4f 20 31 | 38 37 30 0d 0a 31 38 36 | GOTO 1|870..186|
|00006310| 30 20 20 20 57 47 48 54 | 28 49 29 3d 31 20 20 3a |0 WGHT|(I)=1 :|
|00006320| 20 47 4f 54 4f 20 31 39 | 33 30 0d 0a 31 38 37 30 | GOTO 19|30..1870|
|00006330| 20 20 49 46 20 4d 4f 44 | 45 3c 3e 2d 31 20 20 47 | IF MOD|E<>-1 G|
|00006340| 4f 54 4f 20 31 39 30 30 | 0d 0a 31 38 38 30 20 20 |OTO 1900|..1880 |
|00006350| 20 49 46 20 59 28 49 29 | 3c 30 20 20 54 48 45 4e | IF Y(I)|<0 THEN|
|00006360| 20 57 47 48 54 28 49 29 | 3d 20 31 2f 28 2d 59 31 | WGHT(I)|= 1/(-Y1|
|00006370| 29 20 20 45 4c 53 45 20 | 20 20 57 47 48 54 28 49 |) ELSE | WGHT(I|
|00006380| 29 3d 20 31 2f 59 31 0d | 0a 31 38 39 30 20 20 20 |)= 1/Y1.|.1890 |
|00006390| 20 47 4f 54 4f 20 31 39 | 33 30 0d 0a 31 39 30 30 | GOTO 19|30..1900|
|000063a0| 20 20 49 46 20 4d 4f 44 | 45 3d 2d 32 20 20 47 4f | IF MOD|E=-2 GO|
|000063b0| 54 4f 20 31 39 32 30 0d | 0a 31 39 31 30 20 20 20 |TO 1920.|.1910 |
|000063c0| 50 52 49 4e 54 22 20 45 | 52 52 4f 52 20 49 4e 20 |PRINT" E|RROR IN |
|000063d0| 4d 4f 44 45 20 53 50 45 | 43 49 46 49 43 41 54 49 |MODE SPE|CIFICATI|
|000063e0| 4f 4e 22 20 20 3a 20 52 | 45 54 55 52 4e 0d 0a 31 |ON" : R|ETURN..1|
|000063f0| 39 32 30 20 20 20 57 47 | 48 54 28 49 29 3d 20 53 |920 WG|HT(I)= S|
+--------+-------------------------+-------------------------+--------+--------+
Only 25.0 KB of data is shown above.
