mirror of
https://github.com/retro-software/B5500-software.git
synced 2026-03-02 17:44:40 +00:00
2c72f7fd1d
1. Commit library tape images, directories, and extracted text files. 2. Commit additional utilities under Unisys-Emode-Tools.
512 lines
40 KiB
Plaintext
512 lines
40 KiB
Plaintext
BEGIN
|
|
COMMENT ED SEIDMAN ;
|
|
BEGIN SRFT 1
|
|
COMMENT MULTIDIMENSIONAL SURFACE FIT PROGRAM. THIS METHOD USES SRFT 2
|
|
ORTHOGONAL POLYNOMIALS WHICH PERMIT TRUNCATION WITHOUT SRFT 3
|
|
RECOMPUTING THE COEFFICIENTS. A NEW VALUE FOR THE ERROR ISSRFT 4
|
|
COMPUTED FOR EACH TRUNCATION PASS. ; SRFT 5
|
|
BOOLEAN SW, ST ; SRFT 6
|
|
FILE IN CARD (2,10) ; SRFT 8
|
|
FILE OUT PAPER 1(3,15), SRFT 9
|
|
PUNCH 0(2,10) ; SRFT 10
|
|
INTEGER KMAX, PMAX, V, N ; SRFT 16
|
|
LIST SIZESI (KMAX, PMAX, V, N, ST, SW) ; SRFT 17
|
|
FORMAT IN CFI (6I10) ; SRFT 18
|
|
READ (CARD, CFI, SIZESI) ; SRFT 19
|
|
BEGIN SRFT 20
|
|
INTEGER I, J, D, M, X, T, T1, T2, T3, KM ; SRFT 21
|
|
INTEGER ARRAY NAME[0:V], K[0:V], P[0:V], DEGREE[0:V], DEG[0:V], SRFT 22
|
|
NEWP[0:V], TOTAL[0:V], TALLY[0:V] ; SRFT 23
|
|
REAL ALPHATEMP, BETATEMP, DIFF, RHOTEMP, SUM, SUMS ; SRFT 24
|
|
REAL ARRAY TT[0:V], DATA[0:KMAX], POWERS[0:KMAX], ALEPH[0:PMAX], SRFT 25
|
|
BETA[0:PMAX], PHI[0:N], RHO[0:V, 0:PMAX+1], SRFT 26
|
|
POLY[0:V | KMAX | (PMAX+1)], COEFF[0:3 | (PMAX+1)] ; SRFT 27
|
|
LIST SIZESF (FOR M ~ 1 STEP 1 UNTIL V DO SRFT 33
|
|
[NAME[M], K[M], P[M]]), SRFT 34
|
|
POINTS (FOR D ~ 1 STEP 1 UNTIL K[M] DO SRFT 35
|
|
DATA[D]), SRFT 36
|
|
HO (M), SRFT 37
|
|
ALPHABETA (ALEPH[I], BETA[I]), SRFT 38
|
|
RHOVALUES (FOR D ~ 1 STEP 1 UNTIL P[M]+1 DO SRFT 39
|
|
RHO[M,D]), SRFT 40
|
|
POLYS (FOR D ~ T3 STEP 1 UNTIL K[M]+T3-1 DO SRFT 41
|
|
POLY[D]), SRFT 42
|
|
PABO (M, P[M], FOR D ~ 1 STEP 1 UNTIL P[M] DO SRFT 43
|
|
[ALEPH[D], BETA[D]]) ; SRFT 44
|
|
FORMAT IN CFF (3I10), SRFT 66
|
|
CF2 (4E20.11) ; SRFT 67
|
|
FORMAT OUT HEAD (X42, "MULTIDIMENSIONAL SURFACE FIT PROGRAM" //), SRFT 68
|
|
HF (// X50, "VARIABLE", I5 //), SRFT 69
|
|
TITLEAB (X61, "ALPHA[I]", X13, "BETA[I]" //), SRFT 70
|
|
OOT (X50, 2E20.11 /), SRFT 71
|
|
TITLER (/ X50, "RHO[I]" //), SRFT 72
|
|
CF3 (6E20.11), SRFT 73
|
|
TITLEP (// X50, "POLYNOMIALS" //), SRFT 74
|
|
PABF (2I5 / (2E20.11)) ; SRFT 75
|
|
READ (CARD, CFF, SIZESF) ; SRFT 78
|
|
WRITE (PAPER[PAGE], OOT) ;
|
|
WRITE (PAPER, HEAD) ; SRFT 79
|
|
D ~ 1 ; SRFT 28
|
|
FOR I ~ 1 STEP 1 UNTIL V DO SRFT 29
|
|
D ~ (P[I] + 1) | D ; SRFT 30
|
|
BEGIN SRFT 31
|
|
REAL ARRAY COEFFICIENTS[0:D] ; SRFT 32
|
|
LIST MAG (ALPHATEMP, BETATEMP, D, DIFF, I, J, M, N, V, SRFT 45
|
|
PMAX, RHOTEMP, SUM, SUMS, T, T1, T2, T3, X, KMAX, KM, SRFT 46
|
|
FOR I ~ 1 STEP 1 UNTIL V DO SRFT 47
|
|
[DEG[I], DEGREE[I], K[I], NAME[I], NEWP[I], P[I], TT[I], SRFT 48
|
|
TALLY[I], TOTAL[I]], SRFT 49
|
|
FOR I ~ 1 STEP 1 UNTIL KMAX DO SRFT 50
|
|
[DATA[I], POWERS [I]], SRFT 51
|
|
FOR I ~ 1 STEP 1 UNTIL PMAX DO SRFT 52
|
|
[ALEPH[I], BETA[I]], SRFT 53
|
|
FOR I ~ 1 STEP 1 UNTIL V DO SRFT 54
|
|
FOR J ~ 1 STEP 1 UNTIL PMAX+1 DO SRFT 55
|
|
RHO[I,J], SRFT 56
|
|
FOR I ~ 1 STEP 1 UNTIL 3 | (PMAX+1) DO SRFT 59
|
|
COEFF[I]), SRFT 60
|
|
MAGPO (FOR I ~ 1 STEP 1 UNTIL X-1 DO SRFT 57
|
|
POLY[I]), SRFT 58
|
|
MAGC (FOR I ~ 1 STEP 1 UNTIL T3 DO SRFT 61
|
|
COEFFICIENTS[I]), SRFT 62
|
|
MAGPH (FOR I ~ 1 STEP 1 UNTIL N DO SRFT 63
|
|
PHI[I]) ; SRFT 64
|
|
LABEL VARIABLE, MAJORLOOP, RESET, LAST, FURTHER, ONE, TRUNCATE, SRFT 76
|
|
COEFFICIENTC, THEEND, REPLACE ; SRFT 77
|
|
IF SW THEN GO TO REPLACE ; SRFT 14
|
|
TOTAL[1] ~ X ~ I ~ T3 ~ 1 ; SRFT 80
|
|
FOR M ~ 2 STEP 1 UNTIL V DO SRFT 81
|
|
TOTAL[M] ~ (P[M-1] + 1) | K[M-1] + TOTAL[M-1] ; SRFT 82
|
|
M ~ 1 ; SRFT 83
|
|
COMMENT THE FOLLOWING SECTION COMPUTES THE RECU-SIVE COEFFICIENTS SRFT 84
|
|
NEEDED FOR THE ORTHOGONAL POLYNOMIALS, AND THEN THE SRFT 85
|
|
POLYNOMIALS ; SRFT 86
|
|
VARIABLE : WRITE (PAPER, HF, HO) ; SRFT 87
|
|
WRITE (PAPER, TITLEAB) ; SRFT 88
|
|
READ (CARD, CF2, POINTS) ; SRFT 89
|
|
BETA[1] ~ 0 ; KM ~ K[M] ; SRFT 90
|
|
FOR J ~ 1 STEP 1 UNTIL KM DO SRFT 91
|
|
BEGIN SRFT 92
|
|
POWERS[J] ~ POLY[X] ~ 1.0 ; X ~ X + 1 SRFT 93
|
|
END ; SRFT 94
|
|
RHO[M,1] ~ RHOTEMP ~ KM ; SUM ~ 0 ; SRFT 95
|
|
FOR J ~ 1 STEP 1 UNTIL KM DO SRFT 96
|
|
SUM ~ DATA[J] + SUM ; SRFT 97
|
|
ALEPH[1] ~ ALPHATEMP ~ -SUM / RHOTEMP ; SRFT 98
|
|
FOR J ~ 1 STEP 1 UNTIL KM DO SRFT 99
|
|
BEGIN SRFT 100
|
|
POLY[X] ~ DATA[J] + ALPHATEMP ; X ~ X + 1 SRFT 101
|
|
END ; SRFT 102
|
|
MAJORLOOP : RHOTEMP ~ 0 ; SRFT 103
|
|
WRITE (PAPER, OOT, ALPHABETA) ; SRFT 104
|
|
I ~ I + 1 ; T ~ X - KM ; SRFT 105
|
|
FOR J ~ 1 STEP 1 UNTIL KM DO SRFT 106
|
|
BEGIN SRFT 107
|
|
POWERS[J] ~ POWERS[J] | DATA[J] ; SRFT 108
|
|
RHOTEMP ~ POWERS[J] | POLY[T] + RHOTEMP ; SRFT 109
|
|
T ~ T + 1 SRFT 110
|
|
END ; SRFT 111
|
|
COMMENT THE FOLLOWING TESTS FOR ZERO VALUES OF RHO[M,I] WHICH SRFT 112
|
|
WOULD CAUSE ARITHMETIC OVERFLOW LATER ON. THE FIT FOR SRFT 113
|
|
VARIABLE M IS TERMINATED AT THIS POINT AND P[M] IS SRFT 114
|
|
ADJUSTED ACCORDINGLY ; SRFT 115
|
|
IF RHOTEMP = 0 THEN SRFT 116
|
|
BEGIN SRFT 117
|
|
P[M] ~ I - 2 ; GO TO RESET SRFT 118
|
|
END ; SRFT 119
|
|
RHO[M,I] ~ RHOTEMP ; SRFT 120
|
|
IF I = P[M] + 1 THEN GO TO RESET ; SRFT 121
|
|
SUM ~ 0 ; T ~ X - KM ; SRFT 122
|
|
FOR J ~ 1 STEP 1 UNTIL KM DO SRFT 123
|
|
BEGIN SRFT 124
|
|
SUM ~ POLY[T]*2 | DATA[J] + SUM ; SRFT 125
|
|
T ~ T + 1 SRFT 126
|
|
END ; SRFT 127
|
|
ALEPH[I] ~ ALPHATEMP ~ -SUM / RHOTEMP ; SRFT 128
|
|
BETA[I] ~ BETATEMP ~ -RHOTEMP / RHO[M,I-1] ; SRFT 129
|
|
T ~ X ; T1 ~ X - KM ; T2 ~ T1 - KM ; SRFT 130
|
|
FOR J ~ 1 STEP 1 UNTIL KM DO SRFT 131
|
|
BEGIN SRFT 132
|
|
POLY[T] ~ (DATA[J] + ALPHATEMP) | POLY[T1] + SRFT 133
|
|
POLY[T2] | BETATEMP ; SRFT 134
|
|
T ~ T + 1 ; T1 ~ T1 + 1 ; T2 ~ T2 + 1 SRFT 135
|
|
END ; SRFT 136
|
|
X ~ X + KM ; GO TO MAJORLOOP ; SRFT 137
|
|
RESET : WRITE (PAPER, TITLER) ; SRFT 138
|
|
WRITE (PAPER, CF3, RHOVALUES) ; SRFT 139
|
|
WRITE (PAPER, TITLEP) ; SRFT 140
|
|
FOR T2 ~ 1 STEP 1 UNTIL P[M]+1 DO SRFT 141
|
|
BEGIN SRFT 142
|
|
WRITE (PAPER, CF3, POLYS) ; SRFT 143
|
|
T3 ~ T3 + KM SRFT 144
|
|
END ; SRFT 145
|
|
WRITE (PUNCH, PABF, PABO) ; SRFT 146
|
|
COEFFICIENTC : SRFT 147
|
|
BEGIN SRFT 148
|
|
INTEGER J2 ; SRFT 149
|
|
LABEL TESTI ; SRFT 150
|
|
LIST ID (NAME[M], P[M]), SRFT 151
|
|
CIJO (FOR T ~ J2+1 STEP 1 UNTIL J2+I DO SRFT 152
|
|
COEFF[T]) ; SRFT 153
|
|
FORMAT OUT IDENT (// X40, "POWER COEFFICIENTS", 2I5 //), SRFT 154
|
|
OUTPOLY (6E20.11) ; SRFT 155
|
|
J ~ PMAX + 1 ; J2 ~ J + J ; SRFT 156
|
|
WRITE (PAPER, IDENT, ID) ; SRFT 157
|
|
FOR D ~ 1 STEP 1 UNTIL J2+J DO SRFT 158
|
|
COEFF[D] ~ 0 ; SRFT 159
|
|
I ~ 1 ; COEFF[J+2] ~ COEFF[J2+1] ~ ALEPH[1] ; SRFT 160
|
|
WRITE (PAPER, OUTPOLY, CIJO) ; SRFT 161
|
|
COEFF[1] ~ COEFF[J+3] ~ 1.0 ; SRFT 162
|
|
TESTI : IF I = P[M] THEN GO TO LAST ; SRFT 163
|
|
I ~ I + 1 ; T ~ J + 1 ; SRFT 164
|
|
FOR D ~ J2+1 STEP 1 UNTIL J2+I DO SRFT 165
|
|
BEGIN SRFT 166
|
|
COEFF[D] ~ ALEPH[I] | COEFF[T+1] + COEFF[T] + SRFT 167
|
|
BETA[I] | COEFF[T-J] ; SRFT 168
|
|
T ~ T + 1 SRFT 169
|
|
END ; SRFT 170
|
|
WRITE (PAPER, OUTPOLY, CIJO) ; SRFT 171
|
|
T ~ J + 1 ; SRFT 172
|
|
FOR D ~ J+2 STEP 1 UNTIL J+I+1 DO SRFT 173
|
|
BEGIN SRFT 174
|
|
COEFF[T-J] ~ COEFF[D] ; SRFT 175
|
|
COEFF[D] ~ COEFF[T+J] ; T ~ T + 1 SRFT 176
|
|
END ; SRFT 177
|
|
COEFF[D] ~ 1.0 ; GO TO TESTI SRFT 178
|
|
END COEFFICIENTC ; SRFT 179
|
|
LAST : IF M = V THEN GO TO FURTHER ; SRFT 180
|
|
M ~ M + 1 ; I ~ 1 ; GO TO VARIABLE ; SRFT 181
|
|
COMMENT THE FOLLOWING SECTION COMPUTES THE COEFFICIENTS OF THE SRFT 182
|
|
ORTHOGONAL POLYNOMIALS ; SRFT 183
|
|
FURTHER : SRFT 184
|
|
BEGIN SRFT 185
|
|
LABEL BEG1, BACK1, HERE1, REPEAT1, CHECK1, TOWRITE ; SRFT 186
|
|
LIST DEPENDENTS (FOR D ~ 1 STEP 1 UNTIL N DO SRFT 187
|
|
PHI[D]), SRFT 188
|
|
OUTCOEFF (FOR J ~ 1 STEP 1 UNTIL T3 DO SRFT 189
|
|
COEFFICIENTS[J]) ; SRFT 190
|
|
FORMAT IN CFD (4E20.11) ; SRFT 191
|
|
FORMAT OUT ACOEFF (// X50, "ORTHOGONAL COEFFICIENTS" //), SRFT 192
|
|
FORMCOEFF (4E20.11) ; SRFT 193
|
|
READ (CARD, CFD, DEPENDENTS) ; SRFT 194
|
|
CLOSE (CARD, RELEASE) ;
|
|
WRITE (PAPER, ACOEFF) ; SRFT 195
|
|
FOR M ~ 1 STEP 1 UNTIL V DO SRFT 196
|
|
BEGIN SRFT 197
|
|
DEGREE[M] ~ TOTAL[M] ; SRFT 198
|
|
TT[M] ~ 0 ; DEG[M] ~ 1 SRFT 199
|
|
END ; SRFT 200
|
|
T ~ 1 ; SRFT 201
|
|
SUM ~ 0 ;
|
|
BEG1 : D ~ 1 ; SRFT 202
|
|
T2 ~ TOTAL[V] + K[V] | DEG[V] - 1 ; SRFT 203
|
|
T1 ~ -K[V] + T2 + 1 ; SRFT 204
|
|
BACK1 : FOR I ~ T1 STEP 1 UNTIL T2 DO SRFT 205
|
|
BEGIN SRFT 206
|
|
SUM ~ SUM + PHI[D] | POLY[I] ; SRFT 207
|
|
D ~ D + 1 SRFT 208
|
|
END ; SRFT 209
|
|
TT[V] ~ SUM ; SUM ~ 0 ;
|
|
IF V = 1 THEN GO TO REPEAT1 ; SRFT 210
|
|
M ~ V - 1 ; SRFT 211
|
|
HERE1 : TT[M] ~ POLY[K[M]|(DEG[M]-1)+DEGREE[M]] | TT[M+1] + TT[M];SRFT 212
|
|
TT[M+1] ~ 0 ; SRFT 213
|
|
IF DEGREE[M] ~ DEGREE[M] + 1 > TOTAL[M] + K[M] - 1 THEN SRFT 214
|
|
BEGIN SRFT 215
|
|
IF M = 1 THEN GO TO REPEAT1 ; SRFT 216
|
|
DEGREE[M] ~ TOTAL[M] ; SRFT 217
|
|
M ~ M - 1 ; GO TO HERE1 SRFT 218
|
|
END ; SRFT 219
|
|
GO TO BACK1 ; SRFT 220
|
|
REPEAT1 : FOR M ~ 1 STEP 1 UNTIL V DO SRFT 221
|
|
TT[1] ~ TT[1] / RHO[M,DEG[M]] ; SRFT 222
|
|
COEFFICIENTS[T] ~ TT[1] ; T ~ T + 1 ; SRFT 223
|
|
DEGREE[1] ~ TOTAL[1] ; M ~ V ; SRFT 224
|
|
CHECK1 : IF DEG[M] ! P[M] + 1 THEN SRFT 225
|
|
BEGIN SRFT 226
|
|
DEG[M] ~ DEG[M] + 1 ; TT[1] ~ 0 ; SRFT 227
|
|
GO TO BEG1 SRFT 228
|
|
END ; SRFT 229
|
|
IF M = 1 THEN GO TO TOWRITE ; SRFT 230
|
|
DEG[M] ~ 1 ; M ~ M - 1 ; SRFT 231
|
|
GO TO CHECK1 ; SRFT 232
|
|
TOWRITE : T3 ~ T - 1 ; SRFT 233
|
|
WRITE (PAPER, FORMCOEFF, OUTCOEFF) ; SRFT 234
|
|
WRITE (PUNCH, FORMCOEFF, OUTCOEFF) ; SRFT 235
|
|
CLOSE (PUNCH, RELEASE) ;
|
|
END FURTHER ; SRFT 236
|
|
COMMENT WRITE MEMORY ON MAGNETIC TAPE FOR LATER TRUNCATION ; SRFT 237
|
|
BEGIN SRFT 238
|
|
FILE OUT TAPE (2, 1023) ; SRFT 239
|
|
WRITE (TAPE, *, MAG) ; SRFT 241
|
|
WRITE (TAPE, *, MAGPO) ;
|
|
WRITE (TAPE, *, MAGC) ;
|
|
WRITE (TAPE, *, MAGPH) ;
|
|
GO TO ONE SRFT 242
|
|
END ; SRFT 243
|
|
ONE : SRFT 244
|
|
BEGIN SRFT 245
|
|
LABEL BEG2, BACK2, HERE2, REPEAT2, CHECK2 ; SRFT 246
|
|
LIST DIFFERENCES (PHI[T], TT[1], DIFF), SRFT 247
|
|
FINALSUM (SUMS, FOR M ~ 1 STEP 1 UNTIL V DO SRFT 248
|
|
P[M]) ; SRFT 249
|
|
FORMAT OUT DF (X35, 3E20.11 /), SRFT 250
|
|
FF (X40, "VARIANCE", E20.11, X5, 10I3 //), SRFT 251
|
|
ERRORTITLE (// X42, "ORIGINAL", X11, "CALCULATED", X10, SRFT 252
|
|
"DIFFERENCE" //) ; SRFT 253
|
|
WRITE (PAPER[PAGE], DF) ;
|
|
WRITE (PAPER, ERRORTITLE) ; SRFT 254
|
|
FOR M ~ 1 STEP 1 UNTIL V DO SRFT 255
|
|
BEGIN SRFT 256
|
|
TALLY[M] ~ K[M] | P[M] ; SRFT 257
|
|
DEGREE[M] ~ TOTAL[M] ; SRFT 258
|
|
DEG[M] ~ 1 ; TT[M] ~ 0 SRFT 259
|
|
END ; SRFT 260
|
|
SUMS ~ 0 ; T ~ 1 ; SRFT 261
|
|
SUM ~ 0 ;
|
|
BEG2 : D ~ 1 ; SRFT 262
|
|
T1 ~ TOTAL[V] + DEG[V] - 1 ; SRFT 263
|
|
T2 ~ TALLY[V] + T1 ; SRFT 264
|
|
BACK2 : FOR I ~ T1 STEP K[V] UNTIL T2 DO SRFT 265
|
|
BEGIN SRFT 266
|
|
SUM ~ SUM + COEFFICIENTS[D] | POLY[I] ; SRFT 267
|
|
D ~ D + 1 SRFT 268
|
|
END ; SRFT 269
|
|
TT[V] ~ SUM ; SUM ~ 0 ;
|
|
IF V = 1 THEN GO TO REPEAT2 ; SRFT 270
|
|
M ~ V - 1 ; SRFT 271
|
|
HERE2 : TT[M] ~ POLY[DEGREE[M]+DEG[M]-1] | TT[M+1] + TT[M] ; SRFT 272
|
|
TT[M+1] ~ 0 ; SRFT 273
|
|
IF DEGREE[M] ~ DEGREE[M] + K[M] > TOTAL[M] + TALLY[M] THENSRFT 274
|
|
BEGIN SRFT 275
|
|
IF M = 1 THEN GO TO REPEAT2 ; SRFT 276
|
|
DEGREE[M] ~ TOTAL[M] ; SRFT 277
|
|
M ~ M - 1 ; GO TO HERE2 SRFT 278
|
|
END ; SRFT 279
|
|
GO TO BACK2 ; SRFT 280
|
|
REPEAT2 : DIFF ~ PHI[T] - TT[1] ; SRFT 281
|
|
WRITE (PAPER, DF, DIFFERENCES) ; SRFT 282
|
|
SUMS ~ DIFF | DIFF + SUMS ; SRFT 283
|
|
DEGREE[1] ~ TOTAL[1] ; SRFT 284
|
|
T ~ T + 1 ; M ~ V ; SRFT 285
|
|
CHECK2 : IF DEG[M] < K[M] THEN SRFT 286
|
|
BEGIN SRFT 287
|
|
DEG[M] ~ DEG[M] + 1 ; TT[1] ~ 0 ; SRFT 288
|
|
GO TO BEG2 SRFT 289
|
|
END ; SRFT 290
|
|
IF M = 1 THEN SRFT 291
|
|
BEGIN SRFT 292
|
|
SUMS ~ SUMS / N ; SRFT 293
|
|
WRITE (PAPER, FF, FINALSUM) ; SRFT 294
|
|
IF ST THEN GO TO TRUNCATE ; SRFT 295
|
|
GO TO THEEND SRFT 296
|
|
END ; SRFT 297
|
|
DEG[M] ~ 1 ; M ~ M - 1 ; SRFT 298
|
|
GO TO CHECK2 SRFT 299
|
|
END ONE ; SRFT 300
|
|
COMMENT REPLACE MEMORY FROM PREVIOUSLY WRITTEN MAGNETIC TAPE ; SRFT 301
|
|
REPLACE : SRFT 302
|
|
BEGIN SRFT 303
|
|
FILE IN TAPE (2, 1023) ; SRFT 304
|
|
READ (TAPE, *, MAG) ; SRFT 306
|
|
READ (TAPE, *, MAGPO) ;
|
|
READ (TAPE, *, MAGC) ;
|
|
READ (TAPE, *, MAGPH) ;
|
|
END REPLACE ; SRFT 307
|
|
COMMENT THE FOLLOWING BLOCK TRUNCATES THE ORTHOGONAL COEFFICIENTS SRFT 308
|
|
BY OVERLAYING AND THEN THE ERROR IS RECOMPUTED. THIS IS SRFT 309
|
|
DONE ACCORDING TO NEW VALUES OF P[M] THAT ARE INPUT. ; SRFT 310
|
|
TRUNCATE : SRFT 311
|
|
BEGIN SRFT 312
|
|
LABEL TRUNC, DEGE, REPEAT ; SRFT 313
|
|
LIST NEW (FOR M ~ 1 STEP 1 UNTIL V DO SRFT 314
|
|
NEWP[M]) ; SRFT 315
|
|
READ (CARD, CFI, NEW) [THEEND] ; SRFT 316
|
|
FOR M ~ 1 STEP 1 UNTIL V DO SRFT 317
|
|
DEGREE[M] ~ 0 ; SRFT 318
|
|
T ~ NEWP[V] + 1 ; I ~ 1 ; SRFT 319
|
|
DEG[V] ~ P[V] + 1 ; SRFT 320
|
|
IF V = 1 THEN GO TO REPEAT ; SRFT 321
|
|
FOR M ~ V-1 STEP -1 UNTIL 1 DO SRFT 322
|
|
DEG[M] ~ DEG[M+1] | (P[M] + 1) ; SRFT 323
|
|
COMMENT COMPUTE THE NUMBER OF UNWANTED COEFFICIENTS TO BE JUMPED SRFT 324
|
|
OVER ; SRFT 325
|
|
TRUNC : J ~ 0 ; SRFT 326
|
|
FOR M ~ 1 STEP 1 UNTIL V-1 DO SRFT 327
|
|
J ~ DEG[M+1] | DEGREE[M] + J ; SRFT 328
|
|
FOR D ~ 1 STEP 1 UNTIL T DO SRFT 329
|
|
BEGIN SRFT 330
|
|
T1 ~ J + D ; SRFT 331
|
|
COEFFICIENTS[I] ~ COEFFICIENTS[T1] ; SRFT 332
|
|
I ~ I + 1 SRFT 333
|
|
END ; SRFT 334
|
|
M ~ V - 1 ; SRFT 335
|
|
DEGE : DEGREE[M] ~ DEGREE[M] + 1 ; SRFT 336
|
|
IF DEGREE[M] ! NEWP[M] + 1 THEN GO TO TRUNC ; SRFT 337
|
|
IF M = 1 THEN SRFT 338
|
|
BEGIN SRFT 339
|
|
REPEAT : FOR M ~ 1 STEP 1 UNTIL V DO SRFT 340
|
|
P[M] ~ NEWP[M] ; SRFT 341
|
|
GO TO ONE SRFT 342
|
|
END ; SRFT 343
|
|
DEGREE[M] ~ 0 ; M ~ M - 1 ; GO TO DEGE SRFT 344
|
|
END TRUNCATE ; SRFT 345
|
|
THEEND : END SRFT 346
|
|
END SRFT 347
|
|
END ; SRFT 349
|
|
END .
|
|
9 6 3 567 0 0 INITIAL DATA
|
|
1111 9 6 VARIABLE 1
|
|
2222 7 6 VARIABLE 2
|
|
3333 9 6 VARIABLE 3
|
|
0.1 , 3 0.2 , 3 0.3 , 3 0.4 , 3
|
|
0.5 , 3 0.75 , 3 1.0 , 3 1.6 , 3
|
|
2.0 , 3 VARIABLE 1
|
|
0.16 , 2 0.3 , 2 0.4 , 2 0.5 , 2
|
|
0.6 , 2 0.8 , 2 1.0 , 2 VARIABLE 2
|
|
0.0 , 0.5 , 2 1.0 , 2 1.5 , 2
|
|
2.0 , 2 3.0 , 2 4.0 , 2 5.0 , 2
|
|
6.0 , 2 VARIABLE 3
|
|
.48000000000, 04 .52000000000, 04 .56300000000, 04 .61000000000, 04
|
|
.65600000000, 04 .75500000000, 04 .86100000000, 04 .97600000000, 04
|
|
.10960000000, 05 .84900000000, 04 .91500000000, 04 .98400000000, 04
|
|
.10530000000, 05 .11250000000, 05 .12780000000, 05 .14550000000, 05
|
|
.16470000000, 05 .18570000000, 05 .10800000000, 05 .11640000000, 05
|
|
.12560000000, 05 .13440000000, 05 .14360000000, 05 .16280000000, 05
|
|
.18560000000, 05 .21000000000, 05 .23640000000, 05 .13150000000, 05
|
|
.14150000000, 05 .15200000000, 05 .16250000000, 05 .17400000000, 05
|
|
.19850000000, 05 .22500000000, 05 .25500000000, 05 .28700000000, 05
|
|
.15600000000, 05 .16740000000, 05 .18000000000, 05 .19200000000, 05
|
|
.20520000000, 05 .23400000000, 05 .26520000000, 05 .29940000000, 05
|
|
.33660000000, 05 .20000000000, 05 .21600000000, 05 .23200000000, 05
|
|
.24880000000, 05 .26620000000, 05 .30080000000, 05 .34160000000, 05
|
|
.38480000000, 05 .43760000000, 05 .24500000000, 05 .26500000000, 05
|
|
.28500000000, 05 .30700000000, 05 .32900000000, 05 .37500000000, 05
|
|
.42700000000, 05 .48000000000, 05 .53900000000, 05 .47200000000, 04
|
|
.51500000000, 04 .56000000000, 04 .60600000000, 04 .65300000000, 04
|
|
.75200000000, 04 .85900000000, 04 .96800000000, 04 .10930000000, 05
|
|
.82900000000, 04 .89700000000, 04 .96600000000, 04 .10410000000, 05
|
|
.11160000000, 05 .12780000000, 05 .14550000000, 05 .16410000000, 05
|
|
.18450000000, 05 .10560000000, 05 .11400000000, 05 .12360000000, 05
|
|
.13240000000, 05 .14200000000, 05 .16240000000, 05 .18480000000, 05
|
|
.20840000000, 05 .23480000000, 05 .12900000000, 05 .13850000000, 05
|
|
.14900000000, 05 .15950000000, 05 .17050000000, 05 .19600000000, 05
|
|
.22250000000, 05 .25250000000, 05 .28500000000, 05 .15180000000, 05
|
|
.16280000000, 05 .17460000000, 05 .18780000000, 05 .20100000000, 05
|
|
.23040000000, 05 .26160000000, 05 .29700000000, 05 .33600000000, 05
|
|
.19760000000, 05 .21200000000, 05 .22560000000, 05 .24160000000, 05
|
|
.25920000000, 05 .29820000000, 05 .33920000000, 05 .38480000000, 05
|
|
.43520000000, 05 .24000000000, 05 .25700000000, 05 .27610000000, 05
|
|
.29500000000, 05 .31600000000, 05 .36100000000, 05 .41360000000, 05
|
|
.47300000000, 05 .53400000000, 05 .46400000000, 04 .50700000000, 04
|
|
.55200000000, 04 .59800000000, 04 .64600000000, 04 .74600000000, 04
|
|
.85000000000, 04 .96200000000, 04 .10880000000, 05 .81000000000, 04
|
|
.88200000000, 04 .95200000000, 04 .10320000000, 05 .11100000000, 05
|
|
.12720000000, 05 .14430000000, 05 .16350000000, 05 .18450000000, 05
|
|
.10250000000, 05 .11200000000, 05 .12120000000, 05 .13000000000, 05
|
|
.14000000000, 05 .16080000000, 05 .18400000000, 05 .20800000000, 05
|
|
.23400000000, 05 .12380000000, 05 .13500000000, 05 .14500000000, 05
|
|
.15550000000, 05 .16750000000, 05 .19300000000, 05 .22200000000, 05
|
|
.25100000000, 05 .28500000000, 05 .14420000000, 05 .15520000000, 05
|
|
.16800000000, 05 .18000000000, 05 .19380000000, 05 .22380000000, 05
|
|
.25800000000, 05 .29400000000, 05 .33420000000, 05 .18400000000, 05
|
|
.20000000000, 05 .21390000000, 05 .23050000000, 05 .24930000000, 05
|
|
.29050000000, 05 .33200000000, 05 .38000000000, 05 .43120000000, 05
|
|
.22260000000, 05 .24200000000, 05 .26110000000, 05 .28100000000, 05
|
|
.30220000000, 05 .35000000000, 05 .40360000000, 05 .46380000000, 05
|
|
.52590000000, 05 .45900000000, 04 .50400000000, 04 .54900000000, 04
|
|
.59500000000, 04 .64200000000, 04 .74200000000, 04 .84800000000, 04
|
|
.96200000000, 04 .10860000000, 05 .77400000000, 04 .85200000000, 04
|
|
.93000000000, 04 .10110000000, 05 .10950000000, 05 .12630000000, 05
|
|
.14400000000, 05 .16350000000, 05 .18420000000, 05 .98400000000, 04
|
|
.10760000000, 05 .11720000000, 05 .12680000000, 05 .13720000000, 05
|
|
.15840000000, 05 .18200000000, 05 .20720000000, 05 .23400000000, 05
|
|
.11750000000, 05 .12900000000, 05 .13950000000, 05 .15050000000, 05
|
|
.16300000000, 05 .19000000000, 05 .21950000000, 05 .25000000000, 05
|
|
.28250000000, 05 .13640000000, 05 .14800000000, 05 .16020000000, 05
|
|
.17400000000, 05 .18900000000, 05 .21900000000, 05 .25380000000, 05
|
|
.29100000000, 05 .33130000000, 05 .17090000000, 05 .18740000000, 05
|
|
.20360000000, 05 .22140000000, 05 .24000000000, 05 .28230000000, 05
|
|
.32480000000, 05 .37380000000, 05 .42560000000, 05 .20570000000, 05
|
|
.22690000000, 05 .24730000000, 05 .26830000000, 05 .29060000000, 05
|
|
.34010000000, 05 .39500000000, 05 .45390000000, 05 .51800000000, 05
|
|
.44800000000, 04 .49600000000, 04 .54100000000, 04 .59000000000, 04
|
|
.63800000000, 04 .73900000000, 04 .84300000000, 04 .96000000000, 04
|
|
.10850000000, 05 .75900000000, 04 .83700000000, 04 .91800000000, 04
|
|
.10020000000, 05 .10860000000, 05 .12600000000, 05 .14400000000, 05
|
|
.16260000000, 05 .18300000000, 05 .96800000000, 04 .10520000000, 05
|
|
.11480000000, 05 .12480000000, 05 .13600000000, 05 .15800000000, 05
|
|
.18160000000, 05 .20640000000, 05 .23360000000, 05 .11500000000, 05
|
|
.12600000000, 05 .13750000000, 05 .14900000000, 05 .16100000000, 05
|
|
.18750000000, 05 .21600000000, 05 .24800000000, 05 .28200000000, 05
|
|
.13140000000, 05 .14400000000, 05 .15660000000, 05 .17100000000, 05
|
|
.18600000000, 05 .21600000000, 05 .25020000000, 05 .28800000000, 05
|
|
.32840000000, 05 .16420000000, 05 .17980000000, 05 .19700000000, 05
|
|
.21600000000, 05 .23480000000, 05 .27520000000, 05 .31890000000, 05
|
|
.36690000000, 05 .41920000000, 05 .19370000000, 05 .21540000000, 05
|
|
.23710000000, 05 .25920000000, 05 .28230000000, 05 .33260000000, 05
|
|
.38800000000, 05 .44500000000, 05 .50800000000, 05 .44000000000, 04
|
|
.48600000000, 04 .53300000000, 04 .57900000000, 04 .62600000000, 04
|
|
.72800000000, 04 .84200000000, 04 .96000000000, 04 .10800000000, 05
|
|
.75000000000, 04 .82500000000, 04 .90600000000, 04 .98700000000, 04
|
|
.10680000000, 05 .12330000000, 05 .14160000000, 05 .16110000000, 05
|
|
.18210000000, 05 .95600000000, 04 .10400000000, 05 .11360000000, 05
|
|
.12320000000, 05 .13360000000, 05 .15440000000, 05 .17800000000, 05
|
|
.20280000000, 05 .23000000000, 05 .11200000000, 05 .12250000000, 05
|
|
.13450000000, 05 .14650000000, 05 .15850000000, 05 .18400000000, 05
|
|
.21200000000, 05 .24250000000, 05 .27700000000, 05 .12600000000, 05
|
|
.14100000000, 05 .15300000000, 05 .16620000000, 05 .18060000000, 05
|
|
.21120000000, 05 .24300000000, 05 .27900000000, 05 .32100000000, 05
|
|
.16100000000, 05 .17360000000, 05 .18980000000, 05 .20780000000, 05
|
|
.22660000000, 05 .26540000000, 05 .30720000000, 05 .35200000000, 05
|
|
.40480000000, 05 .18100000000, 05 .20000000 00, 05 .22200000000, 05
|
|
.24610000000, 05 .27000000000, 05 .31700000000, 05 .37000000000, 05
|
|
.42600000000, 05 .49000000000, 05 .42700000000, 04 .47000000000, 04
|
|
.51700000000, 04 .56600000000, 04 .61400000000, 04 .71800000000, 04
|
|
.83200000000, 04 .95000000000, 04 .10770000000, 05 .72000000000, 04
|
|
.79800000000, 04 .87900000000, 04 .96000000000, 04 .10440000000, 05
|
|
.12150000000, 05 .13980000000, 05 .15960000000, 05 .18120000000, 05
|
|
.90000000000, 04 .10000000000, 05 .10960000000, 05 .11960000000, 05
|
|
.13000000000, 05 .15160000000, 05 .17480000000, 05 .20000000000, 05
|
|
.22800000000, 05 .10600000000, 05 .11750000000, 05 .12900000000, 05
|
|
.14100000000, 05 .15350000000, 05 .18000000000, 05 .20800000000, 05
|
|
.23900000000, 05 .27250000000, 05 .11640000000, 05 .13020000000, 05
|
|
.14400000000, 05 .15840000000, 05 .17280000000, 05 .20580000000, 05
|
|
.24000000000, 05 .27480000000, 05 .31500000000, 05 .14000000000, 05
|
|
.15760000000, 05 .17600000000, 05 .19380000000, 05 .21440000000, 05
|
|
.25600000000, 05 .30000000000, 05 .34720000000, 05 .40000000000, 05
|
|
.16300000000, 05 .18300000000, 05 .20410000000, 05 .22710000000, 05
|
|
.25300000000, 05 .30500000000, 05 .36000000000, 05 .42000000000, 05
|
|
.48200000000, 05 .41100000000, 04 .45600000000, 04 .50200000000, 04
|
|
.55200000000, 04 .60300000000, 04 .71200000000, 04 .82400000000, 04
|
|
.93800000000, 04 .10640000000, 05 .66900000000, 04 .75000000000, 04
|
|
.84000000000, 04 .92400000000, 04 .10140000000, 05 .11940000000, 05
|
|
.13800000000, 05 .15720000000, 05 .17820000000, 05 .83600000000, 04
|
|
.94000000000, 04 .10400000000, 05 .11480000000, 05 .12640000000, 05
|
|
.15000000000, 05 .17360000000, 05 .19800000000, 05 .22400000000, 05
|
|
.10000000000, 05 .11150000000, 05 .12300000000, 05 .13650000000, 05
|
|
.14900000000, 05 .17550000000, 05 .20450000000, 05 .23500000000, 05
|
|
.26850000000, 05 .11160000000, 05 .12240000000, 05 .13800000000, 05
|
|
.15240000000, 05 .16680000000, 05 .19920000000, 05 .23340000000, 05
|
|
.27000000000, 05 .30900000000, 05 .13520000000, 05 .15200000000, 05
|
|
.17030000000, 05 .18800000000, 05 .20800000000, 05 .24800000000, 05
|
|
.29040000000, 05 .33680000000, 05 .38800000000, 05 .15700000000, 05
|
|
.17370000000, 05 .19800000000, 05 .22390000000, 05 .24790000000, 05
|
|
.29280000000, 05 .34370000000, 05 .40400000000, 05 .46600000000, 05
|
|
.40000000000, 04 .44500000000, 04 .49100000000, 04 .54200000000, 04
|
|
.59200000000, 04 .70200000000, 04 .81600000000, 04 .93600000000, 04
|
|
.10620000000, 05 .65100000000, 04 .72000000000, 04 .79800000000, 04
|
|
.88200000000, 04 .96900000000, 04 .11460000000, 05 .13350000000, 05
|
|
.15450000000, 05 .17580000000, 05 .76000000000, 04 .85600000000, 04
|
|
.95600000000, 04 .10640000000, 05 .11760000000, 05 .14080000000, 05
|
|
.16560000000, 05 .19080000000, 05 .21880000000, 05 .85500000000, 04
|
|
.97000000000, 04 .10850000000, 05 .12150000000, 05 .13500000000, 05
|
|
.16250000000, 05 .19250000000, 05 .22500000000, 05 .26050000000, 05
|
|
.93000000000, 04 .10680000000, 05 .12120000000, 05 .13500000000, 05
|
|
.15060000000, 05 .18420000000, 05 .21900000000, 05 .25800000000, 05
|
|
.30180000000, 05 .10800000000, 05 .12800000000, 05 .14640000000, 05
|
|
.16560000000, 05 .18670000000, 05 .22720000000, 05 .27200000000, 05
|
|
.32320000000, 05 .38080000000, 05 .12400000000, 05 .14500000000, 05
|
|
.16860000000, 05 .19190000000, 05 .21600000000, 05 .26800000000, 05
|
|
.32500000000, 05 .38900000000, 05 .46100000000, 05
|