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

Diff of /gsl/cdf/gammainv.c

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

revision 1.1.2.6 by bjg, Mon Jul 14 12:07:33 2003 UTC revision 1.1.2.7 by bjg, Tue Jul 15 17:56:47 2003 UTC
# Line 51  gsl_cdf_gamma_Pinv (double p, double a, Line 51  gsl_cdf_gamma_Pinv (double p, double a,
51      }      }
52    else if (p > 0.95)    else if (p > 0.95)
53      {      {
54        double x0 = -log (1 - p) + lgamma (a);        double x0 = -log1p (-p) + lgamma (a);
55        x = x0;        x = x0;
56      }      }
57    else    else

Legend:
Removed from v.1.1.2.6  
changed lines
  Added in v.1.1.2.7

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