[sage-support] Re: Integration issue

2012-05-17 Thread kcrisman
On Thursday, May 17, 2012 1:52:38 AM UTC-4, ketchers wrote: I don't know how to get sage to understand domain : complex so I tried with assume and here is what happened. Does it make sense? Yes, it does. Our assumptions go through Maxima, and apparently assuming a variable is

[sage-support] Re: Integration issue

2012-05-16 Thread ketchers
https://lh6.googleusercontent.com/-BDfI2b3v4HA/T7SRx5zaqBI/AGo/Q6NEGFw8xQw/s1600/Screenshot+from+2012-05-17+00%3A50%3A25.png I don't know how to get sage to understand domain : complex so I tried with assume and here is what happened. Does it make sense? On Sunday, May 13, 2012

[sage-support] Re: Integration issue

2012-05-15 Thread Robert Dodier
On 2012-05-14, JamesHDavenport j.h.davenp...@bath.ac.uk wrote: It may be branch cut strangeness, but if so it is very strange. The integrand is clearly well-behaved, and the integral, while in terms of the incomplete gamma function, seems to be off the usual branch cut (negative real axis).

[sage-support] Re: Integration issue

2012-05-15 Thread kcrisman
It may be branch cut strangeness, but if so it is very strange. The integrand is clearly well-behaved, and the integral, while in terms of the incomplete gamma function, seems to be off the usual branch cut (negative real axis). Try domain:complex before calling integrate; that

[sage-support] Re: Integration issue

2012-05-15 Thread kcrisman
On Tuesday, May 15, 2012 2:36:34 AM UTC-4, Keshav Kini wrote: John H Palmieri jhpalmier...@gmail.com writes: This works for me: sage: numerical_integral(x*cos(x^3), 0, 0.5) (0.1247560409610376, 1.3850702913602309e-15) Interesting... sage:

[sage-support] Re: Integration issue

2012-05-15 Thread Robert Dodier
On 2012-05-15, kcrisman kcris...@gmail.com wrote: (%i3) domain:complex; (%o3) complex (%i4) integrate(x*cos(x^3),x,0,1/2); (%o4) gamma_incomplete(2/3,%i/8)/6+gamma_incomplete(2/3,-%i/8)/6-gamma(2/3)/3 Hmm. I get a different result. I am using the current Git version. domain : complex;

[sage-support] Re: Integration issue

2012-05-15 Thread Jason Grout
On 5/15/12 8:33 PM, Keshav Kini wrote: And maybe that's why plot3d(), unlike plot(), does seem to generate the deprecation warning. Sorry---what plot command doesn't generate a deprecation warning? Thanks, Jason -- To post to this group, send email to sage-support@googlegroups.com To

[sage-support] Re: Integration issue

2012-05-14 Thread kcrisman
This is now http://trac.sagemath.org/sage_trac/ticket/12947. We've had some issues with incomplete gamma functions translating properly in the past, and/or errors in Maxima, but I didn't have time to either look into that or whether there was another ticket open for this, apologies if there

[sage-support] Re: Integration issue

2012-05-14 Thread Robert Dodier
On 2012-05-14, kcrisman kcris...@gmail.com wrote: This is now http://trac.sagemath.org/sage_trac/ticket/12947. We've had some issues with incomplete gamma functions translating properly in the past, and/or errors in Maxima, but I didn't have time to either look into that or whether there

[sage-support] Re: Integration issue

2012-05-14 Thread John H Palmieri
On Sunday, May 13, 2012 8:46:41 PM UTC-7, ketchers wrote: Sage returns negative value for the integral of a positive function x*cos(x^3) on (0,0.5), if I use abs(cos(x^3))*x, then it gets it correct? This works for me: sage: numerical_integral(x*cos(x^3), 0, 0.5)

[sage-support] Re: Integration issue

2012-05-14 Thread JamesHDavenport
It may be branch cut strangeness, but if so it is very strange. The integrand is clearly well-behaved, and the integral, while in terms of the incomplete gamma function, seems to be off the usual branch cut (negative real axis). On Monday, 14 May 2012 15:35:01 UTC+1, Robert Dodier wrote: On