Thanks. The Sage trac ticket is https://trac.sagemath.org/ticket/21801 
(thanks to kcrisman).

  John


On Wednesday, November 2, 2016 at 4:17:00 PM UTC-7, Dima Pasechnik wrote:
>
> This is now https://sourceforge.net/p/maxima/bugs/3236/
>
>
> On Wednesday, November 2, 2016 at 11:06:45 PM UTC, Dima Pasechnik wrote:
>>
>> It is a bug in Maxima (I also checked with Maxima compiled with SBCL, 
>> another lisp compiler, more popular than ECL we use). I think it does not 
>> pay attention
>> to the 0th term being negative.
>> Indeed, if you sum from 1, not from 0, it correctly
>> outputs (2 log(2)+1)/3.
>>
>> Dima
>>
>> On Wednesday, November 2, 2016 at 7:09:13 PM UTC, John H Palmieri wrote:
>>>
>>> Sage is not computing a particular infinite sum correctly.
>>>
>>> sage: sum(1/((n+1)*(2*n-1)), n, 0, 1000).n()
>>> -0.205068171626375
>>> sage: sum(1/((n+1)*(2*n-1)), n, 0, 10000).n()
>>> -0.204618542543703
>>> sage: sum(1/((n+1)*(2*n-1)), n, 0, 100000).n()  # seems to be converging
>>> -0.204573546255870
>>> sage: sum(1/((n+1)*(2*n-1)), n, 0, oo).n()      # but not to this number
>>> -1.09345743518226
>>>
>>> sage: sum(1/((n+1)*(2*n-1)), n, 0, oo)
>>> 2/3*log(2) - 14/9
>>>
>>> I think the answer should be 2/3*log(2) - 2/3 -- that's what Mathematica 
>>> says, and it is also consistent with the partial sums. (See 
>>> https://ask.sagemath.org/question/35354/sage-seems-to-be-improperly-computing-an-infinite-sum-and-giving-an-incorrect-answer/.)
>>>  
>>> Is this a bug in Maxima? I don't know enough about Maxima's syntax to 
>>> evaluate it directly there. Just for kicks, I also installed "giac", since 
>>> that was listed as an option for the "algorithm" keyword for "sum":
>>>
>>> sage: sum(1/((n+1)*(2*n-1)), n, 0, oo, algorithm='giac')
>>> +Infinity
>>>
>>> Not an improvement.
>>>
>>> -- 
>>> John
>>>
>>>

-- 
You received this message because you are subscribed to the Google Groups 
"sage-support" group.
To unsubscribe from this group and stop receiving emails from it, send an email 
to sage-support+unsubscr...@googlegroups.com.
To post to this group, send email to sage-support@googlegroups.com.
Visit this group at https://groups.google.com/group/sage-support.
For more options, visit https://groups.google.com/d/optout.

Reply via email to