30  COM B[300],C[300],D[300],E[300],F[300],G[300]
40  COM H[25],N0,M0,I0,I1,I2,I3,I4,I5,I6,I8
50  COM I[100],J[100],M1,M2,A0
60  COM K[20],L[20],M[20],N[20]
61  COM R$[7],R9
63  COM A
65  COM G1,G2,G3,G4,G5,G6,G7,G8,G9
66  COM H1,H2,H3,H4,H5,H6,H7,H8
70  COM A$[4],P[4],B$[4],Q[4]
75  COM Q$[7],Q9
80  COM F1,F2,F3,F4
81  COM I9
82  COM R[9],S[9],T[9],U[9]
83  COM V[200,10]
84  COM W[9,9],X[9]
85  COM F5,F6,F7,F8,F9,L9
86  COM L8
87  COM S6
109  DIM O[3]
110  N1=N0+1
120  M3=M0
130  M4=M3+101
630  FOR I7=M4 TO N0
644  Z3=I7
645  GOSUB 1175
646  E[I7]=D[I7]-Q0
649  NEXT I7
650  I9=101
660  Z0=100-M3
670  J0=1
680  Z1=0
690  T0=0
700  FOR I7=1 TO 3
710  O[I7]=0
711  NEXT I7
720  FOR I7=1 TO Z0
730  Z2=N0+I7
744  Z3=Z2
745  GOSUB 1175
746  X0=Q0
750  T0=T0-O[J0]
760  O[J0]=ABS(X0-A0)
770  T0=T0+O[J0]
780  J0=J0+1
790  IF (J0 <= 3) THEN 820
800  J0=1
810  Z1=1
820  IF (Z1=0) THEN 860
830  IF (T0>(.0001)) THEN 860
840  I9=I9+2
850  GOTO 890
860  I9=101-I7
870  D[Z2]=X0
880  B[I9]=X0
885  NEXT I7
889  I9=I9+1
890  CHAIN "$IDA53C"
1175  S1=0
1185  IF (M1=0) THEN 1235
1195  FOR J1=1 TO M1
1205  K1=Z3-J1
1215  S1=S1+J[J1+1]*(D[K1]-A0)
1225  NEXT J1
1235  IF (M2=0) THEN 1285
1245  FOR J1=1 TO M2
1255  K1=Z3-J1
1265  S1=S1-I[J1+1]*E[K1]
1275  NEXT J1
1285  Q0=S1+A0
1295  RETURN 
1305  END 
