Quote:
Originally Posted by paul0
The estimate of the probability that a random number below X is Bsmooth is given by u^{u}, where u=ln(X)/ln(B).

No need for Newton method, because it will immediately exceed the precision of the
estimate that you started with. There are
better estimates (and then see Dickman's and/or Buchstab's).