Научная работа

PASCAL:Задача решения систем уравнений методом Гаусса
Богданова М.В.

{Метод Гаусса}
uses crt;
const n=3;

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

begin
clrscr;
for i:=1 to n do begin
writeln('Уравнение:',i);
for j:=1 to n do
begin
write('Введите A(',i,',',j,'):');
readln(a[i,j]);
end;
write('Свободный член: B(',i,'):');
readln(b[i]);
end;
for i:=1 to n do
begin
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;
x[n]:=b[n]/a[n,n];
for i:=n-1 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;
end.
writeln('Решение системы:');
for i:=1 to n do
writeln('X(',i,')=',x[i]:3:5);
readln;
end.