Changed algorithm for const-mu, const-bound case.

The implemented algorithm is based on Navarro & Fuss (2009) and switches
between two different series expansions, depending on which one is more
accurate.