Аппроксимация табличной функции с помощью интерполяционных полиномов Лагранжа (доработать программу)

Uses crt;
Var n,k,m,i,j,mi,ii : integer;
    ct,u,zi,zn:real;
    xn: array [1..10] of real;
    x: array [1..20] of real;
    y: array [1..20] of real;
    r,pt:array[1..20] of real;
begin
y[1]:= -1.84; y[11]:=-1.33;
y[2]:=-1.98;  y[12]:=-1.47;
y[3]:=-1.72;  y[13]:=-1.50;
y[4]:=-1.58;  y[14]:=-2.65;
y[5]:=-1.69;  y[15]:=-1.65;
y[6]:=-1.59;  y[16]:=-1.87;
y[7]:=-1.58;  y[17]:=-1.61;
y[8]:=-1.64;  y[18]:=-1.86;
y[9]:=-1.55;  y[19]:=-1.84;
y[10]:=-1.35; y[20]:=-1.91;
n:=20;
writeln('Введите массив узлов интерполирования:');
for i:=1 to n do
 begin
  x[i]:=i/10;
 end;
clrscr;
for i:=1 to n do R[i]:=0;
for i:=1 to n do
 begin
  ct:=1;pt[1]:=1;k:=1;
  for j:=1 to n do
  if i<>j then
   begin
    ct:=ct*(x[i]-x[j]);
    pt[k+1]:=1;
    if k>1 then
    for m:=k downto 2 do
    pt[m]:=pt[m-1]-pt[m]*x[j];
    pt[1]:=-pt[1]*x[j];
    k:=k+1;
   end;
  u:=y[i]/ct;
  for m:=1 to n do
  r[m]:=r[m]+pt[m]*u;
 end;
writeln('Массив коэффициентов:');
for i:=1 to n do begin
writeln('R[', i, ']=', r[i]:12:2,' ');
end;
readln;
clrscr;
 
{for ii:=1 to 10 do
for i:=1 to n do
xn[ii]:=x[i]+x[i]/2;
for mi:=n-1 downto 0 do begin
zn:=zn+r[i]*(exp(ln(xn[ii])*mi));
writeln(zn);
end;} в чем тут ошибка?
 
end.