Приближенное значение в точке с использованием сплайн-функции

Program program7;
Const
  N = 16; {число переменных в таблице}
Var
  x, y : Array [1..N] Of Real;
 
 
Function f(x : Real) : Real;
Var
  A, B, C : Real;
  K, N : Integer;
Begin
  N := 7;
  K := N + 10;
  A := 0.35 + 0.05 * Sin(N);
  B := 0.45 + 0.05 * Sin(N);
  C := 0.50 + 0.30 * Sin(N);
  f := A * Sqr(x) + B * x + C + 0.3 * Sin(K*(x + 1));
End;
  
Procedure CreateTable;
Var
  i : Integer;
  xx : Real;
Begin
  WriteLn('====================================================================');
  WriteLn('Таблица экперементальных данных');
  WriteLn('    x    |    y    ');
  xx := 0.0;
  For i := 1 To N Do
  Begin
    x[i] := xx;
    y[i] := f(xx);
    WriteLn(x[i]:10:5, ' ', y[i]:10:5);
    xx := xx + 0.2;
  End;
End;
 
Procedure Langrange(LnX : Real);   {LnX это х для которого нужно найти L}
Var
  K, i : Integer;
  L, Ln : Real;
Begin
  WriteLn('====================================================================');
  WriteLn('Интерполяционный многочлен Лангранжа');
  WriteLn('X -> ', LnX);
 
  L := 0;
  For i := 1 To N Do
  Begin
    Ln := y[i];
    For k := 1 To N Do
      If i <> k Then
        Ln := Ln*((Lnx - x[k])/(x[i] - x[k]));  {Рекурентная формула для вычисления}
    L := Ln + L;
  End;
  WriteLn('Корень -> ',L:4:10);
End;
Begin
  CreateTable;
  Langrange(0.55);
  Langrange(1.33);
  Langrange(2.75);
End.