Построение модели движения краевых дислокаций в кристалле

uses graphABC;
var
  parametr,r,a,b,c,d,e,f,dx:integer;
  x,y:array[1..10,1..14] of integer;
  Lmax,L:real;
begin
  parametr:=60;
  Lmax:=sqrt(5/4)*parametr;
  r:=15;
  For c:=1 to 10 do
    For d:=1 to 14 do
    begin
      x[c,d]:=50+(d-1)*parametr;
      y[c,d]:=50+(c-1)*parametr;
    end;
  SetWindowWidth(900);
  SetWindowHeight(700);
  SetBrushColor(ClPink);
 
  dx:=1;
  For f:=1 to 14 do
  begin
    For e:=0 to parametr do
    begin
      clearwindow;
      if e>0 then
        For a:=1 to 5 do x[a,f]:=x[a,f]+dx;
      if (f<14) then
        if ((x[1,f+1]-x[1,f])=2*r) then
          For a:=1 to 5 do x[a,f+1]:=x[a,f+1]+dx;
      For a:=1 to 10 do
        For b:=1 to 14 do
          For c:=1 to 10 do
            For d:=1 to 14 do
            begin
              L:=sqrt(sqr(x[c,d]-x[a,b])+sqr(y[c,d]-y[a,b]));
              If L<Lmax Then
              begin
                MoveTo(x[a,b],y[a,b]);
                LineTo(x[c,d],y[c,d]);
              end;
            end;
      For c:=1 to 10 do
        For d:=1 to 14 do circle(x[c,d],y[c,d],r);
      sleep(200);
    end;
  end;
end.

uses graphABC;
var
  parametr,r,dx:integer;
  x,y:array[1..10,1..14] of integer;
  Lmax,L:real;
begin
  parametr:=60;
  Lmax:=sqrt(5/4)*parametr;
  r:=15;
  For var c:=1 to 10 do
    For var d:=1 to 14 do
    begin
      x[c,d]:=50+(d-1)*parametr;
      y[c,d]:=50+(c-1)*parametr;
    end;
  SetWindowWidth(900);
  SetWindowHeight(700);
  Centerwindow;
  SetBrushColor(ClPink);
  LockDrawing;
  dx:=1;
  For var f:=1 to 14 do
  begin
    For var e:=0 to parametr do
    begin
      clearwindow;
      if e>0 then
        For var a:=1 to 5 do x[a,f]:=x[a,f]+dx;
      if (f<14) then
        if ((x[1,f+1]-x[1,f])=2*r) then
          For var a:=1 to 5 do x[a,f+1]:=x[a,f+1]+dx;
      For var a:=1 to 10 do
        For var b:=1 to 14 do
          For var c:=1 to 10 do
            For var d:=1 to 14 do
            begin
              L:=sqrt(sqr(x[c,d]-x[a,b])+sqr(y[c,d]-y[a,b]));
              If L<Lmax Then
              begin
                MoveTo(x[a,b],y[a,b]);
                LineTo(x[c,d],y[c,d]);
              end;
            end;
      For var c:=1 to 10 do
        For var d:=1 to 14 do circle(x[c,d],y[c,d],r);
      sleep(200);
      redraw
    end;
  end;
end.

uses
  GraphABC;
 
const
  x0 = 20;
  y0 = 20;
  dX = 40;
  dY = 40;
  R = 10;
  Rows = 10;
  Cols = 10;
  yTerm = 200;
  xMove = 2;
 
var
  maxDist := (dY+dX)*Sqrt(2)/2;
  
type
  Atom = class
    x, y : Real;
    con : List<Atom>;
    
    constructor (x, y : Real);
    begin
      Self.x := x; Self.y := y;
      con := New List<Atom>;
    end;
    
    procedure drawCon;
    begin
      pen.Style := psDash;
      pen.Color := clBlack;
      foreach var other in con do
        Line(Round(x), Round(y), Round(other.x), Round(other.y));
    end;
    
    procedure drawAtom;
    begin
      pen.Style := psClear;
      brush.Color := clPink;
      circle(Round(x), Round(y), r);
    end;
    
    procedure moveRight;
    begin
      if y < yTerm then
        x += xMove;
    end;
    
    function distance(other : Atom) := Sqrt((x-other.x)**2 + (y-other.y)**2);
    
  end;
 
function InitAtoms : List<Atom>;
begin
  Result := New List<Atom>;
  // сами атомы
  for var row := 0 to Rows-1 do
    for var col := 0 to Cols-1 do
      Result.Add(New Atom(x0 + dx*col, y0 + dy*row) );
end;
 
procedure Reconnection(var atoms : List<Atom>);
begin
  for var i := 0 to atoms.Count-1 do
    begin
      atoms[i].con.Clear;
      for var j := 0 to atoms.Count-1 do
        if i <> j then
          if atoms[i].distance(atoms[j]) < maxDist then
            atoms[i].con.Add( atoms[j] );
    end;
end;
 
  
begin
  var atoms := InitAtoms;
  LockDrawing;
  loop 100 do
    begin
      Reconnection(atoms);
      
      Window.Clear;
      foreach var curAtom in atoms do curAtom.drawCon;
      foreach var curAtom in atoms do curAtom.drawAtom;
      Redraw;
      Sleep(20);
  
      for var i := 0 to atoms.Count-1 do
        atoms[i].moveRight;
    end;
end.