Use Axler (2013) nth_prime bounds
Add Harmonic numbers
nth_prime_lower/upper adjustments