Аппроксимация / Polinom_Hytona

unit Polinom_Hytona;

interface
type znach=array[1..50] of extended;
function Hyton(t:extended; x,y:znach; n:integer):extended;

implementation
function Hyton;
var i,k:integer;
p,z:extended;
rr:znach;
begin
for i:=1 to n do rr[i]:=y[i];
p:=1;
z:=y[1];
for k:=1 to n-1 do
begin
for i:=1 to n-k do
rr[i]:=(rr[i]-rr[i+1])/(x[i]-x[i+k]);
p:=p*(t-x[k]);
z:=z+p*rr[1];
end;
result:=z;
end;
end.