10 REM H. DORFMAN, PIB, 8/1/69
20 REM STATISTICAL ANALYSIS OF LABORATORY DATA
30 REM STAT - COPYRIGHT 1971, POLYTECHNIC INSTITUTE OF BROOKLYN
99 INPUT'DO YOU WANT A BRIEF SUMMARY';A$:IFA$='YES'THEN GOSUB 1750
160 DIM A(50),R(50),Y(51),W(50)
170 DIM Q(50)
310 PRINT'WHAT IS THE TOTAL NO. OF STUDENTS AND THE CALC. VALUE';
320 INPUT B,K
330 PRINT
335  PRINT'ENTER THE STUDENT NUMBER, FOLLOWED BY ITS VALUE, SEPARATED BY A'
336 PRINT"COMMA, WHEN A '?' IS ENTERED.'
340 FOR I=1 TO B
350 INPUT R(I),A(I)
360 NEXT I
370 PRINT'DATA LISTED BY STUDENT NUMBER'
380 PRINT'---- ------ -- ------- ------'
390 PRINT
400 PRINT' STUDENT NO.','VALUE','ERROR','PERCENT ERROR'
410 PRINT
420 FOR I=1 TO B
430 Z=A(I)-K
440 A8=A8+ABS(Z)/B
450 P=P+A(I)/B
460 W(I)=(ABS(Z/K))*100
470 A9=A9+W(I)/B
480 PRINT R(I),A(I),Z,W(I)
490 NEXT I
500 PRINT
510 PRINT'THE ARITHMETIC MEAN (AVERAGE) IS   ';P
520 PRINT'THE AVERAGE ERROR (ABSOLUTE) IS    ';A8
530 PRINT'THE AVERAGE PERCENT ERROR IS       ';A9
540 PRINT
560 PRINT' DO YOU DESIRE ADDITIONAL INFORMATION? IF SO, TYPE'
570 PRINT' IN THE NUMBER OF YOUR CHOICE.';
580 INPUT N
590 IF N=2 THEN 1010
600 IF N=3 THEN 1200
620 IF N=6 THEN 1730
630 PRINT
640 PRINT'FOR THE DITRIBUTION, WHAT LOWER LIMIT, UPPER LIMIT,'
650 PRINT'AND STEP DO YOU DESIRE';
660 INPUT E,F,G
670 PRINT
680 PRINT
690 PRINT'EXPERIMENTAL VALUE DISTRIBUTION'
700 PRINT'------------ ----- ------------'
710 PRINT
720 FOR I=1 TO B
730 M=0
740 IF A(I)>=F THEN 850
750 IF A(I)>=E THEN 770
760 Y(51)=Y(51)+1
770 FOR J=E TO (F-G) STEP G
780 M=M+1
790 IF A(I)>=J THEN 810
800 GOTO 860
810 IF A(I)>=(J+G) THEN 830
820 Y(M)=Y(M)+1
830 NEXT J
840 GOTO 860
850 Y(25)=Y(25)+1
860NEXT I
870 M=1
880 PRINT'FROM','TO LESS THAN',' NO.',' BAR GRAPH'
881 PRINT
882 PRINTO,E,Y(51),
890 A5=Y(51)
900 GOSUB 1996
910 FOR J=E TO (F-G) STEP G
920 PRINTJ,J+G,Y(M),
930 A5=Y(M)
940 GOSUB 1996
950 M=M+1
960 NEXT J
970 PRINTF,'INFINITY',Y(25),
980 A5=Y(25)
990 GOSUB 1996
1000 IF N=1 THEN 540
1010 PRINT
1020 PRINT
1030 PRINT'RANKING BY PERCENT ERROR'
1040 PRINT'------- -- ------- -----'
1050 PRINT
1060 PRINT'RANK','STUDENT NO.','PERCENT ERROR'
1070 PRINT
1080 FOR S=1 TO B
1090 T=1.E+25
1100 FOR I=1 TO B
1110 IF W(I)>=T THEN 1140
1120 T=W(I)
1130 V=I
1140 NEXT I
1150 PRINTS,V,W(V)
1160 W(V)=1.E+25
1170 NEXT S
1180 PRINT
1190 IF N=2 THEN 540
1200 PRINT
1210 PRINT
1220 PRINT
1230 PRINT'RANKING BY EXPERIMENTAL VALUE'
1240 PRINT'------- -- ------------ -----'
1250 PRINT
1260 PRINT' RANK','STUDENT NO.','EXPER. VALUE'
1261 PRINT
1270 FOR S=1 TO B
1280 T=1.E+25
1290 FOR I=1 TO B
1300 IF A(I)>=T THEN 1330
1310 T=A(I)
1320 V=I
1330 NEXT I
1340 PRINTS,V,A(V)
1350 Q(S)=A(V)
1360 A(V)=1.E+25
1370 NEXT S
1380 F=0
1390 G7=B/2
1400 FOR I=1 TO 31
1410 IF G7=I THEN 1470
1420 NEXT I
1430 Z2=INT(G7)
1440 Z4=INT(G7+1)
1450 F=(Q(Z2)+Q(Z4))/2
1460 GOTO 1540
1470 F=(Q(G7))
1480 IF N=3 THEN 540
1490 PRINT
1500 PRINT
1510  PRINT'OTHER INFORMATION'
1520 PRINT'----- -----------'
1530 PRINT
1540 PRINT'THE MEDIAN VALUE IS ';F
1550 PRINT'THE LOWEST VALUE IS ';Q(1)
1560 PRINT'THE HIGHEST VALUE IS  ';Q(B)
1570 M=0
1580 P9=F
1590 FOR I=1 TO B
1600 P3=Q(I)-P9
1610 M=M+P3
1620 NEXT I
1630 M7=M/B
1640 PRINT'THE MEAN DEVIATION (AVERAGE DEVIATION) IS'M7
1650 M=0
1660 FOR I=1 TO B
1670 P3=(Q(I)-P2)^2
1680 M=M+P3
1690 NEXT I
1700 M7=SQR(M/B)
1710 PRINT'THE STANDARD DEVIATION IS 'M7
1720 IF N=4 THEN 540
1730 PRINT'THIS CONCLUDES THE RUN.'
1740 GOTO 2001
1750 PRINT'THIS PROGRAM WILL DO THE FOLLOWING:'
1760 PRINT'	1.  PRINT OUT A TABLE OF EXPERIMENTAL VALUES, ERRORS,'
1770 PRINT'	    AND PERCENT ERRORS INCLUDING AVERAGES FOR ALL.'
1780 PRINT'	2.  GIVE YOU A CHOICE OF ALL OF THE FOLLOWING:'
1790 PRINT'	      CHOICE 1 - EXPERIMENTAL VALUE DISTRIBUTION'
1800 PRINT'	      		 INCLUDING A BAR GRAPH'
1810 PRINT'	      CHOICE 2 - RANKING BY PERCENT ERROR'
1820 PRINT'	      CHOICE 3 - RANKING BY EXPERIMENTAL VALUE'
1840 PRINT'           CHOICE 4 - ALL OF THE ABOVE CHOICES'
1990 PRINT
1995 RETURN
1996 FOR I=1 TO A5
1997 PRINT'*';
1998 NEXT I
1999 PRINT
2000 RETURN
2001 END