Интерполяция по схеме Эйткена

[Gonzo]

По известным в узлах x0,x1,..,x5 значениям функции F(x) построить интерполяционный многочлен по схеме Эйткена. Вычислить значения этого многочлена в равномерной сетке точек на отрезке [a,b] с шагом h и сравнить их со значениями функции F(x) в этих точках.
F(x)=x*x*ln(x)
x0=10, x1=10.5, x2=11, x3=12, x4=12.5, x5=13
a=11, b=15, h=0.1

Есть такая процедурка:

PROCEDURE AITKEN(X:REAL;Y:MAS1;N:INTEGER;XP:REAL; VAR L : REAL);
VAR I,J,N1: INTEGER;
BEGIN 	
	N1:=N-1; 
	FOR J:=0 TO N1 DO 
		FOR I:=J+1 TO N DO
		Y[I]:=((XP-X[J])*Y[I]-(XP-X[I])*Y[J])/(X[I]-X[J]); 	L:=Y[N];
END.

Входные: x[0:n], y[0:n] (тип real) - массивы абсцисс и ординат таблично заданной функции y = [yi(xi)]; n (тип integer) - количество значений табличной функции (количество узлов интерполяции); xp (тип real) - абсцисса, для которой вычисляется интерполированноe значение функции. Выходные: l (тип real) - интерполированное значение функции y (значение полинома Лагранжа в точке xp).

Подскажите что мне еще нужно. Буду признателен за любую помощь.

[Gonzo]

Разобрался сам. Тему можно закрывать.

[anytka29]

кинь код программы пожалуйста)

[ZX Spectrum-128]

Сообщение от anytka29

кинь код программы пожалуйста)

Спустя 5 лет?

[Gonzo]

Так выше же есть исходник и описание. Полного кода с использованием данной подпрограммы и не найду сейчас.