Revision 1258 trunk/src/df1b2-separable/df1b2invcumdbeta.cpp

df1b2invcumdbeta.cpp (revision 1258)
26 26
  df3_three_variable betai(const df3_three_variable& a,const df3_three_variable& b,
27 27
    double x,int maxit);
28 28

  
29
/*
29 30
  static double lnbeta(double a,double b)
30 31
  {
31 32
    return gammln(a)+gammln(b)-gammln(a+b);
......
38 39
    else
39 40
      return -1;
40 41
  }
42
*/
41 43

  
42 44
  double inv_cumd_beta(double _a,double _b,double _y);
43 45
  double inv_cumd_beta_stable(double _a,double _b,double _y,double eps);
......
331 333
static df3_three_variable gammlnguts(const df3_three_variable& _z)
332 334
{
333 335
  df3_three_variable x;
334
  const double lpi =1.1447298858494001741434272;
335
  const double pi =3.1415926535897932384626432;
336 336
  const double lpp =0.9189385332046727417803297;
337 337
  int n=7;
338 338
  const double c[9]={0.99999999999980993,

Also available in: Unified diff