Вычисление интеграла методом средних прямоугольников - QBasic

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

{ Вычисление интеграла методом средних прямоугольников }

var
  n : integer;
  s: real;
 
function f(x:real):real;
begin
  f:=1-0.25*sin(x)*sin(x);
end;
 
procedure rect (a,b:real; n:integer; var s:real);
var
  i : integer;
  h,x: real;
begin
  h:=(b-a)/n;
  x:=a+h/2;
  s:=0.0;
  for i:=1 to n do
    begin
      s:=s+f(x);
      x:=x+h;
    end;
  s:=s*h;
end;
 
begin
  rect (0,pi/2,40, s);
  write('Интеграл = ',s:16:4);
end.

{ Ответ: } { Интеграл = 1.3744 }

Решение задачи: «Вычисление интеграла методом средних прямоугольников»

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

a = 0
b = 3.14159 / 2
n = 40
h = (b - a) / n
s = 0
FOR x = a + h / 2 TO b STEP h
s = s + (1 - .25 * SIN(x) ^ 2) * h
NEXT
PRINT s

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

В данном коде на языке QBasic выполняется вычисление интеграла методом средних прямоугольников. Список действий, выполняемых в коде, можно оформить в виде следующих элементов:

Установка начального значения переменной a равным 0.
Установка конечного значения переменной b равным 3.14159 / 2.
Установка количества прямоугольников (переменная n) равным 40.
Вычисление шагового значения переменной h как (b - a) / n.
Инициализация переменной s равной 0.
Выполнение цикла от a + h / 2 до b с шагом h.
Внутри цикла вычисление значения переменной s как сумма (1 - .25 SIN(x) ^ 2) h.
После выполнения цикла выводится значение переменной s.