Зацикливание в программе после операторов for для i и j

program Kursovaja;
 
{$APPTYPE CONSOLE}
 
uses
  SysUtils;
 
var
  c,en,kn,a,emean,anstep,prop,prop1,anconf,ancon,nstep,nconfi,nprint,ran,ncon,eold,temp,kt,u,ann:real;nall,new,im,k,l,i,j:integer;
  xold:array [1..6] of real;
  inn:array [1..6] of integer;
  xhqold:array[1..4,1..3] of real;
  xhq:array [1..1000,1..4,1..3] of real;
  xx:array [1..1000,1..6] of real;
  eij,enew,ener:array [1..1000] of real;
  f1,f2,f3,f4,f5:file of real;
//       common /a/xx(1000,6)/m/delta,psi,per/d/r(4,4)
//       common r2oo/c/xhq(1000,4,3)/r/im,inn,in(6)
const
       mprint=25000;
       per=31.0734;
       delta=0.15;
       psi=0.1;
       nn=1000;
                FUNCTION urand(ISEED:integer):real;
          var IA,IC,IM:integer;
          begin
          IA:=7141;IC:=54773;IM:=259200;
          ISEED:=(ISEED*IA+IC) mod IM;
          urand:=frac(ISEED)/frac(IM);
          end;
                procedure distan(m,n,nall:integer);
        var a,b,c,d:real;r2,cut{60.347}:real;k,g,h:integer;
            a1:array [1..3] of real;
            xx:array [1..1000,1..3] of real;
            r:array [1..4,1..4] of real;
            xhq:array [1..1000,1..4,1..3] of real;
//      common /a/xx(500,6)/m/delta,psi,per/d/r(4,4),r2oo
//        common /c/xhq(500,4,3)
        begin
        cut:=60.347;
        r2:=0;
        nall:=0;
        for g:=1 to 4 do begin
        for h:=1 to 4 do
        r[g,h]:=0;
        end;
        for k:=1 to 3 do begin
        a:=xx[m,k]-xx[n,k];
        b:=abs(a);
        c:=abs(a+per);
        d:=abs(a-per);
        a1[k]:=a;
        if c<b then if c<d then a1[k]:=a+per;
        if d<b then if d<c then a1[k]:=a-per;
        r2:=r2+sqr(a1[k]);
        end;
        if r2>cut then
        nall:=1;
        for k:=1 to 3 do begin
        a:=a1[k];
        for g:=1 to 4 do begin
        for h:=1 to 4 do
        r[g,h]:=r[g,h]+sqr((a+xhq[m,g,k]-xhq[n,h,k]));
        end;
        end;
        for g:=1 to 4 do begin
        for h:=1 to 4 do
        r[g,h]:=sqrt(r[g,h]);
        end;
        end;
        function potent(l,m:integer):real;
        var potent,pot,r6,aoo{238661.8},coo{268.914},rr,s,r12,r1{2.016},r2,ru2{9.7888} :real;
        q:array [1..4] of real;
        i,j:integer;
        r:array [1..4,1..4] of real;
        begin
             aoo:=238661.8;
             coo:=268.914;
             q[1]:=0.2357;
             q[2]:=0.2357;
             q[3]:=-0.2357;
             q[4]:=-0.2357;
             r1:=2.016;
             r2:=1;
             r12:=4.064256;
             ru2:=9.7888;
//             common /d/r(4,4),r2oo
             r6:=r2*r2*r2;
             r12:=sqr(r6);
             pot:=aoo/r12-coo/r6;
             u:=0;
             if(r2>r12) then begin
             for i:=1 to 4 do begin
             for j:=1 to 4 do
             u:=u+332*q[i]*q[j]/r[i,j];
               if(r2<ru2) then begin
               rr:=sqrt(r2);
               s:=(sqr(rr-r1))*(7.3701-2*rr)/1.37764;
              end
               else
               s:=1;
               end;
               u:=u*s;
               end;
               potent:=pot+u;
               end;
            procedure cohh(n:integer);
//            common /a/xx(500,6)/c/xhq(500,4,3)
            var a,b:array [1..3] of real;r2{0.46184},r1{0.5773},c,d:real; i:integer;
            xhq:array [1..1000,1..4,1..3] of real;
            begin
            r1:=0.5773;
            r2:=0.46184;
            for i:=1 to 3 do begin
            u:=xx[n,i+3];
            a[i]:=cos(u);
            b[i]:=sin(u);
             c:=a[3]+b[3];
             d:=a[3]-b[3];
            end;
      xhq[n,1,1]:=((a[1]*a[2]*d-b[1]*c)*r1+a[1]*b[2]*r1);
      xhq[n,1,2]:=((b[1]*a[2]*d+a[1]*c)*r1+b[1]*b[2]*r1);
      xhq[n,1,3]:=((-b[2])*d*r1+a[2]*r1);
      xhq[n,2,1]:=((-a[1]*a[2]*d-b[1]*c)*r1+a[1]*b[2]*r1);
      xhq[n,2,2]:=((-b[1]*a[2]*d+a[1]*c)*r1+b[1]*b[2]*r1);
      xhq[n,2,3]:=(b[2]*d*r1+a[2]*r1);
      xhq[n,3,1]:=(a[1]*a[2]*c+b[1]*d-a[1]*b[2])*r2;
      xhq[n,3,2]:=(b[1]*a[2]*c-a[1]*d-b[1]*b[2])*r2;
      xhq[n,3,3]:=(-b[2]*c-a[2])*r2;
      xhq[n,4,1]:=(-a[1]*a[2]*c-b[1]*d-a[1]*b[2])*r2;
      xhq[n,4,2]:=(-b[1]*a[2]*c+a[1]*d-b[1]*b[2])*r2;
      xhq[n,4,3]:=(b[2]*c-a[2])*r2;
      end;
      procedure move(new:integer;xold:array of real);
