* Данный текст распознан в автоматическом режиме, поэтому может содержать ошибки
490
Приближение функций и прогноз
6.4. Выбор аппроксимации для сложной функции
6.4.1. Задание исходной функции и построение ее графика
Покажем возможности аппроксимации сложной функции средствами систем Maple на одном из комплексных примеров, давно помещенных в библиотеку пользователей системы Maple V R2, и переработанном для Maple 8 и более позд них версий. Воспользуемся для этого ранее описанными возможностями пакета numapprox, для чего прежде всего подключим его:
> restart:with(numapprox):
Будем искать приемлемую аппроксимацию для следующей, отнюдь не прос той тестовой функции:
> f := x -> int(1/GAMMA(t), t=0..x ) / x^2;
> plot(f,0..4,color=black);
График этой функции представлен на рис. 6.20. С первого взгляда это простой график, но тут как раз тот случай, когда простота обманчива. Вы сразу заметите, что график строится медленно, поскольку в каждой из множества его точек систе ме Maple приходится вычислять значение интеграла с подынтегральной функци ей, содержащей гамма функцию. Maple делает это по сложному и медленному ал горитму адаптивного численного интегрирования. Итак, вычисление f(x) по ее интегральному представлению совершенно не эф фективно. Наша цель состоит в разработке процедуры вычислений, которая дала
Рис. 6.20. График аппроксимируемой функции