BUG: special: scrub corner cases in zeta, exp1

Add asymptotic expansion to zeta, so it overflows later. Improve exp1/expi accuracy on complex plane near negative real axis.

Compares OK to mpmath now. Fixes gh-1967 I believe.