1  REM  ****  HP BASIC PROGRAM LIBRARY  *******************************
2  REM
3  REM        TESTUD:   TEST UNKNOWN POPULATION MEAN
4  REM
5  REM        36722 REV  A   10/73
6  REM
7  REM  ****  CONTRIBUTED PROGRAM  ************************************
12  DATA 5.E+06,5.39828E+06,5.7926E+06,6.17911E+06,6.55422E+06,6.91462E+06,7.25747E+06
13  DATA 7.58036E+06,7.88145E+06,8.1594E+06,8.41345E+06,8.64334E+06,8.8493E+06,9.032E+06
14  DATA 9.19243E+06,9.33193E+06,9.45201E+06,9.55435E+06,9.6407E+06,9.71283E+06,9.7725E+06
15  DATA 9.82136E+06,9.86097E+06,9.89276E+06,9.91803E+06,9.9379E+06,9.95339E+06,9.96533E+06
16  DATA 9.97445E+06,9.98134E+06,9.9865E+06,9.99032E+06,9.99313E+06,9.99517E+06,9.99663E+06
17  DATA 9.99767E+06,9.99841E+06,9.99892E+06,9.99928E+06,9.99952E+06,9.99968E+06,9.99979E+06
18  DATA 9.99987E+06,9.99992E+06,9.99995E+06,9.99997E+06,9.99998E+06,9.99999E+06,9.99999E+06
19  DIM X[49]
180  DEF FND()=X[]-X[-1]
190  DEF FNT()=1-((^2)+1)/(4*D)+(13*(^2)^2+8*(^2)+3)/(96*D^2)
199  FOR I=1 TO 49
200  READ X[I]
201  NEXT I
210  PRINT 
220  PRINT "THIS PROGRAM PERFORMS CALCULATIONS NECESSARY FOR"
232  PRINT "TESTING AN UNKNOWN POPULATION MEAN USING SAMPLE"
240  PRINT "STATISTICS. WHAT ARE N (THE SAMPLE SIZE), M (THE"
250  PRINT "SAMPLE MEAN), S (THE SAMPLE STANDARD DEVIATION),"
260  PRINT "W (POPULATION SIZE, ZERO IF INFINITE), AND X (THE"
270  PRINT "POPULATION MEAN TO BE TESTED)";
280  INPUT N,M,S,W,X
290  LET D=N-1
300  PRINT 
310  PRINT "  BASED ON THE STUDENT'S T-DISTRIBUTION WITH";D
320  PRINT "  DEGREES OF FREEDOM, THE PROBABILITY OF FINDING"
330  PRINT "  A SAMPLE MEAN THIS MUCH";
340  IF M<X THEN 370
350  PRINT "GREATER";
360  GOTO 380
370  PRINT "LESS";
380  PRINT " THAN THE POPUL-"
390  PRINT "  ATION MEAN IS ";
400  IF W>0 THEN 420
410  LET W=1.E+25
420  LET S=S*SQR((W-1)/(W*D))
430  LET B1=(M-X)/S
440  LET B1=B1*FNT(B1)
450  GOSUB 530
460  IF B2<.5 THEN 480
470  LET B2=1-B2
480  IF B2<.00001 THEN 510
490  PRINT .00001*INT(.5+100000.*B2)
500  STOP 
510  PRINT "LESS THAN 1 IN 100,000"
520  STOP 
530  IF B1<-4.5 THEN 640
540  IF B1<0 THEN 610
550  IF B1<4.5 THEN 580
560  LET B2=1
570  GOTO 650
580  GOSUB 660
590  LET B2=Q
600  GOTO 650
610  GOSUB 660
620  LET B2=1-Q
630  GOTO 650
640  LET B2=0
650  RETURN 
660  LET Z=10*ABS(B1)
670  LET K=INT(Z)
680  LET D1=Z-K
690  LET Q=X[K]+D1*FND(K+1)+(D1*(D1-1)/2)*(FND(K+2)-FND(K+1))
700  LET Q=.000001*INT(.5+.1*Q)
710  RETURN 
720  END 
