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
920  S0=0
930  FOR I7=I9 TO N0
944  Z3=I7
945  GOSUB 1035
946  C[I7]=Q0
950  E[I7]=B[I7]-C[I7]
960  S0=S0+E[I7]*E[I7]
961  NEXT I7
970  A=S0
980  CHAIN Q$,Q9
1035  REM
1045  S1=0
1055  IF (M1=0) THEN 1105
1065  FOR J1=1 TO M1
1075  K1=Z3-J1
1085  S1=S1+J[J1+1]*(B[K1]-A0)
1095  NEXT J1
1105  IF (M2=0) THEN 1155
1115  FOR J1=1 TO M2
1125  K1=Z3-J1
1135  S1=S1-I[J1+1]*E[K1]
1145  NEXT J1
1155  Q0=S1+A0
1165  RETURN 
1305  END 
