интерполирование по формулам Ньютона и Лагранжа Опции

[Aleks]

Помогите выполнить работу
или дайте ссылку, если знаете где есть информация по этой теме

Входные параметры:
N:Integer - количество узлов в исходной сетке;
A1,B1:Real - границы сетки интерполирования;
N1:Integer - количество узлов сетки интерполирования ([0..N1]);
X:IntervalType - исходная сетка узлов
( X[1,i] - значения узлов;
X[2,i] - значения функции в этих узлах )
Выходные параметры, возвращаемые процедурой:
Y:IntervalType - сетка интерполяции
( Y[1,i] - значения узлов сетки интерполяции;
Y[2,i] - значения функции в этих узлах;
Y[3,i] - значения первой производной;
Y[4,i] - значения второй производной ) *)

по формуле (см.рис.)