ДЕЛФИ_ГОТОВО / 02 вар / 02

program Z2V2;
uses crt;
type
  mas = array[1..10] of real;
  mat = array[1..10] of mas;
var
  A,C: mat;
  b,x,b1: mas;
  d,s: real;
  i,n,j,k: integer;
  f: boolean;

procedure vvod(var n: integer; var b: mas; var A: mat);
var i,j: integer;
begin
     write('‚ўҐ¤ЁвҐ а §¬Ґа­®бвм n='); readln(n);
     writeln;
     writeln('‚ў®¤ ¬ ваЁжл A:');
     for i:=1 to n do
         for j:=1 to n do
         begin
              write('A[',i,',',j,']='); readln(A[i,j]);
         end;
     writeln;
     writeln('‚ў®¤ ўҐЄв®а  b:');
     for i:=1 to n do
     begin write('b[',i,']='); readln(b[i]); end;
     writeln;
     writeln('‚лў®¤ Ёб室­ле ¤ ­­ле:');
     writeln('Њ ваЁж  A:');
     for i:=1 to n do
     begin
          for j:=1 to n do
              write(A[i,j]:8:3);
          writeln;
     end;
     writeln('‚ҐЄв®а B:');
     for i:=1 to n do
         write(b[i]:8:3);
     writeln;
     readln;
end;

procedure vyvod(n: integer; x: mas);
var i,j: integer;
    r: mas;
begin
     writeln;
     writeln('‚лў®¤ ўҐЄв®а  x:');
     for i:=1 to n do
         write(x[i]:8:3);
     writeln;
     for i:=1 to n do
     begin
          r[i]:=0;
          for j:=1 to n do
              r[i]:=r[i]+A[i,j]*x[j];
     end;
     for i:=1 to n do
         r[i]:=r[i]-b[i];
     writeln;
     writeln('‚лў®¤ ўҐЄв®а -­Ґўп§ЄЁ r:');
     for i:=1 to n do
         write(r[i]:8:3);
     writeln;
     readln;
end;

begin
  clrscr;
  vvod(n,b,A);
  f:=true; k:=1; C:=A; b1:=b;
  while f do
  begin
       d:=A[k,k];
       for i:=k to n do
       begin
            for j:=k to n do
                A[i,j]:=A[i,j]/d;
            b[i]:=b[i]/d;
       end;
       for i:=1 to n do
         if i<>k then
         begin
            d:=-A[i,k];
            for j:=1 to n do
              A[i,j]:=A[i,j]+d*A[k,j];
            b[i]:=b[i]+d*b[k];
         end;
       k:=k+1; f:=false;
       for j:=1 to n do
           f:=f or (A[k,j]<>0)
  end;
  for i:=n downto 1 do
  begin
       s:=0;
       for j:=i+1 to n do
           s:=s+x[j]*A[i,j];
       x[i]:=(b[i]-s)/A[i,i];
  end;
  A:=C; b:=b1;
  vyvod(n,x);
end.