Re: [sage-devel] Bug in integration on a Riemann surface

2024-01-22 Thread Dima Pasechnik
On Mon, Jan 22, 2024 at 11:03 AM Linden Disney wrote: > > In Sage 9.8 I ran the following code: > > R. = QQ[] > f = x*(1^5+z^5) + (x*1*z)^2 - x^4*1*z - 2*1^3*z^3 > S = Curve(f).riemann_surface() > S.riemann_matrix() > > and got a ValueError occurring inside rigorous_line_integral.

[sage-devel] Bug in integration on a Riemann surface

2024-01-22 Thread Linden Disney
In Sage 9.8 I ran the following code: R. = QQ[] f = x*(1^5+z^5) + (x*1*z)^2 - x^4*1*z - 2*1^3*z^3 S = Curve(f).riemann_surface() S.riemann_matrix() and got a ValueError occurring inside rigorous_line_integral. This error doesn't occur if instead for the first line I used "R. = QQ

Re: [sage-devel] bug in integration

2012-04-17 Thread Andrew Fleckenstein
I reproduced this in sage 4.8, but with different results: sage: f(x) = sin(x)^2/x^2 sage: f.integral(x, -5, 5).n() -0.20071537570 sage: f.integral(x, -500, 500).n() -2.0008410959e-7 sage: f.integral(x, -5, 5).n() -2.109169e-9 The answer seems t

Re: [sage-devel] bug in integration

2012-04-17 Thread Andrew Fleckenstein
I reproduced this in sage 4.8, but with different results: sage: f(x) = sin(x)^2/x^2 sage: f.integral(x, -5, 5).n() -0.20071537570 sage: f.integral(x, -500, 500).n() -2.0008410959e-7 sage: f.integral(x, -5, 5).n() -2.109169e-9 The answer seems t

[sage-devel] bug in integration

2012-04-17 Thread ggrafendorfer
Hi, sage 5.0 beta11, on Fedora 15, AMD Phenom II X4: sage: f(x) = sin(x)^2/x^2 sage: f.integral(x, -infinity, infinity) x |--> pi sage: f.integral(x, -infinity, infinity).n() 3.14159265358979 sage: f.integral(x, -5, 0).n() + f.integral(x, 0, 5).n() 3.1415769097886317 everthing fine so fa