Аппроксимация / unit 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.