`integrate()' function on non-elementary integrals outputs sometimes a
limit function, ie:
In: k = var("m", domain='positive', latex_name=r"\hat k" );
In: N_nr = integrate(k^2/(exp(k)+1), (k, 0, oo))
Out: `limit(1/3*k^3 - k^2*log(e^k + 1) - 2*k*dilog(-e^k) + 2*polylog(3,
-e^k), k, +Infinity, minus) + 3/2*zeta(3)`
Analytically, the part under the limit is equal to 0 (compare to
mathematica's N_nr = Integrate[k^2/(Exp[k] + 1), {k, 0, Infinity}]
On simplify, the engine tries to evaluate the limit using maxima and fails
In: N_nr.simplify()
Out: TypeError: ECL says: Error executing code in Maxima: limit: direction
must be either 'plus' or 'minus'; found: _SAGE_VAR_minus
Something seems to be off from putting this sage expression to maxima.
Version: Sage 9.3
OS: Windows 10 Pro, 21H2, x64
--
You received this message because you are subscribed to the Google Groups
"sage-devel" group.
To unsubscribe from this group and stop receiving emails from it, send an email
to [email protected].
To view this discussion on the web visit
https://groups.google.com/d/msgid/sage-devel/3a020049-3903-4bd0-b81e-133f936edaean%40googlegroups.com.
