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

unit Gauss_Krist_3;

interface
type fun=function(x:extended):extended;
function gauss_3(a,b:extended; n:integer; f:fun):extended;

implementation
function gauss_3;
const c=0.7745966692;
var x,h,sum:extended;
k:integer;
begin
h:=(b-a)/n;
x:=a+h/2; sum:=0;
for k:=1 to n do
begin
sum:=sum+f(x-h/2*c)+f(x+h/2*c)+1.6*f(x);
x:=x+h;
end;
result:=5*h*sum/18;
end;
end.