[an error occurred while processing this directive]
|
Для вычисления таких интегралов по большому счету годятся почти все известные методы. В самом деле: аналитическая функция (единственная проблема в нуле решается); монотонно убывает и имеет ассимптоту равную нулю... Тут и метод Эйлера сгодится или трапеций, если не нужна высокая точность. Если же надо поточнее, то можно и методами Рунге-Кутты воспользоваться - либо в точности выиграете при том же количестве шагов, либо в скорости при сохранении той же точности.
Полиномами я бы не стал аппроксимировать ввиду наличия горизонтальной ассимптоты, которая под концепцию полиномов не очень подходит.
Добавлю, что вот над чем бы я еще подумал ради ускорения вычислений (если таковое ускорение требуется) - это принудительно ограничил бы интервал интегрирования значением, полученным на основе анализа функции ошибки усечения:
E(b)=инт(f(x), x=от a до беск) - инт(f(x), x= от а до b)
В этом случае нужно получить заранее хорошую аппроксимацию функции ошибки.
E-mail: info@telesys.ru