bdtrc

Complemented binomial distribution.

SYNOPSIS:

int k, n;
double p, y, bdtrc();
y = bdtrc(k, n, p);

DESCRIPTION:

Returns the sum of the terms k+1 through n of the Binomial
probability density:

   n
   -- ( n )   j      n-j
   >   (   ) p (1-p)
   -- ( j )
j=k+1

The terms are not summed directly; instead the incomplete beta integral is employed, according to the formula:

y = bdtrc(k, n, p) = incbet(k+1, n-k, p).

The arguments must be positive, with p ranging from 0 to 1.

ACCURACY:

Tested at random points (a,b,p).

               a,b                     Relative error:
arithmetic domain     # trials      peak         rms
For p between 0.001 and 1:
    IEEE     0,100       100000      6.7e-15     8.2e-16
For p between 0 and .001:
    IEEE     0,100       100000      1.5e-13     2.7e-15

ERROR MESSAGES:

message       condition      value returned
domain       x<0, x>1, n<k         0.0