Интегрирование методом прямоугольников с оценкой точности

Собственное значение определенного интеграла

находится методом прямоугольников. Отрезок [a, b] разбивается на n частей x₀ =a, x₁ =a+h, ..., x_n =b с шагом h=(b-a)/n. На каждом отрезке [x_i , x_i+1 ] вычисляется значения интеграла по формуле прямоугольников, таким образом:

Затем количество точек разбиения удваивается и производится оценка точности вычислений:

Если R_n  > 3*e, то количество точек разбиения удваивается.

Если нашли ошибку в алгоритме - сообщите!

Реализации алгоритма на различных языках:

Блоксхемы:

Реализация алгоритма на AlgoPascal:

unit IntegralRectUnit;

declaration
(*
интегрируемая функция
*)
function F(X:Real):Real;

interface
    IntegralRect;
implementation

(*********************************************************
Интегрирование методом прямоугольников с оценкой точности.

Считается интеграл функции F на отрезке [a,b] с погрешностью
порядка Epsilon.

function IntegralRect(a:Real;b:Real;Epsilon:real):real;
*********************************************************)
function IntegralRect(
    const   a   :   Real;
    const   b   :   Real;
    const   Epsilon:Real):Real;
var
    i   :   Integer;
    n   :   Integer;
    h   :   Real;
    s1  :   Real;
    s2  :   Real;
begin
    n:=1;
    h:=b-a;
    s2:=h*F((a+b)/2);
    repeat
        n:=2*n;
        s1:=s2;
        h:=h/2;
        s2:=0;
        i:=1;
        repeat
            s2:=s2+F(a+h/2+h*(i-1));
            i:=i+1;
        until not(i<=n);
        s2:=s2*h;
    until not(AbsReal(s2-s1)>3*Epsilon);
    Result:=s2;
end;

end.

Бочканов Сергей, Быстрицкий Владимир
Copyright © 1999-2004
При поддержке проекта MANUAL.RU