Revision 1267 trunk/src/linad99/ccumdbetainv.cpp

ccumdbetainv.cpp (revision 1267)
10 10
  return gammln(a)+gammln(b)-gammln(a+b);
11 11
}
12 12

  
13
/*
13 14
static int sgn(double z)
14 15
{
15 16
  if (z>=0)
......
17 18
  else
18 19
    return -1;
19 20
}
21
*/
20 22

  
21

  
22 23
double inv_cumd_beta_stable(double a,double b,double y,double eps)
23 24
{
24 25
  double eps1=1.0-eps;

Also available in: Unified diff