On 2014-07-04, Chris Maness <
ch...@chrismaness.com> wrote:
> integrate(exp(-k^2/(4*a))*exp(i*k*x-i*hbar*k^2*t/(2*m)),k,-oo,oo)
Well, I see the integrand is a quadratic (with constant term = 0) in k,
so you can try integrate(exp(A*k^2 + B*k), k, -oo, oo) (Maxima, and
therefore Sage, can solve that) and then substitute A and B into the
result. You can get A and B via ratcoef(foo, k^2) and ratcoef(foo, k)
where foo is the quadratic (I assume there's a way to call ratcoef).
It's certainly reasonable to expect that Maxima should carry out that
analysis on its own; it's not a bug although it is certainly a drawback.
I don't know if e.g. SymPy could solve it; I didn't try.
best
Robert Dodier