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

This value is consistent with this reference from the docstring:
http://functions.wolfram.com/GammaBetaErf/HarmonicNumber/
Merge pull request #7420 from skirpichev/4668-harmonic-at-oo

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