Function Intg(p, b: real; n: integer; eps: real): real; {Метод прямоугольников, рекурсия} var a, x, h, s: real; i: integer; begin a := 0; h := (b - a) / n; x := a; s := 0; for i := 1 to n do begin s := s + sin(x); x := x + h; end; if abs(s * h - p) > eps then Intg := Intg(s * h, b, 2 * n, eps) else Intg := s * h; end;
var b, eps: real;
begin write('Задайте верхний предел интегрирования и точность: '); Readln(b, eps); writeln('Значение интеграла: ', Intg(0, b, 2, eps):0:7); end.
Тестовое решение:
Задайте верхний предел интегрирования и точность: 2 0.000001 Значение интеграла: 1.4161460
Точное значение интеграла равно 1-cos(2) ≈1.416146836, т.е. заданная точность обеспечивается.
var x: array[1..m, 1..n] of integer; i, j, a, b, k: integer;
begin write('Введите числа a,b: '); readln(a, b); Randomize; writeln(#13#10, 'Исходный массив'); for i := 1 to m do begin for j := 1 to n do begin x[i, j] := Random(100) - 50; write(x[i, j]:5) end; writeln end; writeln('Количество элементов матрицы, кратных ', a); for j := 1 to n do begin k := 0; for i := 1 to m do if x[i, j] mod a = 0 then k := k + 1; write(k:5) end; writeln; writeln('Количество элементов матрицы, кратных ', b); for j := 1 to n do begin k := 0; for i := 1 to m do if x[i, j] mod b = 0 then k := k + 1; write(k:5) end; writeln end.
{Метод прямоугольников, рекурсия}
var
a, x, h, s: real;
i: integer;
begin
a := 0;
h := (b - a) / n;
x := a;
s := 0;
for i := 1 to n do
begin
s := s + sin(x);
x := x + h;
end;
if abs(s * h - p) > eps then
Intg := Intg(s * h, b, 2 * n, eps)
else
Intg := s * h;
end;
var
b, eps: real;
begin
write('Задайте верхний предел интегрирования и точность: ');
Readln(b, eps);
writeln('Значение интеграла: ', Intg(0, b, 2, eps):0:7);
end.
Тестовое решение:
Задайте верхний предел интегрирования и точность: 2 0.000001
Значение интеграла: 1.4161460
Точное значение интеграла равно 1-cos(2) ≈1.416146836, т.е. заданная точность обеспечивается.
m = 5;
n = 8;
var
x: array[1..m, 1..n] of integer;
i, j, a, b, k: integer;
begin
write('Введите числа a,b: ');
readln(a, b);
Randomize;
writeln(#13#10, 'Исходный массив');
for i := 1 to m do
begin
for j := 1 to n do
begin
x[i, j] := Random(100) - 50;
write(x[i, j]:5)
end;
writeln
end;
writeln('Количество элементов матрицы, кратных ', a);
for j := 1 to n do
begin
k := 0;
for i := 1 to m do
if x[i, j] mod a = 0 then k := k + 1;
write(k:5)
end;
writeln;
writeln('Количество элементов матрицы, кратных ', b);
for j := 1 to n do
begin
k := 0;
for i := 1 to m do
if x[i, j] mod b = 0 then k := k + 1;
write(k:5)
end;
writeln
end.
Тестовое решение:
Введите числа a,b: 4 7
Исходный массив
-16 3 26 4 46 -6 -31 21
-35 -34 46 42 15 -26 -9 4
-19 -47 8 -48 13 14 42 -46
41 -10 13 45 -20 15 13 -44
18 -11 23 -30 44 -36 -8 -42
Количество элементов матрицы, кратных 4
1 0 1 2 2 1 1 2
Количество элементов матрицы, кратных 7
1 0 0 1 0 1 1 2