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