From: SourceForge.net <no...@so...> - 2012-11-15 13:53:56
|
Bugs item #3587514, was opened at 2012-11-15 05:53 Message generated for change (Tracker Item Submitted) made by You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3587514&group_id=4933 Please note that this message will contain a full copy of the comment thread, including the initial issue submission, for this request, not just the latest update. Category: Share Libraries Group: None Status: Open Resolution: None Priority: 5 Private: No Submitted By: https://www.google.com/accounts () Assigned to: Nobody/Anonymous (nobody) Summary: simplify_sum simplifies incorrectly Initial Comment: Maxima 5.28.0 http://maxima.sourceforge.net using Lisp SBCL 1.0.55.0-abb03f9 Distributed under the GNU Public License. See the file COPYING. Dedicated to the memory of William Schelter. The function bug_report() provides bug reporting information. (%i1) load(simplify_sum); (%o1) /Applications/MathApps/Maxima.app/Contents/Resources/maxima/share/maxima\ /5.28.0/share/solve_rec/simplify_sum.mac (%i2) 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 (%o2) ---- + -------- 128 128 But this is not right. The extra -1/32 is spurious. (%i4) float(simplify_sum(sum(1/((2*n-1)^2*(2*n+1)^2*(2*n+3)^2), n, 0, 100))); (%o4) .1156594265749686 (%i5) float (simplify_sum(sum(1/((2*n-1)^2*(2*n+1)^2*(2*n+3)^2), n, 0, inf))); (%o5) .08440942657526591 (%i6) float(1/32); (%o6) 0.03125 This was first reported at http://trac.sagemath.org/sage_trac/ticket/13712 and http://ask.sagemath.org/question/1985/unexpected-result-for-the-sum-of-a-series ---------------------------------------------------------------------- You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3587514&group_id=4933 |