htswanalysis/MACS/lib/gsl/gsl-1.11/cdf/beta.c

   1 /* cdf/cdf_beta.c
   2  *
   3  * Copyright (C) 2003, 2007 Brian Gough.
   4  *
   5  * This program is free software; you can redistribute it and/or modify
   6  * it under the terms of the GNU General Public License as published by
   7  * the Free Software Foundation; either version 3 of the License, or (at
   8  * your option) any later version.
   9  *
  10  * This program is distributed in the hope that it will be useful, but
  11  * WITHOUT ANY WARRANTY; without even the implied warranty of
  12  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  13  * General Public License for more details.
  14  *
  15  * You should have received a copy of the GNU General Public License
  16  * along with this program; if not, write to the Free Software
  17  * Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307, USA.
  18  */
  19
  20 #include <config.h>
  21 #include <math.h>
  22 #include <gsl/gsl_cdf.h>
  23 #include <gsl/gsl_sf_gamma.h>
  24 #include <gsl/gsl_math.h>
  25
  26 #include "beta_inc.c"
  27
  28 double
  29 gsl_cdf_beta_P (const double x, const double a, const double b)
  30 {
  31   double P;
  32
  33   if (x <= 0.0 )
  34     {
  35       return 0.0;
  36     }
  37
  38   if ( x >= 1.0 )
  39     {
  40       return 1.0;
  41     }
  42
  43   P = beta_inc_AXPY (1.0, 0.0, a, b, x);
  44
  45   return P;
  46 }
  47
  48 double
  49 gsl_cdf_beta_Q (const double x, const double a, const double b)
  50 {
  51   double Q;
  52
  53   if ( x >= 1.0)
  54     {
  55       return 0.0;
  56     }
  57
  58   if ( x <= 0.0 )
  59     {
  60       return 1.0;
  61     }
  62
  63   Q = beta_inc_AXPY (-1.0, 1.0, a, b, x);
  64
  65   return Q;
  66 }