https://github.com/JuliaMath/SpecialFunctions.jl/pull/236 should replace the version in `src/expint.jl` (once a new version of SpecialFunctions is tagged).