//        common /a/xx(500,6)/m/delta,psi,per/r/im,inn,in(6)
        var pi2{6.28318}:real;
        i,iy:integer;
        per1,per2,per3,d,a:real;
        begin
        pi2:=6.28318;
        per1:=per/2;
        per2:=-per1;
        per3:=per;
        d:=delta;
        new:=round(1000*random(1000))+1;
        for i:=1 to 6 do begin
        if i=3 then begin
        xold[3]:=xx[new,3];
        iy:=inn[3];
        a:=xx[new,3]+psi*(2*urand(iy)-1);
        inn[3]:=iy;
        if(a>=pi2) then a:=a-per;
        if(a<0) then a:=a+per;
        xx[new,3]:=a;
        end else begin
        xold[i]:=xx[new,i];
        iy:=inn[i];
        a:=xx[new,i]+d*(2*urand(iy)-1);
        inn[i]:=iy;
        if(a>=per1) then a:=a-per3;
        if(a<per2) then a:=a+per3;
        xx[new,i]:=a;
        end;
        end;
        end;
 
begin
       temp:=298.15;
       kt:=1.985e-3*temp;
       nprint:=0;
       ncon:=0;
       nstep:=0;
       nconfi:=0;
       assign(f1,'mcst1000.txt');
       reset(f1);
       for i:=1 to nn do begin
       for j:=1 to 6 do
       read(f1,xx[i,j]);
       end;
       write(xx[nn,6]);
       readln;
       close(f1);
//       format(2(5i10/),i10)
         for i:=1 to nn do
         cohh(i);
         l:=0;
         for i:=1 to nn do begin
         for j:=1 to nn do begin
         if(j<i) then begin
         u:=0;
         l:=l+1;
        distan(i,j,nall);
         if nall<>1 then begin
         u:=potent(i,j);
         eij[l]:=u;
         end
         else
         eij[l]:=0;
         end;
         end;
         end;
{!}
       repeat
       repeat
           nstep:=nstep+1;
           nprint:=nprint+1;
       {процедура не происходит из-за отсутствия new} move(new,xold);
             eold:=ener[new];
       for i:=1 to 4 do begin
        for j:=1 to 3 do
        xhqold[i,j]:=xhq[new,i,j];
        end;
         cohh(new);
         en:=0;
        for i:=1 to nn do begin
        if i<>new then begin
        u:=0;
        distan(i,new,nall);
        end;
        if nall<>1 then begin
         u:=potent(i,new);
         en:=en+u;
         end;
         enew[i]:=u;
         end;
{!}
           if en>eold then begin
           ran:=urand(im);
           a:=(eold-en)/kt;
           a:=exp(a);
           if a<=ran then begin
           for i:=1 to 6 do
           xx[new,i]:=xold[i];
           for i:=1 to 4 do begin
           for j:=1 to 3 do
           xhq[new,i,j]:=xhqold[i,j];
           end;
           end;
           end;
        until false;
{!}
               nconfi:=nconfi+1;
               ncon:=ncon+1;
               {for i:=1 to nn do
//               if(i-new) then
               kn:=((new-1)*(new-2))/(2+i);
//               else
               kn:=((i-1)*(i-2))/(2+new);
               ener[i]:=ener[i]-eij[kn]+enew[i];
//               (eij[kn]):=enew[i];
}               ener[new]:=en;
               if ncon<mprint then c:=0;
               until false;
            emean:=0;
            for i:=1 to nn do
       emean:=emean+ener[i];
//       ann:=frac(nn);
       emean:=emean/nn;
       prop:= {frac}(nconfi)/nstep;
       prop1:= {frac}(ncon)/nprint;
       anstep:={frac}(nstep)/nn;
       anconf:={frac}(nconfi)/nn;
       ancon:={frac}(ncon)/nn;
       assignfile(f2,'dannye.txt');
      rewrite(f2);
      write (f2,nstep,anstep,nconfi,anconf,ncon,ancon,{prop,prop1,}emean,nstep,nconfi);
//    format(2,'полное число',10x,'в среднем на частицу'/1x,'рассмотр. конфиг.=',i8,5X,f10.2/1x,'принятых конфиг.=',i8,5x,f10.2/1x,'принято на отрезке-',i8,5x,f10.2/))
//    format(1x,'пропорция принятых конфиг.=',f7.3/1x,'на отрезке=',f7.3,/1x,'средняя энергия частицы=',f8.4/)
       close(f2);
       assignfile(f3,'coordsOiUgly.txt');
       rewrite(f3);
       for i:=1 to nn do begin
       for j:=1 to 6 do
       write (f3,xx[i,j]);
       end;
       close(f3);
       assignfile(f4,'coordsHiQ.txt');
       rewrite(f4);
       for i:=1 to nn do begin
       for j:=1 to 4 do begin
       for k:=1 to 3 do
       write(f4,xhq[i,j,k]);
       end;
       end;
       assign(f5,'energiya.txt');
       rewrite(f5);
       for i:=1 to nn do
       write(f5,ener[i]);
       close(f5)
end.

repeat
repeat

repeat
repeat
until ncon<mprint 
until a>ran