#4730 ei-asymptotic-expansion

[Barton Willis]

This is Problem 96 (line 352) of rtest_limit_gruntz:

(%i6) limit(exp(2*expintegral_ei(-x))/x^2, x, 0,plus);

(%o6) 'limit(%e^(2*expintegral_ei(-x))/x^2,x,0,plus)

Why a limit nounform? This limit problem makes its way to limfact, and then applies stirling0. But the function ei-asymptotic-expansion (defined in expintegral.lisp) doesn't handle this case.

It's a bit surprising that limfact handles limits of expintegral expressions, but it is supposed to.

[Dan Gildea]

Fixed by [06f26693d2a59da6bae8c4c701c9532ec1884360]