program PR_New;

{$APPTYPE CONSOLE}

uses
  Fil_XYZ_11;
 { LIZA_uvelich; }
const
  SaveFileName1='101u_x.prn';
  SaveFileName2='1u_y.prn';
  SaveFileName3='1u_z.prn';
  FileName1='1u_1.prn';
  FileName2='1u_2.prn';
  FileName3='1u_3.prn';
  FileName4='1u_4.prn';
  SaveFileName4='c:\program Files\Borland\Delphi6\Projects\10-й сем-р\моя\1u_gamma.prn';
  SaveFileName6='c:\program Files\Borland\Delphi6\Projects\10-й сем-р\моя\1u_Uk.prn';
  SaveFileName5='c:\program Files\Borland\Delphi6\Projects\10-й сем-р\моя\1u_U1.prn';

procedure SetInitial;

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

begin
 for k:=50 to Mk do
  begin
    t:=(k-1)/(Mk-1);
    uu:=1*(1-t)-1*t;
   for i:=1 to Mr do
    for j:=1 to Mz 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 Mr  do
  begin
    for j:=1 to Mz  do
    begin
      write(f,u[i,j,120]:16);
      if i<Mr  then write(f,Tab);
    end;
    writeln(f,'');
  end;

  close(f);


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

  for i:=1 to Mr  do
  begin
    for j:=1 to Mz  do
    begin
      write(f,u[i,j,50]:16);
      if i<Mr  then write(f,Tab);
    end;
    writeln(f,'');
  end;

  close(f);


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

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

  close(f);

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

  for i:=1 to Mr  do
  begin
    for k:=1 to Mk  do
    begin
      write(f,gamma[i,20,k]:16);
      if i<Mr  then write(f,Tab);
    end;
    writeln(f,'');
  end;
  close(f);

   assign(f,SaveFileName5);
  rewrite(f);
   for i:=1 to Mr  do
  begin
    for k:=1 to Mk  do
    begin
      write(f,u[i,21,k]:16);
      if i<Mr  then write(f,Tab);
    end;
    writeln(f,'');
  end;

  close(f);


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


{procedure SaveBoundaryToFile;

const
  Tab=chr(9);

var
  i, j : word;
  f : text;

begin
  assign(f,SaveFileName1); {U
  rewrite(f);
  for i:=32 to 53 do
    writeln(f,u[i,1,1]);
  close(f);

  assign(f,SaveFileName2);  {B
  rewrite(f);
  for i:=32 to 53 do
    writeln(f,u[i,16,1]);
  close(f);

  assign(f,SaveFileName3);  {L
  rewrite(f);
  for i:=1 to 16 do
    writeln(f,u[32,i,1]);
  close(f);

  assign(f,SaveFileName4);  {R
  rewrite(f);
  for i:=1 to 16 do
    writeln(f,u[53,i]);
  close(f);
end;{--------------------SaveBoundaryToFile--------------------}


Begin

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