10X$=SYS(CHR$(7)):R$=LEFT(X$,7):L$=RIGHT(X$,8)
13  REM *****************************************************************
14  REM    CMOD7   CMOD7     CMOD7      CMOD7      CMOD7     CMOD7
16  REM *****************************************************************
20  REM
30  REM         ENTRY OF TWO T DISTRIBUTIONS FOR COMP. OF MEANS
40  REM
50  REM****************************************************************
60  DIM #1,A1(111),A2(111),A3(111)
70  OPEN R$ AS FILE 1
110A1(0)=7
120  PRINT L$
130  PRINT "                 COMPARISON OF NORMAL MEANS"
140  PRINT
150  PRINT "THIS PROGRAM ALLOWS YOU TO EXAMINE THE CHARACTERISTICS OF THE"
160  PRINT "POSTERIOR BEHRENS-FISHER DISTRIBUTION ON THE DIFFERENCE OF TWO"
170  PRINT "NORMAL MEANS.  YOU NEED TO KNOW THE PARAMETERS OF THE MARGINAL"
180  PRINT "POSTERIOR T-DISTRIBUTIONS ON THE TWO MEANS.  IF YOU DO NOT KNOW"
190  PRINT "THESE THEY CAN BE OBTAINED BY DOING TWO SEPARATE TWO-PARAMETER"
200  PRINT "NORMAL ANALYSES."
210  PRINT
220  PRINT "IF YOU KNOW THE PARAMETERS AND WANT TO PROCEED TYPE '1'."
230  PRINT "IF YOU WANT TO EXIT MODULE TYPE '0'."
240  GOSUB 9000
250  IF O1 <> 0 THEN 280
270  CLOSE 1: CHAIN "RSTRT"
280  IF O1=1 THEN 360
300  PRINT
310  PRINT "REENTER.  INPUT MUST BE 0 OR 1."
320  PRINT
330  GOTO 240
360  PRINT L$
370  PRINT "INPUT THE PARAMETERS OF THE FIRST T-DISTRIBUTION."
380  PRINT
390  PRINT "DEGREES OF FREEDOM";
400  GOSUB 9000
410  N0=O1
420  IF N0 >= 1 THEN 440
430  GOTO 550
440  PRINT
450  PRINT "MEAN";
460  GOSUB 9000
470  M0=O1
480  PRINT
490  PRINT "SCALE PARAMETER";
500  GOSUB 9000
510  IF O1>0 THEN 590
520  PRINT
530  PRINT "REENTER.  SCALE PARAMETER MUST BE POSITIVE AND NON-ZERO."
540  GOTO 480
550  PRINT
560  PRINT "REENTER.  DEGREES OF FREEDOM MUST BE AT LEAST 1."
570  PRINT
580  GOTO 390
590  PRINT
600  Q0=O1
610  PRINT "INPUT THE PARAMETERS OF THE SECOND T-DISTRIBUTION."
620  PRINT
630  PRINT "DEGREES OF FREEDOM";
640  GOSUB 9000
650  IF O1 >= 1 THEN 690
660  PRINT
670  PRINT "REENTER.  DEGREES OF FREEDOM MUST BE AT LEAST 1."
680  GOTO 620
690  PRINT
700  N1=O1
710  PRINT "MEAN";
720  GOSUB 9000
730  M1=O1
740  PRINT
750  PRINT "SCALE PARAMETER";
760  GOSUB 9000
770  IF O1>0 THEN 810
780  PRINT
790  PRINT "REENTER.  SCALE PARAMETER MUST BE POSITIVE AND NON-ZERO."
800  GOTO 740
810  Q1=O1
820  M=M0-M1
830  E0=SQR(Q0/N0+Q1/N1)
840  T0=ATN(SQR((Q0/N0)/(Q1/N1)))
850  T0=T0*180/3.14159
860  PRINT L$
870  PRINT "THE DISTRIBUTION OF THE DIFFERENCE OF THE TWO T-DISTRIBUTIONS IS A"
880  PRINT "BEHRENS-FISHER DISTRIBUTION WITH THE FOLLOWING 5 PARAMETERS:"
890  PRINT
910A$="NU1 (DEGREES OF FREEDOM OF 1ST DIST.)     ######.###"
911  PRINT  USING  A$,N0
930A$="NU2 (DEGREES OF FREEDOM OF 2ND DIST.)     ######.###"
931  PRINT  USING  A$,N1
950A$="PSI (IN DEGREES)                           #####.###"
951  PRINT  USING  A$,T0
970A$="ZETA (MEAN OF 1ST MINUS MEAN OF 2ND)     #######.###"
971  PRINT  USING  A$,M
990A$="EPSILON (SCALE PARAMETER)                #######.###"
991  PRINT  USING  A$,E0
1000  PRINT
1010  PRINT "WHEN YOU ARE READY TO CONTINUE TYPE '1'.";
1020  GOSUB 9000
1030A2(0)=N0:A2(1)=N1:A2(2)=T0:A2(3)=E0:A2(4)=M
1040  CLOSE 1: CHAIN "CMODF"
9000  REM--SUBROUTINE THAT DETERMINES IF RESTART HAS BEEN REQUESTED.
9005  INPUT O1
9015  IF O1=-9999 THEN 9025
9020  RETURN
9025  CLOSE 1: CHAIN "RSTRT"
9035  REM*************END ROUTINE
9999  END