|
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: 43008 (0xa800)
Types: Rc489kSubCat
Names: »rymac«
└─⟦00964e8f7⟧ Bits:30007478 RC8000 Dump tape fra HCØ.
└─⟦b2ec5d50f⟧
└─⟦this⟧ »rymac«
0x0000…0033 (0,) Rc489kSubCatEnt {f00=0x0031c2, f01=0x000078, f02=0x000081, filename=»atall «, entry_tail=0-4-0-1-1982-04-28T16:32-disc2-0x10028-0x0-0x0-0x0}
0x0033…0066 Rc489kSubCatEnt {f00=0x004dda, f01=0x000078, f02=0x000081, filename=»atprint «, entry_tail=0-4-0-1-1982-04-28T16:17-disc2-0xe033-0x0-0x0-0x0}
0x0066…0099 Rc489kSubCatEnt {f00=0x005b72, f01=0x000078, f02=0x000081, filename=»attau «, entry_tail=0-4-0-1-1982-04-28T16:13-disc2-0xe032-0x0-0x0-0x0}
0x0099…00cc Rc489kSubCatEnt {f00=0x006052, f01=0x000078, f02=0x000081, filename=»atcompz «, entry_tail=0-4-0-1-1982-04-27T12:58-disc2-0x0-0x0-0x0-0x0}
0x00cc…00ff Rc489kSubCatEnt {f00=0x007e42, f01=0x000078, f02=0x000081, filename=»atcompute «, entry_tail=0-4-0-1-1982-04-20T12:47-disc2-0xb023-0x0-0x0-0x0}
0x00ff…0132 Rc489kSubCatEnt {f00=0x0088a2, f01=0x000078, f02=0x000081, filename=»atinit «, entry_tail=0-4-0-1-1982-04-20T12:34-disc2-0x0-0x0-0x0-0x0}
0x0132…0165 Rc489kSubCatEnt {f00=0x0098ba, f01=0x000078, f02=0x000081, filename=»atsort «, entry_tail=0-4-0-1-1981-02-03T13:35-disc2-0x0-0x0-0x0-0x0}
0x0165…0198 Rc489kSubCatEnt {f00=0x00a9f2, f01=0x000078, f02=0x000081, filename=»atlinelist «, entry_tail=0-4-0-1-1980-11-25T09:29-disc2-0x0-0x0-0x0-0x0}
0x0198…01cb Rc489kSubCatEnt {f00=0x00b63a, f01=0x000078, f02=0x000081, filename=»atinitfit «, entry_tail=0-4-0-1-1980-11-21T12:46-disc2-0x0-0x0-0x0-0x0}
0x01cb…01fe Rc489kSubCatEnt {f00=0x00cd62, f01=0x000078, f02=0x000081, filename=»atlist «, entry_tail=0-4-0-1-1980-11-20T14:24-disc2-0x0-0x0-0x0-0x0}
0x01fe…0231 Rc489kSubCatEnt {f00=0x00d41a, f01=0x000078, f02=0x000081, filename=»atline «, entry_tail=0-4-0-1-1980-10-15T16:10-disc3-0x0-0x0-0x0-0x0}
0x0231…0264 Rc489kSubCatEnt {f00=0x00e0c8, f01=0x000079, f02=0x000079, filename=»atresort «, entry_tail=0-4-0-1-1980-10-13T13:15-disc2-0x0-0x0-0x0-0x0}
0x0264…0297 Rc489kSubCatEnt {f00=0x00f003, f01=0x000078, f02=0x000081, filename=»gensimmax «, entry_tail=0-4-0-1-1980-05-07T14:44--0x11009-0x0-0x0-0x0}
0x0297…02ca Rc489kSubCatEnt {f00=0x010003, f01=0x000078, f02=0x000081, filename=»gensim «, entry_tail=0-4-0-2-1980-05-07T14:44--0x9016-0x0-0x0-0x0}
0x02ca…02fd Rc489kSubCatEnt {f00=0x012003, f01=0x000078, f02=0x000081, filename=»atcomps4 «, entry_tail=0-4-0-1-1980-05-07T14:44--0xe00c-0x0-0x0-0x0}
0x02fd…0300 0x000029
0x0300…0333 Rc489kSubCatEnt {f00=0x013003, f01=0x000078, f02=0x000081, filename=»atcompsdr «, entry_tail=0-4-0-1-1980-05-07T14:44--0xb020-0x0-0x0-0x0}
0x0333…0366 Rc489kSubCatEnt {f00=0x014003, f01=0x000078, f02=0x000081, filename=»atpub «, entry_tail=0-4-0-2-1980-05-07T14:44--0x9027-0x0-0x0-0x0}
0x0366…0399 Rc489kSubCatEnt {f00=0x016003, f01=0x000078, f02=0x000081, filename=»atpunch «, entry_tail=0-4-0-1-1980-05-07T14:44--0xa038-0x0-0x0-0x0}
0x0399…03cc Rc489kSubCatEnt {f00=0x017003, f01=0x000078, f02=0x000081, filename=»attau1 «, entry_tail=0-4-0-1-1980-05-07T14:44--0xa00e-0x0-0x0-0x0}
0x03cc…03ff Rc489kSubCatEnt {f00=0x018003, f01=0x000078, f02=0x000081, filename=»attau2 «, entry_tail=0-4-0-1-1980-05-07T14:44--0xd013-0x0-0x0-0x0}
0x03ff…0432 Rc489kSubCatEnt {f00=0x019003, f01=0x000078, f02=0x000081, filename=»attusch «, entry_tail=0-4-0-1-1980-05-07T14:44--0xc032-0x0-0x0-0x0}
0x0432…0465 Rc489kSubCatEnt {f00=0x01a003, f01=0x000078, f02=0x000081, filename=»atinit3 «, entry_tail=0-4-0-1-1980-05-07T14:44--0xe028-0x0-0x0-0x0}
0x0465…0498 Rc489kSubCatEnt {f00=0x01b003, f01=0x000078, f02=0x000081, filename=»atinittst «, entry_tail=0-4-0-1-1980-05-07T14:44--0x0-0x0-0x0-0x0}
0x0498…04cb Rc489kSubCatEnt {f00=0x01c003, f01=0x000078, f02=0x000081, filename=»atnjob «, entry_tail=0-4-0-1-1980-05-07T14:44--0xe023-0x0-0x0-0x0}
0x04cb…04fe Rc489kSubCatEnt {f00=0x01d003, f01=0x000078, f02=0x000081, filename=»atpjob «, entry_tail=0-4-0-1-1980-05-07T14:44--0xe00a-0x0-0x0-0x0}
0x04fe…0531 Rc489kSubCatEnt {f00=0x01e003, f01=0x000078, f02=0x000081, filename=»atplot «, entry_tail=0-4-0-1-1980-05-07T14:44--0xb012-0x0-0x0-0x0}
0x0531…0564 Rc489kSubCatEnt {f00=0x01f003, f01=0x000078, f02=0x000081, filename=»atplot1 «, entry_tail=0-4-0-1-1980-05-07T14:44--0xd025-0x0-0x0-0x0}
0x0564…0597 Rc489kSubCatEnt {f00=0x020003, f01=0x000078, f02=0x000081, filename=»atplott «, entry_tail=0-4-0-2-1980-05-07T14:44--0x9022-0x0-0x0-0x0}
0x0597…05ca Rc489kSubCatEnt {f00=0x022003, f01=0x000078, f02=0x000081, filename=»atplotz «, entry_tail=0-4-0-1-1980-05-07T14:44--0xd01d-0x0-0x0-0x0}
0x05ca…05fd Rc489kSubCatEnt {f00=0x023003, f01=0x000078, f02=0x000081, filename=»atcompss «, entry_tail=0-4-0-1-1980-05-07T14:44--0xd02a-0x0-0x0-0x0}
0x05fd…0600 00 00 00 ┆ ┆
0x0600…0633 Rc489kSubCatEnt {f00=0x024003, f01=0x000078, f02=0x000081, filename=»atcuts «, entry_tail=0-4-0-1-1980-05-07T14:44--0xa012-0x0-0x0-0x0}
0x0633…0666 Rc489kSubCatEnt {f00=0x025003, f01=0x000078, f02=0x000081, filename=»atinit1 «, entry_tail=0-4-0-1-1980-05-07T14:44--0xe026-0x0-0x0-0x0}
0x0666…0699 Rc489kSubCatEnt {f00=0x026003, f01=0x000078, f02=0x000081, filename=»atwf «, entry_tail=0-4-0-1-1980-05-07T14:44--0xc007-0x0-0x0-0x0}
0x0699…06cc Rc489kSubCatEnt {f00=0x027003, f01=0x000078, f02=0x000081, filename=»atalls «, entry_tail=0-4-0-1-1980-05-07T14:44--0x10029-0x0-0x0-0x0}
0x06cc…06ff Rc489kSubCatEnt {f00=0x028003, f01=0x000078, f02=0x000081, filename=»atcomphl «, entry_tail=0-4-0-1-1980-05-07T14:44--0x900c-0x0-0x0-0x0}
0x06ff…0732 Rc489kSubCatEnt {f00=0x029003, f01=0x000078, f02=0x000081, filename=»atcomps «, entry_tail=0-4-0-1-1980-05-07T14:44--0x10001-0x0-0x0-0x0}
0x0732…0765 Rc489kSubCatEnt {f00=0x02a003, f01=0x000078, f02=0x000081, filename=»atcompsl «, entry_tail=0-4-0-1-1980-05-07T14:44--0xb02f-0x0-0x0-0x0}
0x0765…0798 Rc489kSubCatEnt {f00=0x02b003, f01=0x000078, f02=0x000081, filename=»atcompsz «, entry_tail=0-4-0-1-1980-05-07T14:44--0xc037-0x0-0x0-0x0}
0x0798…07cb Rc489kSubCatEnt {f00=0x02c003, f01=0x000078, f02=0x000081, filename=»atcompwf «, entry_tail=0-4-0-1-1980-05-07T14:44--0xd017-0x0-0x0-0x0}
0x07cb…07fe Rc489kSubCatEnt {f00=0x02d003, f01=0x000078, f02=0x000081, filename=»atcompszdr «, entry_tail=0-4-0-1-1980-05-07T14:44--0x0-0x0-0x0-0x0}
0x07fe…0831 Rc489kSubCatEnt {f00=0x02e003, f01=0x000078, f02=0x000081, filename=»atdiagram «, entry_tail=0-4-0-1-1980-05-07T14:44--0xc032-0x0-0x0-0x0}
0x0831…0840 01 f0 03 00 00 78 00 00 81 61 74 70 6c 6f 74 ┆ x atplot┆
0x0840…0860 31 00 00 00 00 00 00 00 01 00 00 00 00 00 00 00 00 00 00 00 00 71 6c 2f 00 d0 25 00 00 00 00 00 ┆1 ql/ % ┆
0x0860…0880 00 00 00 00 02 00 03 00 00 78 00 00 81 61 74 70 6c 6f 74 74 00 00 00 00 00 00 00 02 00 00 00 00 ┆ x atplott ┆
0x0880…08a0 00 00 00 00 00 00 00 00 71 6c 2f 00 90 22 00 00 00 00 00 00 00 00 00 02 20 03 00 00 78 00 00 81 ┆ ql/ " x ┆
0x08a0…08c0 61 74 70 6c 6f 74 7a 00 00 00 00 00 00 00 01 00 00 00 00 00 00 00 00 00 00 00 00 71 6c 2f 00 d0 ┆atplotz ql/ ┆
0x08c0…08e0 1d 00 00 00 00 00 00 00 00 00 02 30 03 00 00 78 00 00 81 61 74 63 6f 6d 70 73 73 00 00 00 00 00 ┆ 0 x atcompss ┆
0x08e0…0900 00 01 00 00 00 00 00 00 00 00 00 00 00 00 71 6c 2f 00 d0 2a 00 00 00 00 00 00 00 00 00 00 00 00 ┆ ql/ * ┆
0x0900…0c00 (1,) Data for atall
0x0c00…0f00 Data for atprint
0x0f00…1200 Data for attau
0x1200…1500 (2,) Data for atcompz
0x1500…1800 Data for atcompute
0x1800…1b00 Data for atinit
0x1b00…1e00 (3,) Data for atsort
0x1e00…2100 Data for atlinelist
0x2100…2400 Data for atinitfit
0x2400…2700 (4,) Data for atlist
0x2700…2a00 Data for atline
0x2a00…2d00 Data for atresort
0x2d00…3000 (5,) Data for gensimmax
0x3000…3600 Data for gensim
0x3600…3900 (6,) Data for atcomps4
0x3900…3c00 Data for atcompsdr
0x3c00…4200 Data for atpub
0x4200…4500 Data for atpunch
0x4500…4800 Data for attau1
0x4800…4b00 (8,) Data for attau2
0x4b00…4e00 Data for attusch
0x4e00…5100 Data for atinit3
0x5100…5400 (9,) Data for atinittst
0x5400…5700 Data for atnjob
0x5700…5a00 Data for atpjob
0x5a00…5d00 (10,) Data for atplot
0x5d00…6000 Data for atplot1
0x6000…6600 Data for atplott
0x6600…6900 Data for atplotz
0x6900…6c00 Data for atcompss
0x6c00…6f00 (12,) Data for atcuts
0x6f00…7200 Data for atinit1
0x7200…7500 Data for atwf
0x7500…7800 (13,) Data for atalls
0x7800…7b00 Data for atcomphl
0x7b00…7e00 Data for atcomps
0x7e00…8100 (14,) Data for atcompsl
0x8100…8400 Data for atcompsz
0x8400…8700 Data for atcompwf
0x8700…8a00 (15,) Data for atcompszdr
0x8a00…8d00 Data for atdiagram
0x8d00…8d20 64 2c 20 61 6e 64 20 73 6f 20 69 74 20 73 68 6f 75 6c 64 20 62 65 20 0a 75 73 65 64 20 77 68 65 ┆d, and so it should be used whe┆
0x8d20…8d40 6e 20 73 6f 6d 65 20 6f 66 20 74 68 65 20 7a 65 72 6f 73 20 61 72 65 20 65 78 70 65 63 74 65 64 ┆n some of the zeros are expected┆
0x8d40…8d60 20 6f 72 20 6b 6e 6f 77 6e 20 74 6f 20 62 65 20 72 65 61 6c 2e 20 28 63 66 72 2e 20 73 70 65 72 ┆ or known to be real. (cfr. sper┆
0x8d60…8d80 63 29 2e 0a 22 6e 6c 32 22 0a 63 61 6c 6c 3a 20 73 70 65 72 63 72 65 28 41 52 2c 6e 2c 41 49 29 ┆c). "nl2" call: spercre(AR,n,AI)┆
0x8d80…8da0 3b 0a 22 6c 6d 20 32 35 22 0a 22 6e 6c 32 22 0a 22 6d 74 20 31 2c 20 20 20 20 20 20 41 52 20 28 ┆; "lm 25" "nl2" "mt 1, AR (┆
0x8da0…8dc0 22 0a 63 61 6c 6c 20 61 6e 64 20 72 65 74 75 72 6e 20 76 61 6c 75 65 2c 20 61 72 72 61 79 28 30 ┆" call and return value, array(0┆
0x8dc0…8de0 3a 6e 29 29 2e 22 6e 6c 22 0a 41 74 20 63 61 6c 6c 3a 20 41 52 28 6a 29 3d 20 63 6f 65 66 66 69 ┆:n))."nl" At call: AR(j)= coeffi┆
0x8de0…8e00 63 69 65 6e 74 20 74 6f 20 7a 2a 2a 28 6e 2d 6a 29 20 6a 3d 30 2c 2e 2e 2c 6e 22 6e 6c 22 0a 61 ┆cient to z**(n-j) j=0,..,n"nl" a┆
0x8e00…8e20 74 20 72 65 74 75 72 6e 3a 20 54 68 65 20 72 65 61 6c 20 70 61 72 74 20 6f 66 20 74 68 65 20 7a ┆t return: The real part of the z┆
0x8e20…8e40 65 72 6f 73 20 3d 20 41 52 28 6a 29 20 6a 3d 31 2c 2e 2e 2c 6e 2e 0a 22 6e 6c 32 22 0a 22 6d 74 ┆eros = AR(j) j=1,..,n. "nl2" "mt┆
0x8e40…8e60 20 31 2c 20 20 20 20 20 20 6e 20 20 28 22 0a 63 61 6c 6c 20 76 61 6c 75 65 2c 20 69 6e 74 65 67 ┆ 1, n (" call value, integ┆
0x8e60…8e80 65 72 29 2e 20 54 68 65 20 64 65 67 72 65 65 20 6f 66 20 74 68 65 20 65 71 75 61 74 69 6f 6e 2e ┆er). The degree of the equation.┆
0x8e80…8ea0 0a 22 6e 6c 32 22 0a 22 6d 74 20 31 2c 20 20 20 20 20 20 41 49 20 28 22 0a 20 72 65 74 75 72 6e ┆ "nl2" "mt 1, AI (" return┆
0x8ea0…8ec0 20 76 61 6c 75 65 2c 20 61 72 72 61 79 28 30 3a 6e 29 29 2e 22 6e 6c 22 0a 41 74 20 72 65 74 75 ┆ value, array(0:n))."nl" At retu┆
0x8ec0…8ee0 72 6e 3a 20 54 68 65 20 69 6d 61 67 69 6e 61 72 79 20 70 61 72 74 20 6f 66 20 74 68 65 20 7a 65 ┆rn: The imaginary part of the ze┆
0x8ee0…8f00 72 6f 73 20 20 41 49 28 6a 29 20 6a 3d 31 2c 2e 2e 2c 6e 2e 0a 22 6c 6d 30 22 22 6e 6c 33 22 0a ┆ros AI(j) j=1,..,n. "lm0""nl3" ┆
0x8f00…8f20 6c 61 6e 67 75 61 67 65 3a 20 41 4c 47 4f 4c 3b 20 74 65 78 74 20 69 6e 20 73 70 65 72 63 74 78 ┆language: ALGOL; text in sperctx┆
0x8f20…8f40 74 22 6e 6c 32 22 0a 4d 61 69 6e 74 65 6e 61 6e 63 65 20 6c 65 76 65 6c 3a 20 75 73 65 72 0a 22 ┆t"nl2" Maintenance level: user "┆
0x8f40…8f60 6e 6c 32 22 52 65 66 65 72 65 6e 63 65 73 3a 0a 4f 2e 4a 2e 48 65 69 6c 6d 61 6e 6e 20 26 20 52 ┆nl2"References: O.J.Heilmann & R┆
0x8f60…8f80 2e 4d 6f 73 73 2c 20 72 61 70 70 6f 72 74 20 33 39 2c 20 31 39 37 33 2c 20 4b 65 6d 2e 20 4c 61 ┆.Moss, rapport 39, 1973, Kem. La┆
0x8f80…8fa0 62 2e 20 49 49 49 2e 0a 22 70 73 30 22 0a 22 71 72 22 31 39 37 38 2d 39 2d 32 32 22 6e 6c 22 52 ┆b. III. "ps0" "qr"1978-9-22"nl"R┆
0x8fa0…8fc0 65 6e 65 20 4d 6f 73 73 22 72 6a 22 22 63 74 22 22 6e 6c 34 22 0a 70 72 6f 63 65 64 75 72 65 20 ┆ene Moss"rj""ct""nl4" procedure ┆
0x8fc0…8fe0 73 70 6c 6e 33 22 72 6a 22 22 6e 6c 32 22 0a 64 65 74 65 72 6d 69 6e 65 73 20 61 20 63 75 62 69 ┆spln3"rj""nl2" determines a cubi┆
0x8fe0…9000 63 20 73 70 6c 69 6e 65 20 74 68 72 6f 75 67 68 20 61 20 73 65 74 20 6f 66 20 70 6f 69 6e 74 73 ┆c spline through a set of points┆
0x9000…9020 (16,) 3a 22 6e 70 31 30 22 0a 20 58 28 69 29 2c 20 59 28 69 29 20 69 3d 31 2c 2e 2e 2e 2c 6e 3b 0a 58 ┆:"np10" X(i), Y(i) i=1,...,n; X┆
0x9020…9040 28 69 29 20 3c 20 58 28 6a 29 20 66 6f 72 20 69 3c 6a 2e 22 6e 6c 22 0a 20 73 70 6c 6e 33 20 63 ┆(i) < X(j) for i<j."nl" spln3 c┆
0x9040…9060 61 6c 63 75 6c 61 74 65 73 20 74 68 65 20 32 2e 20 64 65 72 69 76 61 74 69 76 65 20 6f 66 20 74 ┆alculates the 2. derivative of t┆
0x9060…9080 68 65 0a 73 70 6c 69 6e 65 20 4d 28 69 29 20 61 74 20 74 68 65 20 78 2d 63 6f 6f 72 64 69 6e 61 ┆he spline M(i) at the x-coordina┆
0x9080…90a0 74 65 73 20 58 28 69 29 3b 20 58 2c 59 2c 4d 20 64 65 66 69 6e 65 20 74 68 65 20 73 70 6c 69 6e ┆tes X(i); X,Y,M define the splin┆
0x90a0…90c0 65 0a 63 6f 6d 70 6c 65 74 65 6c 79 2e 20 44 69 66 66 65 72 65 6e 74 20 62 6f 75 6e 64 61 72 79 ┆e completely. Different boundary┆
0x90c0…90e0 2d 63 6f 6e 64 69 74 69 6f 6e 73 20 61 72 65 20 70 6f 73 73 69 62 6c 65 20 61 73 20 73 70 65 63 ┆-conditions are possible as spec┆
0x90e0…9100 69 66 69 65 64 20 69 6e 0a 54 59 50 45 2e 22 6e 6c 22 54 68 65 20 66 75 6e 63 74 69 6f 6e 20 61 ┆ified in TYPE."nl"The function a┆
0x9100…9120 6e 64 20 74 68 65 20 64 65 72 69 76 61 74 69 76 65 73 20 61 72 65 20 65 76 61 6c 75 61 74 65 64 ┆nd the derivatives are evaluated┆
0x9120…9140 0a 62 79 20 73 70 6c 6e 33 76 61 6c 2c 20 73 70 6c 6e 33 64 69 66 2c 20 73 70 6c 6e 33 64 69 66 ┆ by spln3val, spln3dif, spln3dif┆
0x9140…9160 32 20 61 6e 64 20 73 70 6c 6e 33 69 6e 74 20 61 6e 64 20 74 68 65 20 69 6e 76 65 72 73 65 0a 66 ┆2 and spln3int and the inverse f┆
0x9160…9180 75 6e 63 74 69 6f 6e 20 69 73 20 67 69 76 65 6e 20 62 79 20 73 70 6c 6e 33 69 6e 76 2e 22 6e 6c ┆unction is given by spln3inv."nl┆
0x9180…91a0 32 22 0a 63 61 6c 6c 3a 20 73 70 6c 6e 33 28 54 59 50 45 2c 58 2c 59 2c 4d 2c 6e 29 0a 22 6e 6c ┆2" call: spln3(TYPE,X,Y,M,n) "nl┆
0x91a0…91c0 32 22 22 6c 6d 33 31 22 0a 22 6d 74 20 31 2c 20 20 20 20 20 20 54 59 50 45 20 28 22 0a 63 61 6c ┆2""lm31" "mt 1, TYPE (" cal┆
0x91c0…91e0 6c 20 76 61 6c 75 65 2c 20 61 72 72 61 79 28 30 3a 34 29 29 20 54 59 50 45 28 30 29 20 73 65 6c ┆l value, array(0:4)) TYPE(0) sel┆
0x91e0…9200 65 63 74 73 20 74 68 65 20 74 79 70 65 20 6f 66 20 65 6e 64 2d 63 6f 6e 64 69 74 69 6f 6e 2c 0a ┆ects the type of end-condition, ┆
0x9200…9220 77 68 69 6c 65 20 74 68 65 20 65 6e 64 2d 63 6f 6e 64 69 74 69 6f 6e 20 69 73 20 73 70 65 63 69 ┆while the end-condition is speci┆
0x9220…9240 66 69 65 64 20 69 6e 20 54 59 50 45 28 31 3a 34 29 20 61 73 20 69 6e 64 69 63 61 74 65 64 3a 0a ┆fied in TYPE(1:4) as indicated: ┆
0x9240…9260 22 6c 6d 33 36 22 22 6e 6c 32 22 0a 22 6d 74 20 31 2c 54 59 50 45 28 30 29 20 3d 20 30 20 3a 20 ┆"lm36""nl2" "mt 1,TYPE(0) = 0 : ┆
0x9260…9280 22 0a 20 67 65 6e 65 72 61 6c 20 65 6e 64 2d 63 6f 6e 64 69 74 69 6f 6e 73 3a 0a 22 6e 6c 22 0a ┆" general end-conditions: "nl" ┆
0x9280…92a0 4d 28 31 29 2b 54 59 50 45 28 32 29 2a 4d 28 32 29 23 23 23 23 23 23 23 23 3d 23 54 59 50 45 28 ┆M(1)+TYPE(2)*M(2)########=#TYPE(┆
0x92a0…92c0 31 29 22 6e 6c 22 0a 23 23 23 23 23 54 59 50 45 28 33 29 2a 4d 28 6e 2d 31 29 2b 4d 28 6e 29 23 ┆1)"nl" #####TYPE(3)*M(n-1)+M(n)#┆
0x92c0…92e0 3d 23 54 59 50 45 28 34 29 22 6e 6c 31 22 0a 22 6d 74 20 31 2c 20 20 20 20 20 20 20 20 20 20 31 ┆=#TYPE(4)"nl1" "mt 1, 1┆
0x92e0…9300 20 3a 20 22 0a 20 32 2e 20 64 65 72 69 76 61 74 69 76 65 20 3d 20 30 20 61 74 20 58 28 31 29 20 ┆ : " 2. derivative = 0 at X(1) ┆
0x9300…9320 61 6e 64 20 58 28 6e 29 2e 22 6e 6c 31 22 0a 22 6d 74 20 31 2c 20 20 20 20 20 20 20 20 20 20 32 ┆and X(n)."nl1" "mt 1, 2┆
0x9320…9340 20 3a 20 22 0a 20 54 68 65 20 65 6e 64 73 6c 6f 70 65 73 20 61 74 20 58 28 31 29 20 61 6e 64 20 ┆ : " The endslopes at X(1) and ┆
0x9340…9360 58 28 6e 29 20 61 72 65 0a 73 70 65 63 69 66 69 65 64 20 69 6e 20 54 59 50 45 28 31 29 20 72 65 ┆X(n) are specified in TYPE(1) re┆
0x9360…9380 73 70 65 63 74 69 76 65 6c 79 20 54 59 50 45 28 34 29 2e 22 6e 6c 31 22 0a 22 6d 74 20 31 2c 20 ┆spectively TYPE(4)."nl1" "mt 1, ┆
0x9380…93a0 20 20 20 20 20 20 20 20 20 33 20 3a 20 22 0a 20 32 2e 20 64 65 72 69 76 61 74 69 76 65 20 3d 20 ┆ 3 : " 2. derivative = ┆
0x93a0…93c0 30 20 28 69 2e 65 2e 20 68 61 73 20 73 75 70 70 6f 72 74 20 61 74 0a 78 30 2c 78 70 20 77 68 65 ┆0 (i.e. has support at x0,xp whe┆
0x93c0…93e0 72 65 0a 58 28 31 29 2d 78 30 3d 33 28 58 28 32 29 2d 58 28 31 29 29 2c 20 78 30 3c 58 28 31 29 ┆re X(1)-x0=3(X(2)-X(1)), x0<X(1)┆
0x93e0…9400 22 6e 6c 22 0a 23 23 23 23 23 23 78 70 2d 58 28 6e 29 3d 33 28 58 28 6e 29 2d 58 28 6e 2d 31 29 ┆"nl" ######xp-X(n)=3(X(n)-X(n-1)┆
0x9400…9420 29 2c 20 78 70 3e 58 28 6e 29 2e 0a 54 68 69 73 20 63 6f 6e 64 69 74 69 6f 6e 20 69 73 20 72 65 ┆), xp>X(n). This condition is re┆
0x9420…9440 63 6f 6d 6d 65 6e 64 65 64 20 74 6f 20 67 69 76 65 20 61 20 73 6d 6f 6f 74 68 20 72 75 6e 2d 6f ┆commended to give a smooth run-o┆
0x9440…9460 75 74 2e 22 6e 6c 32 22 0a 22 6c 6d 33 31 22 22 6d 74 20 31 2c 20 20 20 20 20 20 58 20 20 20 20 ┆ut."nl2" "lm31""mt 1, X ┆
0x9460…9480 28 22 0a 63 61 6c 6c 20 76 61 6c 75 65 2c 20 61 72 72 61 79 20 28 31 3a 6e 29 29 20 74 68 65 20 ┆(" call value, array (1:n)) the ┆
0x9480…94a0 78 2d 63 6f 6f 72 64 69 6e 61 74 65 73 20 69 6e 20 69 6e 63 72 65 61 73 69 6e 67 20 6f 72 64 65 ┆x-coordinates in increasing orde┆
0x94a0…94c0 72 2e 22 6e 6c 32 22 0a 22 6d 74 20 31 2c 20 20 20 20 20 20 59 20 20 20 20 28 22 0a 63 61 6c 6c ┆r."nl2" "mt 1, Y (" call┆
0x94c0…94e0 20 76 61 6c 75 65 2c 20 61 72 72 61 79 20 28 31 3a 6e 29 29 20 74 68 65 20 63 6f 72 72 65 73 70 ┆ value, array (1:n)) the corresp┆
0x94e0…9500 6f 6e 64 69 6e 67 20 79 2d 63 6f 6f 72 64 69 6e 61 74 65 73 2e 22 6e 6c 32 22 0a 22 6d 74 20 31 ┆onding y-coordinates."nl2" "mt 1┆
0x9500…9520 2c 20 20 20 20 20 20 4d 20 20 20 20 28 22 0a 72 65 74 75 72 6e 20 76 61 6c 75 65 2c 20 61 72 72 ┆, M (" return value, arr┆
0x9520…9540 61 79 20 28 31 3a 6e 29 29 20 74 68 65 20 63 61 6c 63 75 6c 61 74 65 64 20 73 70 6c 69 6e 65 20 ┆ay (1:n)) the calculated spline ┆
0x9540…9560 70 61 72 61 6d 65 74 65 72 73 2c 20 74 68 65 0a 73 65 63 6f 6e 64 20 64 65 72 69 76 61 74 69 76 ┆parameters, the second derivativ┆
0x9560…9580 65 73 2e 22 6e 6c 32 22 0a 22 6d 74 20 31 2c 20 20 20 20 20 20 6e 20 20 20 20 28 22 0a 63 61 6c ┆es."nl2" "mt 1, n (" cal┆
0x9580…95a0 6c 20 76 61 6c 75 65 2c 20 69 6e 74 65 67 65 72 29 20 74 68 65 20 6e 75 6d 62 65 72 20 6f 66 20 ┆l value, integer) the number of ┆
0x95a0…95c0 70 6f 69 6e 74 73 20 61 62 73 28 6e 29 3e 3d 32 2e 20 41 20 63 61 6c 6c 20 6f 66 20 73 70 6c 6e ┆points abs(n)>=2. A call of spln┆
0x95c0…95e0 33 20 77 69 74 68 0a 6e 3e 3d 32 20 69 6e 69 74 69 61 74 69 7a 65 20 73 70 6c 6e 33 69 6e 76 20 ┆3 with n>=2 initiatize spln3inv ┆
0x95e0…9600 61 6e 64 20 73 70 6c 6e 33 69 6e 74 20 62 75 74 20 74 68 69 73 20 69 73 20 73 75 70 70 72 65 73 ┆and spln3int but this is suppres┆
0x9600…9620 73 65 64 20 62 79 0a 61 20 63 61 6c 6c 20 77 69 74 68 20 6e 3c 3d 2d 32 2e 22 6c 6d 30 22 22 6e ┆sed by a call with n<=-2."lm0""n┆
0x9620…9640 6c 33 22 0a 22 73 6a 22 65 72 72 6f 72 20 6d 65 73 73 61 67 65 3a 0a 2a 2a 2a 73 70 6c 6e 33 20 ┆l3" "sj"error message: ***spln3 ┆
0x9640…9660 6e 3c 32 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 2c 20 63 61 6c 6c 65 64 20 77 69 74 68 ┆n<2 , called with┆
0x9660…9680 20 61 62 73 28 6e 29 3c 32 2e 0a 2a 2a 2a 73 70 6c 6e 33 20 69 6c 6c 65 67 61 6c 20 58 28 69 29 ┆ abs(n)<2. ***spln3 illegal X(i)┆
0x9680…96a0 20 20 20 20 20 20 20 2c 20 58 28 69 29 3e 3d 58 28 6a 29 20 66 6f 72 20 73 6f 6d 65 20 69 3c 6a ┆ , X(i)>=X(j) for some i<j┆
0x96a0…96c0 2e 0a 22 72 6a 22 0a 22 6e 6c 31 22 0a 4c 61 6e 67 75 61 67 65 3a 20 41 4c 47 4f 4c 3b 20 74 65 ┆. "rj" "nl1" Language: ALGOL; te┆
0x96c0…96e0 78 74 20 20 69 6e 20 73 70 6c 6e 33 74 78 74 0a 22 6e 6c 32 22 0a 4d 61 69 6e 74 65 6e 61 6e 63 ┆xt in spln3txt "nl2" Maintenanc┆
0x96e0…9700 65 20 6c 65 76 65 6c 3a 20 75 73 65 72 0a 22 70 73 30 22 0a 22 71 72 22 31 39 37 38 2d 39 2d 32 ┆e level: user "ps0" "qr"1978-9-2┆
0x9700…9720 32 0a 52 65 6e 65 20 4d 6f 73 73 22 72 6a 22 22 63 74 22 22 6e 6c 34 22 0a 72 65 61 6c 20 70 72 ┆2 Rene Moss"rj""ct""nl4" real pr┆
0x9720…9740 6f 63 65 64 75 72 65 20 73 70 6c 6e 33 69 6e 74 22 72 6a 22 22 6e 6c 32 22 0a 66 69 6e 64 73 20 ┆ocedure spln3int"rj""nl2" finds ┆
0x9740…9760 74 68 65 20 69 6e 74 65 67 72 61 6c 20 6f 66 20 61 20 63 75 62 69 63 20 73 70 6c 69 6e 65 20 66 ┆the integral of a cubic spline f┆
0x9760…9780 72 6f 6d 20 58 28 31 29 20 74 6f 20 78 2e 20 54 68 65 20 69 6e 74 65 67 72 61 6c 20 66 72 6f 6d ┆rom X(1) to x. The integral from┆
0x9780…97a0 0a 58 28 31 29 20 74 6f 20 58 28 69 29 20 77 68 65 72 65 20 58 28 69 29 3c 78 3c 58 28 69 2b 31 ┆ X(1) to X(i) where X(i)<x<X(i+1┆
0x97a0…97c0 29 20 69 73 20 73 74 6f 72 65 64 20 69 6e 20 61 6e 20 6f 77 6e 20 76 61 72 69 61 62 6c 65 2c 20 ┆) is stored in an own variable, ┆
0x97c0…97e0 61 6e 64 20 69 6e 69 74 69 61 6c 69 7a 65 64 0a 62 79 20 61 20 63 61 6c 6c 20 6f 66 20 73 70 6c ┆and initialized by a call of spl┆
0x97e0…9800 6e 33 2e 22 6e 6c 32 22 0a 63 61 6c 6c 3a 20 73 70 6c 6e 33 69 6e 74 28 78 2c 58 2c 59 2c 4d 2c ┆n3."nl2" call: spln3int(x,X,Y,M,┆
0x9800…9820 6e 29 0a 22 6e 6c 32 22 22 6c 6d 34 31 22 0a 22 6d 74 20 31 2c 20 20 20 20 20 20 73 70 6c 6e 33 ┆n) "nl2""lm41" "mt 1, spln3┆
0x9820…9840 69 6e 74 20 28 22 0a 72 65 74 75 72 6e 20 76 61 6c 75 65 2c 72 65 61 6c 29 20 74 68 65 20 69 6e ┆int (" return value,real) the in┆
0x9840…9860 74 65 67 72 61 6c 20 66 72 6f 6d 20 58 28 31 29 20 74 6f 20 78 2e 22 6e 6c 32 22 0a 22 6d 74 20 ┆tegral from X(1) to x."nl2" "mt ┆
0x9860…9880 31 2c 20 20 20 20 20 20 78 20 20 20 20 20 20 20 20 28 22 0a 63 61 6c 6c 20 76 61 6c 75 65 2c 20 ┆1, x (" call value, ┆
0x9880…98a0 72 65 61 6c 29 20 74 68 65 20 75 70 70 65 72 20 6c 69 6d 69 74 2e 22 6e 6c 32 22 0a 22 6d 74 20 ┆real) the upper limit."nl2" "mt ┆
0x98a0…98c0 31 2c 20 20 20 20 20 20 58 20 20 20 20 20 20 20 20 28 22 0a 63 61 6c 6c 20 76 61 6c 75 65 2c 20 ┆1, X (" call value, ┆
0x98c0…98e0 61 72 72 61 79 29 20 73 65 65 20 73 70 6c 6e 33 22 6e 6c 32 22 0a 22 6d 74 20 31 2c 20 20 20 20 ┆array) see spln3"nl2" "mt 1, ┆
0x98e0…9900 20 20 59 20 20 20 20 20 20 20 20 28 22 0a 63 61 6c 6c 20 76 61 6c 75 65 2c 20 61 72 72 61 79 29 ┆ Y (" call value, array)┆
0x9900…9920 (17,) 20 73 65 65 20 73 70 6c 6e 33 22 6e 6c 32 22 0a 22 6d 74 20 31 2c 20 20 20 20 20 20 4d 20 20 20 ┆ see spln3"nl2" "mt 1, M ┆
0x9920…9940 20 20 20 20 20 28 22 0a 63 61 6c 6c 20 76 61 6c 75 65 2c 20 61 72 72 61 79 29 20 73 65 65 20 73 ┆ (" call value, array) see s┆
0x9940…9960 70 6c 6e 33 22 6e 6c 32 22 0a 22 6d 74 20 31 2c 20 20 20 20 20 20 6e 20 20 20 20 20 20 20 20 28 ┆pln3"nl2" "mt 1, n (┆
0x9960…9980 22 0a 63 61 6c 6c 20 76 61 6c 75 65 2c 20 69 6e 74 65 67 65 72 29 20 73 65 65 20 73 70 6c 6e 33 ┆" call value, integer) see spln3┆
0x9980…99a0 22 6c 6d 30 22 22 6e 6c 33 22 0a 22 73 6a 22 65 72 72 6f 72 20 6d 65 73 73 61 67 65 3a 0a 2a 2a ┆"lm0""nl3" "sj"error message: **┆
0x99a0…99c0 2a 73 70 6c 6e 33 69 6e 74 20 69 6c 6c 65 67 61 6c 20 78 20 3d 20 3c 76 61 6c 75 65 3e 20 20 20 ┆*spln3int illegal x = <value> ┆
0x99c0…99e0 20 20 20 20 2c 20 78 3c 58 28 31 29 20 6f 72 20 78 3e 58 28 6e 29 0a 22 72 6a 22 22 6e 6c 31 22 ┆ , x<X(1) or x>X(n) "rj""nl1"┆
0x99e0…9a00 0a 4c 61 6e 67 75 61 67 65 3a 20 41 4c 47 4f 4c 3b 20 74 65 78 74 20 69 6e 20 73 70 6c 6e 33 74 ┆ Language: ALGOL; text in spln3t┆
0x9a00…9a20 78 74 22 6e 6c 32 22 0a 4d 61 69 6e 74 65 6e 61 6e 63 65 20 6c 65 76 65 6c 3a 20 75 73 65 72 0a ┆xt"nl2" Maintenance level: user ┆
0x9a20…9a40 22 70 73 30 22 0a 22 71 72 22 31 39 37 38 2d 39 2d 32 30 22 6e 6c 22 52 65 6e 65 20 4d 6f 73 73 ┆"ps0" "qr"1978-9-20"nl"Rene Moss┆
0x9a40…9a60 22 72 6a 22 22 6e 6c 34 22 22 63 74 22 0a 62 6f 6f 6c 65 61 6e 20 70 72 6f 63 65 64 75 72 65 20 ┆"rj""nl4""ct" boolean procedure ┆
0x9a60…9a80 73 70 6c 6e 33 69 6e 76 22 72 6a 22 22 6e 6c 32 22 0a 63 61 6c 63 75 6c 61 74 65 73 20 74 68 65 ┆spln3inv"rj""nl2" calculates the┆
0x9a80…9aa0 20 69 6e 76 65 72 73 65 20 66 75 6e 63 74 69 6f 6e 20 74 6f 20 74 68 65 20 63 75 62 69 63 20 73 ┆ inverse function to the cubic s┆
0x9aa0…9ac0 70 6c 69 6e 65 2e 20 47 69 76 65 6e 0a 74 68 65 20 6f 72 64 69 6e 61 74 65 20 79 20 74 68 65 20 ┆pline. Given the ordinate y the ┆
0x9ac0…9ae0 70 72 6f 63 65 64 75 72 65 20 6c 6f 6f 6b 73 20 66 6f 72 20 74 68 65 20 73 6d 61 6c 6c 65 73 74 ┆procedure looks for the smallest┆
0x9ae0…9b00 0a 69 6e 76 65 72 73 65 20 76 61 6c 75 65 20 69 6e 20 74 68 65 20 69 6e 74 65 72 76 61 6c 20 66 ┆ inverse value in the interval f┆
0x9b00…9b20 72 6f 6d 20 78 20 74 6f 20 58 28 6e 29 2e 20 54 68 65 20 69 6e 74 65 72 76 61 6c 0a 69 73 20 6f ┆rom x to X(n). The interval is o┆
0x9b20…9b40 70 65 6e 20 74 6f 20 74 68 65 20 6c 65 66 74 20 69 66 20 78 20 65 71 75 61 6c 73 20 74 68 65 20 ┆pen to the left if x equals the ┆
0x9b40…9b60 6c 61 74 65 73 74 20 64 65 74 65 72 6d 69 6e 65 64 0a 69 6e 76 65 72 73 65 20 76 61 6c 75 65 20 ┆latest determined inverse value ┆
0x9b60…9b80 74 6f 20 79 20 65 6c 73 65 20 74 68 65 20 69 6e 74 65 72 76 61 6c 20 69 73 20 63 6c 6f 73 65 64 ┆to y else the interval is closed┆
0x9b80…9ba0 2e 20 49 66 20 74 68 65 72 65 20 65 78 69 73 74 20 61 20 73 6f 6c 75 74 69 6f 6e 2c 0a 74 68 69 ┆. If there exist a solution, thi┆
0x9ba0…9bc0 73 20 69 73 20 64 65 6c 69 76 65 72 65 64 20 69 6e 20 78 20 61 6e 64 20 74 68 65 20 70 72 6f 63 ┆s is delivered in x and the proc┆
0x9bc0…9be0 65 64 75 72 65 20 69 73 20 73 65 74 20 74 6f 20 74 72 75 65 2e 22 6e 6c 32 22 0a 63 61 6c 6c 3a ┆edure is set to true."nl2" call:┆
0x9be0…9c00 20 73 70 6c 6e 33 69 6e 76 28 78 2c 79 2c 58 2c 59 2c 4d 2c 6e 29 22 6e 6c 32 22 22 6c 6d 34 31 ┆ spln3inv(x,y,X,Y,M,n)"nl2""lm41┆
0x9c00…9c20 22 0a 22 6d 74 20 31 2c 20 20 20 20 20 20 73 70 6c 6e 33 69 6e 76 20 28 22 0a 72 65 74 75 72 6e ┆" "mt 1, spln3inv (" return┆
0x9c20…9c40 20 76 61 6c 75 65 2c 20 62 6f 6f 6c 65 61 6e 29 20 74 72 75 65 2c 20 69 66 20 61 6e 20 69 6e 76 ┆ value, boolean) true, if an inv┆
0x9c40…9c60 65 72 73 65 20 76 61 6c 75 65 20 65 78 69 73 74 73 20 69 6e 0a 74 68 65 20 69 6e 74 65 72 76 61 ┆erse value exists in the interva┆
0x9c60…9c80 6c 20 66 72 6f 6d 20 78 20 74 6f 20 58 28 6e 29 2e 22 6e 6c 32 22 0a 22 6d 74 20 31 2c 20 20 20 ┆l from x to X(n)."nl2" "mt 1, ┆
0x9c80…9ca0 20 20 20 78 20 20 20 20 20 20 20 20 28 22 0a 63 61 6c 6c 20 61 6e 64 20 72 65 74 75 72 6e 20 76 ┆ x (" call and return v┆
0x9ca0…9cc0 61 6c 75 65 2c 20 72 65 61 6c 29 20 61 74 20 63 61 6c 6c 20 74 68 65 20 6c 6f 77 65 72 20 69 6e ┆alue, real) at call the lower in┆
0x9cc0…9ce0 74 65 72 76 61 6c 2c 0a 61 74 20 72 65 74 75 72 6e 20 74 68 65 20 69 6e 76 65 72 73 65 20 76 61 ┆terval, at return the inverse va┆
0x9ce0…9d00 6c 75 65 20 69 66 20 69 74 20 65 78 69 73 74 73 2e 20 58 28 31 29 3c 78 3c 58 28 6e 29 0a 6f 74 ┆lue if it exists. X(1)<x<X(n) ot┆
0x9d00…9d20 68 65 72 77 69 73 65 20 61 6e 20 61 6c 61 72 6d 20 69 73 20 67 69 76 65 6e 2e 22 6e 6c 32 22 0a ┆herwise an alarm is given."nl2" ┆
0x9d20…9d40 22 6d 74 20 31 2c 20 20 20 20 20 20 79 20 20 20 20 20 20 20 20 28 22 0a 63 61 6c 6c 20 76 61 6c ┆"mt 1, y (" call val┆
0x9d40…9d60 75 65 2c 20 72 65 61 6c 29 20 74 68 65 20 6f 72 64 69 6e 61 74 65 2e 22 6e 6c 32 22 0a 22 6d 74 ┆ue, real) the ordinate."nl2" "mt┆
0x9d60…9d80 20 31 2c 20 20 20 20 20 20 58 20 20 20 20 20 20 20 20 28 22 0a 63 61 6c 6c 20 76 61 6c 75 65 2c ┆ 1, X (" call value,┆
0x9d80…9da0 20 61 72 72 61 79 29 20 73 65 65 20 73 70 6c 6e 33 2e 22 6e 6c 32 22 0a 22 6d 74 20 31 2c 20 20 ┆ array) see spln3."nl2" "mt 1, ┆
0x9da0…9dc0 20 20 20 20 59 20 20 20 20 20 20 20 20 28 22 0a 63 61 6c 6c 20 76 61 6c 75 65 2c 20 61 72 72 61 ┆ Y (" call value, arra┆
0x9dc0…9de0 79 29 20 73 65 65 20 73 70 6c 6e 33 2e 22 6e 6c 32 22 0a 22 6d 74 20 31 2c 20 20 20 20 20 20 4d ┆y) see spln3."nl2" "mt 1, M┆
0x9de0…9e00 20 20 20 20 20 20 20 20 28 22 0a 63 61 6c 6c 20 76 61 6c 75 65 2c 20 61 72 72 61 79 29 20 73 65 ┆ (" call value, array) se┆
0x9e00…9e20 65 20 73 70 6c 6e 33 2e 22 6e 6c 32 22 0a 22 6d 74 20 31 2c 20 20 20 20 20 20 6e 20 20 20 20 20 ┆e spln3."nl2" "mt 1, n ┆
0x9e20…9e40 20 20 20 28 22 0a 63 61 6c 6c 20 76 61 6c 75 65 2c 20 69 6e 74 65 67 65 72 29 20 74 68 65 20 6e ┆ (" call value, integer) the n┆
0x9e40…9e60 75 6d 62 65 72 20 6f 66 20 73 70 6c 69 6e 65 70 6f 69 6e 74 73 2e 20 46 6f 72 20 20 6e 3c 30 0a ┆umber of splinepoints. For n<0 ┆
0x9e60…9e80 74 68 65 20 70 72 6f 63 65 64 75 72 65 20 69 73 20 69 6e 69 74 69 61 6c 69 7a 65 64 20 62 65 66 ┆the procedure is initialized bef┆
0x9e80…9ea0 6f 72 65 20 69 74 20 69 73 20 63 61 6c 6c 65 64 20 77 69 74 68 20 6e 3d 61 62 73 28 6e 29 3b 0a ┆ore it is called with n=abs(n); ┆
0x9ea0…9ec0 6e 3d 30 20 6a 75 73 74 20 69 6e 69 74 69 61 6c 69 7a 65 20 69 2e 65 2e 20 74 68 65 20 6b 6e 6f ┆n=0 just initialize i.e. the kno┆
0x9ec0…9ee0 77 6c 65 64 67 65 20 6f 66 20 74 68 65 20 6c 61 74 65 73 74 0a 73 6f 6c 75 74 69 6f 6e 20 69 73 ┆wledge of the latest solution is┆
0x9ee0…9f00 20 6c 6f 73 74 2e 20 41 66 74 65 72 20 61 20 6e 6f 72 6d 61 6c 20 63 61 6c 6c 20 6f 66 20 73 70 ┆ lost. After a normal call of sp┆
0x9f00…9f20 6c 6e 33 2c 20 73 70 6c 6e 33 69 6e 76 0a 69 73 20 69 6e 69 74 69 61 6c 69 7a 65 64 2e 22 6c 6d ┆ln3, spln3inv is initialized."lm┆
0x9f20…9f40 30 22 22 6e 6c 33 22 0a 22 73 6a 22 65 72 72 6f 72 20 6d 65 73 73 61 67 65 3a 0a 2a 2a 2a 73 70 ┆0""nl3" "sj"error message: ***sp┆
0x9f40…9f60 6c 6e 33 69 6e 76 20 69 6c 6c 65 67 61 6c 20 78 20 3d 20 3c 76 61 6c 75 65 3e 20 20 20 20 20 20 ┆ln3inv illegal x = <value> ┆
0x9f60…9f80 20 2c 20 78 3c 58 28 31 29 20 6f 72 20 78 3e 58 28 6e 29 0a 22 72 6a 22 0a 22 6e 6c 31 22 0a 4c ┆ , x<X(1) or x>X(n) "rj" "nl1" L┆
0x9f80…9fa0 61 6e 67 75 61 67 65 3a 20 41 4c 47 4f 4c 3b 20 74 65 78 74 20 69 6e 20 73 70 6c 6e 33 74 78 74 ┆anguage: ALGOL; text in spln3txt┆
0x9fa0…9fc0 0a 22 6e 6c 32 22 0a 4d 61 69 6e 74 65 6e 61 6e 63 65 20 6c 65 76 65 6c 3a 20 75 73 65 72 0a 22 ┆ "nl2" Maintenance level: user "┆
0x9fc0…9fe0 70 73 30 22 0a 22 71 72 22 31 39 37 38 2d 39 2d 32 32 0a 52 65 6e 65 20 4d 6f 73 73 22 72 6a 22 ┆ps0" "qr"1978-9-22 Rene Moss"rj"┆
0x9fe0…a000 22 6e 6c 34 22 22 63 74 22 0a 62 6f 6f 6c 65 61 6e 20 70 72 6f 63 65 64 75 72 65 20 73 70 6c 6e ┆"nl4""ct" boolean procedure spln┆
0xa000…a020 33 73 6d 74 68 22 72 6a 22 22 6e 6c 32 22 0a 63 68 61 6e 67 65 73 20 74 68 65 20 59 2d 63 6f 6f ┆3smth"rj""nl2" changes the Y-coo┆
0xa020…a040 72 64 69 6e 61 74 65 73 20 6f 66 20 61 20 73 65 74 20 6f 66 20 70 6f 69 6e 74 73 20 73 75 63 68 ┆rdinates of a set of points such┆
0xa040…a060 20 74 68 61 74 20 0a 61 20 73 70 6c 69 6e 65 20 64 72 61 77 6e 20 74 68 72 6f 75 67 68 20 74 68 ┆ that a spline drawn through th┆
0xa060…a080 65 73 65 20 6e 65 77 20 70 6f 69 6e 74 73 20 69 73 20 73 6d 6f 6f 74 68 65 72 2c 0a 69 2e 65 2e ┆ese new points is smoother, i.e.┆
0xa080…a0a0 20 74 68 65 20 73 71 75 61 72 65 20 64 65 76 69 61 74 69 6f 6e 20 64 69 76 69 64 65 64 20 62 79 ┆ the square deviation divided by┆
0xa0a0…a0c0 20 74 68 65 20 76 61 72 69 61 6e 63 65 28 53 44 28 69 29 2a 2a 32 29 0a 23 2b 23 74 68 65 20 63 ┆ the variance(SD(i)**2) #+#the c┆
0xa0c0…a0e0 75 72 76 61 74 75 72 65 20 69 73 20 6d 69 6e 69 6d 69 7a 65 64 2e 22 6e 6c 32 22 0a 63 61 6c 6c ┆urvature is minimized."nl2" call┆
0xa0e0…a100 3a 20 73 70 6c 6e 33 73 6d 74 68 28 58 2c 59 53 2c 53 44 2c 6e 29 22 6e 6c 32 22 0a 22 6c 6d 34 ┆: spln3smth(X,YS,SD,n)"nl2" "lm4┆
0xa100…a120 33 22 0a 22 6d 74 20 31 2c 20 20 20 20 20 20 73 70 6c 6e 33 73 6d 74 68 20 28 22 0a 72 65 74 75 ┆3" "mt 1, spln3smth (" retu┆
0xa120…a140 72 6e 20 76 61 6c 75 65 2c 20 62 6f 6f 6c 65 61 6e 29 20 74 72 75 65 20 69 66 20 73 6d 6f 6f 74 ┆rn value, boolean) true if smoot┆
0xa140…a160 68 69 6e 67 20 77 61 73 20 70 6f 73 73 69 62 6c 65 0a 61 6e 64 20 66 61 6c 73 65 20 6f 74 68 65 ┆hing was possible and false othe┆
0xa160…a180 72 77 69 73 65 2e 20 49 6e 20 74 68 69 73 20 63 61 73 65 20 59 53 20 69 73 20 75 6e 63 68 61 6e ┆rwise. In this case YS is unchan┆
0xa180…a1a0 67 65 64 2e 22 6e 6c 32 22 0a 22 6d 74 20 31 2c 20 20 20 20 20 20 58 20 20 20 20 20 20 20 20 20 ┆ged."nl2" "mt 1, X ┆
0xa1a0…a1c0 28 22 0a 63 61 6c 6c 20 76 61 6c 75 65 2c 20 61 72 72 61 79 28 31 3a 6e 29 29 2e 20 78 2d 63 6f ┆(" call value, array(1:n)). x-co┆
0xa1c0…a1e0 6f 72 64 69 6e 61 74 65 73 2c 20 58 28 69 29 2c 0a 69 6e 20 69 6e 63 72 65 61 73 69 6e 67 20 6f ┆ordinates, X(i), in increasing o┆
0xa1e0…a200 72 64 65 72 2e 22 6e 6c 32 22 0a 22 6d 74 20 31 2c 20 20 20 20 20 20 59 20 20 20 20 20 20 20 20 ┆rder."nl2" "mt 1, Y ┆
0xa200…a220 (18,) 20 28 22 0a 63 61 6c 6c 20 61 6e 64 20 72 65 74 75 72 6e 20 76 61 6c 75 65 2c 20 61 72 72 61 79 ┆ (" call and return value, array┆
0xa220…a240 28 31 3a 6e 29 29 20 74 68 65 20 79 2d 63 6f 6f 72 64 69 6e 61 74 65 73 2c 0a 59 28 69 29 2e 22 ┆(1:n)) the y-coordinates, Y(i)."┆
0xa240…a260 6e 6c 32 22 0a 22 6d 74 20 31 2c 20 20 20 20 20 20 53 44 20 20 20 20 20 20 20 20 28 22 0a 63 61 ┆nl2" "mt 1, SD (" ca┆
0xa260…a280 6c 6c 20 76 61 6c 75 65 2c 20 61 72 72 61 79 28 31 3a 6e 29 29 20 74 68 65 20 73 74 61 6e 64 61 ┆ll value, array(1:n)) the standa┆
0xa280…a2a0 72 64 20 64 65 76 69 61 74 69 6f 6e 0a 6f 66 20 74 68 65 20 70 6f 69 6e 74 73 2c 20 53 44 28 69 ┆rd deviation of the points, SD(i┆
0xa2a0…a2c0 29 2e 22 6e 6c 32 22 0a 22 6d 74 20 31 2c 20 20 20 20 20 20 6e 20 20 20 20 20 20 20 20 20 28 22 ┆)."nl2" "mt 1, n ("┆
0xa2c0…a2e0 0a 63 61 6c 6c 20 76 61 6c 75 65 20 69 6e 74 65 67 65 72 29 20 74 68 65 20 6e 75 6d 62 65 72 20 ┆ call value integer) the number ┆
0xa2e0…a300 6f 66 20 70 6f 69 6e 74 73 2e 20 33 3c 6e 3c 32 39 0a 6f 6e 20 52 43 34 30 30 30 20 61 74 20 74 ┆of points. 3<n<29 on RC4000 at t┆
0xa300…a320 68 65 20 6d 6f 6d 65 6e 74 2e 20 44 65 70 65 6e 64 73 20 6f 6e 20 73 69 7a 65 20 6f 66 20 63 6f ┆he moment. Depends on size of co┆
0xa320…a340 72 65 73 74 6f 72 65 2e 22 6c 6d 30 22 0a 22 6e 6c 33 22 0a 22 73 6a 22 65 72 72 6f 72 20 6d 65 ┆restore."lm0" "nl3" "sj"error me┆
0xa340…a360 73 73 61 67 65 3a 0a 2a 2a 2a 73 70 6c 6e 33 73 6d 74 68 20 6e 3c 34 20 20 20 20 20 20 20 20 20 ┆ssage: ***spln3smth n<4 ┆
0xa360…a380 20 20 20 20 2c 20 6e 3c 34 0a 2a 2a 2a 73 70 6c 6e 33 20 69 6c 6c 65 67 61 6c 20 58 28 69 29 20 ┆ , n<4 ***spln3 illegal X(i) ┆
0xa380…a3a0 20 20 20 20 20 20 20 2c 20 58 28 69 29 3e 3d 58 28 6a 29 20 66 6f 72 20 73 6f 6d 65 20 69 3c 6a ┆ , X(i)>=X(j) for some i<j┆
0xa3a0…a3c0 2e 0a 22 72 6a 22 22 6e 6c 31 22 0a 4c 61 6e 67 75 61 67 65 3a 20 41 4c 47 4f 4c 3b 20 74 65 78 ┆. "rj""nl1" Language: ALGOL; tex┆
0xa3c0…a3e0 74 20 69 6e 20 73 70 6c 6e 33 74 78 74 22 6e 6c 32 22 0a 4d 61 69 6e 74 65 6e 61 6e 63 65 20 6c ┆t in spln3txt"nl2" Maintenance l┆
0xa3e0…a400 65 76 65 6c 3a 20 75 73 65 72 22 70 73 30 22 0a 22 71 72 22 31 39 37 38 2d 39 2d 32 32 0a 52 65 ┆evel: user"ps0" "qr"1978-9-22 Re┆
0xa400…a420 6e 65 20 4d 6f 73 73 22 72 6a 22 22 6e 6c 34 22 22 63 74 22 0a 72 65 61 6c 20 70 72 6f 63 65 64 ┆ne Moss"rj""nl4""ct" real proced┆
0xa420…a440 75 72 65 20 73 70 6c 6e 33 76 61 6c 20 0a 72 65 61 6c 20 70 72 6f 63 65 64 75 72 65 20 73 70 6c ┆ure spln3val real procedure spl┆
0xa440…a460 6e 33 64 69 66 20 0a 72 65 61 6c 20 70 72 6f 63 65 64 75 72 65 20 73 70 6c 6e 33 64 69 66 32 22 ┆n3dif real procedure spln3dif2"┆
0xa460…a480 72 6a 22 22 6e 6c 32 22 0a 66 69 6e 64 20 74 68 65 20 76 61 6c 75 65 2c 20 31 73 74 20 64 65 72 ┆rj""nl2" find the value, 1st der┆
0xa480…a4a0 69 76 61 74 69 76 65 2c 20 32 6e 64 20 64 65 72 69 76 61 74 69 76 65 20 6f 66 20 61 20 63 75 62 ┆ivative, 2nd derivative of a cub┆
0xa4a0…a4c0 69 63 0a 73 70 6c 69 6e 65 20 61 74 20 78 2e 22 6e 6c 32 22 0a 63 61 6c 6c 3a 20 73 70 6c 6e 33 ┆ic spline at x."nl2" call: spln3┆
0xa4c0…a4e0 76 61 6c 28 78 2c 58 2c 59 2c 4d 2c 6e 29 22 6e 6c 22 0a 22 6c 6d 31 35 22 0a 73 70 6c 6e 33 64 ┆val(x,X,Y,M,n)"nl" "lm15" spln3d┆
0xa4e0…a500 69 66 28 78 2c 58 2c 59 2c 4d 2c 6e 29 22 6e 6c 22 0a 73 70 6c 6e 33 64 69 66 32 28 78 2c 58 2c ┆if(x,X,Y,M,n)"nl" spln3dif2(x,X,┆
0xa500…a520 59 2c 4d 2c 6e 29 22 6e 6c 32 22 0a 22 6c 6d 34 33 22 0a 22 6d 74 20 31 2c 20 20 20 20 20 20 73 ┆Y,M,n)"nl2" "lm43" "mt 1, s┆
0xa520…a540 70 6c 6e 33 76 61 6c 20 20 28 22 0a 72 65 74 75 72 6e 20 76 61 6c 75 65 2c 20 72 65 61 6c 29 20 ┆pln3val (" return value, real) ┆
0xa540…a560 74 68 65 20 76 61 6c 75 65 20 61 74 20 78 2e 22 6e 6c 32 22 0a 22 6d 74 20 31 2c 20 20 20 20 20 ┆the value at x."nl2" "mt 1, ┆
0xa560…a580 20 73 70 6c 6e 33 64 69 66 20 20 28 22 0a 72 65 74 75 72 6e 20 76 61 6c 75 65 2c 20 72 65 61 6c ┆ spln3dif (" return value, real┆
0xa580…a5a0 29 20 31 73 74 20 64 65 72 69 76 61 74 69 76 65 20 61 74 20 78 2e 22 6e 6c 32 22 0a 22 6d 74 20 ┆) 1st derivative at x."nl2" "mt ┆
0xa5a0…a5c0 31 2c 20 20 20 20 20 20 73 70 6c 6e 33 64 69 66 32 20 28 22 0a 72 65 74 75 72 6e 20 76 61 6c 75 ┆1, spln3dif2 (" return valu┆
0xa5c0…a5e0 65 2c 20 72 65 61 6c 29 20 32 6e 64 20 64 65 72 69 76 61 74 69 76 65 20 61 74 20 78 2e 22 6e 6c ┆e, real) 2nd derivative at x."nl┆
0xa5e0…a600 32 22 0a 22 6d 74 20 31 2c 20 20 20 20 20 20 78 20 20 20 20 20 20 20 20 20 28 22 0a 63 61 6c 6c ┆2" "mt 1, x (" call┆
0xa600…a620 20 76 61 6c 75 65 2c 20 72 65 61 6c 29 20 58 28 31 29 3c 3d 78 3c 3d 58 28 6e 29 2e 22 6e 6c 32 ┆ value, real) X(1)<=x<=X(n)."nl2┆
0xa620…a640 22 0a 22 6d 74 20 31 2c 20 20 20 20 20 20 58 20 20 20 20 20 20 20 20 20 28 22 0a 63 61 6c 6c 20 ┆" "mt 1, X (" call ┆
0xa640…a660 76 61 6c 75 65 2c 20 61 72 72 61 79 29 20 73 65 65 20 73 70 6c 6e 33 2e 22 6e 6c 32 22 0a 22 6d ┆value, array) see spln3."nl2" "m┆
0xa660…a680 74 20 31 2c 20 20 20 20 20 20 59 20 20 20 20 20 20 20 20 20 28 22 0a 63 61 6c 6c 20 76 61 6c 75 ┆t 1, Y (" call valu┆
0xa680…a6a0 65 2c 20 61 72 72 61 79 29 20 73 65 65 20 73 70 6c 6e 33 2e 22 6e 6c 32 22 0a 22 6d 74 20 31 2c ┆e, array) see spln3."nl2" "mt 1,┆
0xa6a0…a6c0 20 20 20 20 20 20 4d 20 20 20 20 20 20 20 20 20 28 22 0a 63 61 6c 6c 20 76 61 6c 75 65 2c 20 61 ┆ M (" call value, a┆
0xa6c0…a6e0 72 72 61 79 29 20 73 65 65 20 73 70 6c 6e 33 2e 22 6e 6c 32 22 0a 22 6d 74 20 31 2c 20 20 20 20 ┆rray) see spln3."nl2" "mt 1, ┆
0xa6e0…a700 20 20 6e 20 20 20 20 20 20 20 20 20 28 22 0a 63 61 6c 6c 20 76 61 6c 75 65 2c 20 69 6e 74 65 67 ┆ n (" call value, integ┆
0xa700…a720 65 72 29 20 73 65 65 20 73 70 6c 6e 33 2e 22 6c 6d 30 22 22 6e 6c 33 22 0a 22 73 6a 22 65 72 72 ┆er) see spln3."lm0""nl3" "sj"err┆
0xa720…a740 6f 72 20 6d 65 73 73 61 67 65 3a 0a 2a 2a 2a 73 70 6c 6e 33 76 61 6c 20 69 6c 6c 65 67 61 6c 20 ┆or message: ***spln3val illegal ┆
0xa740…a760 78 20 3d 20 3c 76 61 6c 75 65 3e 20 20 20 20 20 20 20 20 2c 20 78 3c 58 28 31 29 20 6f 72 20 78 ┆x = <value> , x<X(1) or x┆
0xa760…a780 3e 58 28 6e 29 0a 2a 2a 2a 73 70 6c 6e 33 64 69 66 20 69 6c 6c 65 67 61 6c 20 78 20 3d 20 3c 76 ┆>X(n) ***spln3dif illegal x = <v┆
0xa780…a7a0 61 6c 75 65 3e 20 20 20 20 20 20 20 20 2c 20 20 20 20 20 20 20 20 2d 20 20 20 0a 2a 2a 2a 73 70 ┆alue> , - ***sp┆
0xa7a0…a7c0 6c 6e 33 64 69 66 32 20 69 6c 6c 65 67 61 6c 20 78 20 3d 20 3c 76 61 6c 75 65 3e 20 20 20 20 20 ┆ln3dif2 illegal x = <value> ┆
0xa7c0…a7e0 20 20 2c 20 20 20 20 20 20 20 20 2d 0a 22 72 6a 22 22 6e 6c 31 22 0a 4c 61 6e 67 75 61 67 65 3a ┆ , - "rj""nl1" Language:┆
0xa7e0…a800 20 41 4c 47 4f 4c 22 6e 6c 32 22 0a 4d 61 69 6e 74 65 6e 61 6e 63 65 20 6c 65 76 65 6c 3a 20 75 ┆ ALGOL"nl2" Maintenance level: u┆