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?