Розв’язати систему алгебраїчних рівнянь методом Жордана-Гаусса

8.6x1+1.8x2-4.7x3=3.8

2.7x1-3.6x2+1.9x3=2.6

1.5x1+4.5x2+9.3x3=-1.6

Результати розв’язку в середовищі Turbo Pascal

Program gayss;

uses crt;

const n=3;

var i,j,k,m:integer;

a:array[1..n,1..n+1];

x:array[1..n+1];

c,r:real;

Begin

clrscr;

writeln('Введіть значення коефіцієнтів: ');

for i:=1 to n do

for j:=1 to n+1 do

Begin

writeln('a[',i,',',j,']=');

readln(a[i,j]);

for k=1 to n=1 do

Begin

c:=a[k,k];

for m:=1 to n+1 do

a[k,m]:=a[k,m]/c;

for i:=k+1 to n do

Begin

c:=a[i,k];

for j:=1 to n+1 do

a[i,j]:=a[i,j]-c*a[k,j]/a[k,k];

end;

end;

c:=a[n,n];

for m:=1 to n+1 do

a[n,m]:=a[n,m]/c;

writeln('Результати:');

writeln('Трикутна матриця А:');

for i:=1 to n do

Begin

for j:=1 to n+1 do

writeln(' ',a[i,j]:4:3);

writeln;

end;

x[n]:=a[n,n+1];

i:=n-1;

for i:=n-1 downto 1 do

Begin

s:=0;

for j:=i+1 to n do

s:=s+a[i,j]*x[j];

x[i]:=(a[i,n+1]-s);

end;

writeln('Корені рівняння:');

for i:=1 to n do

writeln('x[',i,']=',x:5:3);

readln;

end.