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

procedure CountViaBar;
var
  xx1,xx2:real;
  c:longint;
begin
  writeln('------------------------------------------------');
  writeln('metod srednih pryamougolnikov.');
  writeln('vsego iteracij:',round(abs(x2-x1)/e));
  i:=0;
  for c:=1 to round(abs(x2-x1)/e) do begin
    write('iteraciya ',c,chr(13));
    xx1:=Fx(x1+c*e);
    xx2:=Fx(x1+c*e+e);
    i:=i+abs(xx1+xx2)/2*e;
  end;
  writeln('------------------------------------------------');
  writeln('integral=',i);
end;

type
TIntFunc = function(X: Double): Double;
 
function LeftPramInt(a, b: Double; eps: Double; IntF: TIntFunc): Double;
var
//S - площадь на предыдущей итерации,
//step - "толщина" прямоугольника
//gran - передвигаемая от a до b граница
//n - число прямоугольников, удваивается на каждой итерации
S, step, gran: Double;
n: integer;
begin
//Сначала приближение одного прямоугольника
step := b - a;
Result := IntF(a) * step;
n := 1;
repeat
S := Result;
n := n * 2;
step := (b - a) / n;
Gran := a;
Result := 0;
//Ниже - просто вычисляем площади новых прямоугольников
while gran < b do
begin
Result := Result + IntF(gran) * step;
gran := gran + step;
end;
until abs(S - Result) <= eps;
end;
 
function RightPramInt(a, b: Double; eps: Double; IntF: TIntFunc): Double;
var
//S - площадь на предыдущей итерации,
//step - "толщина" прямоугольника
//gran - передвигаемая от a до b граница
//n - число прямоугольников, удваивается на каждой итерации
S, step, gran: Double;
n: integer;
begin
//Сначала приближение одного прямоугольника
step := b - a;
Result := IntF(b) * step;
n := 1;
repeat
S := Result;
n := n * 2;
step := (b - a) / n;
Gran := b;
Result := 0;
//Ниже - просто вычисляем площади новых прямоугольников
while a < gran do
begin
Result := Result + IntF(gran) * step;
gran := gran - step;
end;
until abs(S - Result) <= eps;
end;
 
end.