Program metod_iteracia;
uses crt;
const
maxn=10;
type
matrix=array [1..maxn,1..maxn] of real;
vector=array [1..maxn] of real;
var
a:matrix;
b,x:vector;
i,n,k:integer;
procedure DiaG(a:matrix; b:vector; var x:vector);
{процедура пошуку початкового наближення}
var
i:integer;
begin
for i:=1 to n do
x[i]:=b[i]/abs(a[i,i]);
end;
procedure SymmA(a:matrix; b:vector; var t:vector);
{процедура пошуку суми елементів рівняння }
var
i,j:integer;
begin
for i:=1 to n do begin
t[i]:=0;
for j:=1 to n do
t[i]:=t[i]+(a[i,j]*b[j]);
end;
end;
procedure MinusveC(a:vector; b:vector; var t:vector);
{процедура пошуку різниці вектора і знайденого значення}
var
i:integer;
begin
for i:=1 to n do
t[i]:=a[i]-b[i]
end;
procedure IteraciA(a:matrix; b:vector; var y:vector);
{процедура розрахунків методом ітерації}
const
m=100;
var
i,j:integer;
d:matrix;
buf,buf1 : vector;
begin
for i:= 1 to n do
begin
for j:=1 to n do
d[i,j]:=0;
d[i,i]:=a[i,i];
y[i]:=0
end;
for i:= 1 to m do
begin
SymmA(a,y,buf);
MinusveC(buf,b,buf);
SymmA(d,y,buf1);
MinusveC(buf1,buf,buf);
DiaG(d,buf,y);
end;
end;
procedure RivN(A: Matrix; b: Vector);
{процедура виводу на екран системи заданих рівнянь}
var
i,j:integer;
begin
writeln('zadana systema rivnan');
writeln('-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=');
for i:=1 to N do
begin
for j:=1 to N do
begin
write('+','(',A[i,j]:2:2,')',' * X[',j,'] ');
end;
writeln(' = ',b[i]:2:2);
end;
writeln('-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=');
writeln;
end;
procedure VvodMatriX(var A: Matrix; var b: Vector);
{процедура введення коефіцієнтів системи рівнянь}
var
i,j:integer;
begin
for i:=1 to N do
begin
for j:=1 to N do
begin
write('A[',i,'.',j,'] = ');
read(A[i,j]);
end;
write('b[',i,'] = ');
readln(b[i]);
end;
end;
{основна програма}
begin
clrscr;
writeln('vvedit n-rozmirnist matruci');
readln(n);
VvodMatriX(a,b);
RivN(a,b);
writeln('pochatkove nabligina');
for i:=1 to n do
begin
x[i]:=b[i]/abs(a[i,i]);
writeln(' x[',i,']=',x[i]:3:4);
x[i]:=0;
end;
writeln;
IteraciA(a,b,x);
writeln('Rishenie metodom iteracii');
for i:=1 to n do
writeln(' X[',I,']=',x[i]:3:3,' ');
writeln;
readln;
end.
Смотри, что выдает твоя программа:
