/[gsl]/gsl/cdf/f.c
ViewVC logotype

Diff of /gsl/cdf/f.c

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 1.1.2.8 by bjg, Sun Jun 22 21:20:36 2003 UTC revision 1.1.2.9 by bjg, Mon Jun 23 14:05:03 2003 UTC
# Line 20  Line 20 
20  #include <config.h>  #include <config.h>
21  #include <gsl/gsl_cdf.h>  #include <gsl/gsl_cdf.h>
22  #include <gsl/gsl_sf_gamma.h>  #include <gsl/gsl_sf_gamma.h>
23    #include <gsl/gsl_math.h>
24  #include <gsl/gsl_errno.h>  #include <gsl/gsl_errno.h>
25    
26    #include "beta_inc.c"
27    
28  /*  /*
29   * Lower tail.   * Lower tail.
30   */   */
# Line 35  gsl_cdf_F_ratio_P (double x, double nu1, Line 38  gsl_cdf_F_ratio_P (double x, double nu1,
38      {      {
39        double u = x / (r + x);        double u = x / (r + x);
40    
41        P = gsl_sf_beta_inc (nu1 / 2.0, nu2 / 2.0, u);        P = beta_inc_AXPY (1.0, 0.0, nu1 / 2.0, nu2 / 2.0, u);
42      }      }
43    else    else
44      {      {
45        double u = r / (r + x);        double u = r / (r + x);
46    
47        P = 1 - gsl_sf_beta_inc (nu2 / 2.0, nu1 / 2.0, u);        P = beta_inc_AXPY (-1.0, 1.0, nu2 / 2.0, nu1 / 2.0, u);
48      }      }
49    
50    return P;    return P;

Legend:
Removed from v.1.1.2.8  
changed lines
  Added in v.1.1.2.9

savannah-hackers-public@gnu.org
ViewVC Help
Powered by ViewVC 1.1.26