BUG: Improved accuracy of scipy.stats.truncexpon

Fixes the bug related on #4035
IIRC numpy log1p and/or expm1 have problems in corner cases on Windows in some numpy versions --- maybe better use versions from `scipy.special`