program Project;

{$APPTYPE CONSOLE}

uses
  Maxv in 'MAXV.PAS';

const
  SaveFileName={'gamma.prn''u_matr3.prn'}'U1.prn';
  SaveFileName1='U3.prn';
  SaveFileName2='U5.prn';
  SaveFileName3='U7.prn';
  SaveFileName4='U9.prn';
  El_Potential=-1;
  FileName1='u_i.prn'{'Un_b.prn'};
  FileName2='u_j.prn'{'U_jq.prn'};
  File1='u_ii.prn'{'J_iq.prn'};
  File2='u_jj.prn'{'J_jq.prn'};

procedure SetInitial;

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

begin
  for n:=0 to nt do
  begin
    for i:=1 to Mr do
      for j:=1 to Mz do
      {if n=0 then}
        u[i,j,n]:= 0; {else  }

        for i:=40 to 210 do
  begin
    t:=(i-40)/(210-1);
    uu:=El_potential*(1-t)-El_potential*t;
    for j:=15 to Mz-20 do u[i,j,n]:=uu;
  end;


  {  for i:=Start_A to Last_A do
      for j:=15 to Mz do
      u[i,j,n]:= - El_Potential;
      for i:=Start_P to Last_P do
       for j:=15 to Mz do
      u[i,j,n]:=  El_Potential; }
  end;
end;{--------------------SetInitial--------------------}


procedure SaveMatrixToFile;

const
  Tab=chr(9);

var
  i, j : word;
  f : text;

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

  for i:=1 to Mr  do
  begin
    for j:=1 to Mz  do
    begin
      write(f,{gamma[i*2,(Mz+1-j)*2]*sqrt(sqr(u[(i-1)*2,(Mz+1-j)*2]-u[i*2,(Mz+1-j)*2])+
                        sqr(u[i*2,(Mz-j)*2]-u[i*2,(Mz+1-j)*2]))/0.02}
                        u[i,Mz+1-j,3]:16);
      if i<Mr  then write(f,Tab);
    end;
    writeln(f,'');
  end;

  close(f);

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

  for i:=1 to Mr  do
  begin
    for j:=1 to Mz  do
    begin
      write(f,{gamma[i*2,(Mz+1-j)*2]*sqrt(sqr(u[(i-1)*2,(Mz+1-j)*2]-u[i*2,(Mz+1-j)*2])+
                        sqr(u[i*2,(Mz-j)*2]-u[i*2,(Mz+1-j)*2]))/0.02}
                        u[i,Mz+1-j,7]: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,{gamma[i*2,(Mz+1-j)*2]*sqrt(sqr(u[(i-1)*2,(Mz+1-j)*2]-u[i*2,(Mz+1-j)*2])+
                        sqr(u[i*2,(Mz-j)*2]-u[i*2,(Mz+1-j)*2]))/0.02}
                        u[i,Mz+1-j,10]: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 j:=1 to Mz  do
    begin
      write(f,{gamma[i*2,(Mz+1-j)*2]*sqrt(sqr(u[(i-1)*2,(Mz+1-j)*2]-u[i*2,(Mz+1-j)*2])+
                        sqr(u[i*2,(Mz-j)*2]-u[i*2,(Mz+1-j)*2]))/0.02}
                        u[i,Mz+1-j,10]: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 j:=1 to Mz  do
    begin
      write(f,{gamma[i*2,(Mz+1-j)*2]*sqrt(sqr(u[(i-1)*2,(Mz+1-j)*2]-u[i*2,(Mz+1-j)*2])+
                        sqr(u[i*2,(Mz-j)*2]-u[i*2,(Mz+1-j)*2]))/0.02}
                        u[i,Mz+1-j,10]:16);
      if i<Mr  then write(f,Tab);
    end;
    writeln(f,'');
  end;

  close(f);

 { assign(f,FileName1);
  rewrite(f);
  for n:=0 to Nt-1 do
  begin
    write(f, u[10,2,n]:16, u[11,2,n]:16, u[10,3,n]:16, u[11,3,n]:16, tab);
    writeln(f,'');
  end;
  close(f);  }


  assign(f,FileName1);
  rewrite(f);
 { for i:=0 to 24 do
  begin
    for j:=25 downto 0 do
    begin
      write(f,(-u[1+i*2,j*2+1]+u[1+(i+1)*2,j*2+1]):16, tab);
      if i<40 then write(f,Tab);  }
       for i:=2 to 125 do
  begin
    for j:=1 to 50 do
    begin
      write(f,(u[(i-1)*2,(Mz+1-j)*2,9]-u[i*2,(Mz+1-j)*2,9]):16);
      if i<50 then write(f,Tab);
    end;
    writeln(f,'');
  end;
  close(f);

  assign(f,FileName2);
  rewrite(f);
 { for i:=0 to 25 do
  begin
    for j:=24 downto 0 do
    begin
      write(f,(u[1+i*2,j*2+1]-u[1+i*2,(j+1)*2+1]):16, tab);
      if i<49 then write(f,Tab); }
  for i:=1 to 125 do
  begin
    for j:=1 to 50 do
    begin
      write(f,(-u[i*2,(Mz+1-(j-1))*2,9]+u[i*2,(Mz+1-j)*2,9]):16);
      if i<25 then write(f,Tab);
    end;
    writeln(f,'');
  end;
  close(f);

 {----------Current-----------------}

   assign(f,File1);
  rewrite(f);
  for i:=1 to Mr-1 do
  begin
    for j:=1 to Mz-1 do
    begin
      write(f,(S[i,(Mz+1-j)]*(-u[i,(Mz+1-j),9]+u[i+1,(Mz+1-j),9])/hh):16);
      if i<Mr-1 then write(f,Tab);
    end;
    writeln(f,'');
  end;
  close(f);

  assign(f,File2);
  rewrite(f);
  for i:=1 to Mr-1 do
  begin
    for j:=1 to Mz-1 do
    begin
      write(f,(S[i,(Mz+1-j)]*(u[i,(Mz+1-j),9]-u[i,(Mz+1-(j+1)),9])/hh):16);
      if i<Mr-1 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]);
  close(f);

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

  assign(f,SaveFileName3);  {L
  rewrite(f);
  for i:=1 to 16 do
    writeln(f,u[32,i]);
  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;
  
 { SaveBoundarytoFile;  }
  SaveMatrixToFile;
   Curr;
End.