/ ARM Club 1 / ARM_CLUB_CD.iso / contents / apps / maths / progs / maths1 / A_Pole / Matrix

next >

MacOS 8.1

Win98

DOS

view JSON data

view as text

+--------+-------------------------+-------------------------+--------+--------+
|00000000| 0d 00 0a 72 f4 3e 4d 41 | 54 52 49 58 20 4d 61 74 |...r.>MA|TRIX Mat|
|00000010| 72 69 78 20 4c 55 20 64 | 65 63 6f 6d 70 2e 2c 20 |rix LU d|ecomp., |
|00000020| 43 68 6f 6c 65 73 6b 79 | 20 64 65 63 6f 6d 70 2e |Cholesky| decomp.|
|00000030| 2c 20 69 6e 76 65 72 73 | 65 20 28 4c 55 20 61 6e |, invers|e (LU an|
|00000040| 64 20 43 68 6f 6c 65 73 | 6b 79 29 2c 20 64 65 74 |d Choles|ky), det|
|00000050| 2e 3b 20 63 75 62 69 63 | 20 73 70 6c 69 6e 65 73 |.; cubic| splines|
|00000060| 3b 20 6c 69 6e 65 61 72 | 20 65 71 75 61 74 69 6f |; linear| equatio|
|00000070| 6e 73 0d 00 14 05 f4 0d | 00 1e 23 f4 20 4d 61 74 |ns......|..#. Mat|
|00000080| 68 65 6d 61 74 69 63 61 | 6c 20 6c 69 62 72 61 72 |hematica|l librar|
|00000090| 79 20 72 6f 75 74 69 6e | 65 73 0d 00 28 04 0d 00 |y routin|es..(...|
|000000a0| 32 04 0d 00 78 31 f4 20 | 20 20 20 20 20 20 20 20 |2...x1. |        |
|000000b0| 20 20 20 20 20 20 20 20 | 20 20 2a 2a 2a 20 4d 41 |        |  *** MA|
|000000c0| 54 52 49 58 20 4f 50 45 | 52 41 54 49 4f 4e 53 20 |TRIX OPE|RATIONS |
|000000d0| 2a 2a 2a 0d 00 82 05 3a | 0d 00 8c 32 3a f4 20 44 |***....:|...2:. D|
|000000e0| 65 74 65 72 6d 69 6e 61 | 6e 74 20 6f 66 20 6d 61 |etermina|nt of ma|
|000000f0| 74 72 69 78 2d 2d 76 69 | 61 20 4c 55 20 64 65 63 |trix--vi|a LU dec|
|00000100| 6f 6d 70 6f 73 69 74 69 | 6f 6e 0d 00 96 05 3a 0d |ompositi|on....:.|
|00000110| 00 a0 3c dd a4 44 45 54 | 28 58 28 29 29 3a ea 20 |..<..DET|(X()):. |
|00000120| 4e 25 2c 44 25 2c 49 4e | 44 58 25 28 29 2c 4a 25 |N%,D%,IN|DX%(),J%|
|00000130| 2c 4c 55 28 29 2c 44 45 | 54 3a 4e 25 3d de 28 58 |,LU(),DE|T:N%=.(X|
|00000140| 28 29 2c 31 29 3a 44 45 | 54 3d 31 0d 00 aa 21 de |(),1):DE|T=1...!.|
|00000150| 49 4e 44 58 25 28 4e 25 | 29 2c 4c 55 28 4e 25 2c |INDX%(N%|),LU(N%,|
|00000160| 4e 25 29 3a 4c 55 28 29 | 3d 58 28 29 0d 00 b4 1f |N%):LU()|=X()....|
|00000170| f2 4c 55 44 43 4d 50 28 | 4c 55 28 29 2c 4e 25 2c |.LUDCMP(|LU(),N%,|
|00000180| 49 4e 44 58 25 28 29 2c | 44 25 29 0d 00 be 25 e3 |INDX%(),|D%)...%.|
|00000190| 4a 25 3d 30 b8 4e 25 3a | 44 45 54 3d 44 45 54 2a |J%=0.N%:|DET=DET*|
|000001a0| 4c 55 28 4a 25 2c 4a 25 | 29 3a ed 3a 3d 44 45 54 |LU(J%,J%|):.:=DET|
|000001b0| 0d 00 bf 05 3a 0d 00 c0 | 3c 3a f4 20 44 65 74 65 |....:...|<:. Dete|
|000001c0| 72 6d 69 6e 61 6e 74 20 | 6f 66 20 73 70 64 20 6d |rminant |of spd m|
|000001d0| 61 74 72 69 78 2d 2d 76 | 69 61 20 43 68 6f 6c 65 |atrix--v|ia Chole|
|000001e0| 73 6b 79 20 64 65 63 6f | 6d 70 6f 73 69 74 69 6f |sky deco|mpositio|
|000001f0| 6e 0d 00 c1 05 3a 0d 00 | c2 3c dd 20 a4 43 44 45 |n....:..|.<. .CDE|
|00000200| 54 28 58 28 29 29 3a ea | 4e 25 2c 44 45 54 2c 49 |T(X()):.|N%,DET,I|
|00000210| 25 2c 55 28 29 3a 4e 25 | 3d de 28 58 28 29 2c 31 |%,U():N%|=.(X(),1|
|00000220| 29 3a de 20 55 28 4e 25 | 2c 4e 25 29 3a 44 45 54 |):. U(N%|,N%):DET|
|00000230| 3d 31 0d 00 c3 39 f2 43 | 48 4f 4c 45 53 4b 59 28 |=1...9.C|HOLESKY(|
|00000240| 58 28 29 2c 55 28 29 29 | 3a e3 49 25 3d 30 b8 4e |X(),U())|:.I%=0.N|
|00000250| 25 3a 44 45 54 3d 44 45 | 54 2a 55 28 49 25 2c 49 |%:DET=DE|T*U(I%,I|
|00000260| 25 29 5e 32 3a ed 3a 3d | 44 45 54 0d 00 c8 05 3a |%)^2:.:=|DET....:|
|00000270| 0d 00 d2 16 3a f4 20 49 | 64 65 6e 74 69 74 79 20 |....:. I|dentity |
|00000280| 6d 61 74 72 69 78 0d 00 | dc 05 3a 0d 00 e6 2c dd |matrix..|..:...,.|
|00000290| f2 49 44 4d 41 54 28 f8 | 20 49 28 29 29 3a ea 20 |.IDMAT(.| I()):. |
|000002a0| 4e 25 2c 4a 25 3a 49 28 | 29 3d 30 3a 4e 25 3d de |N%,J%:I(|)=0:N%=.|
|000002b0| 28 49 28 29 2c 31 29 0d | 00 f0 1e e3 20 4a 25 3d |(I(),1).|.... J%=|
|000002c0| 30 20 b8 20 4e 25 3a 49 | 28 4a 25 2c 4a 25 29 3d |0 . N%:I|(J%,J%)=|
|000002d0| 31 3a ed 3a e1 0d 00 fa | 05 3a 0d 01 04 3f 3a f4 |1:.:....|.:...?:.|
|000002e0| 20 4d 61 74 72 69 78 20 | 69 6e 76 65 72 73 69 6f | Matrix |inversio|
|000002f0| 6e 20 28 6d 6f 64 69 66 | 69 65 64 20 66 72 6f 6d |n (modif|ied from|
|00000300| 20 4e 75 6d 65 72 69 63 | 61 6c 20 52 65 63 69 70 | Numeric|al Recip|
|00000310| 65 73 20 70 2e 20 33 38 | 29 0d 01 0e 05 3a 0d 01 |es p. 38|)....:..|
|00000320| 18 15 dd f2 49 4d 41 54 | 28 4d 28 29 2c f8 20 59 |....IMAT|(M(),. Y|
|00000330| 28 29 29 0d 01 22 2e ea | 20 4e 25 2c 20 44 25 2c |()).."..| N%, D%,|
|00000340| 4a 25 2c 4b 25 2c 41 28 | 29 2c 42 28 29 2c 49 4e |J%,K%,A(|),B(),IN|
|00000350| 44 58 25 28 29 3a 4e 25 | 3d de 28 4d 28 29 2c 31 |DX%():N%|=.(M(),1|
|00000360| 29 0d 01 2c 26 de 20 49 | 4e 44 58 25 28 4e 25 29 |)..,&. I|NDX%(N%)|
|00000370| 2c 42 28 4e 25 29 2c 41 | 28 4e 25 2c 4e 25 29 3a |,B(N%),A|(N%,N%):|
|00000380| 41 28 29 3d 4d 28 29 0d | 01 36 2a f2 49 44 4d 41 |A()=M().|.6*.IDMA|
|00000390| 54 28 59 28 29 29 3a f2 | 4c 55 44 43 4d 50 28 41 |T(Y()):.|LUDCMP(A|
|000003a0| 28 29 2c 4e 25 2c 49 4e | 44 58 25 28 29 2c 44 25 |(),N%,IN|DX%(),D%|
|000003b0| 29 0d 01 40 04 0d 01 4a | 20 e7 20 44 25 3d 30 3a |)..@...J| . D%=0:|
|000003c0| e1 20 20 3a f4 20 6d 61 | 74 72 69 78 20 73 69 6e |.  :. ma|trix sin|
|000003d0| 67 75 6c 61 72 0d 01 54 | 04 0d 01 5e 3e e3 20 4a |gular..T|...^>. J|
|000003e0| 25 3d 30 20 b8 4e 25 3a | e3 20 4b 25 3d 30 20 b8 |%=0 .N%:|. K%=0 .|
|000003f0| 4e 25 3a 42 28 4b 25 29 | 3d 59 28 4b 25 2c 4a 25 |N%:B(K%)|=Y(K%,J%|
|00000400| 29 3a ed 3a 20 f4 20 73 | 65 65 20 63 6f 6d 6d 6d |):.: . s|ee commm|
|00000410| 65 6e 74 73 20 69 6e 0d | 01 68 46 f2 4c 55 42 4b |ents in.|.hF.LUBK|
|00000420| 53 42 28 41 28 29 2c 4e | 25 2c 49 4e 44 58 25 28 |SB(A(),N|%,INDX%(|
|00000430| 29 2c 42 28 29 29 20 20 | 20 20 20 20 20 20 20 20 |),B())  |        |
|00000440| 20 20 20 20 20 20 20 3a | 20 f4 20 61 62 6f 76 65 |       :| . above|
|00000450| 20 72 65 66 65 72 65 6e | 63 65 2e 2e 2e 0d 01 72 | referen|ce.....r|
|00000460| 1f e3 20 4b 25 3d 30 b8 | 20 4e 25 3a 59 28 4b 25 |.. K%=0.| N%:Y(K%|
|00000470| 2c 4a 25 29 3d 42 28 4b | 25 29 3a ed 0d 01 7c 07 |,J%)=B(K|%):...|.|
|00000480| ed 3a e1 0d 01 86 04 0d | 01 90 05 3a 0d 01 9a 3e |.:......|...:...>|
|00000490| 3a f4 20 4c 55 20 64 65 | 63 6f 6d 70 6f 73 69 74 |:. LU de|composit|
|000004a0| 69 6f 6e 20 28 6d 6f 64 | 69 66 69 65 64 20 66 72 |ion (mod|ified fr|
|000004b0| 6f 6d 20 4e 75 6d 65 72 | 69 63 61 6c 20 52 65 63 |om Numer|ical Rec|
|000004c0| 69 70 65 73 20 70 2e 33 | 35 29 0d 01 a4 2d 3a f4 |ipes p.3|5)...-:.|
|000004d0| 20 52 65 73 75 6c 74 20 | 69 73 20 72 65 74 75 72 | Result |is retur|
|000004e0| 6e 65 64 20 69 6e 20 6f | 72 69 67 69 6e 61 6c 20 |ned in o|riginal |
|000004f0| 6d 61 74 72 69 78 3a 0d | 01 ae 34 3a f4 20 55 20 |matrix:.|..4:. U |
|00000500| 69 73 20 74 68 65 20 75 | 70 70 65 72 20 74 72 69 |is the u|pper tri|
|00000510| 61 6e 67 75 6c 61 72 20 | 70 61 72 74 20 6f 66 20 |angular |part of |
|00000520| 74 68 65 20 6c 61 74 74 | 65 72 3b 0d 01 b8 3d 3a |the latt|er;...=:|
|00000530| f4 20 4c 20 69 73 20 74 | 68 65 20 73 74 72 69 63 |. L is t|he stric|
|00000540| 74 20 6c 6f 77 65 72 20 | 70 61 72 74 20 74 6f 67 |t lower |part tog|
|00000550| 65 74 68 65 72 20 77 69 | 74 68 20 75 6e 69 74 20 |ether wi|th unit |
|00000560| 64 69 61 67 6f 6e 61 6c | 0d 01 c2 41 3a f4 20 4e |diagonal|...A:. N|
|00000570| 6f 74 65 20 74 68 61 74 | 20 74 68 65 20 72 65 73 |ote that| the res|
|00000580| 75 6c 74 20 69 73 20 74 | 68 65 20 4c 55 20 64 65 |ult is t|he LU de|
|00000590| 63 6f 6d 70 6f 73 69 74 | 69 6f 6e 20 6f 66 20 61 |composit|ion of a|
|000005a0| 20 72 6f 77 2d 77 69 73 | 65 0d 01 cc 49 3a f4 20 | row-wis|e...I:. |
|000005b0| 70 65 72 6d 75 74 61 74 | 69 6f 6e 20 6f 66 20 74 |permutat|ion of t|
|000005c0| 68 65 20 6f 72 69 67 69 | 6e 61 6c 20 6d 61 74 72 |he origi|nal matr|
|000005d0| 69 78 3a 20 49 4e 44 58 | 25 20 72 65 63 6f 72 64 |ix: INDX|% record|
|000005e0| 73 20 74 68 65 20 70 65 | 72 6d 75 74 61 74 69 6f |s the pe|rmutatio|
|000005f0| 6e 2e 0d 01 d6 05 3a 0d | 01 e0 25 dd f2 4c 55 44 |n.....:.|..%..LUD|
|00000600| 43 4d 50 28 f8 20 41 28 | 29 2c 4e 25 2c f8 20 49 |CMP(. A(|),N%,. I|
|00000610| 4e 44 58 25 28 29 2c f8 | 20 44 25 29 0d 01 ea 34 |NDX%(),.| D%)...4|
|00000620| ea 20 4e 4d 41 58 25 2c | 54 49 4e 59 2c 49 25 2c |. NMAX%,|TINY,I%,|
|00000630| 4a 25 2c 4b 25 2c 41 41 | 4d 41 58 2c 58 2c 73 55 |J%,K%,AA|MAX,X,sU|
|00000640| 4d 2c 44 55 4d 2c 49 4d | 41 58 25 2c 56 56 28 29 |M,DUM,IM|AX%,VV()|
|00000650| 0d 01 f4 25 4e 4d 41 58 | 25 3d 31 30 30 3a 54 49 |...%NMAX|%=100:TI|
|00000660| 4e 59 3d 30 3a 44 25 3d | 31 3a de 20 56 56 28 4e |NY=0:D%=|1:. VV(N|
|00000670| 4d 41 58 25 29 0d 01 fe | 16 e3 20 49 25 3d 30 b8 |MAX%)...|.. I%=0.|
|00000680| 20 4e 25 3a 41 41 4d 41 | 58 3d 30 0d 02 08 2e e3 | N%:AAMA|X=0.....|
|00000690| 20 4a 25 3d 30 b8 4e 25 | 3a 58 3d 94 28 41 28 49 | J%=0.N%|:X=.(A(I|
|000006a0| 25 2c 4a 25 29 29 3a e7 | 20 58 3e 41 41 4d 41 58 |%,J%)):.| X>AAMAX|
|000006b0| 20 3a 41 41 4d 41 58 3d | 58 0d 02 12 05 ed 0d 02 | :AAMAX=|X.......|
|000006c0| 1c 31 e7 20 41 41 4d 41 | 58 3d 30 3a 20 f1 27 22 |.1. AAMA|X=0: .'"|
|000006d0| 2a 2a 2a 20 53 49 4e 47 | 55 4c 41 52 20 4d 41 54 |*** SING|ULAR MAT|
|000006e0| 52 49 58 20 2a 2a 2a 22 | 3a 44 25 3d 30 3a e1 0d |RIX ***"|:D%=0:..|
|000006f0| 02 26 15 56 56 28 49 25 | 29 3d 31 2f 41 41 4d 41 |.&.VV(I%|)=1/AAMA|
|00000700| 58 20 3a ed 0d 02 30 0e | e3 20 4a 25 3d 30 b8 20 |X :...0.|. J%=0. |
|00000710| 4e 25 0d 02 3a 0c e7 20 | 4a 25 3e 31 20 8c 0d 02 |N%..:.. |J%>1 ...|
|00000720| 44 1d e3 20 49 25 3d 30 | 20 b8 4a 25 2d 31 3a 73 |D.. I%=0| .J%-1:s|
|00000730| 55 4d 3d 41 28 49 25 2c | 4a 25 29 0d 02 4e 0c e7 |UM=A(I%,|J%)..N..|
|00000740| 20 49 25 3e 30 20 8c 0d | 02 58 37 e3 4b 25 3d 30 | I%>0 ..|.X7.K%=0|
|00000750| b8 49 25 2d 31 3a 73 55 | 4d 3d 73 55 4d 2d 41 28 |.I%-1:sU|M=sUM-A(|
|00000760| 49 25 2c 4b 25 29 2a 41 | 28 4b 25 2c 4a 25 29 3a |I%,K%)*A|(K%,J%):|
|00000770| ed 3a 41 28 49 25 2c 4a | 25 29 3d 73 55 4d 0d 02 |.:A(I%,J|%)=sUM..|
|00000780| 62 07 cd 3a ed 0d 02 6c | 05 cd 0d 02 76 0b 41 41 |b..:...l|....v.AA|
|00000790| 4d 41 58 3d 30 0d 02 80 | 1a e3 49 25 3d 4a 25 b8 |MAX=0...|..I%=J%.|
|000007a0| 4e 25 3a 73 55 4d 3d 41 | 28 49 25 2c 4a 25 29 0d |N%:sUM=A|(I%,J%).|
|000007b0| 02 8a 0c e7 20 4a 25 3e | 30 20 8c 0d 02 94 37 e3 |.... J%>|0 ....7.|
|000007c0| 4b 25 3d 30 b8 4a 25 2d | 31 3a 73 55 4d 3d 73 55 |K%=0.J%-|1:sUM=sU|
|000007d0| 4d 2d 41 28 49 25 2c 4b | 25 29 2a 41 28 4b 25 2c |M-A(I%,K|%)*A(K%,|
|000007e0| 4a 25 29 3a ed 3a 41 28 | 49 25 2c 4a 25 29 3d 73 |J%):.:A(|I%,J%)=s|
|000007f0| 55 4d 0d 02 9e 05 cd 0d | 02 a8 15 44 55 4d 3d 56 |UM......|...DUM=V|
|00000800| 56 28 49 25 29 2a 94 28 | 73 55 4d 29 0d 02 b2 24 |V(I%)*.(|sUM)...$|
|00000810| e7 20 44 55 4d 3e 3d 41 | 41 4d 41 58 20 3a 49 4d |. DUM>=A|AMAX :IM|
|00000820| 41 58 25 3d 49 25 3a 41 | 41 4d 41 58 3d 44 55 4d |AX%=I%:A|AMAX=DUM|
|00000830| 0d 02 bc 05 ed 0d 02 c6 | 11 e7 20 4a 25 3c 3e 49 |........|.. J%<>I|
|00000840| 4d 41 58 25 20 8c 0d 02 | d0 26 e3 4b 25 3d 30 b8 |MAX% ...|.&.K%=0.|
|00000850| 4e 25 3a c8 94 20 41 28 | 49 4d 41 58 25 2c 4b 25 |N%:.. A(|IMAX%,K%|
|00000860| 29 2c 41 28 4a 25 2c 4b | 25 29 3a ed 0d 02 da 1b |),A(J%,K|%):.....|
|00000870| 44 25 3d 2d 44 25 3a 56 | 56 28 49 4d 41 58 25 29 |D%=-D%:V|V(IMAX%)|
|00000880| 3d 56 56 28 4a 25 29 0d | 02 e4 05 cd 0d 02 ee 13 |=VV(J%).|........|
|00000890| 49 4e 44 58 25 28 4a 25 | 29 3d 49 4d 41 58 25 0d |INDX%(J%|)=IMAX%.|
|000008a0| 02 f8 0e e7 20 4a 25 3c | 3e 4e 25 20 8c 0d 03 02 |.... J%<|>N% ....|
|000008b0| 1e e7 20 41 28 4a 25 2c | 4a 25 29 3d 30 3a 41 28 |.. A(J%,|J%)=0:A(|
|000008c0| 4a 25 2c 4a 25 29 3d 54 | 49 4e 59 0d 03 0c 12 44 |J%,J%)=T|INY....D|
|000008d0| 55 4d 3d 31 2f 41 28 4a | 25 2c 4a 25 29 0d 03 16 |UM=1/A(J|%,J%)...|
|000008e0| 27 e3 49 25 3d 4a 25 2b | 31 b8 4e 25 3a 41 28 49 |'.I%=J%+|1.N%:A(I|
|000008f0| 25 2c 4a 25 29 3d 41 28 | 49 25 2c 4a 25 29 2a 44 |%,J%)=A(|I%,J%)*D|
|00000900| 55 4d 3a ed 0d 03 20 05 | cd 0d 03 2a 08 ed 20 4a |UM:... .|...*.. J|
|00000910| 25 0d 03 34 1e e7 20 41 | 28 4e 25 2c 4e 25 29 3d |%..4.. A|(N%,N%)=|
|00000920| 30 3a 41 28 4e 25 2c 4e | 25 29 3d 54 49 4e 59 0d |0:A(N%,N|%)=TINY.|
|00000930| 03 3e 05 e1 0d 03 48 05 | 3a 0d 03 52 47 3a f4 20 |.>....H.|:..RG:. |
|00000940| 42 61 63 6b 73 75 62 73 | 74 69 74 75 74 69 6f 6e |Backsubs|titution|
|00000950| 20 61 66 74 65 72 20 4c | 55 20 64 65 63 6f 6d 70 | after L|U decomp|
|00000960| 6f 73 69 74 69 6f 6e 20 | 62 79 20 4c 55 44 43 4d |osition |by LUDCM|
|00000970| 50 20 28 6d 6f 64 69 66 | 69 65 64 20 66 72 6f 6d |P (modif|ied from|
|00000980| 0d 03 5c 1e 3a f4 20 4e | 75 6d 65 72 69 63 61 6c |..\.:. N|umerical|
|00000990| 20 52 65 63 69 70 65 73 | 20 70 2e 33 36 29 0d 03 | Recipes| p.36)..|
|000009a0| 66 05 3a 0d 03 70 20 dd | f2 4c 55 42 4b 53 42 28 |f.:..p .|.LUBKSB(|
|000009b0| 41 28 29 2c 4e 25 2c 49 | 4e 44 58 25 28 29 2c 42 |A(),N%,I|NDX%(),B|
|000009c0| 28 29 29 0d 03 7a 25 ea | 20 49 49 25 2c 4c 4c 25 |())..z%.| II%,LL%|
|000009d0| 2c 73 55 4d 2c 49 25 2c | 4a 25 2c 64 6f 73 75 6d |,sUM,I%,|J%,dosum|
|000009e0| 3a 64 6f 73 75 6d 3d a3 | 0d 03 84 32 e3 49 25 3d |:dosum=.|...2.I%=|
|000009f0| 30 b8 4e 25 3a 4c 4c 25 | 3d 49 4e 44 58 25 28 49 |0.N%:LL%|=INDX%(I|
|00000a00| 25 29 3a 73 55 4d 3d 42 | 28 4c 4c 25 29 3a 42 28 |%):sUM=B|(LL%):B(|
|00000a10| 4c 4c 25 29 3d 42 28 49 | 25 29 0d 03 8e 0d e7 20 |LL%)=B(I|%)..... |
|00000a20| 64 6f 73 75 6d 20 8c 0d | 03 98 2a e3 4a 25 3d 49 |dosum ..|..*.J%=I|
|00000a30| 49 25 b8 49 25 2d 31 3a | 3a 73 55 4d 3d 73 55 4d |I%.I%-1:|:sUM=sUM|
|00000a40| 2d 41 28 49 25 2c 4a 25 | 29 2a 42 28 4a 25 29 3a |-A(I%,J%|)*B(J%):|
|00000a50| ed 0d 03 a2 1d cc 20 e7 | 20 73 55 4d 3c 3e 30 3a |...... .| sUM<>0:|
|00000a60| 64 6f 73 75 6d 3d b9 3a | 49 49 25 3d 49 25 0d 03 |dosum=.:|II%=I%..|
|00000a70| ac 05 cd 0d 03 b6 0f 42 | 28 49 25 29 3d 73 55 4d |.......B|(I%)=sUM|
|00000a80| 3a ed 0d 03 c0 19 e3 49 | 25 3d 4e 25 b8 30 88 2d |:......I|%=N%.0.-|
|00000a90| 31 3a 73 55 4d 3d 42 28 | 49 25 29 0d 03 ca 0d e7 |1:sUM=B(|I%).....|
|00000aa0| 20 49 25 3c 4e 25 20 8c | 0d 03 d4 28 e3 4a 25 3d | I%<N% .|...(.J%=|
|00000ab0| 49 25 2b 31 b8 4e 25 3a | 73 55 4d 3d 73 55 4d 2d |I%+1.N%:|sUM=sUM-|
|00000ac0| 41 28 49 25 2c 4a 25 29 | 2a 42 28 4a 25 29 3a ed |A(I%,J%)|*B(J%):.|
|00000ad0| 0d 03 de 05 cd 0d 03 e8 | 1a 42 28 49 25 29 3d 73 |........|.B(I%)=s|
|00000ae0| 55 4d 2f 41 28 49 25 2c | 49 25 29 3a ed 3a e1 0d |UM/A(I%,|I%):.:..|
|00000af0| 03 f2 05 3a 0d 03 fc 45 | 3a f4 20 43 68 6f 6c 65 |...:...E|:. Chole|
|00000b00| 73 6b 69 20 64 65 63 6f | 6d 70 6f 73 69 74 69 6f |ski deco|mpositio|
|00000b10| 6e 20 28 6f 66 20 73 79 | 6d 6d 65 74 72 69 63 20 |n (of sy|mmetric |
|00000b20| 70 6f 73 69 74 69 76 65 | 20 64 65 66 69 6e 69 74 |positive| definit|
|00000b30| 65 20 6d 61 74 72 69 78 | 29 0d 04 06 17 3a f4 20 |e matrix|)....:. |
|00000b40| 28 73 65 65 20 53 65 62 | 65 72 20 70 33 30 34 29 |(see Seb|er p304)|
|00000b50| 0d 04 10 05 3a 0d 04 1a | 1a dd 20 f2 43 48 4f 4c |....:...|.. .CHOL|
|00000b60| 45 53 4b 59 28 58 28 29 | 2c f8 20 55 28 29 29 0d |ESKY(X()|,. U()).|
|00000b70| 04 24 32 ea 4e 25 2c 49 | 25 2c 4a 25 2c 4b 25 2c |.$2.N%,I|%,J%,K%,|
|00000b80| 73 75 6d 3a 20 4e 25 3d | de 28 58 28 29 2c 31 29 |sum: N%=|.(X(),1)|
|00000b90| 3a 55 28 30 2c 30 29 3d | b6 28 58 28 30 2c 30 29 |:U(0,0)=|.(X(0,0)|
|00000ba0| 29 0d 04 2e 25 e3 4a 25 | 3d 31 b8 4e 25 3a 55 28 |)...%.J%|=1.N%:U(|
|00000bb0| 30 2c 4a 25 29 3d 58 28 | 30 2c 4a 25 29 2f 55 28 |0,J%)=X(|0,J%)/U(|
|00000bc0| 30 2c 30 29 3a ed 0d 04 | 38 32 e3 49 25 3d 31 b8 |0,0):...|82.I%=1.|
|00000bd0| 4e 25 3a 73 75 6d 3d 30 | 3a e3 4b 25 3d 30 b8 49 |N%:sum=0|:.K%=0.I|
|00000be0| 25 2d 31 3a 73 75 6d 3d | 73 75 6d 2b 55 28 4b 25 |%-1:sum=|sum+U(K%|
|00000bf0| 2c 49 25 29 5e 32 3a ed | 0d 04 42 1c 55 28 49 25 |,I%)^2:.|..B.U(I%|
|00000c00| 2c 49 25 29 3d b6 28 58 | 28 49 25 2c 49 25 29 2d |,I%)=.(X|(I%,I%)-|
|00000c10| 73 75 6d 29 0d 04 4c 2b | e7 55 28 49 25 2c 49 25 |sum)..L+|.U(I%,I%|
|00000c20| 29 3c 3d 30 3a f1 22 4e | 4f 54 20 50 4f 53 49 54 |)<=0:."N|OT POSIT|
|00000c30| 49 56 45 20 44 45 46 49 | 4e 49 54 45 22 3a e1 0d |IVE DEFI|NITE":..|
|00000c40| 04 56 0d e7 20 49 25 3d | 4e 25 3a e1 0d 04 60 3c |.V.. I%=|N%:...`<|
|00000c50| e3 4a 25 3d 49 25 2b 31 | b8 4e 25 3a 73 75 6d 3d |.J%=I%+1|.N%:sum=|
|00000c60| 30 3a e3 4b 25 3d 30 b8 | 49 25 2d 31 3a 73 75 6d |0:.K%=0.|I%-1:sum|
|00000c70| 3d 73 75 6d 2b 55 28 4b | 25 2c 49 25 29 2a 55 28 |=sum+U(K|%,I%)*U(|
|00000c80| 4b 25 2c 4a 25 29 3a ed | 0d 04 6a 2a 55 28 49 25 |K%,J%):.|..j*U(I%|
|00000c90| 2c 4a 25 29 3d 28 58 28 | 49 25 2c 4a 25 29 2d 73 |,J%)=(X(|I%,J%)-s|
|00000ca0| 75 6d 29 2f 55 28 49 25 | 2c 49 25 29 3a ed 3a ed |um)/U(I%|,I%):.:.|
|00000cb0| 3a e1 0d 04 74 05 3a 0d | 04 7e 49 3a f4 20 49 6e |:...t.:.|.~I:. In|
|00000cc0| 76 65 72 73 65 20 6f 66 | 20 73 79 6d 6d 65 74 72 |verse of| symmetr|
|00000cd0| 69 63 20 70 6f 73 69 74 | 69 76 65 20 64 65 66 69 |ic posit|ive defi|
|00000ce0| 6e 69 74 65 20 6d 61 74 | 72 69 78 20 76 69 61 20 |nite mat|rix via |
|00000cf0| 43 68 6f 6c 65 73 6b 79 | 20 64 65 63 6f 6d 70 2e |Cholesky| decomp.|
|00000d00| 0d 04 88 05 3a 0d 04 92 | 1a dd 20 f2 49 53 50 44 |....:...|.. .ISPD|
|00000d10| 4d 41 54 28 58 28 29 2c | f8 20 49 58 28 29 29 0d |MAT(X(),|. IX()).|
|00000d20| 04 9c 3c ea 49 25 2c 4a | 25 2c 4b 25 2c 4e 25 2c |..<.I%,J|%,K%,N%,|
|00000d30| 54 28 29 2c 55 28 29 2c | 73 75 6d 3a 4e 25 3d de |T(),U(),|sum:N%=.|
|00000d40| 28 58 28 29 2c 31 29 3a | de 20 54 28 4e 25 2c 4e |(X(),1):|. T(N%,N|
|00000d50| 25 29 2c 55 28 4e 25 2c | 4e 25 29 0d 04 a6 16 f2 |%),U(N%,|N%).....|
|00000d60| 43 48 4f 4c 45 53 4b 59 | 28 58 28 29 2c 55 28 29 |CHOLESKY|(X(),U()|
|00000d70| 29 0d 04 b0 22 e3 49 25 | 3d 30 b8 4e 25 3a 54 28 |)...".I%|=0.N%:T(|
|00000d80| 49 25 2c 49 25 29 3d 31 | 2f 55 28 49 25 2c 49 25 |I%,I%)=1|/U(I%,I%|
|00000d90| 29 3a ed 0d 04 ba 11 e3 | 49 25 3d 4e 25 2d 31 b8 |):......|I%=N%-1.|
|00000da0| 30 88 2d 31 0d 04 c4 3e | 20 e3 4a 25 3d 49 25 2b |0.-1...>| .J%=I%+|
|00000db0| 31 b8 4e 25 3a 73 75 6d | 3d 30 3a e3 4b 25 3d 49 |1.N%:sum|=0:.K%=I|
|00000dc0| 25 2b 31 b8 4a 25 3a 73 | 75 6d 3d 73 75 6d 2b 55 |%+1.J%:s|um=sum+U|
|00000dd0| 28 49 25 2c 4b 25 29 2a | 54 28 4b 25 2c 4a 25 29 |(I%,K%)*|T(K%,J%)|
|00000de0| 3a ed 0d 04 ce 1f 20 54 | 28 49 25 2c 4a 25 29 3d |:..... T|(I%,J%)=|
|00000df0| 2d 73 75 6d 2f 55 28 49 | 25 2c 49 25 29 3a ed 3a |-sum/U(I|%,I%):.:|
|00000e00| ed 0d 04 d8 26 f2 54 52 | 41 4e 53 50 4f 53 45 28 |....&.TR|ANSPOSE(|
|00000e10| 54 28 29 2c 55 28 29 29 | 3a 49 58 28 29 3d 54 28 |T(),U())|:IX()=T(|
|00000e20| 29 2e 55 28 29 3a e1 0d | 04 e2 05 3a 0d 04 ec 17 |).U():..|...:....|
|00000e30| 3a f4 20 54 72 61 6e 73 | 70 6f 73 65 20 6d 61 74 |:. Trans|pose mat|
|00000e40| 72 69 78 0d 04 ed 05 3a | 0d 04 f6 1c dd 20 f2 54 |rix....:|..... .T|
|00000e50| 52 41 4e 53 50 4f 53 45 | 28 58 28 29 2c f8 20 58 |RANSPOSE|(X(),. X|
|00000e60| 54 28 29 29 0d 05 00 1a | ea 20 49 25 2c 4a 25 2c |T())....|. I%,J%,|
|00000e70| 4e 25 3a 4e 25 3d de 28 | 58 28 29 2c 31 29 0d 05 |N%:N%=.(|X(),1)..|
|00000e80| 0a 2e e3 49 25 3d 30 b8 | 4e 25 3a e3 4a 25 3d 30 |...I%=0.|N%:.J%=0|
|00000e90| b8 4e 25 3a 58 54 28 49 | 25 2c 4a 25 29 3d 58 28 |.N%:XT(I|%,J%)=X(|
|00000ea0| 4a 25 2c 49 25 29 3a ed | 3a ed 3a e1 0d 05 0b 04 |J%,I%):.|:.:.....|
|00000eb0| 0d 05 14 33 f4 20 20 20 | 20 20 20 20 20 20 20 20 |...3.   |        |
|00000ec0| 20 20 20 2a 2a 2a 20 45 | 4e 44 20 4f 46 20 4d 41 |   *** E|ND OF MA|
|00000ed0| 54 52 49 58 20 4f 50 45 | 52 41 54 49 4f 4e 53 20 |TRIX OPE|RATIONS |
|00000ee0| 2a 2a 2a 0d 05 78 04 0d | 05 79 04 0d 05 82 3b f4 |***..x..|.y....;.|
|00000ef0| 20 20 20 20 20 20 20 20 | 2a 2a 2a 20 43 55 42 49 |        |*** CUBI|
|00000f00| 43 20 53 50 4c 49 4e 45 | 20 46 49 54 54 49 4e 47 |C SPLINE| FITTING|
|00000f10| 20 41 4e 44 20 49 4e 54 | 45 52 50 4f 4c 41 54 49 | AND INT|ERPOLATI|
|00000f20| 4f 4e 20 2a 2a 2a 0d 05 | 8c 26 f4 20 53 65 65 3a |ON ***..|.&. See:|
|00000f30| 20 4e 75 6d 65 72 69 63 | 61 6c 20 52 65 63 69 70 | Numeric|al Recip|
|00000f40| 65 73 20 43 68 61 70 74 | 65 72 20 33 0d 05 96 05 |es Chapt|er 3....|
|00000f50| 3a 0d 05 a0 3c 3a f4 20 | 43 61 6c 63 75 6c 61 74 |:...<:. |Calculat|
|00000f60| 65 20 73 65 63 6f 6e 64 | 20 64 65 72 69 76 61 74 |e second| derivat|
|00000f70| 69 76 65 73 20 67 69 76 | 65 6e 20 74 61 62 75 6c |ives giv|en tabul|
|00000f80| 61 74 65 64 20 66 75 6e | 63 74 69 6f 6e 0d 05 aa |ated fun|ction...|
|00000f90| 1f 3a 20 20 20 20 28 6d | 6f 64 69 66 69 65 64 20 |.:    (m|odified |
|00000fa0| 66 72 6f 6d 20 4e 52 20 | 70 38 38 29 0d 05 b4 1f |from NR |p88)....|
|00000fb0| dd 20 f2 43 5f 53 50 4c | 49 4e 45 28 58 28 29 2c |. .C_SPL|INE(X(),|
|00000fc0| 59 28 29 2c f8 20 59 32 | 28 29 29 0d 05 be 35 ea |Y(),. Y2|())...5.|
|00000fd0| 20 4e 25 2c 4e 4d 41 58 | 25 2c 55 28 29 2c 49 25 | N%,NMAX|%,U(),I%|
|00000fe0| 2c 53 49 47 2c 50 2c 51 | 4e 2c 55 4e 3a 4e 25 3d |,SIG,P,Q|N,UN:N%=|
|00000ff0| de 28 58 28 29 2c 31 29 | 3a de 20 55 28 4e 25 29 |.(X(),1)|:. U(N%)|
|00001000| 0d 05 c8 04 0d 05 d2 20 | f4 20 6e 6f 74 2d 61 2d |....... |. not-a-|
|00001010| 6b 6e 6f 74 20 65 6e 64 | 20 63 6f 6e 64 69 74 69 |knot end| conditi|
|00001020| 6f 6e 73 3a 0d 05 dc 43 | 59 50 31 3d 28 59 28 31 |ons:...C|YP1=(Y(1|
|00001030| 29 2d 59 28 30 29 29 2f | 28 58 28 31 29 2d 58 28 |)-Y(0))/|(X(1)-X(|
|00001040| 30 29 29 3a 59 50 4e 3d | 28 59 28 4e 25 29 2d 59 |0)):YPN=|(Y(N%)-Y|
|00001050| 28 4e 25 2d 31 29 29 2f | 28 58 28 4e 25 29 2d 58 |(N%-1))/|(X(N%)-X|
|00001060| 28 4e 25 2d 31 29 29 0d | 05 e6 04 0d 05 f0 41 59 |(N%-1)).|......AY|
|00001070| 32 28 30 29 3d 2d 30 2e | 35 3a 55 28 30 29 3d 28 |2(0)=-0.|5:U(0)=(|
|00001080| 33 2f 28 58 28 31 29 2d | 58 28 30 29 29 29 2a 28 |3/(X(1)-|X(0)))*(|
|00001090| 28 59 28 31 29 2d 59 28 | 30 29 29 2f 28 58 28 31 |(Y(1)-Y(|0))/(X(1|
|000010a0| 29 2d 58 28 30 29 29 2d | 59 50 31 29 0d 05 fa 46 |)-X(0))-|YP1)...F|
|000010b0| e3 20 49 25 3d 31 b8 4e | 25 2d 31 3a 53 49 47 3d |. I%=1.N|%-1:SIG=|
|000010c0| 28 58 28 49 25 29 2d 58 | 28 49 25 2d 31 29 29 2f |(X(I%)-X|(I%-1))/|
|000010d0| 28 58 28 49 25 2b 31 29 | 2d 58 28 49 25 2d 31 29 |(X(I%+1)|-X(I%-1)|
|000010e0| 29 3a 50 3d 53 49 47 2a | 59 32 28 49 25 2d 31 29 |):P=SIG*|Y2(I%-1)|
|000010f0| 2b 32 0d 06 04 14 59 32 | 28 49 25 29 3d 28 53 49 |+2....Y2|(I%)=(SI|
|00001100| 47 2d 31 29 2f 50 0d 06 | 0e 6f 55 28 49 25 29 3d |G-1)/P..|.oU(I%)=|
|00001110| 28 36 2a 28 28 59 28 49 | 25 2b 31 29 2d 59 28 49 |(6*((Y(I|%+1)-Y(I|
|00001120| 25 29 29 2f 28 58 28 49 | 25 2b 31 29 2d 58 28 49 |%))/(X(I|%+1)-X(I|
|00001130| 25 29 29 2d 28 59 28 49 | 25 29 2d 59 28 49 25 2d |%))-(Y(I|%)-Y(I%-|
|00001140| 31 29 29 2f 28 58 28 49 | 25 29 2d 58 28 49 25 2d |1))/(X(I|%)-X(I%-|
|00001150| 31 29 29 29 2f 28 58 28 | 49 25 2b 31 29 2d 58 28 |1)))/(X(|I%+1)-X(|
|00001160| 49 25 2d 31 29 29 2d 53 | 49 47 2a 55 28 49 25 2d |I%-1))-S|IG*U(I%-|
|00001170| 31 29 29 2f 50 0d 06 18 | 05 ed 0d 06 22 47 51 4e |1))/P...|...."GQN|
|00001180| 3d 30 2e 35 3a 55 4e 3d | 28 33 2f 28 58 28 4e 25 |=0.5:UN=|(3/(X(N%|
|00001190| 29 2d 58 28 4e 25 2d 31 | 29 29 29 2a 28 59 50 4e |)-X(N%-1|)))*(YPN|
|000011a0| 2d 28 59 28 4e 25 29 2d | 59 28 4e 25 2d 31 29 29 |-(Y(N%)-|Y(N%-1))|
|000011b0| 2f 28 58 28 4e 25 29 2d | 58 28 4e 25 2d 31 29 29 |/(X(N%)-|X(N%-1))|
|000011c0| 29 0d 06 2c 2a 59 32 28 | 4e 25 29 3d 28 55 4e 2d |)..,*Y2(|N%)=(UN-|
|000011d0| 51 4e 2a 55 28 4e 25 2d | 31 29 29 2f 28 51 4e 2a |QN*U(N%-|1))/(QN*|
|000011e0| 59 32 28 4e 25 2d 31 29 | 2b 31 29 0d 06 36 34 e3 |Y2(N%-1)|+1)..64.|
|000011f0| 49 25 3d 4e 25 2d 31 b8 | 30 20 88 20 2d 31 3a 59 |I%=N%-1.|0 . -1:Y|
|00001200| 32 28 49 25 29 3d 59 32 | 28 49 25 29 2a 59 32 28 |2(I%)=Y2|(I%)*Y2(|
|00001210| 49 25 2b 31 29 2b 55 28 | 49 25 29 3a ed 3a e1 0d |I%+1)+U(|I%):.:..|
|00001220| 06 40 05 3a 0d 06 4a 43 | 3a f4 20 43 61 6c 63 75 |.@.:..JC|:. Calcu|
|00001230| 6c 61 74 65 20 69 6e 74 | 65 72 70 6f 6c 61 74 65 |late int|erpolate|
|00001240| 64 20 76 61 6c 75 65 20 | 67 69 76 65 6e 20 72 65 |d value |given re|
|00001250| 73 75 6c 74 20 6f 66 20 | 53 50 4c 49 4e 45 20 28 |sult of |SPLINE (|
|00001260| 4e 52 20 70 38 39 29 0d | 06 54 05 3a 0d 06 5e 23 |NR p89).|.T.:..^#|
|00001270| dd 20 f2 43 5f 53 50 4c | 49 4e 54 28 58 28 29 2c |. .C_SPL|INT(X(),|
|00001280| 59 28 29 2c 59 32 28 29 | 2c 58 2c f8 20 59 29 0d |Y(),Y2()|,X,. Y).|
|00001290| 06 68 37 ea 20 4b 4c 4f | 25 2c 4b 48 49 25 2c 4e |.h7. KLO|%,KHI%,N|
|000012a0| 25 2c 42 49 47 3a 4e 25 | 3d de 28 58 28 29 2c 31 |%,BIG:N%|=.(X(),1|
|000012b0| 29 3a 4b 4c 4f 25 3d 30 | 3a 4b 48 49 25 3d 4e 25 |):KLO%=0|:KHI%=N%|
|000012c0| 3a 42 49 47 3d b9 0d 06 | 72 0a c8 95 20 42 49 47 |:BIG=...|r... BIG|
|000012d0| 0d 06 7c 13 e7 20 4b 48 | 49 25 2d 4b 4c 4f 25 3e |..|.. KH|I%-KLO%>|
|000012e0| 31 20 8c 0d 06 86 2e 4b | 25 3d 28 4b 48 49 25 2b |1 .....K|%=(KHI%+|
|000012f0| 4b 4c 4f 25 29 2f 32 3a | e7 58 28 4b 25 29 3e 58 |KLO%)/2:|.X(K%)>X|
|00001300| 3a 4b 48 49 25 3d 4b 25 | 3a 8b 4b 4c 4f 25 3d 4b |:KHI%=K%|:.KLO%=K|
|00001310| 25 0d 06 90 0b cc 20 42 | 49 47 3d a3 0d 06 9a 05 |%..... B|IG=.....|
|00001320| cd 0d 06 a4 25 ce 20 20 | 3a f4 20 4b 4c 4f 25 20 |....%.  |:. KLO% |
|00001330| 61 6e 64 20 4b 48 49 25 | 20 6e 6f 77 20 62 72 61 |and KHI%| now bra|
|00001340| 63 6b 65 74 20 58 0d 06 | ae 15 48 3d 58 28 4b 48 |cket X..|..H=X(KH|
|00001350| 49 25 29 2d 58 28 4b 4c | 4f 25 29 0d 06 b8 2d e7 |I%)-X(KL|O%)...-.|
|00001360| 20 48 3d 30 3a 20 f1 20 | 22 58 20 76 61 6c 75 65 | H=0: . |"X value|
|00001370| 73 20 6d 75 73 74 20 62 | 65 20 64 69 73 74 69 6e |s must b|e distin|
|00001380| 63 74 2e 2e 2e 22 3a e1 | 0d 06 c2 23 41 3d 28 58 |ct...":.|...#A=(X|
|00001390| 28 4b 48 49 25 29 2d 58 | 29 2f 48 3a 42 3d 28 58 |(KHI%)-X|)/H:B=(X|
|000013a0| 2d 58 28 4b 4c 4f 25 29 | 29 2f 48 0d 06 cc 45 59 |-X(KLO%)|)/H...EY|
|000013b0| 3d 41 2a 59 28 4b 4c 4f | 25 29 2b 42 2a 59 28 4b |=A*Y(KLO|%)+B*Y(K|
|000013c0| 48 49 25 29 2b 28 28 41 | 5e 33 2d 41 29 2a 59 32 |HI%)+((A|^3-A)*Y2|
|000013d0| 28 4b 4c 4f 25 29 2b 28 | 42 5e 33 2d 42 29 2a 59 |(KLO%)+(|B^3-B)*Y|
|000013e0| 32 28 4b 48 49 25 29 29 | 2a 28 48 5e 32 29 2f 36 |2(KHI%))|*(H^2)/6|
|000013f0| 0d 06 d6 05 e1 0d 06 e0 | 05 3a 0d 06 f4 27 3a f4 |........|.:...':.|
|00001400| 20 32 2d 64 69 6d 65 6e | 73 69 6f 6e 61 6c 20 63 | 2-dimen|sional c|
|00001410| 61 73 65 20 28 4e 52 20 | 70 31 30 30 2d 31 30 31 |ase (NR |p100-101|
|00001420| 29 0d 06 fe 05 3a 0d 07 | 08 34 3a f4 20 43 6f 6d |)....:..|.4:. Com|
|00001430| 70 75 74 65 20 74 68 65 | 20 61 75 78 69 6c 69 61 |pute the| auxilia|
|00001440| 72 79 20 73 65 63 6f 6e | 64 2d 64 65 72 69 76 61 |ry secon|d-deriva|
|00001450| 74 69 76 65 20 74 61 62 | 6c 65 0d 07 12 05 3a 0d |tive tab|le....:.|
|00001460| 07 1c 26 dd 20 f2 43 32 | 5f 53 50 4c 49 4e 45 28 |..&. .C2|_SPLINE(|
|00001470| 58 31 28 29 2c 58 32 28 | 29 2c 59 28 29 2c f8 20 |X1(),X2(|),Y(),. |
|00001480| 59 32 28 29 29 0d 07 26 | 20 ea 20 4e 25 2c 4d 25 |Y2())..&| . N%,M%|
|00001490| 2c 59 54 4d 50 28 29 2c | 59 54 4d 50 32 28 29 2c |,YTMP(),|YTMP2(),|
|000014a0| 4a 25 2c 4b 25 0d 07 30 | 2e 4d 25 3d de 28 58 31 |J%,K%..0|.M%=.(X1|
|000014b0| 28 29 29 3a 4e 25 3d de | 28 58 32 28 29 29 3a de |()):N%=.|(X2()):.|
|000014c0| 20 59 54 4d 50 28 4d 25 | 29 2c 59 32 54 4d 50 28 | YTMP(M%|),Y2TMP(|
|000014d0| 4e 25 29 0d 07 3a 2a e3 | 20 4a 25 3d 30 b8 4d 25 |N%)..:*.| J%=0.M%|
|000014e0| 3a e3 4b 25 3d 30 b8 4e | 25 3a 59 54 4d 50 28 4b |:.K%=0.N|%:YTMP(K|
|000014f0| 25 29 3d 59 28 4a 25 2c | 4b 25 29 3a ed 0d 07 44 |%)=Y(J%,|K%):...D|
|00001500| 23 20 f2 43 5f 53 50 4c | 49 4e 45 28 58 32 28 29 |# .C_SPL|INE(X2()|
|00001510| 2c 59 54 4d 50 28 29 2c | 59 32 54 4d 50 28 29 29 |,YTMP(),|Y2TMP())|
|00001520| 0d 07 4e 23 20 e3 4b 25 | 3d 30 b8 4e 25 3a 59 32 |..N# .K%|=0.N%:Y2|
|00001530| 28 4b 25 2c 4a 25 29 3d | 59 32 54 4d 50 28 4b 25 |(K%,J%)=|Y2TMP(K%|
|00001540| 29 3a ed 0d 07 58 07 ed | 3a e1 0d 07 62 05 3a 0d |):...X..|:...b.:.|
|00001550| 07 6c 3e 3a f4 20 42 69 | 63 75 62 69 63 20 73 70 |.l>:. Bi|cubic sp|
|00001560| 6c 69 6e 65 20 69 6e 74 | 65 72 70 6f 6c 61 74 69 |line int|erpolati|
|00001570| 6f 6e 20 75 73 69 6e 67 | 20 72 65 73 75 6c 74 73 |on using| results|
|00001580| 20 6f 66 20 43 32 5f 53 | 50 4c 49 4e 45 0d 07 76 | of C2_S|PLINE..v|
|00001590| 05 3a 0d 07 80 2e dd 20 | f2 43 32 5f 53 50 4c 49 |.:..... |.C2_SPLI|
|000015a0| 4e 54 28 58 31 28 29 2c | 58 32 28 29 2c 59 28 29 |NT(X1(),|X2(),Y()|
|000015b0| 2c 59 32 28 29 2c 78 31 | 2c 78 32 2c f8 20 79 29 |,Y2(),x1|,x2,. y)|
|000015c0| 0d 07 8a 38 ea 20 4d 25 | 2c 4e 25 2c 4a 25 2c 4b |...8. M%|,N%,J%,K|
|000015d0| 25 2c 59 54 4d 50 28 29 | 2c 59 32 54 4d 50 28 29 |%,YTMP()|,Y2TMP()|
|000015e0| 2c 5a 3a 4d 25 3d de 28 | 58 31 28 29 29 3a 4e 25 |,Z:M%=.(|X1()):N%|
|000015f0| 3d de 28 58 32 28 29 29 | 0d 07 94 22 de 20 59 54 |=.(X2())|...". YT|
|00001600| 4d 50 28 4e 25 29 2c 59 | 32 54 4d 50 28 4e 25 29 |MP(N%),Y|2TMP(N%)|
|00001610| 2c 59 59 54 4d 50 28 4d | 25 29 0d 07 9e 3d e3 4a |,YYTMP(M|%)...=.J|
|00001620| 25 3d 30 b8 4d 25 3a e3 | 4b 25 3d 30 b8 4e 25 3a |%=0.M%:.|K%=0.N%:|
|00001630| 59 54 4d 50 28 4b 25 29 | 3d 59 28 4a 25 2c 4b 25 |YTMP(K%)|=Y(J%,K%|
|00001640| 29 3a 59 32 54 4d 50 28 | 4b 25 29 3d 59 32 28 4a |):Y2TMP(|K%)=Y2(J|
|00001650| 25 2c 4b 25 29 3a ed 0d | 07 a8 36 20 f2 43 5f 53 |%,K%):..|..6 .C_S|
|00001660| 50 4c 49 4e 54 28 58 32 | 28 29 2c 59 54 4d 50 28 |PLINT(X2|(),YTMP(|
|00001670| 29 2c 59 32 54 4d 50 28 | 29 2c 78 32 2c 5a 29 3a |),Y2TMP(|),x2,Z):|
|00001680| 59 59 54 4d 50 28 4a 25 | 29 3d 5a 3a ed 0d 07 b2 |YYTMP(J%|)=Z:....|
|00001690| 23 f2 43 5f 53 50 4c 49 | 4e 45 28 58 31 28 29 2c |#.C_SPLI|NE(X1(),|
|000016a0| 59 59 54 4d 50 28 29 2c | 59 32 54 4d 50 28 29 29 |YYTMP(),|Y2TMP())|
|000016b0| 0d 07 bc 2a f2 43 5f 53 | 50 4c 49 4e 54 28 58 31 |...*.C_S|PLINT(X1|
|000016c0| 28 29 2c 59 59 54 4d 50 | 28 29 2c 59 32 54 4d 50 |(),YYTMP|(),Y2TMP|
|000016d0| 28 29 2c 78 31 2c 79 29 | 3a e1 0d 07 c6 04 0d 07 |(),x1,y)|:.......|
|000016e0| d0 32 f4 20 20 20 20 20 | 20 20 20 20 20 20 20 2a |.2.     |       *|
|000016f0| 2a 2a 20 45 4e 44 20 4f | 46 20 43 55 42 49 43 20 |** END O|F CUBIC |
|00001700| 53 50 4c 49 4e 45 20 47 | 52 4f 55 50 20 2a 2a 2a |SPLINE G|ROUP ***|
|00001710| 0d 07 d1 04 0d 07 d2 30 | f4 20 20 20 20 20 20 20 |.......0|.       |
|00001720| 20 20 20 20 20 20 2a 2a | 2a 20 41 52 52 41 59 20 |      **|* ARRAY |
|00001730| 44 49 53 50 4c 41 59 20 | 52 4f 55 54 49 4e 45 53 |DISPLAY |ROUTINES|
|00001740| 20 2a 2a 2a 0d 07 da 04 | 0d 07 e4 05 3a 0d 07 e5 | ***....|....:...|
|00001750| 18 3a f4 20 70 72 69 6e | 74 20 72 65 61 6c 20 76 |.:. prin|t real v|
|00001760| 65 63 74 6f 72 0d 07 e6 | 05 3a 0d 08 34 27 dd f2 |ector...|.:..4'..|
|00001770| 50 52 49 4e 54 56 45 43 | 28 56 28 29 29 20 ea 20 |PRINTVEC|(V()) . |
|00001780| 49 25 2c 4a 25 3a 49 25 | 3d de 28 56 28 29 2c 31 |I%,J%:I%|=.(V(),1|
|00001790| 29 0d 08 3e 1a e3 4a 25 | 3d 30 b8 49 25 3a f1 56 |)..>..J%|=0.I%:.V|
|000017a0| 28 4a 25 29 3b 3a ed 3a | f1 3a e1 0d 08 48 05 3a |(J%);:.:|.:...H.:|
|000017b0| 0d 08 49 1b 3a f4 20 70 | 72 69 6e 74 20 69 6e 74 |..I.:. p|rint int|
|000017c0| 65 67 65 72 20 76 65 63 | 74 6f 72 0d 08 4a 05 3a |eger vec|tor..J.:|
|000017d0| 0d 08 52 2a dd f2 50 52 | 49 4e 54 56 45 43 49 28 |..R*..PR|INTVECI(|
|000017e0| 56 25 28 29 29 20 ea 20 | 49 25 2c 4a 25 3a 49 25 |V%()) . |I%,J%:I%|
|000017f0| 3d de 28 56 25 28 29 2c | 31 29 0d 08 5c 1b e3 4a |=.(V%(),|1)..\..J|
|00001800| 25 3d 30 b8 49 25 3a f1 | 56 25 28 4a 25 29 3b 3a |%=0.I%:.|V%(J%);:|
|00001810| ed 3a f1 3a e1 0d 08 66 | 05 3a 0d 08 67 18 3a f4 |.:.:...f|.:..g.:.|
|00001820| 20 70 72 69 6e 74 20 72 | 65 61 6c 20 6d 61 74 72 | print r|eal matr|
|00001830| 69 78 0d 08 68 05 3a 0d | 08 70 39 dd f2 50 52 49 |ix..h.:.|.p9..PRI|
|00001840| 4e 54 4d 41 54 28 4d 28 | 29 29 20 ea 20 49 25 2c |NTMAT(M(|)) . I%,|
|00001850| 4a 25 2c 4b 25 2c 4c 25 | 3a 49 25 3d de 28 4d 28 |J%,K%,L%|:I%=.(M(|
|00001860| 29 2c 31 29 3a 4a 25 3d | de 28 4d 28 29 2c 32 29 |),1):J%=|.(M(),2)|
|00001870| 0d 08 7a 28 e3 4b 25 3d | 30 b8 49 25 3a e3 4c 25 |..z(.K%=|0.I%:.L%|
|00001880| 3d 30 b8 4a 25 3a f1 4d | 28 4b 25 2c 4c 25 29 3b |=0.J%:.M|(K%,L%);|
|00001890| 3a ed 3a f1 3a ed 3a e1 | 0d 08 84 05 3a 0d 08 85 |:.:.:.:.|....:...|
|000018a0| 1b 3a f4 20 70 72 69 6e | 74 20 69 6e 74 65 67 65 |.:. prin|t intege|
|000018b0| 72 20 6d 61 74 72 69 78 | 0d 08 86 05 3a 0d 08 8e |r matrix|....:...|
|000018c0| 3d dd f2 50 52 49 4e 54 | 4d 41 54 49 28 4d 25 28 |=..PRINT|MATI(M%(|
|000018d0| 29 29 20 ea 20 49 25 2c | 4a 25 2c 4b 25 2c 4c 25 |)) . I%,|J%,K%,L%|
|000018e0| 3a 49 25 3d de 28 4d 25 | 28 29 2c 31 29 3a 4a 25 |:I%=.(M%|(),1):J%|
|000018f0| 3d de 28 4d 25 28 29 2c | 32 29 0d 08 98 29 e3 4b |=.(M%(),|2)...).K|
|00001900| 25 3d 30 b8 49 25 3a e3 | 4c 25 3d 30 b8 4a 25 3a |%=0.I%:.|L%=0.J%:|
|00001910| f1 4d 25 28 4b 25 2c 4c | 25 29 3b 3a ed 3a f1 3a |.M%(K%,L|%);:.:.:|
|00001920| ed 3a e1 0d 08 a2 27 f4 | 20 20 20 20 20 20 20 20 |.:....'.|        |
|00001930| 20 20 20 20 20 20 20 20 | 2a 2a 2a 20 41 4c 4c 20 |        |*** ALL |
|00001940| 44 4f 4e 45 20 21 20 2a | 2a 2a 0d ff             |DONE ! *|**..    |
+--------+-------------------------+-------------------------+--------+--------+