|
DataMuseum.dkPresents historical artifacts from the history of: RC4000/8000/9000 |
This is an automatic "excavation" of a thematic subset of
See our Wiki for more about RC4000/8000/9000 Excavated with: AutoArchaeologist - Free & Open Source Software. |
top - metrics - download
Length: 344064 (0x54000)
Types: Rc489kSubCat
Names: »algbib«
└─⟦621cfb9a2⟧ Bits:30002817 RC8000 Dump tape fra HCØ. Detaljer om "HC8000" projekt.
└─⟦0364f57e3⟧
└─⟦this⟧ »algbib«
0x00000…00033 (0,) Rc489kSubCatEnt {f00=0x0042b2, f01=0x00c382, f02=0x00c38b, filename=»algfuncfit «, entry_tail=0-4-0-15-1983-01-27T14:08-disc-0x0-0x0-0x0-0x0}
0x00033…00066 Rc489kSubCatEnt {f00=0x0134f2, f01=0x00c382, f02=0x00c38b, filename=»algfuncfit1 «, entry_tail=0-4-0-1-1983-01-27T15:32-disc-0x0-0x0-0x0-0x0}
0x00066…00099 Rc489kSubCatEnt {f00=0x01431a, f01=0x00c382, f02=0x00c38b, filename=»algupdate «, entry_tail=0-4-0-5-1983-01-27T13:37-disc-0x0-0x0-0x0-0x0}
0x00099…000cc Rc489kSubCatEnt {f00=0x019762, f01=0x00c382, f02=0x00c38b, filename=»alglinda «, entry_tail=0-4-0-5-1983-01-12T14:13-disc-0x0-0x0-0x0-0x0}
0x000cc…000ff Rc489kSubCatEnt {f00=0x01e000, f01=0x000000, f02=0x000000, filename=»algstruktur «, entry_tail=0-4-0-7-1982-11-01T14:11--0x0-0x0-0x0-0x0}
0x000ff…00132 Rc489kSubCatEnt {f00=0x025000, f01=0x000000, f02=0x000000, filename=»algsymosc «, entry_tail=0-4-0-6-1982-11-01T14:11--0x0-0x0-0x0-0x0}
0x00132…00165 Rc489kSubCatEnt {f00=0x02b000, f01=0x000000, f02=0x000000, filename=»algtdosc «, entry_tail=0-4-0-9-1982-11-01T14:11--0x0-0x0-0x0-0x0}
0x00165…00198 Rc489kSubCatEnt {f00=0x034000, f01=0x000000, f02=0x000000, filename=»algtextprg «, entry_tail=0-4-0-3-1982-11-01T14:11--0x0-0x0-0x0-0x0}
0x00198…001cb Rc489kSubCatEnt {f00=0x037000, f01=0x000000, f02=0x000000, filename=»algthermo «, entry_tail=0-4-0-6-1982-11-01T14:11--0x0-0x0-0x0-0x0}
0x001cb…001fe Rc489kSubCatEnt {f00=0x03d000, f01=0x000000, f02=0x000000, filename=»algvibrot «, entry_tail=0-4-0-35-1982-11-01T14:11--0x0-0x0-0x0-0x0}
0x001fe…00231 Rc489kSubCatEnt {f00=0x060000, f01=0x000000, f02=0x000000, filename=»asm «, entry_tail=0-4-0-59-1982-11-01T14:11--0x0-0x0-0x0-0x0}
0x00231…00264 Rc489kSubCatEnt {f00=0x09b000, f01=0x000000, f02=0x000000, filename=»algprsep «, entry_tail=0-4-0-2-1982-11-01T14:11--0x0-0x0-0x0-0x0}
0x00264…00297 Rc489kSubCatEnt {f00=0x09d000, f01=0x000000, f02=0x000000, filename=»algqrotfit «, entry_tail=0-4-0-13-1982-11-01T14:11--0x0-0x0-0x0-0x0}
0x00297…002ca Rc489kSubCatEnt {f00=0x0aa000, f01=0x000000, f02=0x000000, filename=»algrotcor «, entry_tail=0-4-0-4-1982-11-01T14:11--0x0-0x0-0x0-0x0}
0x002ca…002fd Rc489kSubCatEnt {f00=0x0ae000, f01=0x000000, f02=0x000000, filename=»algroteqq «, entry_tail=0-4-0-2-1982-11-01T14:11--0x0-0x0-0x0-0x0}
0x002fd…00300 0x000030
0x00300…00333 Rc489kSubCatEnt {f00=0x0b0000, f01=0x000000, f02=0x000000, filename=»algroteqq2 «, entry_tail=0-4-0-2-1982-11-01T14:11--0x0-0x0-0x0-0x0}
0x00333…00366 Rc489kSubCatEnt {f00=0x0b2000, f01=0x000000, f02=0x000000, filename=»algrotfit «, entry_tail=0-4-0-7-1982-11-01T14:11--0x0-0x0-0x0-0x0}
0x00366…00399 Rc489kSubCatEnt {f00=0x0b9000, f01=0x000000, f02=0x000000, filename=»algrotfit0 «, entry_tail=0-4-0-35-1982-11-01T14:11--0x0-0x0-0x0-0x0}
0x00399…003cc Rc489kSubCatEnt {f00=0x0dc000, f01=0x000000, f02=0x000000, filename=»algrotfrekv «, entry_tail=0-4-0-36-1982-11-01T14:11--0x0-0x0-0x0-0x0}
0x003cc…003ff Rc489kSubCatEnt {f00=0x100000, f01=0x000000, f02=0x000000, filename=»algrottabel «, entry_tail=0-4-0-2-1982-11-01T14:11--0x0-0x0-0x0-0x0}
0x003ff…00432 Rc489kSubCatEnt {f00=0x102000, f01=0x000000, f02=0x000000, filename=»algsetmove «, entry_tail=0-4-0-2-1982-11-01T14:11--0x0-0x0-0x0-0x0}
0x00432…00465 Rc489kSubCatEnt {f00=0x104000, f01=0x000000, f02=0x000000, filename=»algstark «, entry_tail=0-4-0-7-1982-11-01T14:11--0x0-0x0-0x0-0x0}
0x00465…00498 Rc489kSubCatEnt {f00=0x10b000, f01=0x000000, f02=0x000000, filename=»algstarkkoe «, entry_tail=0-4-0-1-1982-11-01T14:11--0x0-0x0-0x0-0x0}
0x00498…004cb Rc489kSubCatEnt {f00=0x10c000, f01=0x000000, f02=0x000000, filename=»algstarkplo «, entry_tail=0-4-0-8-1982-11-01T14:11--0x0-0x0-0x0-0x0}
0x004cb…004fe Rc489kSubCatEnt {f00=0x114000, f01=0x000000, f02=0x000000, filename=»alggaussbar «, entry_tail=0-4-0-7-1982-11-01T14:11--0x0-0x0-0x0-0x0}
0x004fe…00531 Rc489kSubCatEnt {f00=0x11b000, f01=0x000000, f02=0x000000, filename=»algimom «, entry_tail=0-4-0-2-1982-11-01T14:11--0x0-0x0-0x0-0x0}
0x00531…00564 Rc489kSubCatEnt {f00=0x11d000, f01=0x000000, f02=0x000000, filename=»algimominp «, entry_tail=0-4-0-5-1982-11-01T14:11--0x0-0x0-0x0-0x0}
0x00564…00597 Rc489kSubCatEnt {f00=0x122000, f01=0x000000, f02=0x000000, filename=»algisocomp «, entry_tail=0-4-0-5-1982-11-01T14:11--0x0-0x0-0x0-0x0}
0x00597…005ca Rc489kSubCatEnt {f00=0x127000, f01=0x000000, f02=0x000000, filename=»algkrmn «, entry_tail=0-4-0-6-1982-11-01T14:11--0x0-0x0-0x0-0x0}
0x005ca…005fd Rc489kSubCatEnt {f00=0x12d000, f01=0x000000, f02=0x000000, filename=»alglavib «, entry_tail=0-4-0-24-1982-11-01T14:11--0x0-0x0-0x0-0x0}
0x005fd…00600 00 00 00 ┆ ┆
0x00600…00633 Rc489kSubCatEnt {f00=0x145000, f01=0x000000, f02=0x000000, filename=»alglavib2 «, entry_tail=0-4-0-2-1982-11-01T14:11--0x0-0x0-0x0-0x0}
0x00633…00666 Rc489kSubCatEnt {f00=0x147000, f01=0x000000, f02=0x000000, filename=»alglavib3 «, entry_tail=0-4-0-2-1982-11-01T14:11--0x0-0x0-0x0-0x0}
0x00666…00699 Rc489kSubCatEnt {f00=0x149000, f01=0x000000, f02=0x000000, filename=»alglinbend «, entry_tail=0-4-0-3-1982-11-01T14:11--0x0-0x0-0x0-0x0}
0x00699…006cc Rc489kSubCatEnt {f00=0x14c000, f01=0x000000, f02=0x000000, filename=»algmathieu «, entry_tail=0-4-0-7-1982-11-01T14:11--0x0-0x0-0x0-0x0}
0x006cc…006ff Rc489kSubCatEnt {f00=0x153000, f01=0x000000, f02=0x000000, filename=»algmatpr «, entry_tail=0-4-0-2-1982-11-01T14:11--0x0-0x0-0x0-0x0}
0x006ff…00732 Rc489kSubCatEnt {f00=0x155000, f01=0x000000, f02=0x000000, filename=»algmsa «, entry_tail=0-4-0-3-1982-11-01T14:11--0x0-0x0-0x0-0x0}
0x00732…00765 Rc489kSubCatEnt {f00=0x158000, f01=0x000000, f02=0x000000, filename=»algplankrmn «, entry_tail=0-4-0-7-1982-11-01T14:11--0x0-0x0-0x0-0x0}
0x00765…00798 Rc489kSubCatEnt {f00=0x15f000, f01=0x000000, f02=0x000000, filename=»algbrinvert «, entry_tail=0-4-0-2-1982-11-01T14:10--0x0-0x0-0x0-0x0}
0x00798…007cb Rc489kSubCatEnt {f00=0x161000, f01=0x000000, f02=0x000000, filename=»algcoriol2 «, entry_tail=0-4-0-2-1982-11-01T14:10--0x0-0x0-0x0-0x0}
0x007cb…007fe Rc489kSubCatEnt {f00=0x163000, f01=0x000000, f02=0x000000, filename=»algcoriolis «, entry_tail=0-4-0-8-1982-11-01T14:10--0x0-0x0-0x0-0x0}
0x007fe…00831 Rc489kSubCatEnt {f00=0x16b000, f01=0x000000, f02=0x000000, filename=»algdiagtest «, entry_tail=0-4-0-3-1982-11-01T14:10--0x0-0x0-0x0-0x0}
0x00831…00864 Rc489kSubCatEnt {f00=0x16e000, f01=0x000000, f02=0x000000, filename=»algfjerdegr «, entry_tail=0-4-0-2-1982-11-01T14:10--0x0-0x0-0x0-0x0}
0x00864…00897 Rc489kSubCatEnt {f00=0x170000, f01=0x000000, f02=0x000000, filename=»algflinda «, entry_tail=0-4-0-47-1982-11-01T14:10--0x0-0x0-0x0-0x0}
0x00897…008ca Rc489kSubCatEnt {f00=0x19f000, f01=0x000000, f02=0x000000, filename=»algflinda2 «, entry_tail=0-4-0-2-1982-11-01T14:10--0x0-0x0-0x0-0x0}
0x008ca…008fd Rc489kSubCatEnt {f00=0x1a1000, f01=0x000000, f02=0x000000, filename=»algflinda3 «, entry_tail=0-4-0-2-1982-11-01T14:10--0x0-0x0-0x0-0x0}
0x008fd…00900 00 00 00 ┆ ┆
0x00900…00933 Rc489kSubCatEnt {f00=0x1a3000, f01=0x000000, f02=0x000000, filename=»algfortrat «, entry_tail=0-4-0-12-1982-11-01T14:10--0x0-0x0-0x0-0x0}
0x00933…00966 Rc489kSubCatEnt {f00=0x1af000, f01=0x000000, f02=0x000000, filename=»algfrekvfor «, entry_tail=0-4-0-1-1982-11-01T14:10--0x0-0x0-0x0-0x0}
0x00966…00999 Rc489kSubCatEnt {f00=0x1b0000, f01=0x000000, f02=0x000000, filename=»algfuncdraw «, entry_tail=0-4-0-8-1982-11-01T14:10--0x0-0x0-0x0-0x0}
0x00999…009a0 14 c0 00 00 00 00 00 ┆ ┆
0x009a0…009c0 00 00 61 6c 67 6d 61 74 68 69 65 75 00 00 00 00 07 00 00 00 00 00 00 00 00 00 00 00 00 88 41 1b ┆ algmathieu A ┆
0x009c0…009e0 00 00 00 00 00 00 00 00 00 00 00 00 15 30 00 00 00 00 00 00 00 61 6c 67 6d 61 74 70 72 00 00 00 ┆ 0 algmatpr ┆
0x009e0…00a00 00 00 00 02 00 00 00 00 00 00 00 00 00 00 00 00 88 41 1b 00 00 00 00 00 00 00 00 00 00 00 00 15 ┆ A ┆
0x00a00…00a20 50 00 00 00 00 00 00 00 61 6c 67 6d 73 61 00 00 00 00 00 00 00 00 03 00 00 00 00 00 00 00 00 00 ┆P algmsa ┆
0x00a20…00a40 00 00 00 88 41 1b 00 00 00 00 00 00 00 00 00 00 00 00 15 80 00 00 00 00 00 00 00 61 6c 67 70 6c ┆ A algpl┆
0x00a40…00a60 61 6e 6b 72 6d 6e 00 00 00 07 00 00 00 00 00 00 00 00 00 00 00 00 88 41 1b 00 00 00 00 00 00 00 ┆ankrmn A ┆
0x00a60…00a80 00 00 00 00 00 15 f0 00 00 00 00 00 00 00 61 6c 67 62 72 69 6e 76 65 72 74 00 00 00 02 00 00 00 ┆ algbrinvert ┆
0x00a80…00aa0 00 00 00 00 00 00 00 00 00 88 41 1a 00 00 00 00 00 00 00 00 00 00 00 00 16 10 00 00 00 00 00 00 ┆ A ┆
0x00aa0…00ac0 00 61 6c 67 63 6f 72 69 6f 6c 32 00 00 00 00 02 00 00 00 00 00 00 00 00 00 00 00 00 88 41 1a 00 ┆ algcoriol2 A ┆
0x00ac0…00ae0 00 00 00 00 00 00 00 00 00 00 00 16 30 00 00 00 00 00 00 00 61 6c 67 63 6f 72 69 6f 6c 69 73 00 ┆ 0 algcoriolis ┆
0x00ae0…00b00 00 00 08 00 00 00 00 00 00 00 00 00 00 00 00 88 41 1a 00 00 00 00 00 00 00 00 00 00 00 00 16 b0 ┆ A ┆
0x00b00…00b20 00 00 00 00 00 00 00 61 6c 67 64 69 61 67 74 65 73 74 00 00 00 03 00 00 00 00 00 00 00 00 00 00 ┆ algdiagtest ┆
0x00b20…00b40 00 00 88 41 1a 00 00 00 00 00 00 00 00 00 00 00 00 16 e0 00 00 00 00 00 00 00 61 6c 67 66 6a 65 ┆ A algfje┆
0x00b40…00b60 72 64 65 67 72 00 00 00 02 00 00 00 00 00 00 00 00 00 00 00 00 88 41 1a 00 00 00 00 00 00 00 00 ┆rdegr A ┆
0x00b60…00b80 00 00 00 00 17 00 00 00 00 00 00 00 00 61 6c 67 66 6c 69 6e 64 61 00 00 00 00 00 2f 00 00 00 00 ┆ algflinda / ┆
0x00b80…00ba0 00 00 00 00 00 00 00 00 88 41 1a 00 00 00 00 00 00 00 00 00 00 00 00 19 f0 00 00 00 00 00 00 00 ┆ A ┆
0x00ba0…00bc0 61 6c 67 66 6c 69 6e 64 61 32 00 00 00 00 02 00 00 00 00 00 00 00 00 00 00 00 00 88 41 1a 00 00 ┆algflinda2 A ┆
0x00bc0…00be0 00 00 00 00 00 00 00 00 00 00 1a 10 00 00 00 00 00 00 00 61 6c 67 66 6c 69 6e 64 61 33 00 00 00 ┆ algflinda3 ┆
0x00be0…00c00 00 02 00 00 00 00 00 00 00 00 00 00 00 00 88 41 1a 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 ┆ A ┆
0x00c00…03900 (1,) Data for algfuncfit
0x03900…03c00 Data for algfuncfit1
0x03c00…04b00 (5,) Data for algupdate
0x04b00…05a00 Data for alglinda
0x05a00…06f00 Data for algstruktur
0x06f00…08100 Data for algsymosc
0x08100…09c00 Data for algtdosc
0x09c00…0a500 (13,) Data for algtextprg
0x0a500…0b700 Data for algthermo
0x0b700…12000 Data for algvibrot
0x12000…1d100 (24,) Data for asm
0x1d100…1d700 Data for algprsep
0x1d700…1fe00 Data for algqrotfit
0x1fe00…20a00 Data for algrotcor
0x20a00…21000 Data for algroteqq
0x21000…21600 (44,) Data for algroteqq2
0x21600…22b00 Data for algrotfit
0x22b00…29400 Data for algrotfit0
0x29400…30000 (55,) Data for algrotfrekv
0x30000…30600 (64,) Data for algrottabel
0x30600…30c00 Data for algsetmove
0x30c00…32100 (65,) Data for algstark
0x32100…32400 Data for algstarkkoe
0x32400…33c00 (67,) Data for algstarkplo
0x33c00…35100 (69,) Data for alggaussbar
0x35100…35700 Data for algimom
0x35700…36600 Data for algimominp
0x36600…37500 Data for algisocomp
0x37500…38700 Data for algkrmn
0x38700…3cf00 Data for alglavib
0x3cf00…3d500 Data for alglavib2
0x3d500…3db00 Data for alglavib3
0x3db00…3e400 Data for alglinbend
0x3e400…3f900 (83,) Data for algmathieu
0x3f900…3ff00 Data for algmatpr
0x3ff00…40800 Data for algmsa
0x40800…41d00 (86,) Data for algplankrmn
0x41d00…42300 Data for algbrinvert
0x42300…42900 Data for algcoriol2
0x42900…44100 Data for algcoriolis
0x44100…44a00 Data for algdiagtest
0x44a00…45000 Data for algfjerdegr
0x45000…4dd00 (92,) Data for algflinda
0x4dd00…4e300 Data for algflinda2
0x4e300…4e900 Data for algflinda3
0x4e900…50d00 Data for algfortrat
0x50d00…51000 Data for algfrekvfor
0x51000…52800 (108,) Data for algfuncdraw
0x52800…52820 (110,) 2c 72 2f 2b 2f 2d 2f 2c 72 2f 2f 20 20 2f 2c 6c 31 2c 69 2f 0a 20 20 20 20 20 20 65 6e 64 3b 0a ┆,r/+/-/,r// /,l1,i/ end; ┆
0x52820…52840 2f 2c 6c 2e 2f 4a 28 69 2c 6a 29 2f 2c 6c 2d 32 2c 69 2f 0a 69 66 20 4e 54 3d 30 20 74 68 65 6e ┆/,l./J(i,j)/,l-2,i/ if NT=0 then┆
0x52840…52860 20 62 65 67 69 6e 0a 2f 2c 6c 2e 2f 73 79 6d 69 6e 76 65 72 73 65 2f 2c 69 2f 0a 65 6e 64 20 65 ┆ begin /,l./syminverse/,i/ end e┆
0x52860…52880 6c 73 65 20 62 65 67 69 6e 0a 66 6f 72 20 69 3a 3d 31 20 73 74 65 70 20 31 20 75 6e 74 69 6c 20 ┆lse begin for i:=1 step 1 until ┆
0x52880…528a0 33 20 64 6f 0a 66 6f 72 20 6a 3a 3d 69 20 73 74 65 70 20 31 20 75 6e 74 69 6c 20 33 20 64 6f 0a ┆3 do for j:=i step 1 until 3 do ┆
0x528a0…528c0 4a 28 69 2c 6a 29 3a 3d 20 2d 73 75 6d 28 6d 61 73 73 28 6b 2b 31 29 2a 58 70 28 33 2a 6b 2b 69 ┆J(i,j):= -sum(mass(k+1)*Xp(3*k+i┆
0x528c0…528e0 29 2a 58 28 33 2a 6b 2b 6a 29 2c 6b 2c 30 2c 4e 2d 31 2d 4e 54 29 3b 0a 49 61 3a 3d 20 4a 28 31 ┆)*X(3*k+j),k,0,N-1-NT); Ia:= J(1┆
0x528e0…52900 2c 31 29 3b 20 49 62 3a 3d 20 4a 28 32 2c 32 29 3b 20 49 63 3a 3d 20 4a 28 33 2c 33 29 3b 0a 4a ┆,1); Ib:= J(2,2); Ic:= J(3,3); J┆
0x52900…52920 28 31 2c 31 29 3a 3d 20 2d 49 62 2d 49 63 3b 0a 4a 28 32 2c 32 29 3a 3d 20 2d 49 61 2d 49 63 3b ┆(1,1):= -Ib-Ic; J(2,2):= -Ia-Ic;┆
0x52920…52940 0a 4a 28 33 2c 33 29 3a 3d 20 2d 49 61 2d 49 62 3b 0a 65 6e 64 3b 0a 2f 2c 6c 2e 2f 2d 63 28 6d ┆ J(3,3):= -Ia-Ib; end; /,l./-c(m┆
0x52940…52960 29 2f 2c 72 2f 2d 63 2f 63 2f 2c 66 0a 69 20 6c 61 76 69 62 74 65 78 74 0a 19 00 00 2c 31 29 3b ┆)/,r/-c/c/,f i lavibtext ,1);┆
0x52960…52980 20 49 62 3a 3d 20 4a 28 32 2c 32 29 3b 20 49 63 3a 3d 20 4a 28 33 2c 33 29 3b 0a 4a 28 31 2c 31 ┆ Ib:= J(2,2); Ic:= J(3,3); J(1,1┆
0x52980…529a0 29 3a 3d 20 2d 49 62 2d 49 63 3b 0a 4a 28 32 2c 32 29 3a 3d 20 2d 49 61 2d 49 63 3b 0a 4a 28 33 ┆):= -Ib-Ic; J(2,2):= -Ia-Ic; J(3┆
0x529a0…529c0 2c 33 29 3a 3d 20 2d 49 61 2d 49 62 3b 0a 65 6e 64 3b 0a 2f 2c 6c 2e 2f 2d 63 28 6d 29 2f 2c 72 ┆,3):= -Ia-Ib; end; /,l./-c(m)/,r┆
0x529c0…529e0 2f 2d 63 2f 63 2f 0a 6c 33 2c 72 2f 20 4a 2f 0a 20 20 20 69 66 20 4e 54 3d 30 20 74 68 65 6e 20 ┆/-c/c/ l3,r/ J/ if NT=0 then ┆
0x529e0…52a00 4a 2f 2c 72 2f 3b 2f 20 65 6c 73 65 20 73 75 6d 28 4a 28 6d 2c 69 29 2a 49 52 28 69 2c 6b 29 2c ┆J/,r/;/ else sum(J(m,i)*IR(i,k),┆
0x52a00…52a20 69 2c 31 2c 33 29 3b 2f 2c 66 0a 69 20 6c 61 76 69 62 74 65 78 74 0a 19 00 2e 64 64 64 64 3e 2c ┆i,1,3);/,f i lavibtext .dddd>,┆
0x52a20…52a40 47 46 28 6d 2c 6a 2b 33 29 2a 35 38 33 37 2e 35 38 36 37 29 3b 0a 65 6e 64 3b 0a 77 72 69 74 65 ┆GF(m,j+3)*5837.5867); end; write┆
0x52a40…52a60 28 72 65 73 2c 6e 6c 2c 33 2c 3c 3a 47 2d 4d 41 54 52 49 58 3a 3a 3e 2c 6e 6c 2c 31 29 3b 0a 66 ┆(res,nl,3,<:G-MATRIX::>,nl,1); f┆
0x52a60…52a80 6f 72 20 69 3a 3d 31 20 73 74 65 70 20 31 20 75 6e 74 69 6c 20 4e 46 20 64 6f 20 62 65 67 69 6e ┆or i:=1 step 1 until NF do begin┆
0x52a80…52aa0 20 77 72 69 74 65 28 72 65 73 2c 6e 6c 2c 31 29 3b 0a 66 6f 72 20 6a 3a 3d 31 20 73 74 65 70 20 ┆ write(res,nl,1); for j:=1 step ┆
0x52aa0…52ac0 31 20 75 6e 74 69 6c 20 69 20 64 6f 0a 20 20 20 77 72 69 74 65 28 72 65 73 2c 3c 3c 5f 5f 2d 64 ┆1 until i do write(res,<<__-d┆
0x52ac0…52ae0 2e 64 64 64 64 64 64 3e 2c 47 46 28 69 2b 33 2c 6a 2b 33 29 29 0a 65 6e 64 20 65 6e 64 20 65 6e ┆.dddddd>,GF(i+3,j+3)) end end en┆
0x52ae0…52b00 64 20 47 2c 5a 3b 0a 65 6e 64 20 49 74 3b 0a 0a 69 66 20 70 3d 30 20 74 68 65 6e 20 67 6f 74 6f ┆d G,Z; end It; if p=0 then goto┆
0x52b00…52b20 3b 6b 65 6d 6c 61 62 35 20 31 0a 6c 6f 6f 6b 75 70 20 62 61 6e 64 64 69 61 67 0a 69 66 20 6f 6b ┆;kemlab5 1 lookup banddiag if ok┆
0x52b20…52b40 2e 6e 6f 0a 28 62 61 6e 64 64 69 61 67 3d 20 73 65 74 20 36 0a 62 61 6e 64 64 69 61 67 3d 61 6c ┆.no (banddiag= set 6 banddiag=al┆
0x52b40…52b60 67 6f 6c 20 65 78 74 62 61 6e 64 64 69 61 67 20 69 6e 64 65 78 2e 6e 6f 29 0a 6c 69 6e 62 65 6e ┆gol extbanddiag index.no) linben┆
0x52b60…52b80 64 3d 73 65 74 20 34 30 0a 70 65 72 6d 61 6e 65 6e 74 20 6c 69 6e 62 65 6e 64 2e 31 35 0a 6c 69 ┆d=set 40 permanent linbend.15 li┆
0x52b80…52ba0 6e 62 65 6e 64 3d 61 6c 67 6f 6c 20 69 6e 64 65 78 2e 6e 6f 0a 0c 0a 0a 0a 0a 6c 69 6e 62 65 6e ┆nbend=algol index.no linben┆
0x52ba0…52bc0 64 0a 0a 62 65 67 69 6e 0a 63 6f 6d 6d 65 6e 74 20 69 6e 70 75 74 28 6f 76 65 72 73 6b 72 69 66 ┆d begin comment input(overskrif┆
0x52bc0…52be0 74 2c 67 2c 73 31 2c 73 32 2c 76 6d 61 78 2c 6d 61 74 72 69 78 6f 72 64 65 6e 20 6e 29 3b 0a 69 ┆t,g,s1,s2,vmax,matrixorden n); i┆
0x52be0…52c00 6e 74 65 67 65 72 20 76 2c 6c 2c 76 6d 61 78 2c 6e 2c 69 2c 6b 3b 0a 72 65 61 6c 20 67 2c 73 31 ┆nteger v,l,vmax,n,i,k; real g,s1┆
0x52c00…52c20 2c 73 32 2c 73 71 31 2c 73 71 32 2c 61 31 2c 61 32 2c 62 31 2c 62 32 2c 63 31 2c 63 32 3b 0a 62 ┆,s2,sq1,sq2,a1,a2,b1,b2,c1,c2; b┆
0x52c20…52c40 6f 6f 6c 65 61 6e 20 73 70 2c 6e 6c 3b 0a 61 72 72 61 79 20 68 65 61 64 28 31 3a 31 32 29 3b 0a ┆oolean sp,nl; array head(1:12); ┆
0x52c40…52c60 73 70 3a 3d 66 61 6c 73 65 20 61 64 64 20 33 32 3b 20 6e 6c 3a 3d 20 66 61 6c 73 65 20 61 64 64 ┆sp:=false add 32; nl:= false add┆
0x52c60…52c80 20 31 30 3b 0a 72 65 61 64 68 65 61 64 28 69 6e 2c 68 65 61 64 2c 31 29 3b 0a 72 65 61 64 28 69 ┆ 10; readhead(in,head,1); read(i┆
0x52c80…52ca0 6e 2c 67 2c 73 31 2c 73 32 2c 76 6d 61 78 2c 6e 29 3b 20 69 3a 3d 20 31 3b 0a 77 72 69 74 65 28 ┆n,g,s1,s2,vmax,n); i:= 1; write(┆
0x52ca0…52cc0 6f 75 74 2c 3c 3a 3c 31 32 3e 3a 3e 2c 6e 6c 2c 33 2c 73 74 72 69 6e 67 20 68 65 61 64 28 69 6e ┆out,<:<12>:>,nl,3,string head(in┆
0x52cc0…52ce0 63 72 65 61 73 65 28 69 29 29 2c 6e 6c 2c 31 2c 0a 20 20 20 20 20 20 3c 3a 6c 69 6e 65 61 72 20 ┆crease(i)),nl,1, <:linear ┆
0x52ce0…52d00 62 65 6e 64 65 72 3a 20 20 67 20 3d 3a 3e 2c 3c 3c 64 64 64 64 2e 64 64 64 64 64 3e 2c 67 2c 6e ┆bender: g =:>,<<dddd.ddddd>,g,n┆
0x52d00…52d20 6c 2c 31 2c 0a 20 20 20 20 20 20 3c 3a 72 65 64 2e 20 62 61 72 72 69 65 72 65 2c 20 73 31 20 3d ┆l,1, <:red. barriere, s1 =┆
0x52d20…52d40 20 3a 3e 2c 73 31 2c 3c 3a 20 20 73 32 20 3d 20 3a 3e 2c 73 32 29 3b 0a 76 3a 3d 20 28 76 6d 61 ┆ :>,s1,<: s2 = :>,s2); v:= (vma┆
0x52d40…52d60 78 2b 32 29 2f 2f 32 3b 0a 0a 62 65 67 69 6e 20 72 65 61 6c 20 61 2c 62 2c 63 2c 77 3b 0a 77 3a ┆x+2)//2; begin real a,b,c,w; w:┆
0x52d60…52d80 3d 20 73 31 2b 73 32 3b 20 61 3a 3d 20 28 73 31 2b 34 2a 73 32 29 2f 28 31 32 2a 77 29 3b 20 62 ┆= s1+s2; a:= (s1+4*s2)/(12*w); b┆
0x52d80…52da0 3a 3d 20 28 73 31 2b 31 36 2a 73 32 29 2f 28 32 34 2a 77 29 3b 0a 63 3a 3d 20 28 61 2b 61 2d 67 ┆:= (s1+16*s2)/(24*w); c:= (a+a-g┆
0x52da0…52dc0 29 2f 38 3b 20 77 3a 3d 20 73 71 72 74 28 77 2f 32 29 3b 0a 61 31 3a 3d 20 2d 31 2f 36 2d 61 3b ┆)/8; w:= sqrt(w/2); a1:= -1/6-a;┆
0x52dc0…52de0 20 62 31 3a 3d 20 28 2d 36 2a 61 2d 67 29 2f 38 3b 20 63 31 3a 3d 20 31 2f 36 2b 63 3b 20 61 3a ┆ b1:= (-6*a-g)/8; c1:= 1/6+c; a:┆
0x52de0…52e00 3d 20 61 2a 61 3b 20 63 3a 3d 20 63 2a 63 3b 0a 61 32 3a 3d 20 28 2d 31 2f 31 32 30 2b 67 2f 31 ┆= a*a; c:= c*c; a2:= (-1/120+g/1┆
0x52e00…52e20 32 2b 62 2f 35 2d 61 2d 61 2d 34 2a 63 29 2f 77 20 2b 20 77 3b 0a 62 32 3a 3d 20 28 62 2f 31 32 ┆2+b/5-a-a-4*c)/w + w; b2:= (b/12┆
0x52e20…52e40 2d 61 2d 63 29 2f 77 3b 0a 63 32 3a 3d 20 28 31 2f 33 30 2d 67 2f 31 32 2d 62 2f 32 30 2b 61 2f ┆-a-c)/w; c2:= (1/30-g/12-b/20+a/┆
0x52e40…52e60 32 2b 63 29 2f 77 0a 65 6e 64 3b 0a 0a 62 65 67 69 6e 20 61 72 72 61 79 20 68 28 31 3a 6e 2c 30 ┆2+c)/w end; begin array h(1:n,0┆
0x52e60…52e80 3a 69 66 20 73 32 3d 30 20 74 68 65 6e 20 31 20 65 6c 73 65 20 32 29 2c 20 65 76 28 31 3a 76 29 ┆:if s2=0 then 1 else 2), ev(1:v)┆
0x52e80…52ea0 2c 20 78 28 31 3a 76 2c 31 3a 6e 29 3b 0a 63 6f 6d 6d 65 6e 74 0a 20 20 20 20 20 20 20 20 20 68 ┆, x(1:v,1:n); comment h┆
0x52ea0…52ec0 28 69 2c 30 29 3a 3d 20 3c 76 2c 6c 21 20 68 20 21 76 2c 6c 3e 20 2c 0a 20 20 20 20 20 20 20 20 ┆(i,0):= <v,l! h !v,l> , ┆
0x52ec0…52ee0 20 68 28 69 2c 31 29 3a 3d 20 3c 76 2c 6c 21 20 68 20 21 76 2d 31 2c 6c 3e 2c 0a 20 20 20 20 20 ┆ h(i,1):= <v,l! h !v-1,l>, ┆
0x52ee0…52f00 20 20 20 20 68 28 69 2c 32 29 3a 3d 20 3c 76 2c 6c 21 20 68 20 21 76 2d 32 2c 6c 3e 20 77 69 74 ┆ h(i,2):= <v,l! h !v-2,l> wit┆
0x52f00…52f20 68 20 76 20 3d 20 69 2b 6c 2d 31 3b 0a 0a 66 6f 72 20 6c 3a 3d 30 20 73 74 65 70 20 31 20 75 6e ┆h v = i+l-1; for l:=0 step 1 un┆
0x52f20…52f40 74 69 6c 20 76 6d 61 78 20 64 6f 20 62 65 67 69 6e 0a 73 71 31 3a 3d 20 30 3b 0a 66 6f 72 20 69 ┆til vmax do begin sq1:= 0; for i┆
0x52f40…52f60 3a 3d 20 31 20 73 74 65 70 20 31 20 75 6e 74 69 6c 20 6e 20 64 6f 20 62 65 67 69 6e 0a 20 20 20 ┆:= 1 step 1 until n do begin ┆
0x52f60…52f80 76 3a 3d 20 69 2b 6c 2d 31 3b 20 73 71 32 3a 3d 20 73 71 31 3b 20 73 71 31 3a 3d 20 73 71 72 74 ┆v:= i+l-1; sq2:= sq1; sq1:= sqrt┆
0x52f80…52fa0 28 28 76 2b 6c 29 2a 28 76 2d 6c 29 2f 28 28 76 2b 76 2b 31 29 2a 28 76 2b 76 2d 31 29 29 29 3b ┆((v+l)*(v-l)/((v+v+1)*(v+v-1)));┆
0x52fa0…52fc0 0a 20 20 20 68 28 69 2c 30 29 3a 3d 20 28 28 31 2d 67 29 2a 76 2a 28 76 2b 31 29 20 2b 20 73 31 ┆ h(i,0):= ((1-g)*v*(v+1) + s1┆
0x52fc0…52fe0 0a 20 20 20 20 20 20 20 20 20 20 20 20 2b 20 28 76 2a 28 76 2b 31 29 2b 6c 2a 6c 2d 31 29 2f 28 ┆ + (v*(v+1)+l*l-1)/(┆
0x52fe0…53000 28 76 2b 76 2d 31 29 2a 28 76 2b 76 2b 33 29 29 2a 73 32 29 2f 32 3b 0a 20 20 20 68 28 69 2c 31 ┆(v+v-1)*(v+v+3))*s2)/2; h(i,1┆
0x53000…53020 29 3a 3d 20 73 71 31 2a 28 67 2a 76 2a 76 20 2d 20 73 31 29 2f 32 3b 0a 20 20 20 69 66 20 73 32 ┆):= sq1*(g*v*v - s1)/2; if s2┆
0x53020…53040 3c 3e 30 20 74 68 65 6e 20 68 28 69 2c 32 29 3a 3d 20 2d 20 73 71 31 2a 73 71 32 2a 73 32 2f 34 ┆<>0 then h(i,2):= - sq1*sq2*s2/4┆
0x53040…53060 3b 0a 65 6e 64 3b 0a 6b 3a 3d 20 28 76 6d 61 78 2b 32 2d 6c 29 2f 2f 32 3b 0a 62 61 6e 64 64 69 ┆; end; k:= (vmax+2-l)//2; banddi┆
0x53060…53080 61 67 28 6e 2c 69 66 20 73 32 3d 30 20 74 68 65 6e 20 32 20 65 6c 73 65 20 33 2c 31 2c 6b 2c 65 ┆ag(n,if s2=0 then 2 else 3,1,k,e┆
0x53080…530a0 76 2c 68 2c 78 29 3b 0a 77 72 69 74 65 28 6f 75 74 2c 6e 6c 2c 32 2c 73 70 2c 35 2a 6c 29 3b 0a ┆v,h,x); write(out,nl,2,sp,5*l); ┆
0x530a0…530c0 66 6f 72 20 69 3a 3d 31 20 73 74 65 70 20 31 20 75 6e 74 69 6c 20 6b 20 64 6f 0a 20 20 20 77 72 ┆for i:=1 step 1 until k do wr┆
0x530c0…530e0 69 74 65 28 6f 75 74 2c 3c 3c 2d 64 64 64 2e 64 64 64 64 64 3e 2c 65 76 28 69 29 29 3b 0a 77 72 ┆ite(out,<<-ddd.ddddd>,ev(i)); wr┆
0x530e0…53100 69 74 65 28 6f 75 74 2c 6e 6c 2c 31 2c 73 70 2c 35 2a 6c 29 3b 20 20 76 3a 3d 20 6c 3b 0a 66 6f ┆ite(out,nl,1,sp,5*l); v:= l; fo┆
0x53100…53120 72 20 69 3a 3d 31 20 73 74 65 70 20 31 20 75 6e 74 69 6c 20 6b 20 64 6f 20 62 65 67 69 6e 0a 20 ┆r i:=1 step 1 until k do begin ┆
0x53120…53140 20 20 77 72 69 74 65 28 6f 75 74 2c 3c 3c 2d 64 64 64 2e 64 64 64 64 64 3e 2c 0a 20 20 20 20 20 ┆ write(out,<<-ddd.ddddd>, ┆
0x53140…53160 20 28 76 2b 31 29 2a 28 61 32 2b 62 32 2a 76 2a 28 76 2b 32 29 2b 63 32 2a 6c 2a 6c 29 20 2b 20 ┆ (v+1)*(a2+b2*v*(v+2)+c2*l*l) + ┆
0x53160…53180 61 31 2b 62 31 2a 76 2a 28 76 2b 32 29 2b 63 31 2a 6c 2a 6c 29 3b 0a 20 20 20 76 3a 3d 20 76 2b ┆a1+b1*v*(v+2)+c1*l*l); v:= v+┆
0x53180…531a0 32 20 65 6e 64 0a 65 6e 64 20 6c 0a 65 6e 64 20 68 0a 65 6e 64 0a 19 00 00 29 3b 0a 66 6f 72 20 ┆2 end end l end h end ); for ┆
0x531a0…531c0 70 3a 3d 20 70 2d 31 20 77 68 69 6c 65 20 70 3e 3d 30 20 64 6f 20 62 65 67 69 6e 0a 20 20 20 72 ┆p:= p-1 while p>=0 do begin r┆
0x531c0…531e0 65 61 64 28 69 6e 2c 71 2c 76 29 3b 20 6d 3a 3d 20 74 28 71 2c 31 29 3b 20 67 28 71 29 3a 3d 20 ┆ead(in,q,v); m:= t(q,1); g(q):= ┆
0x531e0…53200 69 66 20 6d 3d 31 20 74 68 65 6e 20 76 20 65 6c 73 65 20 76 2a 70 69 2f 31 38 30 3b 0a 20 20 20 ┆if m=1 then v else v*pi/180; ┆
0x53200…53220 69 3a 3d 20 74 28 71 2c 32 29 3b 20 6a 3a 3d 20 74 28 71 2c 33 29 3b 20 6b 3a 3d 20 74 28 71 2c ┆i:= t(q,2); j:= t(q,3); k:= t(q,┆
0x53220…53240 34 29 3b 20 6c 3a 3d 20 74 28 71 2c 35 29 3b 0a 20 20 20 0a 20 20 20 77 72 69 74 65 28 72 65 73 ┆4); l:= t(q,5); write(res┆
0x53240…53260 2c 6e 6c 2c 31 2c 63 61 73 65 20 6d 20 6f 66 20 28 3c 3a 64 69 73 74 61 6e 63 65 20 3a 3e 2c 3c ┆,nl,1,case m of (<:distance :>,<┆
0x53260…53280 3a 61 6e 67 6c 65 20 3a 3e 2c 0a 20 20 20 20 20 20 20 20 20 20 20 20 3c 3a 6f 75 74 6f 20 3a 3e ┆:angle :>, <:outo :>┆
0x53280…532a0 2c 3c 3a 74 6f 72 73 20 3a 3e 2c 3c 3a 73 75 6d 20 6f 66 20 3a 3e 29 29 3b 0a 20 20 20 77 72 69 ┆,<:tors :>,<:sum of :>)); wri┆
0x532a0…532c0 74 65 28 72 65 73 2c 3c 3c 2d 64 3e 2c 69 2c 2d 6a 29 3b 20 69 66 20 6d 3d 32 20 74 68 65 6e 20 ┆te(res,<<-d>,i,-j); if m=2 then ┆
0x532c0…532e0 77 72 69 74 65 28 72 65 73 2c 3c 3c 2d 64 3e 2c 2d 6b 29 3b 0a 20 20 20 69 66 20 6d 3d 33 20 6f ┆write(res,<<-d>,-k); if m=3 o┆
0x532e0…53300 72 20 6d 3d 34 20 74 68 65 6e 20 77 72 69 74 65 28 72 65 73 2c 3c 3c 2d 64 3e 2c 2d 6b 2c 2d 6c ┆r m=4 then write(res,<<-d>,-k,-l┆
0x53300…53320 29 3b 0a 20 20 20 69 66 20 6d 3c 35 20 74 68 65 6e 20 77 72 69 74 65 28 72 65 73 2c 3c 3a 20 3d ┆); if m<5 then write(res,<: =┆
0x53320…53340 20 3a 3e 2c 3c 3c 2d 64 64 64 2e 64 64 64 30 30 3e 2c 76 29 0a 65 6e 64 3b 0a 63 79 63 6c 3a 3d ┆ :>,<<-ddd.ddd00>,v) end; cycl:=┆
0x53340…53360 20 30 3b 20 73 31 3a 3d 20 31 30 30 3b 20 67 6f 74 6f 20 72 65 70 0a 65 6e 64 20 74 2c 67 20 65 ┆ 0; s1:= 100; goto rep end t,g e┆
0x53360…53380 6e 64 20 58 3b 0a 0a 73 74 6f 70 3a 0a 77 72 69 74 65 28 72 65 73 2c 6e 6c 2c 35 2c 3c 3a 3c 32 ┆nd X; stop: write(res,nl,5,<:<2┆
0x53380…533a0 35 3e 3a 3e 29 3b 20 63 6c 6f 73 65 28 72 65 73 2c 63 6c 6f 73 65 72 65 73 29 0a 65 6e 64 0a 19 ┆5>:>); close(res,closeres) end ┆
0x533a0…533c0 31 20 75 6e 74 69 6c 20 69 20 64 6f 0a 20 20 20 77 72 69 74 65 28 72 65 73 2c 3c 3c 5f 5f 2d 64 ┆1 until i do write(res,<<__-d┆
0x533c0…533e0 2e 64 64 64 64 64 64 3e 2c 47 46 28 69 2b 33 2c 6a 2b 33 29 29 0a 65 6e 64 20 65 6e 64 20 65 6e ┆.dddddd>,GF(i+3,j+3)) end end en┆
0x533e0…53400 64 20 47 2c 5a 3b 0a 65 6e 64 20 49 74 3b 0a 0a 69 66 20 70 3d 30 20 74 68 65 6e 20 67 6f 74 6f ┆d G,Z; end It; if p=0 then goto┆
0x53400…53420 (111,) 3b 6b 65 6d 6c 61 62 35 20 31 0a 72 3d 61 6c 67 6f 6c 20 69 6e 64 65 78 2e 6e 6f 0a 0c 0a 0a 0a ┆;kemlab5 1 r=algol index.no ┆
0x53420…53440 4d 41 54 48 49 45 55 2d 6d 61 74 72 69 78 20 64 69 61 67 6f 6e 61 6c 69 73 65 72 69 6e 67 2e 20 ┆MATHIEU-matrix diagonalisering. ┆
0x53440…53460 20 20 20 20 20 20 20 20 38 2d 33 2d 31 39 37 31 2e 20 20 20 47 4f 53 2e 0a 0a 62 65 67 69 6e 0a ┆ 8-3-1971. GOS. begin ┆
0x53460…53480 63 6f 6d 6d 65 6e 74 20 69 6e 70 75 74 20 72 65 67 6c 65 72 6e 65 20 65 72 3a 0a 20 20 20 31 29 ┆comment input reglerne er: 1)┆
0x53480…534a0 20 42 61 72 72 69 65 72 65 6e 73 20 6d 75 6c 74 69 70 6c 69 63 69 74 65 74 2c 20 20 4e 2e 0a 20 ┆ Barrierens multiplicitet, N. ┆
0x534a0…534c0 20 20 32 29 20 45 6e 20 72 65 64 75 63 65 72 65 74 20 62 61 72 72 69 65 72 65 68 7c 6a 64 65 2c ┆ 2) En reduceret barrierehøjde,┆
0x534c0…534e0 20 68 65 72 20 64 65 66 69 6e 65 72 65 74 0a 20 20 20 20 20 20 20 20 20 42 20 3d 20 56 2a 49 72 ┆ her defineret B = V*Ir┆
0x534e0…53500 65 64 2f 28 32 2a 68 62 61 72 2a 2a 32 29 3d 20 4e 2a 4e 2a 73 2f 31 36 2e 0a 20 20 20 33 29 20 ┆ed/(2*hbar**2)= N*N*s/16. 3) ┆
0x53500…53520 53 79 6d 6d 65 74 72 69 70 61 72 61 6d 65 74 65 72 65 6e 20 73 69 67 6d 61 2c 20 30 3c 3d 73 69 ┆Symmetriparameteren sigma, 0<=si┆
0x53520…53540 67 6d 61 3c 3d 4e 2f 2f 32 2e 0a 20 20 20 34 29 20 4d 69 6e 64 73 74 65 20 6f 67 20 73 74 7c 72 ┆gma<=N//2. 4) Mindste og stør┆
0x53540…53560 73 74 65 20 76 7b 72 64 69 20 66 6f 72 20 76 69 62 72 61 74 69 6f 6e 73 6b 76 61 6e 74 65 74 61 ┆ste værdi for vibrationskvanteta┆
0x53560…53580 6c 6c 65 6e 65 2c 20 76 6d 69 6e 20 6f 67 20 76 6d 61 78 2e 0a 20 20 20 35 29 20 45 74 20 68 65 ┆llene, vmin og vmax. 5) Et he┆
0x53580…535a0 6c 74 61 6c 2c 20 74 20 28 6e 6f 72 6d 61 6c 74 20 3d 20 31 29 2c 20 64 65 72 20 62 65 73 74 65 ┆ltal, t (normalt = 1), der beste┆
0x535a0…535c0 6d 6d 65 72 20 74 72 75 6e 6b 65 72 69 6e 67 65 6e 3b 0a 69 6e 74 65 67 65 72 20 4e 2c 20 6e 2c ┆mmer trunkeringen; integer N, n,┆
0x535c0…535e0 20 69 2c 20 76 2c 20 76 6d 69 6e 2c 20 76 6d 61 78 2c 20 64 2c 20 6b 2c 20 6b 31 2c 20 6c 2c 20 ┆ i, v, vmin, vmax, d, k, k1, l, ┆
0x535e0…53600 6d 2c 20 73 2c 20 74 3b 0a 72 65 61 6c 20 41 31 2c 20 41 32 2c 20 42 2c 20 62 2c 20 79 2c 20 71 ┆m, s, t; real A1, A2, B, b, y, q┆
0x53600…53620 2c 20 71 31 2c 20 70 2c 20 70 31 2c 20 72 2c 20 61 2c 20 61 31 2c 20 67 2c 20 68 2c 20 6c 61 6d ┆, q1, p, p1, r, a, a1, g, h, lam┆
0x53620…53640 62 64 61 2c 20 6e 6f 72 6d 3b 0a 62 6f 6f 6c 65 61 6e 20 65 76 65 6e 2c 20 77 61 6e 67 2c 20 76 ┆bda, norm; boolean even, wang, v┆
0x53640…53660 6c 69 67 65 2c 20 63 6c 6f 73 65 72 65 73 3b 0a 7a 6f 6e 65 20 72 65 73 28 31 32 38 2c 31 2c 73 ┆lige, closeres; zone res(128,1,s┆
0x53660…53680 74 64 65 72 72 6f 72 29 3b 0a 70 72 6f 63 65 64 75 72 65 20 77 72 69 74 65 63 72 3b 0a 20 20 20 ┆tderror); procedure writecr; ┆
0x53680…536a0 62 65 67 69 6e 20 77 72 69 74 65 28 72 65 73 2c 3c 3a 3c 31 30 3e 3a 3e 29 20 65 6e 64 3b 0a 63 ┆begin write(res,<:<10>:>) end; c┆
0x536a0…536c0 6c 6f 73 65 72 65 73 3a 3d 20 6f 75 74 6d 65 64 69 75 6d 28 72 65 73 29 3b 0a 6e 65 77 4e 42 3a ┆loseres:= outmedium(res); newNB:┆
0x536c0…536e0 0a 72 65 61 64 28 69 6e 2c 4e 29 3b 20 69 66 20 4e 3c 31 20 74 68 65 6e 20 67 6f 74 6f 20 73 74 ┆ read(in,N); if N<1 then goto st┆
0x536e0…53700 6f 70 3b 0a 72 65 61 64 28 69 6e 2c 42 29 3b 20 77 72 69 74 65 28 72 65 73 2c 3c 3a 0a 4e 20 3d ┆op; read(in,B); write(res,<: N =┆
0x53700…53720 20 3a 3e 2c 3c 3c 64 64 64 3e 2c 4e 2c 3c 3a 20 20 20 20 42 20 3d 20 3a 3e 2c 3c 3c 64 64 64 2e ┆ :>,<<ddd>,N,<: B = :>,<<ddd.┆
0x53720…53740 64 64 64 64 64 30 30 30 30 30 3e 2c 42 2c 0a 3c 3a 20 20 20 20 20 73 20 3d 20 3a 3e 2c 20 31 36 ┆ddddd00000>,B, <: s = :>, 16┆
0x53740…53760 2a 42 2f 28 4e 2a 4e 29 29 3b 0a 6e 65 77 76 64 3a 0a 72 65 61 64 28 69 6e 2c 64 29 3b 20 77 72 ┆*B/(N*N)); newvd: read(in,d); wr┆
0x53760…53780 69 74 65 63 72 3b 20 77 72 69 74 65 63 72 3b 0a 69 66 20 64 3c 30 20 74 68 65 6e 20 62 65 67 69 ┆itecr; writecr; if d<0 then begi┆
0x53780…537a0 6e 20 77 72 69 74 65 63 72 3b 20 77 72 69 74 65 63 72 3b 20 77 72 69 74 65 63 72 3b 20 67 6f 74 ┆n writecr; writecr; writecr; got┆
0x537a0…537c0 6f 20 6e 65 77 4e 42 20 65 6e 64 3b 0a 72 65 61 64 28 69 6e 2c 76 6d 69 6e 2c 76 6d 61 78 2c 74 ┆o newNB end; read(in,vmin,vmax,t┆
0x537c0…537e0 29 3b 0a 77 72 69 74 65 28 72 65 73 2c 3c 3a 73 69 67 6d 61 2c 20 76 6d 69 6e 2c 20 76 6d 61 78 ┆); write(res,<:sigma, vmin, vmax┆
0x537e0…53800 20 6f 67 20 74 20 3d 20 3a 3e 2c 3c 3c 64 64 64 64 3e 2c 64 2c 76 6d 69 6e 2c 76 6d 61 78 2c 74 ┆ og t = :>,<<dddd>,d,vmin,vmax,t┆
0x53800…53820 29 3b 0a 65 76 65 6e 3a 3d 20 64 3d 30 3b 20 20 77 61 6e 67 3a 3d 20 65 76 65 6e 20 6f 72 20 28 ┆); even:= d=0; wang:= even or (┆
0x53820…53840 64 2b 64 29 3d 4e 3b 0a 62 3a 3d 20 31 2f 28 69 66 20 77 61 6e 67 20 74 68 65 6e 20 34 2a 42 20 ┆d+d)=N; b:= 1/(if wang then 4*B ┆
0x53840…53860 65 6c 73 65 20 42 29 3b 20 20 6c 61 6d 62 64 61 3a 3d 20 2d 32 3b 0a 74 3a 3d 20 74 2a 28 69 66 ┆else B); lambda:= -2; t:= t*(if┆
0x53860…53880 20 77 61 6e 67 20 74 68 65 6e 20 32 2a 4e 20 65 6c 73 65 20 4e 29 3b 0a 0a 66 6f 72 20 76 3a 3d ┆ wang then 2*N else N); for v:=┆
0x53880…538a0 76 6d 69 6e 20 73 74 65 70 20 31 20 75 6e 74 69 6c 20 76 6d 61 78 20 64 6f 20 62 65 67 69 6e 0a ┆vmin step 1 until vmax do begin ┆
0x538a0…538c0 68 3a 3d 20 6c 61 6d 62 64 61 3b 20 20 67 3a 3d 20 34 2a 76 2a 76 2f 42 2b 31 3b 0a 69 66 20 77 ┆h:= lambda; g:= 4*v*v/B+1; if w┆
0x538c0…538e0 61 6e 67 20 74 68 65 6e 20 62 65 67 69 6e 0a 20 20 20 76 6c 69 67 65 3a 3d 20 76 2f 2f 32 2a 32 ┆ang then begin vlige:= v//2*2┆
0x538e0…53900 3d 76 3b 0a 20 20 20 6b 31 3a 3d 20 69 66 20 65 76 65 6e 20 61 6e 64 20 76 6c 69 67 65 20 74 68 ┆=v; k1:= if even and vlige th┆
0x53900…53920 65 6e 20 30 20 65 6c 73 65 20 69 66 20 2d 2c 65 76 65 6e 20 74 68 65 6e 20 4e 20 65 6c 73 65 20 ┆en 0 else if -,even then N else ┆
0x53920…53940 32 2a 4e 3b 0a 20 20 20 41 31 3a 3d 20 62 2a 6b 31 2a 6b 31 2b 20 28 69 66 20 65 76 65 6e 20 74 ┆2*N; A1:= b*k1*k1+ (if even t┆
0x53940…53960 68 65 6e 20 30 20 65 6c 73 65 20 69 66 20 76 6c 69 67 65 20 74 68 65 6e 20 2d 31 20 65 6c 73 65 ┆hen 0 else if vlige then -1 else┆
0x53960…53980 20 31 29 3b 0a 20 20 20 41 32 3a 3d 20 69 66 20 6b 31 3d 30 20 74 68 65 6e 20 32 20 65 6c 73 65 ┆ 1); A2:= if k1=0 then 2 else┆
0x53980…539a0 20 31 3b 20 20 6e 3a 3d 20 20 76 2f 2f 32 3b 20 20 73 3a 3d 20 6b 31 2b 74 2b 4e 2a 34 0a 65 6e ┆ 1; n:= v//2; s:= k1+t+N*4 en┆
0x539a0…539c0 64 20 65 6c 73 65 20 62 65 67 69 6e 0a 20 20 20 73 3a 3d 20 74 2b 4e 2b 4e 3b 20 20 6e 3a 3d 20 ┆d else begin s:= t+N+N; n:= ┆
0x539c0…539e0 76 20 65 6e 64 3b 0a 0a 66 6f 72 20 6c 61 6d 62 64 61 3a 3d 20 28 67 2b 68 29 2a 30 2e 35 20 77 ┆v end; for lambda:= (g+h)*0.5 w┆
0x539e0…53a00 68 69 6c 65 20 67 3c 3e 6c 61 6d 62 64 61 20 64 6f 20 62 65 67 69 6e 0a 69 66 20 77 61 6e 67 20 ┆hile g<>lambda do begin if wang ┆
0x53a00…53a20 74 68 65 6e 20 62 65 67 69 6e 0a 20 20 20 71 3a 3d 20 41 31 2d 6c 61 6d 62 64 61 3b 20 70 3a 3d ┆then begin q:= A1-lambda; p:=┆
0x53a20…53a40 20 41 32 3b 20 20 6b 3a 3d 20 6b 31 3b 20 20 6d 3a 3d 20 69 66 20 71 3c 30 20 74 68 65 6e 20 31 ┆ A2; k:= k1; m:= if q<0 then 1┆
0x53a40…53a60 20 65 6c 73 65 20 30 3b 0a 20 20 20 66 6f 72 20 6b 3a 3d 20 6b 2b 4e 2b 4e 20 77 68 69 6c 65 20 ┆ else 0; for k:= k+N+N while ┆
0x53a60…53a80 6d 3c 3d 6e 20 61 6e 64 20 6b 3c 3d 73 20 64 6f 20 62 65 67 69 6e 0a 20 20 20 20 20 20 79 3a 3d ┆m<=n and k<=s do begin y:=┆
0x53a80…53aa0 20 28 62 2a 6b 2a 6b 2d 6c 61 6d 62 64 61 29 2a 71 2d 20 70 3b 0a 20 20 20 20 20 20 71 31 3a 3d ┆ (b*k*k-lambda)*q- p; q1:=┆
0x53aa0…53ac0 20 70 3b 20 70 3a 3d 20 71 3b 20 71 3a 3d 20 79 3b 0a 20 20 20 20 20 20 69 66 20 71 3e 3d 30 20 ┆ p; p:= q; q:= y; if q>=0 ┆
0x53ac0…53ae0 3d 3d 20 70 3c 30 20 74 68 65 6e 20 62 65 67 69 6e 20 6d 3a 3d 20 6d 2b 31 3b 20 69 66 20 73 3c ┆== p<0 then begin m:= m+1; if s<┆
0x53ae0…53b00 6b 2b 74 20 74 68 65 6e 20 73 3a 3d 20 6b 2b 74 20 65 6e 64 3b 0a 20 20 20 20 20 20 69 66 20 71 ┆k+t then s:= k+t end; if q┆
0x53b00…53b20 2a 71 31 3c 20 70 2a 70 20 61 6e 64 20 73 3c 6b 2b 74 20 74 68 65 6e 20 73 3a 3d 20 6b 2b 74 0a ┆*q1< p*p and s<k+t then s:= k+t ┆
0x53b20…53b40 65 6e 64 20 65 6e 64 20 65 6c 73 65 20 62 65 67 69 6e 0a 20 20 20 71 3a 3d 20 64 2a 64 2a 62 2d ┆end end else begin q:= d*d*b-┆
0x53b40…53b60 6c 61 6d 62 64 61 3b 20 20 70 3a 3d 20 70 31 3a 3d 20 31 3b 20 20 72 3a 3d 20 30 3b 0a 20 20 20 ┆lambda; p:= p1:= 1; r:= 0; ┆
0x53b60…53b80 6d 3a 3d 20 69 66 20 71 3c 30 20 74 68 65 6e 20 31 20 65 6c 73 65 20 30 3b 20 20 20 20 20 6b 3a ┆m:= if q<0 then 1 else 0; k:┆
0x53b80…53ba0 3d 20 30 3b 0a 20 20 20 66 6f 72 20 6b 3a 3d 20 6b 2b 4e 20 77 68 69 6c 65 20 6d 3c 3d 6e 20 61 ┆= 0; for k:= k+N while m<=n a┆
0x53ba0…53bc0 6e 64 20 6b 3c 3d 73 20 64 6f 20 62 65 67 69 6e 20 71 31 3a 3d 20 72 3b 0a 20 20 20 66 6f 72 20 ┆nd k<=s do begin q1:= r; for ┆
0x53bc0…53be0 6c 3a 3d 20 2d 6b 2b 64 2c 6b 2b 64 20 64 6f 20 62 65 67 69 6e 0a 20 20 20 20 20 20 61 3a 3d 20 ┆l:= -k+d,k+d do begin a:= ┆
0x53be0…53c00 6c 2a 6c 2a 62 2d 6c 61 6d 62 64 61 3b 0a 20 20 20 20 20 20 79 3a 3d 20 71 2a 61 2d 70 3b 20 20 ┆l*l*b-lambda; y:= q*a-p; ┆
0x53c00…53c20 70 3a 3d 20 70 31 3b 20 20 70 31 3a 3d 20 71 3b 20 20 71 3a 3d 20 79 3b 0a 20 20 20 20 20 20 79 ┆p:= p1; p1:= q; q:= y; y┆
0x53c20…53c40 3a 3d 20 70 2a 61 2d 72 3b 20 20 72 3a 3d 20 70 3b 20 20 20 70 3a 3d 20 79 3b 0a 20 20 20 20 20 ┆:= p*a-r; r:= p; p:= y; ┆
0x53c40…53c60 20 69 66 20 71 3e 3d 30 20 3d 3d 20 70 31 3c 30 20 74 68 65 6e 20 62 65 67 69 6e 20 6d 3a 3d 20 ┆ if q>=0 == p1<0 then begin m:= ┆
0x53c60…53c80 6d 2b 31 3b 20 69 66 20 73 3c 6b 2b 74 20 74 68 65 6e 20 73 3a 3d 20 6b 2b 74 20 65 6e 64 0a 20 ┆m+1; if s<k+t then s:= k+t end ┆
0x53c80…53ca0 20 20 65 6e 64 3b 0a 20 20 20 69 66 20 28 71 2a 71 31 3c 20 72 2a 72 20 6f 72 20 28 71 3e 3d 30 ┆ end; if (q*q1< r*r or (q>=0┆
0x53ca0…53cc0 20 3d 3d 20 72 3c 30 29 29 20 61 6e 64 20 73 3c 6b 2b 74 20 74 68 65 6e 20 73 3a 3d 20 6b 2b 74 ┆ == r<0)) and s<k+t then s:= k+t┆
0x53cc0…53ce0 0a 65 6e 64 20 65 6e 64 3b 0a 69 66 20 6d 3e 6e 20 74 68 65 6e 20 67 3a 3d 20 6c 61 6d 62 64 61 ┆ end end; if m>n then g:= lambda┆
0x53ce0…53d00 20 65 6c 73 65 20 68 3a 3d 20 6c 61 6d 62 64 61 0a 65 6e 64 20 62 69 73 65 6b 74 69 6f 6e 3b 0a ┆ else h:= lambda end bisektion; ┆
0x53d00…53d20 0a 79 3a 3d 20 6c 61 6d 62 64 61 2a 30 2e 32 35 2b 30 2e 35 3b 20 20 20 77 72 69 74 65 63 72 3b ┆ y:= lambda*0.25+0.5; writecr;┆
0x53d20…53d40 0a 6d 3a 3d 20 69 66 20 77 61 6e 67 20 74 68 65 6e 20 28 28 73 2d 6b 31 29 2f 2f 28 4e 2b 4e 29 ┆ m:= if wang then ((s-k1)//(N+N)┆
0x53d40…53d60 2b 31 29 20 65 6c 73 65 20 73 2f 2f 4e 3b 0a 77 72 69 74 65 28 72 65 73 2c 3c 3a 76 20 3d 20 3a ┆+1) else s//N; write(res,<:v = :┆
0x53d60…53d80 3e 2c 3c 3c 64 64 3e 2c 76 2c 3c 3a 20 20 6f 72 64 65 6e 3a 20 3a 3e 2c 3c 3c 64 64 3e 2c 6d 29 ┆>,<<dd>,v,<: orden: :>,<<dd>,m)┆
0x53d80…53da0 3b 0a 77 72 69 74 65 28 72 65 73 2c 3c 3a 20 20 6c 61 6d 62 64 61 20 3d 20 3a 3e 2c 3c 3c 2d 64 ┆; write(res,<: lambda = :>,<<-d┆
0x53da0…53dc0 64 2e 64 64 64 20 64 64 64 20 64 64 64 20 30 30 30 3e 2c 20 79 29 3b 0a 77 72 69 74 65 28 72 65 ┆d.ddd ddd ddd 000>, y); write(re┆
0x53dc0…53de0 73 2c 3c 3a 20 20 62 20 3d 20 3a 3e 2c 3c 3c 2d 64 64 64 2e 64 64 64 20 64 64 64 20 64 64 30 20 ┆s,<: b = :>,<<-ddd.ddd ddd dd0 ┆
0x53de0…53e00 30 30 3e 2c 20 79 2a 31 36 2a 42 2f 28 4e 2a 4e 29 29 3b 0a 69 66 20 77 61 6e 67 20 74 68 65 6e ┆00>, y*16*B/(N*N)); if wang then┆
0x53e00…53e20 20 62 65 67 69 6e 0a 61 72 72 61 79 20 55 28 31 3a 6d 29 3b 0a 55 28 6d 29 3a 3d 20 31 3b 20 20 ┆ begin array U(1:m); U(m):= 1; ┆
0x53e20…53e40 55 28 6d 2d 31 29 3a 3d 20 62 2a 73 2a 73 2d 6c 61 6d 62 64 61 3b 20 20 6e 3a 3d 20 31 3b 20 20 ┆U(m-1):= b*s*s-lambda; n:= 1; ┆
0x53e40…53e60 6e 6f 72 6d 3a 3d 20 31 3b 0a 6b 3a 3d 20 73 2d 4e 2d 4e 3b 20 20 20 20 20 61 3a 3d 20 62 2a 6b ┆norm:= 1; k:= s-N-N; a:= b*k┆
0x53e60…53e80 2a 6b 2d 6c 61 6d 62 64 61 3b 0a 6b 3a 3d 20 6b 2d 4e 2d 4e 3b 20 20 20 20 20 61 31 3a 3d 20 62 ┆*k-lambda; k:= k-N-N; a1:= b┆
0x53e80…53ea0 2a 6b 2a 6b 2d 6c 61 6d 62 64 61 3b 0a 66 6f 72 20 69 3a 3d 6d 2d 32 20 73 74 65 70 20 2d 31 20 ┆*k*k-lambda; for i:=m-2 step -1 ┆
0x53ea0…53ec0 75 6e 74 69 6c 20 6e 20 64 6f 20 62 65 67 69 6e 0a 20 20 20 55 28 69 29 3a 3d 20 55 28 69 2b 31 ┆until n do begin U(i):= U(i+1┆
0x53ec0…53ee0 29 2a 61 2d 20 55 28 69 2b 32 29 3b 0a 20 20 20 61 3a 3d 20 61 31 3b 20 20 6b 3a 3d 20 6b 2d 4e ┆)*a- U(i+2); a:= a1; k:= k-N┆
0x53ee0…53f00 2d 4e 3b 20 20 61 31 3a 3d 20 62 2a 6b 2a 6b 2d 6c 61 6d 62 64 61 3b 0a 20 20 20 69 66 20 61 62 ┆-N; a1:= b*k*k-lambda; if ab┆
0x53f00…53f20 73 20 61 3c 20 61 62 73 20 61 31 20 74 68 65 6e 20 6e 3a 3d 20 69 0a 65 6e 64 3b 0a 69 66 20 6e ┆s a< abs a1 then n:= i end; if n┆
0x53f20…53f40 3d 31 20 74 68 65 6e 20 62 65 67 69 6e 0a 20 20 20 55 28 31 29 3a 3d 20 55 28 31 29 2f 73 71 72 ┆=1 then begin U(1):= U(1)/sqr┆
0x53f40…53f60 74 28 41 32 29 3b 0a 20 20 20 66 6f 72 20 69 3a 3d 6d 2d 31 20 73 74 65 70 20 2d 31 20 75 6e 74 ┆t(A2); for i:=m-1 step -1 unt┆
0x53f60…53f80 69 6c 20 31 20 64 6f 20 6e 6f 72 6d 3a 3d 20 6e 6f 72 6d 2b 55 28 69 29 2a 2a 32 0a 65 6e 64 20 ┆il 1 do norm:= norm+U(i)**2 end ┆
0x53f80…53fa0 65 6c 73 65 20 62 65 67 69 6e 0a 20 20 20 67 3a 3d 20 31 2f 55 28 6e 29 3b 0a 20 20 20 66 6f 72 ┆else begin g:= 1/U(n); for┆
0x53fa0…53fc0 20 69 3a 3d 6e 2b 31 20 73 74 65 70 20 31 20 75 6e 74 69 6c 20 6d 20 64 6f 20 62 65 67 69 6e 0a ┆ i:=n+1 step 1 until m do begin ┆
0x53fc0…53fe0 20 20 20 20 20 20 79 3a 3d 20 55 28 69 29 3a 3d 20 55 28 69 29 2a 67 3b 20 20 6e 6f 72 6d 3a 3d ┆ y:= U(i):= U(i)*g; norm:=┆
0x53fe0…54000 20 6e 6f 72 6d 2b 79 2a 79 20 65 6e 64 3b 0a 20 20 20 55 28 31 29 3a 3d 20 41 32 3b 20 20 55 28 ┆ norm+y*y end; U(1):= A2; U(┆