Численное интегрирование (по формуле прямоугольника, трапеции, симпсона)

uses crt;
 
function f(x:real):real;
begin
f := sqr(x)/((4*x - 3)*sqrt(4*x - 3));
end;
 
var
  Int : real;
  a, b, n : integer;
  x, h : real;
begin
Write('Нижний предел: ');
Readln(a);
Write('Верхний предел: ');
Readln(b);
Write('Количество узлов: ');
Readln(n);
h := abs(b - a)/n; 
 
//Формула левых прямоугольников
x := a;
Int := 0;
while (x < b) do
  begin
  Int := Int + f(x)*h;
  x := x + h;
  end;
Writeln('Формула прямоугольников с избытком: ', Int:0:4);
 
//Формула правых прямоугольников
x := a + h;
Int := 0;
while (x <= b) do
  begin
  Int := Int + f(x)*h;
  x := x + h;
  end;
Writeln('Формула прямоугольников с недостатком: ', Int:0:4);
 
//Формула трапеции
x := a + h; Int := 0;
while (x < b)do
  begin
  Int := Int + 1/2*h*(f(x - h) + f(x));
  x := x + h;
  end;
Writeln('Формула трапеции: ', Int:0:4); 
 
//Формула Симпсона 
x := a ; Int := 0;
while (x < b) do
  begin
  Int := Int + h*(f(x) + 4*f(x + h/2) + f(x + h))/6;
  x := x + h;
  end;
Writeln('Формула Симпсона: ', Int:0:4);  
end.

function f(x:real):real;
begin
 if 4*x + 1 < 0 then f := 0 else
  f := exp(1/4*ln(4*x + 1));
end;