1 REM:LASTREC@2-28-88@151...TAPCNT=4000-4412...STAT=7973 10 CLS:MEMSET:CLEAR:WIDTH20,60:DEFINTE-L,Z:OPTIONBASE1 20 DIMDEV(24),FRQ(44),NMBR(44),DRWG(44),CMBN(44),NO(6,6),MTRX#(6,6) 30 DIMSUM(4,44),SLOPE(44),YINT(44),Y(44) 40 A$=TIME$:PI#=3.141592653589793:E#=2.718281828459045 50 READFRQ:FORI=1TOFRQ-24:FORJ=1TO6:READK:FRQ(K)=FRQ(K)+1:NEXTJ,I 60 MEAN=((FRQ-24)*6/44):FORI=1TO44:VAR=VAR+((FRQ(I)-MEAN))^2/44:NEXTI 70 DEV=SQR(VAR):SIGMA=VAR*1936+((FRQ-24)*6)^2:FORI=FRQ-23TOFRQ:DELTA=0 80 FORJ=1TO6:READK:FRQ(K)=FRQ(K)+1:DELTA=DELTA+FRQ(K):NEXTJ 90 DELTA=DELTA*88-264:SIGMA=SIGMA+DELTA:MICRO=(I*6)^2 100 VAR=(SIGMA-MICRO)/1936:DEV(I+24-FRQ)=SQR(VAR):NEXTI 110 MICRO=((FRQ+1)*6)^2:MEAN=(FRQ+1)*6/44 120 GOSUB740:DEV1=DEV:VAR1=VAR 130 DELTA=1936*VAR+MICRO-SIGMA:DELTA=(DELTA+264)/88 140 EFRQ1=CINT(DELTA):ESUMSQ1=EFRQ1*2-6 150 GOSUB810:DEV2=DEV:VAR2=VAR 160 DELTA=1936*VAR+MICRO-SIGMA:DELTA=(DELTA+264)/88 170 EFRQ2=CINT(DELTA):ESUMSQ2=EFRQ2*2-6 180 PRINT"LINEAR:":PRINT"MEAN=";MEAN:PRINT"VAR=";VAR1:PRINT"DEV=";DEV1 190 PRINT"SLOPE=";SLOPE:PRINT"YINT=";YINT:PRINT"CORR=";CORR:PRINT"SUM=";EFRQ1 200 PRINT"POLYNOMIAL:":PRINT"DEV=";DEV2:PRINT"VAR=";VAR2:PRINT"SUM=";EFRQ2 210 FORI=1TO6:PRINT"COEF";STR$(I-1);"=";STR$(CSNG(COEF#(I))):NEXTI 220 PRINT"MEAN SUM=";CINT(MEAN*6) 230 RESTORE:READFRQ:FORI=1TO44:FRQ(I)=0:NEXTI 240 FORI=1TOFRQ-24:FORJ=1TO6:READK:FRQ(K)=FRQ(K)+1:NEXTJ,I 250 FORI=FRQ-23TOFRQ:FORJ=1TO6:READK:FRQ(K)=FRQ(K)+1:E(J)=FRQ(K):NEXTJ 260 SUM=0:FORJ=1TO6:SUM=SUM+E(J):NEXTJ:AVE=SUM/6 270 SUM=0:FORJ=1TO6:SUM=SUM+(E(J)-AVE)^2:NEXTJ 280 DEV(I+24-FRQ)=SQR(SUM/6):NEXTI 290 GOSUB740:DEV3=DEV:VAR3=VAR:GOSUB810:DEV4=DEV:VAR4=VAR 300 RESTORE:READFRQ:FORI=1TO44:FRQ(I)=0:NEXTI 310 FORI=1TOFRQ:FORJ=1TO6:READK:FRQ(K)=FRQ(K)+1 320 SUM(1,K)=SUM(1,K)+FRQ(K):SUM(2,K)=SUM(2,K)+I 330 SUM(3,K)=SUM(3,K)+FRQ(K)*I:SUM(4,K)=SUM(4,K)+FRQ(K)^2:NEXTJ,I 340 FORI=1TO44:NUM=(SUM(1,I)*SUM(2,I)/FRQ(I)-SUM(3,I)) 350 DEN=(SUM(1,I)^2/FRQ(I)-SUM(4,I)):SLOPE(I)=NUM/DEN 360 YINT(I)=(SUM(2,I)-SUM(1,I)*SLOPE(I))/FRQ(I) 370 Y(I)=(FRQ(I)+1)*SLOPE(I)+YINT(I):NEXTI 380 FORI=1TO44:PRINTTAB(10*((I-1)MOD2));:PRINTUSING"##";I; 390 PRINTUSING"###";FRQ(I);:PRINTUSING"####";Y(I);:NEXTI 400 PRINT:PRINT:FORI=1TO44:NMBR(I)=I:NEXTI 410 FORL=1TO2:FORI=1TO44:DRWG(I)=ABS(Y(NMBR(I))-FRQ-1):NEXTI:GOSUB940 420 FORJ=1TO6:NO(2,J)=(NMBR(J)):NEXTJ 430 FORI=1TO44:DRWG(I)=-FRQ(NMBR(I)):NEXTI:GOSUB940 440 FORJ=1TO6:NO(1,J)=(NMBR(J)):NEXTJ,L 450 MIN%=ABS(DRWG(44))+1:MAX%=ABS(DRWG(1))+1 460 FORI=1TO44:CMBN(I)=NMBR(I):NEXTI 470 FORZ=0TO1:IFZ THENEFRQ=EFRQ2:DEV=DEV4 ELSEEFRQ=EFRQ1:DEV=DEV3 480 FORI=1TO44:NMBR(I)=CMBN(I):NEXTI 490 AVE=(FRQ+1)*6/44:XINC=DEV/SQR(35/12):XINC0=(EFRQ-6*AVE-21*XINC)/6 500 FORI=1TO6:XINC(I)=AVE+XINC0+I*XINC:NEXTI 510 FORI=1TO6:FORJ=1TO44:DRWG(J)=ABS(FRQ(NMBR(J))-XINC(I)+1):NEXTJ 520 FORJ=43TOI STEP-1:IFDRWG(J+1)=>DRWG(J)THEN540 530 SWAPDRWG(J),DRWG(J+1):SWAPNMBR(J),NMBR(J+1) 540 NEXTJ,I:FORI=1TO6:NO(3+Z,I)=NMBR(I):NEXTI 550 DEV0=FRQ:AVE=FRQ:FORJ=1TO5:ADD%(J)=MIN%:NEXTJ 560 ADD%(6)=EFRQ:FORJ=1TO5:ADD%(6)=ADD%(6)-ADD%(J):NEXTJ 570 IFADD%(6)<=MAX% THEN670ELSEADD%=ADD%(6)-MAX%:ADD%(6)=MAX% 580 FORJ=5TO1STEP-1:ADD%(J)=ADD%(J)+ADD%:ADD%=ADD%(J)-MAX% 590 ADD%=(ABS(ADD%)+ADD%)\2:ADD%(J)=ADD%(J)-ADD%:NEXTJ 600 I=5:J=1:K=6:SUM=0 610 ON(SGN(ADD%(K)-FRQ(CMBN(J))-1)+2)GOTO620,630,660 620 J=J+1:IFJ>44THEN660ELSE610 630 ADD0%(K)=CMBN(J):SUM=SUM+(Y(CMBN(J))-AVE)^2:K=K-1:IFK>0THEN620 640 DEV=SQR(SUM/6):IFDEV>DEV0 THEN660 650 DEV0=DEV:FORJ=1TO6:NO(5+Z,J)=ADD0%(J):NEXTJ 660 ADD%(5)=ADD%(5)+1:ADD%(6)=ADD%(6)-1 670 IFADD%(5)<=ADD%(6)THEN600ELSEI=I-1:IFI=0THEN690 680 ADD%(I)=ADD%(I)+1:FORJ=I TO5:ADD%(J)=ADD%(I):NEXTJ:GOTO560 690 NEXTZ:RANDOMIZE(FRQ):GOSUB970 700 FORI=44TO39STEP-1:PRINTUSING"###";NMBR(I);:NEXTI:PRINT 710 FORI=1TO6:FORJ=1TO6:PRINTUSING"###";NO(I,J);:NEXTJ:PRINT:NEXTI 720 PRINTA$;TAB(10);TIME$:PRINTDATE$:COPY 730 END 740 SUM1=0:SUM2=0:SUM3=0:SUM4=0:FORI=FRQ-23TOFRQ:J=I+24-FRQ 750 SUM1=SUM1+I:SUM2=SUM2+DEV(J):SUM3=SUM3+I*DEV(J):SUM4=SUM4+I^2:NEXTI 760 SLOPE=(SUM1*SUM2/24-SUM3)/(SUM1^2/24-SUM4):YINT=(SUM2-SUM1*SLOPE)/24 770 DEV=(FRQ+1)*SLOPE+YINT:VAR=DEV^2:Y0=1/DEV/SQR(2*PI#) 780 YAVE=0:FORI=1TO24:YAVE=YAVE+DEV(I):NEXTI:YAVE=YAVE/24 790 YVAR=0:FORI=1TO24:YVAR=YVAR+(DEV(I)-YAVE)^2:NEXTI:YVAR=YVAR/24 800 YDEV=SQR(YVAR):XDEV=SQR(575/12):CORR=SLOPE*XDEV/YDEV:RETURN 810 FORI=1TO6:COEF#(I)=0:FORJ=1TO6:MTRX#(I,J)=0:NEXTJ,I 820 FORI=1TO6:FORJ=1TO6:FORK=1TO24 830 MTRX#(I,J)=MTRX#(I,J)+CDBL(K^(I+J-2)):NEXTK,J,I 840 FORI=1TO6:FORJ=1TO24:COEF#(I)=COEF#(I)+CDBL(DEV(J)*J^(I-1)):NEXTJ,I 850 FORI=1TO6:FORJ=I TO6:MTRX#=MTRX#(J,I):FORK=I TO6 860 MTRX#(J,K)=MTRX#(J,K)/MTRX#:NEXTK:COEF#(J)=COEF#(J)/MTRX#:NEXTJ 870 FORJ=I+1TO6:FORK=I TO6:MTRX#(J,K)=MTRX#(J,K)-MTRX#(I,K):NEXTK 880 COEF#(J)=COEF#(J)-COEF#(I):NEXTJ,I 890 FORI=2TO6:FORJ=I TO6:MTRX#=MTRX#(I-1,J):FORK=J TO6 900 MTRX#(I-1,K)=MTRX#(I-1,K)-MTRX#*MTRX#(J,K):NEXTK 910 COEF#(I-1)=COEF#(I-1)-MTRX#*COEF#(J):NEXTJ,I 920 COEF#=0:FORI=1TO6:COEF#=COEF#+COEF#(I)*25^(I-1):NEXTI 930 DEV=CSNG(COEF#):VAR=DEV^2:RETURN 940 FORJ=43TO1STEP-1:FORK=1TOJ:IFDRWG(K+1)=>DRWG(K)THEN960 950 SWAPDRWG(K),DRWG(K+1):SWAPNMBR(K),NMBR(K+1) 960 NEXTK,J:RETURN 970 FORI=1TO44:DRWG(I)=I:NEXTI 980 FORI=44TO1STEP-1:K=INT(RND*I)+1:NMBR(I)=DRWG(K) 990 FORJ=K TOI-1:DRWG(J)=DRWG(J+1):NEXTJ:NEXTI:RETURN 1000 DATA 184 1001 DATA 5,12,18,22,26,44 1002 DATA 3,16,32,34,40,44 1003 DATA 17,21,22,34,39,40 1004 DATA 9,26,28,31,38,42 1005 DATA 7,12,16,20,37,41 1006 DATA 23,25,37,39,41,42 1007 DATA 4,5,18,19,25,43 1008 DATA 6,21,29,31,32,44 1009 DATA 8,23,25,31,34,44 1010 DATA 2,5,23,27,33,43 1011 DATA 19,20,23,31,37,38 1012 DATA 13,24,28,30,32,34 1013 DATA 5,9,18,29,31,39 1014 DATA 13,14,18,27,42,44 1015 DATA 3,8,13,25,34,41 1016 DATA 24,25,27,28,32,42 1017 DATA 9,14,16,19,26,37 1018 DATA 10,13,15,21,29,34 1019 DATA 8,13,16,23,31,40 1020 DATA 3,14,15,18,19,27 1021 DATA 11,14,16,22,23,33 1022 DATA 12,15,19,37,38,40 1023 DATA 1,4,21,23,26,38 1024 DATA 5,14,21,32,40,43 1025 DATA 13,30,32,33,38,40 1026 DATA 7,9,13,19,20,42 1027 DATA 3,5,9,17,35,41 1028 DATA 1,16,24,25,27,38 1029 DATA 2,4,7,16,19,27 1030 DATA 2,10,21,27,31,36 1031 DATA 6,16,17,28,33,43 1032 DATA 1,20,22,24,26,37 1033 DATA 4,6,31,33,36,40 1034 DATA 7,9,11,12,32,38 1035 DATA 2,21,28,35,40,44 1036 DATA 12,20,22,34,35,44 1037 DATA 1,9,16,32,33,36 1038 DATA 2,21,24,35,37,39 1039 DATA 21,31,32,39,41,42 1040 DATA 6,15,19,24,28,36 1041 DATA 5,15,17,27,33,35 1042 DATA 4,7,8,9,12,19 1043 DATA 1,4,10,29,40,42 1044 DATA 3,4,15,27,28,44 1045 DATA 2,12,23,25,34,38 1046 DATA 6,16,19,28,29,40 1047 DATA 3,4,10,25,35,43 1048 DATA 10,16,19,20,29,43 1049 DATA 3,4,6,10,13,19 1050 DATA 6,8,17,25,26,28 1051 DATA 3,5,19,36,37,40 1052 DATA 6,7,8,9,28,30 1053 DATA 9,10,17,20,24,32 1054 DATA 21,23,32,33,35,42 1055 DATA 1,5,19,23,25,28 1056 DATA 7,14,24,31,40,41 1057 DATA 6,10,22,24,31,41 1058 DATA 1,7,10,11,23,33 1059 DATA 7,15,29,32,33,44 1060 DATA 7,9,18,34,41,43 1061 DATA 2,3,14,18,22,31 1062 DATA 1,9,22,24,41,42 1063 DATA 9,17,22,24,31,34 1064 DATA 4,7,18,22,24,38 1065 DATA 1,4,7,8,12,36 1066 DATA 4,15,21,24,35,37 1067 DATA 4,10,19,25,26,43 1068 DATA 19,21,29,30,31,36 1069 DATA 6,10,17,20,25,43 1070 DATA 5,14,18,21,25,38 1071 DATA 3,4,13,21,29,35 1072 DATA 6,18,22,27,33,36 1073 DATA 7,13,26,36,40,42 1074 DATA 15,17,23,29,39,40 1075 DATA 10,12,23,26,27,28 1076 DATA 8,10,13,15,19,30 1077 DATA 6,22,23,28,37,39 1078 DATA 10,14,28,35,42,44 1079 DATA 7,24,30,35,37,44 1080 DATA 2,5,9,10,33,41 1081 DATA 1,7,8,19,23,42 1082 DATA 13,18,19,22,27,28 1083 DATA 1,5,10,22,24,33 1084 DATA 7,10,17,22,38,44 1085 DATA 9,10,22,28,31,33 1086 DATA 1,10,14,28,34,43 1087 DATA 19,28,33,34,35,37 1088 DATA 5,7,25,31,35,36 1089 DATA 21,27,33,35,39,42 1090 DATA 1,5,10,11,35,44 1091 DATA 2,3,12,22,26,41 1092 DATA 2,6,9,21,29,36 1093 DATA 3,15,28,31,34,36 1094 DATA 9,14,16,20,31,34 1095 DATA 7,9,22,23,32,41 1096 DATA 3,20,24,27,29,33 1097 DATA 6,17,28,35,39,42 1098 DATA 18,21,30,32,33,43 1099 DATA 5,8,20,21,24,37 1100 DATA 9,10,21,25,32,38 1101 DATA 6,9,13,21,31,37 1102 DATA 11,16,24,34,39,43 1103 DATA 16,20,27,33,34,42 1104 DATA 5,6,12,28,43,44 1105 DATA 7,12,17,21,22,24 1106 DATA 5,14,16,27,29,34 1107 DATA 5,6,12,17,20,33 1108 DATA 12,18,19,25,32,38 1109 DATA 4,7,15,17,18,35 1110 DATA 2,8,20,25,30,43 1111 DATA 2,11,19,28,32,40 1112 DATA 8,12,14,27,28,39 1113 DATA 14,19,27,38,43,44 1114 DATA 1,9,14,16,28,44 1115 DATA 3,7,16,20,21,34 1116 DATA 22,23,24,31,34,37 1117 DATA 5,8,25,27,30,41 1118 DATA 3,4,8,14,39,44 1119 DATA 4,16,19,22,27,37 1120 DATA 9,12,18,19,22,43 1121 DATA 11,31,36,39,41,44 1122 DATA 4,7,13,20,32,39 1123 DATA 3,8,15,26,40,44 1124 DATA 1,5,6,26,29,40 1125 DATA 1,3,13,14,29,32 1126 DATA 10,20,22,35,39,41 1127 DATA 2,17,24,26,34,39 1128 DATA 16,26,29,36,39,41 1129 DATA 7,19,25,29,35,39 1130 DATA 18,19,25,36,38,42 1131 DATA 1,5,6,25,38,43 1132 DATA 9,10,38,39,41,42 1133 DATA 7,10,27,30,34,43 1134 DATA 15,24,31,37,38,41 1135 DATA 2,21,26,34,41,43 1136 DATA 8,14,18,30,42,43 1137 DATA 16,21,29,33,35,43 1138 DATA 1,5,18,25,39,42 1139 DATA 3,10,30,36,40,44 1140 DATA 10,16,27,28,33,35 1141 DATA 23,25,29,30,42,44 1142 DATA 2,7,14,25,42,43 1143 DATA 5,8,14,21,22,38 1144 DATA 4,8,22,35,40,44 1145 DATA 11,12,29,33,36,39 1146 DATA 7,10,22,26,33,42 1147 DATA 28,29,31,34,35,40 1148 DATA 1,7,12,14,21,42 1149 DATA 6,18,19,31,34,42 1150 DATA 6,8,10,36,40,41 1151 DATA 8,13,27,30,33,44 1152 DATA 1,27,29,39,42,44 1153 DATA 1,3,28,31,38,41 1154 DATA 14,15,16,30,41,42 1155 DATA 8,19,25,27,35,44 1156 DATA 6,9,15,24,25,32 1157 DATA 5,6,7,9,11,13 1158 DATA 1,4,14,21,31,35 1159 DATA 1,2,29,37,40,44 1160 DATA 17,18,22,35,38,44 1161 DATA 5,9,11,20,32,38 1162 DATA 4,18,19,24,35,44 1163 DATA 8,18,27,33,36,43 1164 DATA 2,11,17,23,26,29 1165 DATA 5,8,18,19,34,38 1166 DATA 9,10,12,13,35,40 1167 DATA 1,2,21,29,31,39 1168 DATA 3,14,17,20,26,31 1169 DATA 10,14,15,20,23,34 1170 DATA 8,20,24,25,30,41 1171 DATA 9,31,32,33,40,42 1172 DATA 1,2,10,19,30,37 1173 DATA 6,14,18,19,38,42 1174 DATA 8,12,13,15,30,31 1175 DATA 4,7,8,15,42,43 1176 DATA 9,11,15,23,38,44 1177 DATA 5,9,12,13,40,44 1178 DATA 4,19,23,28,42,44 1179 DATA 8,10,15,24,26,36 1180 DATA 25,26,29,35,38,41 1181 DATA 3,11,32,33,35,37 1182 DATA 9,18,29,30,38,43 1183 DATA 8,11,15,17,26,43 1184 DATA 11,13,15,25,33,37:REM 6-22-88