Hello,
I've got the following problem: Sometimes when I let Maple compute an integral it gives a result which is much smaller than it should be. This happens even for rather simple integrals such as the product of a Gauss function and a 1/x-type function. Here's a minimal working example:
restart;
with(plots):
f := x -> exp(-x^2):
g := (x,a) -> 1/(1+abs(x-a)):
plot([f(x),g(x,3)],x=-3..10):
for a from 1 to 100 do
A[a] := int(f(x)*g(x,a),x=-infinity..a-1)
end do:
B:=[seq(A[a],a=1..100)]:
X:=[seq(a,a=1..100)]:
plot(X,B);
The plot it returns looks like this:
which is clearly wrong between around a=35 and a=60...
Does anyone have an idea how this can be repaired?