program Pr;

{$APPTYPE CONSOLE}

uses
  Fil_XYZ_otlad;

const
  SaveFileName1='u_z_ot.prn';
  SaveFileName2='u_y_ot.prn';
  SaveFileName3='u_x_ot.prn';
  SaveFileName4='gamma_ot.prn';
procedure SetInitial;

var
  i, j, k : word;
  uu : real;

begin
  for k:=1 to Mk do
  begin
    uu:=(k-1)*(abs(-0.8+0.7)/50)-0.8;
   for i:=1 to Mx do
    for j:=1 to Mk do
    u[i,j,k]:=uu;

  end;

  {for k:=1 to 50 do
  begin
     t:=(k-1)/(50-1);
    uu:=0*(1-t)+1*t;
    for i:=1 to Mr do
    for j:=1 to Mz 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}5]:16);
      if i<Mx  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,5,k]:16);
      if i<Mx  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[{35}4,j,k]:16);
      if j<Mx  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[{35}6,j,k]:16);
      if j<Mx  then write(f,Tab);
    end;
    writeln(f,'');
  end;

  close(f);


end;{--------------------SaveMatrixToFile--------------------}



Begin

   SetParameters;
   SetInitial;
   Compute;
   SaveMatrixToFile;
   {Current;}

End.