Уплотнить часть таблицы функции (Delphi)

[FireHawK]

Привет. Есть такое задание: уплотнить часть таблицы заданной функции, пользуясь интерполяционными формулами Ньютона.
Так вот, воспользовался я формулой (в приложении она).

Для нее надо найти конечные разности. нашел.

 
    X              Y                 (разности --->   )
    1,1       0,89121       0,02155       -0,00227       -0,00007       0,00005       -0,00009       0,00018       -0,00035       0,00066       -0,00125       0,00259
       1,15       0,91276       0,01928       -0,00234       -0,00002       -0,00004       0,00009       -0,00017       0,00031       -0,00059       0,00134       0
       1,2       0,93204       0,01694       -0,00236       -0,00006       0,00005       -0,00008       0,00014       -0,00028       0,00075       0       0
       1,25       0,94898       0,01458       -0,00242       -0,00001       -0,00003       0,00006       -0,00014       0,00047       0       0       0
       1,3       0,96356       0,01216       -0,00243       -0,00004       0,00003       -0,00008       0,00033       0       0       0       0
       1,35       0,97572       0,00973       -0,00247       -0,00001       -0,00005       0,00025       0       0       0       0       0
       1,4       0,98545       0,00726       -0,00248       -0,00006       0,0002       0       0       0       0       0       0
       1,45       0,99271       0,00478       -0,00254       0,00014       0       0       0       0       0       0       0
       1,5       0,99749       0,00224       -0,0024       0       0       0       0       0       0       0       0
       1,55       0,99973       -0,00016       0       0       0       0       0       0       0       0       0
       1,6       0,99957       0       0       0       0       0       0       0       0       0       0

мне нужно уплотнить часть таблицы в промежутке a=1.45 b=1.55, h=0.01.
ввел формулу в цикле. t как я понял будет равно просто -10 -9 -8.....
И в итоге выдает мне тотже самый столбец значений Y, которые даны(см. выше). Но я как понял должны вычислится значения Y от 0,99271 до 0,99973. Но судя по формуле он никак не зависит от X . Какже с помощью нее уплотнять?
ЗЫ. в коде f(2).. это факториал.
mas[0..10,0..11] массив с значениями

For i:=0 to 9 do
begin
 t:=((1.45+i*0.01)-1.55)/0.01;
 P:=mas[10,1]+t*mas[9,2]+mas[8,3]*(t*(t+1))/f(2)+mas[7,4]*(t*(t+1)*(t+2))/f(3)+
 mas[6,5]*(t*(t+1)*(t+2)*(t+3))/f(4) + mas[5,6]*(t*(t+1)*(t+2)*(t+3)*(t+4))/f(5)+
  mas[4,7]*(t*(t+1)*(t+2)*(t+3)*(t+4)*(t+5))/f(6)+ mas[3,8]*(t*(t+1)*(t+2)*(t+3)*(t+4)*(t+5)*(t+6))/f(7) +
  mas[2,9]*(t*(t+1)*(t+2)*(t+3)*(t+4)*(t+5)*(t+6)*(t+7))/f(8) +
  mas[1,10]*(t*(t+1)*(t+2)*(t+3)*(t+4)*(t+5)*(t+6)*(t+7)*(t+8))/f(9) +
  mas[0,11]*(t*(t+1)*(t+2)*(t+3)*(t+4)*(t+5)*(t+6)*(t+7)*(t+8)*(t+9))/f(10);
Memo1.Lines.Add(FloatToStr(P));
end;