Вычисление определенного интеграла методом центральных прямоугольников, и методом трапеции

const
a=1;
b=3;
n=76;
 
function f(x : real) : real;
begin
  f :=1/((power(x,0.5))*(3+exp(0.5*x)));
end;
 
var
  
   i : integer;
  h, q : real;
 
begin
 
h := (b - a)/n; 
q := 0;
for i := 1 to n do
 q := q + f(a + h* (i - 0.5));
q := h*q;
Writeln('Ответ: ', q:0:4); 
end.

const a=1;
      b=3;
      eps=0.0001;
function F(x:real):real;
begin
F:=1/(sqrt(x)*(3+exp(0.5*x)));
end;
//метод средних прямоугольников
function sr_prm(n:integer):real; {Вычисление интеграла}
var sum,x,h : real;
    i:integer;
begin
h:=(b-a)/n;
sum:=0;
x:=a+h/2;
for i:=1 to n do
 begin
  sum:=sum+f(x);
  x:=x+h;
 end;
Sr_prm:=sum*h;
end;
//метод трапеций
function Trap(n:integer):real;
var h,sum:real;
    i:integer;
begin
h:=(b-a)/n;
sum:=(f(a)+f(b))/2;
for i:=1 to n-1 do
sum:=sum+f(a+i*h);
sum:=sum*h;
Trap:=sum;
end;
var x1,x2:real;
    si,sn:real;
    n:integer;
begin
n:=1;
sn:=0; //инициируем первое приближение для интеграла
si:=Sr_prm(n); //Второе приближение
while abs(si-sn)>eps do //вычисляем более точно
 begin
  sn:=si;
  n:=2*n;
  si:=Sr_prm(n);
 end;
writeln('Интеграл методом средних прямоугольников=', si:0:4);
n:=1;
sn:=0; //инициируем первое приближение для интеграла
si:=Trap(n); //Второе приближение
while abs(si-sn)>eps do //вычисляем более точно
 begin
  sn:=si;
  n:=2*n;
  si:=Trap(n);
 end;
writeln('Интеграл методом трапеций=', si:0:4);
end.