Интерполяция табл данных через полином Ньютона

procedure Enter(var X,y,n:real); 
var
 i: integer;
label mp;
begin
mp: for i:=0 to n do
 begin
 write('X[',i,'] = '); readln(x[i]);
 write('y[',i,'] = '); readln(y[i]);
 end;
 for i:=0 to n-1 do
 if x[i+1]-x[i]<=0 then
 begin
 writeln ('Ошибка. Повторите ввод.');
 goto mp
 end;
 end;
procedure Print(n: integer; X,y: real); 
var i: integer;
 begin
 for i:=0 to n do
 begin
 write(x[i]:12:6);
 end;
 writeln;
 for i:=0 to n do
 begin
 write(y[i]:12:6);
 end;
 writeln;
end;
Function Polinom(n: integer; d,x,y:real;);
var
l:real;
k,i:integer;
p: real;
begin
L:=y[0];
P:=1;
for k:=1 to n do begin
P:=P*(D-X[k-1]);
for i:=0 to (n-k) do begin
Y[i]:=(y[i+1]-y[i])/(x[i+k]-x[i]);
end;
L:=L+P*y[0];
end;
POlinom:=l;
end;
begin 
writeln('Введите количество узлов n (N=n+1)');
 read(n);
 Enter(X,y);
 Print(n,X,y);
 repeat
 lt:Writeln('Введите X (ot ',x[0]:4:2,' do ',x[n]:4:2,')');
 read(d);
 if d<x[0] then begin
 writeln('Ошибка. x не может быть меньше ',x[0]:4:2);
 goto lt; end;
 if d>x[n] then begin
 writeln('Ошибка. x не может быть больше ',x[n]:4:2);
 goto lt; end;
 writeln(Polinom (n,d,X,y):6:3);
end.