Sorry about the formatting, I'm not sure why it didn't work and I cant edit it.
expintegral_ei incorrect behavior at infinity