home *** CD-ROM | disk | FTP | other *** search
/ Power-Programmierung / CD2.mdf / c / library / dos / math / cephes / cprob / gdtr.c < prev    next >
Encoding:
C/C++ Source or Header  |  1992-11-17  |  1.8 KB  |  128 lines
  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.0:  April, 1987
  95. Copyright 1984, 1987 by Stephen L. Moshier
  96. Direct inquiries to 30 Frost Street, Cambridge, MA 02140
  97. */
  98.  
  99. #include "mconf.h"
  100.  
  101. double gdtr( a, b, x )
  102. double a, b, x;
  103. {
  104. double igam();
  105.  
  106. if( x < 0.0 )
  107.     {
  108.     mtherr( "gdtr", DOMAIN );
  109.     return( 0.0 );
  110.     }
  111. return(  igam( b, a * x )  );
  112. }
  113.  
  114.  
  115.  
  116. double gdtrc( a, b, x )
  117. double a, b, x;
  118. {
  119. double igamc();
  120.  
  121. if( x < 0.0 )
  122.     {
  123.     mtherr( "gdtrc", DOMAIN );
  124.     return( 0.0 );
  125.     }
  126. return(  igamc( b, a * x )  );
  127. }
  128.