Здравствуйте, Посетитель - 381027!

Текст программы:
[code h=200]
const n=5;

var a,aa:array[1..n,1..n] of real;
b,x:array[1..n] of real;
i,j,k:integer;
h:real;
fa,fb:text;

begin
{Считывание матрицы A}
assign(fa,'a.txt');
reset(fa);
for i:=1 to n do
begin
for j:=1 to n do
read(fa,a[i,j]);
readln(fa);
end;
close(fa);
{Копия матрицы для проверки}
aa:=a;

{Считывание вектора B}
assign(fb,'b.txt');
reset(fb);
for i:=1 to n do
readln(fb,b[i]);
close(fb);
bb:=b;

{Прямой ход - исключение переменных}
for i:=1 to n-1 do
begin
j:=i;
{Перестановка строк при нулевом диагональном элементе}
while (a[j,i]=0) and (j<=n) do inc(j);
if j<>i then
begin
for k:=1 to n do
begin
h:=a[i,k];
a[i,k]:=a[j,k];
a[j,k]:=h;
end;
h:=b[i];
b[i]:=b[j];
b[j]:=h;
end;
{Обнуление поддиагональных коэффициентов A}
for j:=i+1 to n do
begin
a[j,i]:=-a[j,i]/a[i,i];
for k:=i+1 to n do
a[j,k]:=a[j,k]+a[j,i]*a[i,k];
b[j]:=b[j]+a[j,i]*b[i]
end;
end;

{Обратный ход - нахождение корней}
for i:=n downto 1 do
begin
h:=b[i];
for j:=i+1 to n do
h:=h-x[j]*a[i,j];
x[i]:=h/a[i,i]
end;

writeln('Ответ: X=');
for i:=1 to n do
writeln(x[i]:0:3);

writeln('Проверка: A*x=');
for i:=1 to n do
begin
b[i]:=0;
for j:=1 to n do
b[i]:=b[i]+aa[i,j]*x[j];
writeln(b[i]:0:3);
end;
readln;

end.[/code]

Исходный файл с матрицей A -- a.txt


Исходный файл с вектором b -- b.txt


Результат выполнения: