← Back to branch summary

~ubuntu-branches/ubuntu/karmic/grace/karmic

« back to all changes in this revision

Viewing changes to cephes/gdtr.c

  • Committer: Bazaar Package Importer
  • Author(s): Torsten Werner
  • Date: 2002-03-19 14:19:58 UTC
  • Revision ID: james.westby@ubuntu.com-20020319141958-5gxna6vo1ek3zjml
Tags: upstream-5.1.7
Import upstream version 5.1.7

Show diffs side-by-side

added added

removed removed

Lines of Context:
cephes/gdtr.c
 
1
/*                                                      gdtr.c
 
2
 *
 
3
 *      Gamma distribution function
 
4
 *
 
5
 *
 
6
 *
 
7
 * SYNOPSIS:
 
8
 *
 
9
 * double a, b, x, y, gdtr();
 
10
 *
 
11
 * y = gdtr( a, b, x );
 
12
 *
 
13
 *
 
14
 *
 
15
 * DESCRIPTION:
 
16
 *
 
17
 * Returns the integral from zero to x of the gamma probability
 
18
 * density function:
 
19
 *
 
20
 *
 
21
 *                x
 
22
 *        b       -
 
23
 *       a       | |   b-1  -at
 
24
 * y =  -----    |    t    e    dt
 
25
 *       -     | |
 
26
 *      | (b)   -
 
27
 *               0
 
28
 *
 
29
 *  The incomplete gamma integral is used, according to the
 
30
 * relation
 
31
 *
 
32
 * y = igam( b, ax ).
 
33
 *
 
34
 *
 
35
 * ACCURACY:
 
36
 *
 
37
 * See igam().
 
38
 *
 
39
 * ERROR MESSAGES:
 
40
 *
 
41
 *   message         condition      value returned
 
42
 * gdtr domain         x < 0            0.0
 
43
 *
 
44
 */
 
45
/*                                                     gdtrc.c
 
46
 *
 
47
 *      Complemented gamma distribution function
 
48
 *
 
49
 *
 
50
 *
 
51
 * SYNOPSIS:
 
52
 *
 
53
 * double a, b, x, y, gdtrc();
 
54
 *
 
55
 * y = gdtrc( a, b, x );
 
56
 *
 
57
 *
 
58
 *
 
59
 * DESCRIPTION:
 
60
 *
 
61
 * Returns the integral from x to infinity of the gamma
 
62
 * probability density function:
 
63
 *
 
64
 *
 
65
 *               inf.
 
66
 *        b       -
 
67
 *       a       | |   b-1  -at
 
68
 * y =  -----    |    t    e    dt
 
69
 *       -     | |
 
70
 *      | (b)   -
 
71
 *               x
 
72
 *
 
73
 *  The incomplete gamma integral is used, according to the
 
74
 * relation
 
75
 *
 
76
 * y = igamc( b, ax ).
 
77
 *
 
78
 *
 
79
 * ACCURACY:
 
80
 *
 
81
 * See igamc().
 
82
 *
 
83
 * ERROR MESSAGES:
 
84
 *
 
85
 *   message         condition      value returned
 
86
 * gdtrc domain         x < 0            0.0
 
87
 *
 
88
 */
 
89
 
 
90
/*                                                      gdtr()  */
 
91
 
 
92
 
 
93
/*
 
94
Cephes Math Library Release 2.3:  March,1995
 
95
Copyright 1984, 1987, 1995 by Stephen L. Moshier
 
96
*/
 
97
 
 
98
#include "mconf.h"
 
99
#include "cephes.h"
 
100
 
 
101
double gdtr( a, b, x )
 
102
double a, b, x;
 
103
{
 
104
 
 
105
if( x < 0.0 )
 
106
        {
 
107
        mtherr( "gdtr", DOMAIN );
 
108
        return( 0.0 );
 
109
        }
 
110
return(  igam( b, a * x )  );
 
111
}
 
112
 
 
113
 
 
114
 
 
115
double gdtrc( a, b, x )
 
116
double a, b, x;
 
117
{
 
118
 
 
119
if( x < 0.0 )
 
120
        {
 
121
        mtherr( "gdtrc", DOMAIN );
 
122
        return( 0.0 );
 
123
        }
 
124
return(  igamc( b, a * x )  );
 
125
}