Вычислить интеграл - Pascal ABC

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

Заранее большое спасибо!

Решение задачи: «Вычислить интеграл»

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

function f(x:real):real;
begin
f:=sin(x);
end;
var a,b,h,s:real;
    n,i:integer;
begin
repeat
writeln('Введите границы интервала b>a');
readln(a,b);
until b>a;
repeat
write('Введите число разбиений от 20 до 1000 n=');
readln(n);
until (n>=20)and(n<=1000);
h:=(b-a)/n;
s:=(f(a)+f(b))/2;
for i:=1 to n-1 do
s:=s+f(a+i*h);
s:=s*h;
write('S=',s:0:5);
end.

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

Введенная функция f(x) вычисляет синус аргумента x и возвращает результат в виде числа с плавающей точкой.
Переменные a, b, h, s и n объявлены как числа с плавающей точкой.
Переменная i объявлена как целочисленный индекс.
Пользователю предлагается ввести границы интервала, в рамках которого будет вычисляться интеграл.
Пользователю предлагается ввести количество разбиений интервала от 20 до 1000.
Переменная h инициализируется как (b-a)/n.
Переменная s инициализируется как (f(a)+f(b))/2.
В цикле от 1 до n-1 происходит вычисление значений функции f(x) для всех точек x, лежащих внутри интервала, и добавление этих значений к переменной s.
Значение переменной s умножается на h.
Результат вычисления интеграла выводится на экран с точностью до 5 знаков после запятой.