Вычисление интеграла методом средних прямоугольников - 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 }

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

textual
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

12   голосов, оценка 3.917 из 5


СОХРАНИТЬ ССЫЛКУ