[Top][All Lists]
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
[Axiom-developer] [#6 wrong integration result] (new)
From: |
Bob McElrath |
Subject: |
[Axiom-developer] [#6 wrong integration result] (new) |
Date: |
Wed, 22 Dec 2004 22:04:17 -0600 |
anonymous address@hidden wrote:
> integrate((x^2)*(%e^(-x^2)),x=0..%plusInfinity) does not give
> sqrt(%pi)/4.
The answer it gives is correct (unless there is a new bug introduced
that is not on my axiom 0.20040831-1). The gamma function has sqrt(pi)
for half-integer values.
\begin{axiom}
integrate(x^2*(exp(-x^2)), x=0..%plusInfinity)
Gamma(3/2)/2::Float
sqrt(%pi)/4::Float
\end{axiom}
What I don't understand is why the integration result is so hard to
convert to a Float. Appending '::Float' to the 'integrate()' command
fails. Why? If I type 'Gamma(3/2)/2' by hand it has no problem.
--
forwarded from http://page.axiom-developer.org/zope/mathaction/address@hidden