wrong evaluation of infinite sum

[burcin]

From this ask.sagemath.org question:

sage: n = var('n')
sage: sum(1/((2*n-1)^2*(2*n+1)^2*(2*n+3)^2), n, 0, oo)
3/256*pi^2 - 1/32

The correct answer is 3/256*pi^2.

We should report this to maxima:

(%i10) load(simplify_sum);
(%o10) /tmpbig/burcin/sage-5.4/local/share/maxima/5.26.0/share/contrib/solve_r\
ec/simplify_sum.mac
(%i11) simplify_sum(sum(1/((2*n-1)^2*(2*n+1)^2*(2*n+3)^2), n, 0, inf));

                                          2
                                       %pi
                                   2   ---- - 4
                                %pi     2
(%o11)                          ---- + --------
                                128      128

Depends on #13973

Component: symbolics

Keywords: maxima

Work Issues: add doctest

Author: Travis Scrimshaw

Branch/Commit: 6e7166a

Reviewer: Jeroen Demeyer, Karl-Dieter Crisman

Issue created by migration from https://trac.sagemath.org/ticket/13712

[kcrisman]

Upstream: Reported upstream. No feedback yet.

[kcrisman]

comment:1

I can confirm this is true in Maxima 5.28.0 as well. This is now Maxima artifact 3587514.

[jdemeyer]

comment:2

Fixed in Maxima 5.30.0: #13973.

[jdemeyer]

Reviewer: Jeroen Demeyer

[jdemeyer]

Changed upstream from Reported upstream. No feedback yet. to Fixed upstream, in a later stable release.