program kate_otlad;

{$APPTYPE CONSOLE}

uses
  unit_kate_otlad;

const
  SaveFileName1='u_z5_ot_MY.prn';
  SaveFileName2='u_y5_ot_MY.prn';
  SaveFileName3='u_x5_ot_MY.prn';
  SaveFileName4='gamma5_ot_MY.prn';
  SaveFileName5='nachalo5_ot_MY.prn';

  procedure SetInitial;
 const
  Tab=chr(9);

var
  i, j, k : word;
  uu : real;
   f:text;
begin
  for k:=1 to Mk do
  begin
       uu:=(k-1)*(abs(1+1)/20)-1;
   { uu:=(k-1)*(abs(-0.8+0.7)/60)-0.8; }
   for i:=1 to Mx do
    for j:=1 to Mk do
    u[i,j,k]:=uu;
end;
end;{--------------------SetInitial--------------------}

procedure SaveMatrixToFile;

const
  Tab=chr(9);
var
  i, j, k : word;
  f : text;

begin
  assign(f,SaveFileName1);
  rewrite(f);
for i:=1 to Mx  do
  begin
    for j:=1 to My  do
    begin
      write(f,u[i,j,20]:16);
      if i<Mx+1  then write(f,Tab);
    end;
    writeln(f,'');
  end;

  close(f);
  assign(f,SaveFileName2);
  rewrite(f);

  for i:=1 to Mx  do
  begin
    for k:=1 to Mk  do
    begin
      write(f,u[i,25,k]:16);
      if i<Mx+1  then write(f,Tab);
    end;
    writeln(f,'');
  end;

  close(f);
  assign(f,SaveFileName3);
  rewrite(f);

  for j:=1 to My  do
  begin
    for k:=1 to Mk  do
    begin
      write(f,u[15,j,k]:16);
      if j<My+1  then write(f,Tab);
    end;
    writeln(f,'');
  end;

  close(f);
  assign(f,SaveFileName4);
  rewrite(f);

  for j:=1 to My  do
  begin
    for k:=1 to Mk  do
    begin
      write(f,gamma[15,j,k]:16);
      if j<Mx  then write(f,Tab);
    end;
    writeln(f,'');
  end;
  close(f);
end;{--------------------SaveMatrixToFile--------------------}
Begin
   SetParameters;
   SetInitial;
   Compute;
   SaveMatrixToFile;
End.