[sage-devel] Re: weird integral behavior

2014-07-22 Thread Emmanuel Charpentier
See answer in text. Le samedi 19 juillet 2014 18:25:42 UTC+2, Anne Schilling a écrit : > > Hi! > > [ Bandwidth savings ] > > Pulling out e^{2\pi i x} to simplify the integral to a one dimensional > integral, > ??? Again, I can't follow you. Care to explain ? > Sage can solve this numerica

[sage-devel] Re: weird integral behavior

2014-07-20 Thread Emmanuel Charpentier
Le samedi 19 juillet 2014 18:25:42 UTC+2, Anne Schilling a écrit : > > Hi! > > David Bailey (http://www.davidhbailey.com/) showed Viviane, Travis, and > myself the > following oddity yesterday. > > Take the integral > > \int_0^1 \int_0^1 |e^{2\pi i x} + e^{2\pi i y}| dx dy > > The answer sh

[sage-devel] Re: weird integral behavior

2014-07-19 Thread Anne Schilling
Dear Nils, Thanks for your answer! Even if the symbolic integration techniques that maxima uses cannot tackle this problem, it would be good if the computer could somehow catch that there is a problem instead of returning the wrong answer. Best wishes, Anne On 7/19/14 6:33 PM, Nils Bruin wrote:

[sage-devel] Re: weird integral behavior

2014-07-19 Thread Nils Bruin
On Saturday, July 19, 2014 9:25:42 AM UTC-7, Anne Schilling wrote: > > Sage can solve this numerically: > > sage: g = lambda x : (1+e^(2*pi*I*x)).abs() > sage: numerical_integral(g,0,1) > (1.2732395447351625, 1.4155343563970746e-14) > sage: n(4/pi) > 1.27323954473516 > > but not symbolically: