This is strange:

sage: (integrate(exp(-x^2)*cos(x),(x,-pi,pi))/pi).n()

complains that "Unable to convert 3.14159265358979 - 0.500000000000000*I
to float". Which is odd, given that I'm integrating a real-valued
function. But it gets weirder: if I try to use .abs(), it *still* complains:

sage: (integrate(exp(-x^2)*cos(x),(x,-pi,pi))/pi).abs().n()

...same error.

Meanwhile:

sage: numerical_integral(exp(-x^2)*cos(x), (-pi,pi))       
(1.3804038617166086, 1.5567746861731315e-14)

I'm guessing this is a Maxima problem. Any ideas what's going on?

Dan

--
---  Dan Drake
-----  http://mathsci.kaist.ac.kr/~drake
-------

Attachment: signature.asc
Description: Digital signature

Reply via email to