Listing of file='3A' on disk='vmedia/glbr_22.wvd.zip'
# Sector 618, program filename = '3A' 0010 COM A(13),Q(7,8),E(8),M,N 0015 PRINT "INPUT M,N" : INPUT M,N : PRINT "INPUT DATA POINTS" 0017 FOR I=2 TO 2*M+1 : A(I)=0 : NEXT I 0018 FOR I=1 TO M+2 : E(I)=0 : NEXT I 0020 A(1)=N 0025 FOR I=1 TO N : INPUT X,Y 0030 FOR J=2 TO 2*M+1 : A(J)=A(J)+X^(J-1) : NEXT J 0035 FOR J=1 TO M+1 : E(J),Q(J,M+2)=E(J)+Y*X^(J-1) : NEXT J 0040 E(M+2)=E(M+2)+Y^2 : NEXT I 0045 FOR I=1 TO M+1 : FOR J=1 TO M+1 : Q(I,J)=A(I+J-1) : NEXT J : NEXT I 0050 FOR S=1 TO M+1 0055 FOR T=S TO M+1 : IF Q(T,S)<>0 THEN 65 : NEXT T 0060 PRINT "NO UNIQUE SOLUTION" : STOP 0065 GOSUB 90 0070 C=1/Q(S,S) : GOSUB 110 0075 FOR T=1 TO M+1 : IF T=S THEN 85 0080 C=-Q(T,S) : GOSUB 120 0085 NEXT T : NEXT S : GOTO 130 0090 FOR J=1 TO M+2 0095 B=Q(S,J) : Q(S,J)=Q(T,J) : Q(T,J)=B 0100 NEXT J : RETURN 0110 FOR J=1 TO M+2 : Q(S,J)=C*Q(S,J) : NEXT J : RETURN 0120 FOR J=1 TO M+2 : Q(T,J)=Q(T,J)+C*Q(S,J) : NEXT J : RETURN 0130 PRINT 0140 FOR I=1 TO M+1 : PRINT I-1;"DEG.COEFF.=";Q(I,M+2) : NEXT I 0145 STOP : PRINT HEX(03) 0200 S=0 0210 FOR I=2 TO M+1 : S=S+Q(I,M+2)*(E(I)-A(I)*E(1)/N) : NEXT I 0220 T=E(M+2)-E(1)^2/N : C=T-S 0230 I=N-M-1 : J=S/M : K=C/I 0240 PRINT : PRINT : PRINT " REGRESSION TABLE" : PRINT 0245 PRINT "SOURCE","SUM OF SQ.","DEG. FREEDOM","MEAN SQ." 0250 PRINT "REGRESSION",S,M,J 0260 PRINT "RESIDUAL",C,I,K 0270 PRINT "TOTAL",T,N-1 : PRINT 0280 PRINT "F=";J/K : PRINT : PRINT 0290 J=S/T : PRINT "COEFF. OF DETERMINATION=";J 0300 PRINT "COEFF. OF CORRELATION=";SQR(J) 0310 PRINT "STANDARD ERROR OF ESTIMATE=";SQR(C/I) : PRINT : PRINT 0320 PRINT "DO YOU WISH TO ESTIMATE VALUES OF Y FROM" 0330 PRINT "THE REGRESSION CURVE? (1=YES,0=NO)" 0340 INPUT I : IF I=0 THEN 999 0350 PRINT : S=Q(1,M+2) 0360 PRINT "INPUT X" : INPUT T 0370 FOR I=1 TO M : S=S+Q(I+1,M+2)*T^I : NEXT I 0380 PRINT "Y=";S : PRINT 0390 PRINT "ANOTHER POINT? (1=YES,0=NO)" : GOTO 340 0999 END