rtoy / maxima

A Clone of Maxima's repo
Other
0 stars 0 forks source link

calculates geometric sum incorrectly inside integrate function #3492

Open rtoy opened 4 months ago

rtoy commented 4 months ago

Imported from SourceForge on 2024-07-08 01:57:01 Created by osmanb on 2020-10-09 15:38:53 Original: https://sourceforge.net/p/maxima/bugs/3661



Maxima version: "5.43.0" Maxima build date: "2019-08-28 00:34:58" Host type: "x86_64-unknown-linux-gnu" Lisp implementation type: "GNU Common Lisp (GCL)" Lisp implementation version: "GCL 2.6.12" User dir: "/home/carol/.maxima" Temp dir: "/tmp" Object dir: "/home/carol/.maxima/binary/5_43_0/gcl/GCL_2_6_12" Frontend: false

simpsum is false by default. Maxima 5.43.0 http://maxima.sourceforge.net using Lisp GNU Common Lisp (GCL) GCL 2.6.12 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) display2d:false;

(%o1) false (%i2) integrate(x^(sum(1/2^k,k,1,inf)),x); Is 'sum(1/2^k,k,1,inf) equal to -1?

no; (%o2) x^('sum(1/2^k,k,1,inf)+1)/('sum(1/2^k,k,1,inf)+1)

This is correct but not in the desired form. simpsum:true prior to invokation results in the correct form. why does it think the sum might be -1? Why ask the user for that simple sum anyways? Same result for SBCL based maxima too.

rtoy commented 4 months ago

Imported from SourceForge on 2024-07-08 01:57:02 Created by robert_dodier on 2020-12-27 02:41:03 Original: https://sourceforge.net/p/maxima/bugs/3661/#4c71


I'm not sure what to think of this. integrate is asking about -1 since if the exponent is -1 it will return a logarithm. Figuring out that the sum is easily calculated isn't really obvious within integrate, and I'm not sure we would want to try all the heuristics known to Maxima (which are collected in the simplify_sum package).

I'm leaning towards thinking that this is a case where it's reasonable to expect the user to help out, by calling simplify_sum for example. If anyone else has an opinion I'll be glad to hear it.

rtoy commented 4 months ago

Imported from SourceForge on 2024-07-08 01:57:06 Created by billingd on 2020-12-27 07:08:57 Original: https://sourceforge.net/p/maxima/bugs/3661/#c873


I also think it is reasonable to expect the user to help out here.

rtoy commented 4 months ago

Imported from SourceForge on 2024-07-08 01:57:10 Created by robert_dodier on 2024-06-13 20:20:00 Original: https://sourceforge.net/p/maxima/bugs/3661/#0fe0