Вычисление значение интеграла по формуле трапеций - Pascal ABC

Формулировка задачи:

Напишите пожалуйста программу Вычислить значение интеграла по формуле трапеций где h=(b-a)/n; x_k=a+kh; f(x_k) - подынтегральная функция. Принять a=1, b=4, n=40.

Решение задачи: «Вычисление значение интеграла по формуле трапеций»

Листинг программы

function F(x:real):real;
begin
 F:=power(x,3)/(power(x,4)+1);
end;
const a=1;
      b=4;
      n=40;
var h,k,s:real;
    i:integer;
begin
h:=(b-a)/n;
s:=(F(a)+F(b))/2;
k:=a;
for i:=1 to n-1 do
 begin
  k:=k+h;
  s:=s+F(k);
 end;
s:=s*h;
write('S=',s:0:4);
end.

Объяснение кода листинга программы

В функции F(x) вычисляется значение интеграла по формуле трапеций для переменной x.
Константы a и b определяют пределы интегрирования.
Переменная n определяет количество точек, через которые будет проходить интегрирование.
Переменные h, k и s используются для промежуточных вычислений.
Переменная i используется для цикла, который проходит n-1 раз.
Внутри цикла переменная k увеличивается на h, а переменная s вычисляется как среднее значение F(k) и F(b).
После завершения цикла значение переменной s умножается на h.
Результат выводится на экран.