Barton Willis - 2026-05-01

A recent build gives

(%i2) display2d : false$

(%i3) qqq : %e^(%e^(1/x)*gamma(x-%e^(-x)))-%e^gamma(x);

(%o3) %e^(%e^(1/x)*gamma(x-%e^-x))-%e^gamma(x)
(%i4) limit(qqq,x,inf);

(%o4) -(%e^-1*(%e^2-sqrt(%e)))

It's nice that we don't get thrown into the debugger, but the answer is wrong--I'm pretty sure that the correct value is infinity.