Turbo Pascal (стр. 3 из 4)

UnitKursUn; {*** Начало модуля KursUn ***}

interface {*** Интерфейсная секция ***}

usescrt;

type

matrix= array [1..10,1..10] of real;

vector= array [1..10] of real;

var

i,j,k:integer;

n:word;

procedure outputmat (n:word; a:matrix; m:char);

procedure inputm (n:word; var a:matrix;m:char);

procedure sred_arifm (n:word; a:matrix;var z:vector);

procedure transpm (n:word; a:matrix; var c:matrix);

procedure sum_m (n:word; a,b:matrix; var c:matrix);

procedure multm (n:word; a,b:matrix; var c:matrix);

procedure multconstm (n:word; r:real; a:matrix;var c:matrix);

implementation {*** ИСПОЛНЯЕМАЯЧАСТЬ ***}

{***************************************************************************}

{*** процедуравводаматриц ***}

procedure inputm;

begin

clrscr;

writeln;

writeln(' Введите матрицу ',m,' размером ',n,'*',n);

for i:=1 to n do

for j:=1 to n do

begin

write(' ',m,'[',i,',',j,']=');

readln(a[i,j]);

end;

end;

{***************************************************************************}

{*** процедуравыводаматриц ***}

procedure outputmat;

begin

writeln;

writeln(' Матрица ',m,'.');

writeln;

for i:=1 to n do

begin

write(' ');

for j:=1 to n do

write(' ',a[i,j]:3:1);

writeln;

end;

end;

{***************************************************************************}

{*** процедура транспонирования матрицы ***}

procedure transpm;

begin

for i:=1 to n do

for j:=1 to n do

c[j,i]:=a[i,j];

end;

{***************************************************************************}

{*** процедура умножения матрицы на число ***}

procedure multconstm;

begin

for i:=1 to n do

for j:=1 to n do

c[i,j]:=a[i,j]*r

end;

{***************************************************************************}

{*** процедура суммирования матриц ***}

proceduresum_m;

begin

for i:=1 to n do

for j:=1 to n do

c[i,j]:=a[i,j]+b[i,j];

end;

{***************************************************************************}

{*** процедураумноженияматриц ***}

procedure multm;

begin

for i:=1 to n do

for j:=1 to n do

begin

c[i,j]:=0;

for k:=1 to n do

c[i,j]:=c[i,j]+a[i,k]*b[k,j];

end;

end;

{***************************************************************************}

{*** процедура формирования вектора из средних ***}

{*** арифметических значений элементов столбцов ***}

proceduresred_arifm;

var

S:real;

begin

S:=0;

for i:=1 to n do

begin

for j:=1 to n do

S:=S+a[j,i];

z[i]:=S/n;

S:=0;

end;

end;

{***************************************************************************}

end. {*** Конец модуля KursUn ***}

Основная программа.

Program Kursach1;

Uses KursUn , Crt;

Var

a,b,c,u : matrix;

vect : vector;

begin

ClrScr; textcolor(LightCyan);

writeln;

writeln(' ╔═══════════════════════════════════════════════════════════════╗');

writeln(' ║ Этапрограммавычисляетматричноевыражение║');

writeln(' ║ ║');

writeln(' ║ T ║');

writeln(' ║ U=C*( A+2*B ) ║');

writeln(' ║ ║');

writeln(' ╚═══════════════════════════════════════════════════════════════╝');

writeln;

write(' Введите размерности матриц: '); readln(n);

if n=0 then {*** проверка размерности матрицы ***}

begin

ClrScr; textcolor(red);

writeln;

writeln(' Такая размерность не допустима!!!');

readkey;

exit;

end;

ClrScr;

inputm(n,a,'A'); {*** вводматрицы A ***}

ClrScr;

inputm(n,b,'B'); {*** вводматрицы B ***}

ClrScr;

inputm(n,c,'C'); {*** вводматрицы C ***}

transpm(n,b,u); {*** транспонирование матрицы B. ***}

multconstm(n,2,u,u); {*** умножения матрицы на 2. ***}

sum_m(n,a,u,u); {*** суммирование матриц A+2*BT. ***} multm(n,c,u,u); {*** умножение матриц С и (A+2*BT). ***}

ClrScr;

writeln;

writeln(' ****************** Исходные значения ********************');

outputmat(n, a, 'A'); {*** вывод матрицы A***}

outputmat(n, b, 'B'); {*** вывод матрицы B***}

outputmat(n, c, 'C'); {*** вывод матрицы C***}

writeln;

writeln(' ***** Для продолжения нажмите любую клавишу *****');

readkey;

outputmat(n, u, 'U'); {*** вывод результата: матрицы U ***}

writeln;

writeln(' ***** Для продолжения нажмите любую клавишу *****');

readkey;

ClrScr;

writeln;

writeln(' *******************************************************');

writeln(' * Вектор из средних арифметических значений элементов *');

writeln(' * столбцов результирующей матрицы. *');

writeln(' *******************************************************');

sred_arifm(n, u, vect);

writeln; write(' ');

for i:=1 to n do

write(' ',vect[i]:5:2);

writeln;

readkey;

end.

Задание 2

1. Анализ задания и математическая постановка задачи

При решении поставленной задачи необходимо выполнить следующие действия:

1. Ввод исходных данных.

2. Нахождение значения определённого интеграла с использованием метода прямоугольников.

3. Вывод результатов.

При численном интегрировании вместо кривой подынтегральной функции используют заменяющие (аппроксимирующие) её кривые или ломаные линии, для которых вычисление ограниченной ими площади производится в соответствии с достаточно несложными формулами.

Принцип метода прямоугольников состоит в том, что исходный отрезок разбивается на достаточно малые части:

a= x₁< x₂< x₃<…< x_n-1< x_n=b; h= x_k-x_k-1;

площадь каждой такой части (прямоугольника): S_k=h*f(x_k);

соответственно площадь всей фигуры, образованной из n-1 таких прямоугольников: S= S₁+S₂+…+ S_n-2+ S_n-1.Величина S является приближённым значением определённого интеграла, она приближается к истинному значению при увеличении числа n.

Погрешность данного метода определяется абсолютным значением разности приближённых значений определённого интеграла при различных n. Если эта разность меньше требуемой погрешности, то необходимая точность достигнута, и дальнейшее увеличение n не требуется.

2. Разработка схемы алгоритма и её описание

По результатам анализа задания можно составить укрупненную схему алгоритма последовательной структуры:

Полный алгоритм:

h:= (b-a)/n; y1:=y2;

y2:=0; x:=a+h;