Площадь криволинейной трапеции приближенно заменяют площадью многоугольника, составленного из n прямоугольников.

Сделаем вывод: вычисление определенного интеграла сводится к нахождению суммы n элементарных прямоугольников, где h=(b-a)/n –ширина прямоугольников.

Формула средних прямоугольников: S _средих = (S _правых + S _левых ) /2

Создадим программу решения интегралов методом левых прямоугольников и результаты выведем на экран, в результате чего получим следующее…

Program levii;{Метод левых прямоугольников}
uses crt;
var i,n:integer; a,b,h,x,xb,s:real;
function f(x:real):real;
begin f:=(1/x)*sin(3.14*x/2); end;
begin
clrscr;
write('Введите нижний предел интегрирования '); readln(a);
write('Введите верхний предел интегрирования '); readln(b);
write('Введите количество отрезков '); readln(n);
h:=(b-a)/n; s:=0; xb:=a;
for i:=0 to n-1 do
begin x:=xb+i*h; s:=s+f(x)*h; end;
writeln('Интеграл равен ',s:12:10); readln;
end.

a=1 b=2 n=10 S= 18,077

a=1 b=2 n=20 S= 18, 208

a=1 b=2 n=100 S= 18, 270

Создадим программу решения интегралов методом средних прямоугольников и результаты выведем на экран, в результате чего получим следующее…

Program srednii; {Метод средних прямоугольников}
uses crt;
var i, n: integer; a, b, dx, x, s, xb : real;
function f(x : real):real;
begin f:=(1/x)*sin(3.14*x/2); end;
begin
clrscr;
write('Введите нижний предел интегрирования '); readln(a);
write('Введите верхний предел интегрирования '); readln(b);
write('Введите количество отрезков '); readln(n);
dx:=(b-a)/n; xb:=a+dx/2;
for i:=0 to n-1 do
begin x:=xb+i*dx; s:=s+f(x)*dx; end;
write('Интеграл равен ',s:15:10); readln;
end.

a=1 b=2 n=10 S=18,07667

a=1 b=2 n=20 S=18,368

a=1 b=2 n=100 S= 18,156

Из рассмотренных выше примеров легко заметить, что при вычислении определенных интегралов методами прямоугольников мы не можем достигнуть точного значения, т.е. чем больше значение n, тем точнее значение интеграла.

Поделитесь этой записью или добавьте в закладки