Интеграл, решенный методом симпсона и прямоугольника, растолковать код

program integral_2;
var
h,a,b,s,s1,x,int:real;
i,m:integer;
 
function f(x:real):real;
begin
f := 1/x;
end;
 
begin
a:=0.5;
b:=1;
read(m);
h:=(b-a)/(2*m);
s:=0;
s1:=0;
x:=a+h;
for i:=1 to (2*m-1) do
begin
if (i mod 2=0) then
s:=s+f(x)
else
s1:=s1+f(x);
x:=x+h;
end;
writeln(s);
writeln(s1);
writeln(f(a));
writeln(f(b));
int:=(h*(f(a)+f(b)+2*s+4*s1))/3;
writeln(int);
 
end.

program integral;
var
s, a, b, h, x: real;
i, n:integer;
 
function f(x:real):real;
begin
f :=1/power(8-x*x*x, 0.2);
end;
 
begin
a:=0;
b:=1;
writeln ('Vvedite chislo razbienii');
readln(n);
writeln;
h:=(b-a)/n;
s:=0;
x:=a+h/2;
for i:=1 to n do
begin
s:=s+f(x);
x:=x+h;
end;
s:=s*h;
writeln('Integral raven');
writeln(s:7:6);
 
end.

program integral_2;
uses
    crt;    // стандартный модуль crt
var
   h,a,b,s,s1,x,int:real;
   i,m:integer;
 
function f(x:real):real;
begin
 f := 1/x;  // нахождение f
end;
 
begin
 a:=0.5; // присвоение значения
 b:=1;  // присвоение значения
 write('m='); // вывод на экран
 readln(m);
 h:=(b-a)/(2*m); //нахождение h по формуле
 s:=0;  // присвоение значения
 s1:=0; // присвоение значения
 x:=a+h;  //нахождение x
 for i:=1 to (2*m-1) do  // цикл от 1 до 2*m-1
 begin
  if (i mod 2=0) then //если i mod 2=0, то вычисляем по 1 формуле, а если нет, то по 2.
   s:=s+f(x)    // 1
  else
   s1:=s1+f(x);    // 2
   x:=x+h;    //нахождение x
 end;
 writeln(s);  // вывод на экран
 writeln(s1);   // вывод на экран
 writeln(f(a));  // вывод на экран
 writeln(f(b)); // вывод на экран
 int:=(h*(f(a)+f(b)+2*s+4*s1))/3;   //нахождение int по формуле
 writeln(int);  // вывод на экран
 
end.

program integral;
uses
    crt;
var
   s, a, b, h, x: real;
   i, n:integer;
 
function f(x:real):real;
begin
 f :=1/power(8-x*x*x, 0.2);  // нахождение f по формуле
end;
 
begin
 a:=0;
 b:=1;
 writeln ('Vvedite chislo razbienii');
 readln(n);
 writeln;
 h:=(b-a)/n;
 s:=0;
 x:=a+h/2;
 for i:=1 to n do  //цикл от 1 до n
 begin
  s:=s+f(x);
  x:=x+h;
 end;
 s:=s*h;
 writeln('Integral raven');
 writeln(s:7:6);
 
end.