Re: [sage-devel] Re: Incorrect double integral result

2019-03-14 Thread Christopher Subich
On 3/14/19 8:54 AM, kcrisman wrote: Following up, does this computation make any sense?  Presumably this is the core issue. That computation makes sense, it's the second integration that breaks: (x,y) = var('x y') z = sqrt(1-x^2-y^2) print(z) assume(abs(x)<1) assume(abs(y)<1) zint_y = z.inte

[sage-devel] Re: Incorrect double integral result

2019-03-14 Thread kcrisman
On Wednesday, March 13, 2019 at 7:11:02 PM UTC-4, Emmanuel Charpentier wrote: > > This integral turns out to be anything but trivial. > > Neither Sympy, Giac nor Fricas can solve it, in its original form nor in > polar coordinates. BTW, neither Mathematica nor Rubi running on Mathematica > can

[sage-devel] Re: Incorrect double integral result

2019-03-13 Thread Emmanuel Charpentier
This integral turns out to be anything but trivial. Neither Sympy, Giac nor Fricas can solve it, in its original form nor in polar coordinates. BTW, neither Mathematica nor Rubi running on Mathematica can solve it in neither form. The analytic answer given (semi-)trivially by Maxima is probably