program PR_Otlad_MY;

{$APPTYPE CONSOLE}

uses
  Fil_XYZ_o_k_s_d_e11;

{Fil_XYZ_otlad_MY;  }
    const
  SaveFileName1='u_z111_ot_kate.prn';
  SaveFileName2='u_y111_ot_kate.prn';
  SaveFileName3='u_x111_ot_kate.prn';
  SaveFileName4='gamma_ot11_kate.prn';
  SaveFileName5='nachalo_ot11_kate.prn';
  SaveFileName6='nachalo2_ot11_kate.prn';

  procedure SetInitial;
 const
  Tab=chr(9);

var
  i, j, k : word;
  uu : real;
   f:text;
begin

   for i:=1 to Mx do
    for j:=1 to Mk do
     for k:=1 to Mk do
     begin
      u1[i,j,k]:=0;
     end;

  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 i:=1 to Mx do
    for j:=1 to Mk do
     for k:=1 to Mk do
     begin
      u[i,j,k]:=1;
     end;     }



  {    assign(f,SaveFileName5);
  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;     }



end;{--------------------SetInitial--------------------}


procedure SaveMatrixToFile;

const
  Tab=chr(9);

var
  i, j, k : word;
  f : text;

begin
  assign(f,SaveFileName1);
  rewrite(f);

  for i:=2 to Mx-1  do
  begin
    for j:=2 to My-1  do
    begin
      write(f,u[i,j,50]:16);
      if i<Mx{+1}  then write(f,Tab);
    end;
    writeln(f,'');
  end;

  close(f);


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

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

  close(f);


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

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

  close(f);

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

  for i:=1 to Mx  do
  begin
    for k:=1 to Mk  do
    begin
      write(f,gamma[i{35},{j}23,k]:16);
      if i{j}<=Mx  then write(f,Tab);
    end;
    writeln(f,'');
  end;

  close(f);

   assign(f,SaveFileName5);
  rewrite(f);

   for j:=4 to My-3  do
  begin
    for k:=Mk-3 downto 4  do

    begin
      write(f,((U[30,j,k]-U[30,j{+1},k+1])*gamma[30,j,k]/0.00005):16);
      if j<Mx-2  then write(f,Tab);
    end;
    writeln(f,'');
  end;
    close(f);
    assign(f,SaveFileName6);
  rewrite(f);

   for i:=3 to Mx-2  do
  begin
    for k:=Mk-2 downto 3  do

    begin
      write(f,((U[i,25,k]-U[i,25,k+1])*gamma[i,25,k]/0.5):16);
      if i<Mx-2 then write(f,Tab);
    end;
    writeln(f,'');
  end;
    close(f);
end;{--------------------SaveMatrixToFile--------------------}



Begin

   SetParameters;
   SetInitial;
   Compute;
   SaveMatrixToFile;
   {Emkost;}
End.