[sage-devel] Re: Fwd: sage-devel "exact" numerical integration

2008-03-05 Thread achrzesz
rzej Chrzęszczyk <[EMAIL PROTECTED]> > > Date: March 5, 2008 6:23:53 PM EST > > To: [EMAIL PROTECTED] > > Subject: sage-devel "exact" numerical integration > > > Dear David > > Try > > > sage: maxima_console() > > (%i1) integrate(%e^(-x^2)

[sage-devel] Re: Fwd: sage-devel "exact" numerical integration

2008-03-05 Thread Jason Grout
Jason Grout wrote: > David Harvey wrote: >> Begin forwarded message: >> >>> *From: *Andrzej Chrzęszczyk <[EMAIL PROTECTED] <mailto:[EMAIL PROTECTED]>> >>> *Date: *March 5, 2008 6:23:53 PM EST >>> *To: [EMAIL PROTECTED] <mailto:[EMAIL

[sage-devel] Re: Fwd: sage-devel "exact" numerical integration

2008-03-05 Thread Jason Grout
David Harvey wrote: > Begin forwarded message: > >> *From: *Andrzej Chrzęszczyk <[EMAIL PROTECTED] <mailto:[EMAIL PROTECTED]>> >> *Date: *March 5, 2008 6:23:53 PM EST >> *To: [EMAIL PROTECTED] <mailto:[EMAIL PROTECTED]> >> *Subject: **sage-devel

[sage-devel] Fwd: sage-devel "exact" numerical integration

2008-03-05 Thread David Harvey
Begin forwarded message: > From: Andrzej Chrzęszczyk <[EMAIL PROTECTED]> > Date: March 5, 2008 6:23:53 PM EST > To: [EMAIL PROTECTED] > Subject: sage-devel "exact" numerical integration > > Dear David > Try > > sage: maxima_co

[sage-devel] "exact" numerical integration

2008-03-05 Thread David Harvey
I tried doing some integrals today and the output doesn't make much sense to me: sage: f = e^(-x^2) sage: f.integrate(x, 0, 0.1) 2066*sqrt(pi)/36741 sage: f.integrate(x, 0, 1/10) sqrt(pi)*erf(1/10)/2 H. Does this mean erf(1/10) is a rational number? That's a little surprising to me. In f