Numerical recipes

#include
6.2 Incomplete Gamma Function 221
float erfcc(float x)
Returns the complementary error function erfc(x) with fractional error everywhere less than
1.2 × 10 −7 .
{
float t,z,ans;
}
z=fabs(x);
t=1.0/(1.0+0.5*z);
ans=t*exp(-z*z-1.26551223+t*(1.00002368+t*(0.37409196+t*(0.09678418+
t*(-0.18628806+t*(0.27886807+t*(-1.13520398+t*(1.48851587+
t*(-0.82215223+t*0.17087277)))))))));
return x >= 0.0 ? ans : 2.0-ans;
There are also some functions of two variables that are special cases of the
incomplete gamma function:
Cumulative Poisson Probability Function
Px(< k), for positive x and integer k ≥ 1, denotes the cumulative Poisson
probability function. It is defined as the probability that the number of Poisson
random events occurring will be between 0 and k − 1 inclusive, if the expected mean
number is x. It has the limiting values
Px(< 1) = e −x
Its relation to the incomplete gamma function is simply
Chi-Square Probability Function
Px(< ∞) =1 (6.2.14)
Px(< k)=Q(k, x) =gammq (k, x) (6.2.15)
P (χ 2 |ν) is defined as the probability that the observed chi-square for a correct
model should be less than a value χ 2 . (We will discuss the use of this function in
Chapter 15.) Its complement Q(χ 2 |ν) is the probability that the observed chi-square
will exceed the value χ 2 by chance even for a correct model. In both cases ν is an
integer, the number of degrees of freedom. The functions have the limiting values
P (0|ν) =0 P (∞|ν) =1 (6.2.16)
Q(0|ν) =1 Q(∞|ν) =0 (6.2.17)
and the following relation to the incomplete gamma functions,
P (χ 2 � �
ν χ2
|ν) =P ,
2 2
�
ν
= gammp
χ2
,
2 2
Q(χ 2 � � � �
ν χ2
ν χ2
|ν) =Q , = gammq ,
2 2
2 2
�
(6.2.18)
(6.2.19)
Sample page from NUMERICAL RECIPES IN C: THE ART OF SCIENTIFIC COMPUTING (ISBN 0-521-43108-5)
Copyright (C) 1988-1992 by Cambridge University Press. Programs Copyright (C) 1988-1992 by Numerical Recipes Software.
Permission is granted for internet users to make one paper copy for their own personal use. Further reproduction, or any copying of machinereadable
files (including this one) to any server computer, is strictly prohibited. To order Numerical Recipes books or CDROMs, visit website
http://www.nr.com or call 1-800-872-7423 (North America only), or send email to directcustserv@cambridge.org (outside North America).