Question: Strange behaviour of the numerical integration using NAG method: method = _d01amc

u := x -> Int(exp(-(xi-x)^2)*exp(xi), xi = -infinity .. infinity, method = _d01amc);
evalf(u(0.));

gives an error

u := x -> Int(exp(-(xi-x)^2)*exp(xi), xi = -infinity .. infinity);
evalf(u(0.));

2.275875794


u := x-> Int(exp(-(xi-x)^2+xi), xi = -infinity .. infinity, method = _d01amc);
evalf(u(0.));

2.275875794


A bug is suspected...

Please Wait...