Задания для самостоятельного выполнения.

с заданной точностью ε.

В программе должны быть следующие подпрограммы:

1. Процедура Proc, вычисляющая интеграл с заданной точностью. Входные параметры: границы интервала [a, b], точность. Выходные: значение интеграла, количество разбиений n Начальное количество разбиений – 100, шаг по количеству разбиений – 100. Вызывается из главной программы

2. Функция Intg, вычисляющая интеграл при фиксированном количестве разбиений заданным методом. Входные параметры для нее: границы интервала [a, b], количество разбиений. Результат: значение интеграла. Вызывается из процедуры Proc.

3. Функции , вычисляющая Вызывается в функции Intg.

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

Пример 1.

Подинтергральная функция:

Интервал [1; 4]. Метод интегрирования – прямоугольников.

,

где .

Program pr6;

uses crt;

var a, b, eps, I2: real; N: longint;

procedure proc(eps, a, b: real; var I2: real; var N: longint);

var I1: real;

function Intg(a, b: real; N: longint): real;

var x, s, dx: real; i: longint;

function f1(x: real): real;

begin

f1:= (1 + sqrt(x))/sqr(x)

end;

begin

s:= f1(a); x:= a;

dx:= (b-a)/N;

for i:=1 to N do begin

x:= x + dx;

s:= s + f1(x)

end;

Intg:= dx*s

end;

begin

I1:= 0;

N:= 100;

repeat

I2:= Intg(a, b, N);

if abs(I1-I2) < eps then break

else begin

I1:= I2; N:= N + 100

end

until false

end;

begin clrscr;

write('Введите a, b, eps -> ');

readln(a, b, eps);

proc(eps, a, b, I2, N);

writeln(' Результаты : ');

writeln('I2 = ', I2:12:6);

writeln('Кол-во итераций = ', N:10);

readkey

end.