Курсовая работа: Приближенное вычисление значений определенного интеграла

function ch:real;

var s, dp, kf, a1, b1:real;

begin

s:=0;

kf:=sqrt (1/3);

for i:=2 to n+1 do

begin

a1:=a+h*(i‑2);

b1:=a1+h;

s:=s+((b1‑a1)/2)*(f((a1+b1)/2‑kf*((b1‑a1)/2))+f((a1+b1)/2+kf*((b1‑a1)/2)));

end;

ch:=s;

end;

function si:real;

var s, x, f1, f2:real;

begin

s:=0;

x:=a;

i:=1;

f1:=0;

repeat

f1:=f1+f (a+h*i);

i:=i+2;

until i>=n;

i:=2;

f2:=0;

repeat

f2:=f2+f (a+h*i);

i:=i+2;