bridge/class-master/depend/pari/src/basemath/trans3.c
/* return L * approximate solution to x log x = B */
return L * approximate solution to x log x = B
return L * approximate solution to x log x = B
static long
bessel_get_lim(double B, double L)
{
  long lim;
  double x = 1 + B; 
  x = (x + B)/(log(x)+1);
  x = (x + B)/(log(x)+1);
  x = (x + B)/(log(x)+1); x = L*x;
  lim = (long)x; if (lim < 2) lim = 2;
  return lim;
}
static long bessel_get_lim ( double B , double L ) { long lim ; double x = 1 + B ; x = ( x + B ) / ( log ( x ) + 1 ) ; x = ( x + B ) / ( log ( x ) + 1 ) ; x = ( x + B ) / ( log ( x ) + 1 ) ; x = L * x ; lim = ( long ) x ; if ( lim < 2 ) lim = 2 ; return lim ; }
return L * approximate solution to x log x = B
return L * approximate solution to x log x = B
/* 3 iterations are enough except in pathological cases */