! fdist.f90
! module contains routines for F-distribution testing
! it uses FDist_K and FDist_BR

module FDist

implicit none

contains

!*******************************************************************************

SUBROUTINE fprob(num, denom, ratio, p, q, ier)
! Calculates the probability p that a variable with the F-distribution
! with num d.f. for the numerator, and denom d.f. for the denominator
! is less than the value ratio.   q gives the probability in the other
! tail, so that p + q = 1.

! Most of the code is translated from the Fortran 66 code for the incomplete
! beta function from the Naval Surface Warfare Center library and was written
! by Alfred Morris.

! This Fortran 90 version by Alan Miller
! e-mail: amiller@bigpond.net.au
! WWW:    http://users.bigpond.net.au/amiller

! Latest revision - 24 June 1997

USE FDist_K
USE FDist_BR
IMPLICIT NONE
INTEGER, INTENT(IN)    :: num, denom
REAL (dp), INTENT(IN)  :: ratio
INTEGER, INTENT(OUT)   :: ier
REAL (dp), INTENT(OUT) :: p, q

!     Local variables
REAL (dp) :: zero = 0.D0, one = 1.D0, two = 2.D0, x

ier = 0
IF (num <= 0 .OR. denom <= 0) THEN
  ier = 1
  RETURN
ELSE IF (ratio < zero) THEN
  ier = 2
  RETURN
END IF

x = DBLE(denom) / (DBLE(denom) + DBLE(num)*ratio)
CALL bratio(DBLE(denom)/two, DBLE(num)/two, x, one-x, q, p, ier)

RETURN
END SUBROUTINE fprob

!*******************************************************************************

end module FDist