Re: [sage-devel] Counterexample of Fubini - Strange Result with Maxima

2020-12-20 Thread Emmanuel Charpentier
It is difficult to know before the fact which library will returnthe expected result. Consider : ``` sage: var("x, a, b") (x, a, b) sage: dgamma(x, a, b)=x^(a-1)*(1-x)^(b-1)/beta(a, b) sage: with assuming(a>0, b>0): dgamma(x, a, b).integrate(x,0,1) 1 This is quasi-immediate. And of course so

Re: [sage-devel] Counterexample of Fubini - Strange Result with Maxima

2020-12-18 Thread Michael Orlitzky
On 12/18/20 9:53 AM, Dima Pasechnik wrote: isn't pynac/ginac used a lot? At least for symbolic integration, the current list in symbolic/integration/integral.py is, self.integrators = [external.maxima_integrator, external.giac_integrator, externa

Re: [sage-devel] Counterexample of Fubini - Strange Result with Maxima

2020-12-18 Thread Dima Pasechnik
isn't pynac/ginac used a lot? On Fri, 18 Dec 2020, 14:35 Michael Orlitzky, wrote: > On 12/18/20 2:55 AM, Sébastien Labbé wrote: > > > > Why do we use maxima first as opposed to giac/sympy? Is it because it is > > faster than giac/sympy? Is it because it returns answers that are > > correct but f

Re: [sage-devel] Counterexample of Fubini - Strange Result with Maxima

2020-12-18 Thread Michael Orlitzky
On 12/18/20 2:55 AM, Sébastien Labbé wrote: Why do we use maxima first as opposed to giac/sympy? Is it because it is faster than giac/sympy? Is it because it returns answers that are correct but for which giac/sympy returns incorrect results? I personally have never tried it because I'm afra

Re: [sage-devel] Counterexample of Fubini - Strange Result with Maxima

2020-12-17 Thread Sébastien Labbé
> There is still a lot of room for improvement. SymPy could be tried first > when integrating expressions containing an absolute value, for one. We > already _fall back_ to giac/sympy if maxima throws an error; but when it > simply returns garbage, the problem goes unnoticed. > Why do we use

Re: [sage-devel] Counterexample of Fubini - Strange Result with Maxima

2020-12-17 Thread Michael Orlitzky
On 3/18/20 11:48 AM, Michael Jung wrote: Dear fellow developers, I've encountered a really strange result in Sage while using Maxima. | sage:f(x,y)=(x^2-y^2)/(x^2+y^2)^2 sage:integrate(integrate(abs(f(x,y)),x,0,1),y,0,1) -1/4*pi ||| This is really weird. At least, the result should be positiv

[sage-devel] Counterexample of Fubini - Strange Result with Maxima

2020-03-18 Thread Michael Jung
Dear fellow developers, I've encountered a really strange result in Sage while using Maxima. sage: f(x,y) = (x^2-y^2)/(x^2+y^2)^2 sage: integrate(integrate(abs(f(x,y)), x, 0, 1), y, 0, 1) -1/4*pi This is really weird. At least, the result should be positive! SymPy however yields th

[sage-devel] Counterexample of Fubini - Strange Result with Maxima

2020-03-18 Thread Michael Jung
Dear fellow developers, I've encountered a really strange result in Sage while using Maxima. | sage:f(x,y)=(x^2-y^2)/(x^2+y^2)^2 sage:integrate(integrate(abs(f(x,y)),x,0,1),y,0,1) -1/4*pi | This is really weird. At least, the result should be positive! SymPy however yields the correct result: