2021-04-09 05:28:59 +0200 received badge ● Famous Question (source) 2020-04-03 20:48:53 +0200 received badge ● Popular Question (source) 2020-04-03 20:48:53 +0200 received badge ● Notable Question (source) 2018-01-07 02:11:46 +0200 asked a question Help summing an infinite series I have been trying to evaluate the following series using SageMath.  sum((x^x)/(factorial(x)*exp(x)) - 1/sqrt(2*pi*x), x, 1, +oo).n()  The answer should be -0.084069508727655... (-2/3 - zeta(1/2)/sqrt(Pi)), but I get the following error when the upper limit is above 96: RuntimeError: maximum recursion depth exceeded while calling a Python object  Is there any way SageMath could give me even the first few digits of the infinite sum? This sum seems to evaluate just fine in Maple and Mathematica, but not in any open source program I have tried, including Axiom, Maxima, PARI, and SymPy. Any help would be greatly appreciated.