@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
C-FOCAL , 8/68

01.01 C F DISTRIBUTION ROUTINE
01.02 C A.S.FIELDS   US NSRDL-ANNAPOLIS  6/29/70
01.10 A !!! "F-RATIO ",RF,!," FIRST DEGREES OF FREEDOM ",D1,!
01.20 A "SECOND DEGREES OF FREEDOM ",D2,!
01.30 D 2; T %5.04,"PROBABILITY OF CHANCE ",P; G 1.1

02.01 C CONTROL
02.10 S A=D2/2; S B=D1/2
02.20 S X=A; D 3; S QF=1/R; S X=B; D 3; S QF=QF/R
02.30 S X=A+B; D 3; S QF=QF*R; S X=D2/[D2+RF*D1]
02.40 S QF=QF*FEXP[FLOG(X)*A]*FEXP[FLOG(1-X)*B]/A
02.50 S Q=QF; S QQ=1; S N=-1
02.60 S N=N+1; D 4; S QF=QF*X; S Q=Q+QF*QQ
02.70 I (QF*QQ-1E-5)2.8,2.6,2.6
02.80 S P=Q

03.01 C GAMMA FUNCTION
03.10 S R=1
03.20 I (X)3.99,3.99; I (X-1)3.3,3.4,3.4
03.30 S R=R/X; S X=X+1
03.40 I (X-2)3.6,3.6,3.5
03.50 S X=X-1; S R=R*X; G 3.4
03.60 S X=X-1; S XX=1-.57486*X+.95123*X^2-.69986*X^3
03.70 S XX=XX+.42455*X^4-.10106*X^5; S R=R*XX; R
03.99 Q

04.01 C RECURSIVE B(A,B)
04.10 S QQ=QQ*(A+B+N)/(A+1+N)
*@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@O