program zeid;
label 10,20;
const       e=0.01;
var i,j,f:integer;
    c:array[1..4,1..4] of real;
    d:array[1..4] of real;
    x,x0,b:array[1..4] of real;
    s_ch_1str,s_ch_2str,s_ch_3str,s_ch_4str:real;
begin
 writeln('    Введи матрицу C');
 for i:=1 to 4 do
  begin
   for j:=1 to 4 do
     readln(c[i,j]);
   writeln;
  end;
 writeln('    Введи матрицу сободных членовD');
 for i:=1 to 4 do
   readln(d[i]);
   writeln;

 for i:=1 to 4 do
     x0[i]:=0;
10:  for i:=1 to 4 do
    begin
     x[i]:=d[i];
     for j:=1 to 4 do
     x[i]:=x[i]+c[i,j]*x0[j];
   end;
   {проверка}
   for i:=1 to 4 do
   b[i]:=abs(x[i]-x0[i]);

      f:=0;
  for i:=1 to 4 do begin
  if b[i]>e then f:=1 end;
  if f=0 then goto 20;
  for i:=1 to 4 do
  x0[i]:=x[i];
  goto 10;
20: for i:=1 to 4 do begin
   write('x[',i,']= ',x[i]);
   writeln;          end;
   s_ch_1str:=1.7*x[1]+0.0003*x[2]+0.0004*x[3]+0.0005*x[4];
   s_ch_2str:=0.8*x[2]+0.0001*x[3]+0.0002*x[4];
   s_ch_3str:=-0.0003*x[1]-0.0002*x[2]-0.0003*x[3]+0*x[4];
   s_ch_4str:=-0.000*x[1]-0.0004*x[2]-0.0003*x[3]-x[4];
   writeln('Свободный член первой строки равен: ',s_ch_1str);
   writeln('Свободный член второй строки равен: ',s_ch_2str);
   writeln('Свободный член третьей строки равен: ',s_ch_3str);
   writeln('Свободный член четвертой строки равен: ',s_ch_4str);
end.