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:
f := x -> exp(-x^2):
g := (x,a) -> 1/(1+abs(x-a)):
for a from 1 to 100 do
A[a] := int(f(x)*g(x,a),x=-infinity..a-1)
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?