BUG: Improved accuracy of scipy.stats.truncexpon

Fixes the bug related on #4035
I'd add a new class, TestTruncExpon, having your edge case tests as a method.
Also, there is an comment in the code, hinting at some numericql problems at large values of a. Can you check these as well?