{ПРОГРАММА ИТЕРАЦИИ}
uses crt;
const n=3;
      e=0.001;
type arr=array[1..n] of real;
var i,j,shetcik:integer;
    a,c:array[1..n,1..n] of real;
    d,x,x0:arr;

procedure osn(y:arr;var kof:arr);
var i,j:integer;
begin
   for i:=1 to n do
     kof[i]:=0;

   for i:=1 to n do
   begin
      for j:=1 to n do
        kof[i]:=kof[i]+c[i,j]*y[j];

        y[i]:=kof[i];
  end;
    {Прибавляем матрицу D}
          for i:=1 to n do kof[i]:=kof[i]+d[i];
end;


begin
clrscr;

{ коэфф.}

{c[1,1]:=3;c[1,2]:=0.0038;c[1,3]:=0.0049; c[1,4]:=0.0059;d[1]:=1.5136;
c[2,1]:=0.0011;c[2,2]:=2.1;c[2,3]:=0.0032; c[2,4]:=0.0043;d[2]:=1.4782;
c[3,1]:=-0.0005;c[3,2]:=0.0005;c[3,3]:=1.2; c[3,4]:=0.0026;d[3]:=1.0830;
c[4,1]:=-0.0022;c[4,2]:=-0.0011;c[4,3]:=-0.0001; c[4,4]:=0.3;d[4]:=0.3280;}
for i:=1 to n do
  begin
      writeln('введи коэф в',i,' строке');
     for j:= 1 to n do
        readln(c[i,j]);
      writeln('введи свободный член');

      readln(d[i]);
 end;

for i:=1 to n do
for j:=1 to n do
 begin
  if (i<>j) and (c[i,i]<>0) then c[i,j]:=-c[i,j]/c[i,i];
 end;
{Делим сводбные члены}
for i:=1 to n do d[i]:=d[i]/c[i,i];
{убираем главную диагональ}
for i:=1 to n do
 begin
 c[i,i]:=0;
   x[i]:=d[i];
   x0[i]:=0;

end;
  shetcik:=0;
 while  shetcik<n    do
 begin
                       begin
                       for i:=1 to n do
                         x0[i]:=x[i];

                         osn(x0,x);
                       end;
    for j:=1 to n do
      if abs(x[j]-x0[j])<e  then shetcik:=shetcik+1;

 end;
for j:=1 to n do
  write(x[j]:2:2,'  ');

readln;
end.