Уважаемые форумчане
Вычисляя на Maple интеграл от exp(-x^2)*exp(x) с использованием NAG-овского метода интегрирования _d01amc (интегрирование по бесконечному промежутку)
Код: w :=x-> exp(-x^2);
f := x-> exp(2*x);
evalf(Int(f(x)*w(x), x = -infinity .. infinity, method = _d01amc)); |
Maple выдает ошибку:
Error, (in evalf/int) NE_QUAD_ROUNDOFF_TOL:
Round-off error prevents the requested tolerance from
being achieved: epsabs = 5.0e-013, epsrel = 5.0e-010.
Если вычислять вот так
Код: f := x-> exp(-x^2+2*x);
evalf(Int(f(x), x = -infinity .. infinity, method = _d01amc)); |
то получается правильный результат: 4.818029095 (что есть e*sqrt(\pi))
У меня подозрение, что причина ошибки в кривости использования Мейплом NAG-овского метода интегрирования _d01amc. Если же с Мейплом всё в порядке, то использование NAG-овской библиотеки в любой другой системе должно дать аналогичный результат. Потому я прошу вас вычислить проблемный интеграл на фортране с использованием метода _d01amc из библиотеки NAG и сообщить мне результат. К сожалению у меня на фортране NAG не установлена.
Заранее спасибо.
Спасибо. А вот о авторах cvode - как-то внимания не обратил. CVODE немножко проблематичный тем что не имеет stiffness-детектора а в lsoda он есть. И вот dlsoda (и например radau) уже протестировал успешно на своей задаче а cvode - еще нет. 
