On Apr 24, 7:03 pm, Chris Seberino <cseber...@gmail.com> wrote:

Hi,

> Notice first and second are very similar except for the cos() and exp
> ().
>
> Why first one ok but second bombs?
>
> sage: numerical_integral(sin(pi*cos(x/2)),0,2)[0]
> 0.85397903781471396

I guess you are using Sage 3.4?


Works for me in my current merge tree:

----------------------------------------------------------------------
| Sage Version 3.4.2.alpha0, Release Date: 2009-04-24                |
| Type notebook() for the GUI, and license() for information.        |
----------------------------------------------------------------------
sage: numerical_integral(sin(pi*cos(x/2)),0,2)[0]
0.85397903781471396
sage: numerical_integral(sin(pi*cos(x/2)),0,2)[0]
0.85397903781471396
sage: numerical_integral(sin(pi*cos(x/2)),0,2)[0]
0.85397903781471396
sage: numerical_integral(sin(pi*cos(x/2)),0,2)[0]
0.85397903781471396
sage: numerical_integral(sin(pi*cos(x/2)),0,2)[0]
0.85397903781471396
sage: numerical_integral(sin(pi*cos(x/2)),0,2)[0]
0.85397903781471396
sage:

3.4.1:

----------------------------------------------------------------------
| Sage Version 3.4.1, Release Date: 2009-04-21                       |
| Type notebook() for the GUI, and license() for information.        |
----------------------------------------------------------------------
sage: numerical_integral(sin(pi*cos(x/2)),0,2)[0]
0.85397903781471396
sage: numerical_integral(sin(pi*cos(x/2)),0,2)[0]
0.85397903781471396
sage: numerical_integral(sin(pi*cos(x/2)),0,2)[0]
0.85397903781471396
sage: numerical_integral(sin(pi*cos(x/2)),0,2)[0]
0.85397903781471396
sage: numerical_integral(sin(pi*cos(x/2)),0,2)[0]
0.85397903781471396
sage: numerical_integral(sin(pi*cos(x/2)),0,2)[0]
0.85397903781471396
sage: numerical_integral(sin(pi*cos(x/2)),0,2)[0]
0.85397903781471396
sage:

And Sage 3.4:

| Sage Version 3.4, Release Date: 2009-03-11                         |
| Type notebook() for the GUI, and license() for information.        |
----------------------------------------------------------------------
sage: numerical_integral(sin(pi*cos(x/2)),0,2)[0]
0.85397903781471396
sage: numerical_integral(sin(pi*cos(x/2)),0,2)[0]
0.85397903781471396
sage: numerical_integral(sin(pi*cos(x/2)),0,2)[0]
0.85397903781471396
sage: numerical_integral(sin(pi*cos(x/2)),0,2)[0]
0.85397903781471396
sage: numerical_integral(sin(pi*cos(x/2)),0,2)[0]
0.85397903781471396
sage: numerical_integral(sin(pi*cos(x/2)),0,2)[0]
0.85397903781471396
sage: numerical_integral(sin(pi*cos(x/2)),0,2)[0]
0.85397903781471396
sage:

So I cannot reproduce this. What platform are you on and are you using
a binary, build from source, etc?

Cheers,

Michael
--~--~---------~--~----~------------~-------~--~----~
To post to this group, send email to sage-support@googlegroups.com
To unsubscribe from this group, send email to 
sage-support-unsubscr...@googlegroups.com
For more options, visit this group at 
http://groups.google.com/group/sage-support
URLs: http://www.sagemath.org
-~----------~----~----~----~------~----~------~--~---

Reply via email to