xren'

program lab_sm1;

{$APPTYPE CONSOLE}

uses
  SysUtils,Math;
// íàõîæäåíèå îáðàòíîå ìàòðèöû
type matr=array [1..50,1..50] of real;
type vect=array [1..50] of real;

var
a,c,b: array [1..1010] of matr;
y,f:array [1..1010] of vect;
i,j,k,n,m:integer;
obr,matr2,matr3,matr0:matr;
matr1:vect;
alfa:array [1..1010] of matr;
betta:array [1..1010] of vect;

function gauss(a:matr;b:vect):vect;
var i,j,k,l:integer;
max,c,z,tmp:real;
y,e:vect;
begin

   // ~~~gauss~~~
   for k:=1 to m-1 do begin
       max:=abs(a[k][k]);
       l:=k;
       for j:=k+1 to m do begin
           if abs(a[j][k])>max then begin
              l:=j;
              max:=abs(a[j][k]);
           end;
       end;
       if l<>k then begin
          for j:=k to m do begin
             tmp:=a[l,j];
             a[l,j]:=a[k,j];
             a[k,j]:=tmp;
          end;
          tmp:=b[k];
          b[k]:=b[l];
          b[l]:=tmp;
       end;

       for i:=k+1 to m do begin
         c:=a[i,k]/a[k,k];
         b[i]:=b[i]-b[k]*c;
         for j:=k to m do a[i][j]:=a[i][j]-a[k][j]*c;
       end;
   end;

   y[m]:=b[m]/a[m,m];

   for k:=m-1 downto 1 do begin
       z:=b[k];
       for i:=k+1 to m do z:=z-a[k][i]*y[i];
       y[k]:=z/a[k][k];
   end;

   // ~~~gauss~~~


      gauss:=y;
end;
function obrat(a:matr):matr;
var
u,l,tmp,obr:matr;
c,sum:real;
i,j,k,f,t:integer;
e:vect;
x1,y1,b1:array [1..50] of vect;
begin

// ðàçëîæåíèå LU

 f:=1;
   for i := f to m do
               Begin
                  for j:=f to m do
                  begin
                     L[i,j]:=0;
                     U[i,j]:=0;
                  end;
               end;
               for i:=f to m do
               Begin
                  for j:=f to m do
                  Begin
                     U[f,i]:=a[f,i];
                     L[i,f]:=a[i,f]/U[f, f];
                     sum:=0;
                     for k:=f to i do sum:=sum+L[i, k]*U[k, j];
                     U[i,j]:=a[i, j]-sum;
                     if i>j then L[j,i]:=0 else
                     Begin
                        sum := 0;
                        for k := f to i do sum := sum + L[j, k] * U[k, i];
                        if U[i, i] <> 0 then L[j, i] :=(a[j, i] - sum)/U[i, i];
                     end;
                  end;
               end;

       { for i:=1 to m do begin
           for j:=1 to m do write(l[i,j]:3:3,' ');
           writeln;
        end;
        for i:=1 to m do begin
           for j:=1 to m do write(u[i,j]:3:3,' ');
           writeln;
        end;  }
   for j:=1 to m do b1[1][j]:=0;
   for k:=2 to m do b1[k]:=b1[1];
   for k:=1 to m do b1[k][k]:=1;

   for j:=1 to m do y1[j]:=gauss(l,b1[j]);
   for j:=1 to m do x1[j]:=gauss(u,y1[j]);

  for j:=1 to m do
  for i:=1 to m do  obr[i][j]:=x1[j][i];

   obrat:=obr;
 {for i:=1 to m do begin
  for j:=1 to m do write(obr[i,j]:3:3,' ');
  writeln;
  end;}
end;
// proizvedenie matrits
function multi(a,b:matr):matr;
var i,j,k:integer;
tmp:real;
c:matr;
begin
for i:=1 to m do
  for j:=1 to m do begin
     tmp:=0;
     for k:=1 to m do begin
         tmp:=tmp+a[i][k]*b[k][j];
     end;
     c[i][j]:=tmp;
   end;
   multi:=c;
end;
// proizvedenie matritsy na stolbets
function pulti(a:matr;b:vect):vect;
var i,j,k:integer;
tmp:real;
c:vect;
begin
  for i:=1 to m do begin
     tmp:=0;
     for k:=1 to m do begin
         tmp:=tmp+a[i][k]*b[k];
     end;
     c[i]:=tmp;
   end;
   pulti:=c;
end;
// slogenie matrits
function plus(a,b:matr):matr;
var i,j,k:integer;
tmp:real;
c:matr;
begin
for i:=1 to m do
  for j:=1 to m do begin
     c[i][j]:=a[i][j]+b[i][j];
  end;
  plus:=c;
end;
//slozhenie stolbtsov
function plus1(a,b:vect):vect;
var
i:integer;
c:vect;
begin
  for i:=1 to m do  c[i]:=a[i]+b[i];
  plus1:=c;
end;

procedure Matrprog(n,m:integer);
var i,j,k:integer;
begin
for i:=1 to n do
for j:=1 to m do f[i][j]:=16;
for k:=1 to n do
for i:=1 to m do  begin
   for j:=1 to m do begin
       if i=j then begin
          a[k][i,i]:=-1;
          b[k][i,i]:=-1;
       end;
   end;
end;

for k:=1 to n do
for i:=1 to m do
   for j:=1 to m do begin
       if i=j then c[k][i,j]:=4;
       if (i-j=1) or (i-j=-1) then c[k][i,j]:=-1
   end;
 { matr0:=obrat(c[1]);
  matr0:=multi(matr0,c[1]);

  for i:=1 to m do begin
  for j:=1 to m do write(matr0[i,j]:3:3,' ');
  writeln;
  end;
  writeln;  }
alfa[1]:=multi(obrat(c[1]),b[1]);
for i:=1 to m do
for j:=1 to m do alfa[1][i][j]:=(-1)*alfa[1][i][j];
betta[1]:=pulti(obrat(c[1]),f[1]);

for k:=2 to n do begin
    matr2:=multi(a[k],alfa[k-1]);
    matr3:=plus(c[k],matr2);
    obr:=obrat(matr3);
    for i:=1 to m do
      for j:=1 to m do obr[i][j]:=(-1)*obr[i][j];
    alfa[k]:=multi(obr,b[k]);
    for i:=1 to m do
      for j:=1 to m do obr[i][j]:=(-1)*obr[i][j];
    matr1:=pulti(a[k],betta[k-1]);
    for i:=1 to m do matr1[i]:=(-1)*matr1[i];
    betta[k]:=pulti(obr,plus1(f[k],matr1));
end;
 writeln('ewrheth');

y[n]:=betta[n];
for i:=n-1 downto 1 do begin
   y[i]:=pulti(alfa[i],y[i+1]);
   y[i]:=plus1(y[i],betta[i]);
end;
  for i:=1 to n do  begin
    for j:=1 to m do write(y[i][j]:3:3,' ');
    writeln;
  end;
end;

begin
assign(output,'output.txt');
rewrite(output);
 n:=7;m:=7;
Matrprog(n,m);

end.