Return oo for harmonic(oo), this fixes #4668

This value is consistent with this reference from the docstring:
http://functions.wolfram.com/GammaBetaErf/HarmonicNumber/