[Axiom-developer] [#6 wrong integration result] (new)

[Bob McElrath]

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