Hallo,
I am having big troubles with an integral: A and B (see below) should clearly
be equal, but they aren't (A is correct, B is wrong).
What's happening?
Do you have any idea of what's going wrong?
Thanks,
Andrea
var('x, p, q')
assume(p,'integer'); assume(p>0)
assume(q,'integer'); assume(q>0)
fun = exp(-x^2) * x^(2*p) * x^q
A = integral(fun, (x,-oo,oo))
B = integral(fun, (x,-oo,0)) + integral(fun, (x,0,oo))
----
print A
sage: 0
print B
sage: 1/2*(-1)^q*gamma(p + 1/2*q + 1/2) + 1/2*gamma(p + 1/2*q + 1/2)
--
You received this message because you are subscribed to the Google Groups
"sage-support" group.
To post to this group, send email to [email protected].
To unsubscribe from this group, send email to
[email protected].
Visit this group at http://groups.google.com/group/sage-support?hl=en.
