Приближение функций. Сплайны.

Проинтерполировать функцию Рунге [Graphics:1.gif] на отрезке [Graphics:2.gif] многочленом Ньютона степени [Graphics:3.gif].

Введем функцию [Graphics:4.gif]

[Graphics:5.gif]

Зададим шаг

[Graphics:6.gif]

Вычислим массив узлов и масств значений функции в узлах

[Graphics:7.gif]

[Graphics:8.gif]

Запишем функцию, возвращающую интерполяционный многочлен в форме Ньютона. Здесь [Graphics:9.gif] - векторы исходных данных, [Graphics:10.gif]- размерность этих векторов.

[Graphics:11.gif]

Например, для [Graphics:12.gif] этот многочлен выглядит так

[Graphics:13.gif]

[Graphics:14.gif]

Теперь нарисуем на графике интерполяционный многочлен, исходную функцию и исходные данные

[Graphics:15.gif]

[Graphics:16.gif]

[Graphics:17.gif]

[Graphics:18.gif]

[Graphics:19.gif]

[Graphics:20.gif]

[Graphics:21.gif]

[Graphics:22.gif]

[Graphics:23.gif]

[Graphics:24.gif]

[Graphics:25.gif]

[Graphics:26.gif]

График погрешности

[Graphics:27.gif]

[Graphics:28.gif]

[Graphics:29.gif]

[Graphics:30.gif]