Revision 1244

trunk/src/linad99/ccumdbetainv.cpp (revision 1244)
21 21

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

  
27 26
  int icount=0;
......
59 58
  double d=0.0;
60 59
  do
61 60
  {
62
    double xm;
63 61
    double f,dx; // der of x wrt s
64 62
    x=1.0/(1.0+exp(-s));  //transform from s to x
65 63

  

Also available in: Unified diff