Программа для графической иллюстрации сортировки массивов слиянием

uses crt, graph;
var drive,mode,x,y,i,t,y0:integer;
begin
 drive:=detect;
 initgraph(drive,mode,'c:\tp\bgi');
 setfillstyle(1,1);
 floodfill(1,4,1);
 t:=–4;
 y0:=10;
 setcolor(16);
 for i:=1 to 150 do
 begin
 t:=t+2;
 y0:=y0+3;
 x:=getmaxx div 2 + trunc(cos(t/10)*i);
 y:=y0 – trunc(sin(t/10)*i);
 setfillstyle(1,14);
 fillellipse(x,y,20,20);
 delay(100);
 end;
 repeat
 until keypressed;
 closegraph;
end.

uses graph,crt;
const n=10; {количество элементов в исходных массивах}
var a,b:array[1..n+1] of integer;{исходные массивы значений}
    c:array[1..2*n] of integer;{итоговый массив}
    x,y1,y2,y3,i,j,k,ia,ib,mx,d,w,h:integer;
    ms:real;
    ma,mb:set of byte;
begin
randomize;
initgraph(x,y1,'');
a[1]:=5+random(5);
b[1]:=5+random(5);
mx:=0;{максимальное значение}
for i:=2 to n do
 begin
  a[i]:=a[i-1]+random(10);
  if a[i]>mx then mx:=a[i];
  b[i]:=b[i-1]+random(10);
  if b[i]>mx then mx:=b[i];
 end;
d:=5;{пропуски между столбцами}
w:=(getmaxX-d*(2*n+1)) div (2*n); {ширина столбцов}
y1:=getmaxY div 3-20;{нижняя граница столбцов 1 массива}
y2:=2*getmaxY div 3-20;{нижняя граница столбцов 2 массива}
y3:=getmaxY-20;
x:=d;{левый край первого столбца}
ms:=(y1-20)/mx;{масштаб для определения высоты столбцов на экране}
for i:=1 to n do{исходное положение столбцов 1 массива}
 begin
  rectangle(x,y1-round(a[i]*ms),x+w,y1);
  x:=x+w+d;
 end;
x:=d;
for i:=1 to n do{исходное положение столбцов 2 массива}
 begin
  rectangle(x,y2-round(b[i]*ms),x+w,y2);
  x:=x+w+d;
 end;
delay(1000);
a[n+1]:=maxint;
b[n+1]:=maxint;
ia:=1;
ib:=1;
ma:=[];
mb:=[];
for i:=1 to 2*n do
 begin
  cleardevice; {очищаем экран}
  if a[ia]>b[ib] then
   begin
    c[i]:=b[ib];
    mb:=mb+[ib];
    ib:=ib+1;
    x:=d;
    for j:=1 to n do {новое положение столбцов 1 массива}
     begin
      if not(j in ma) then rectangle(x,y1-round(a[j]*ms),x+w,y1);
      x:=x+w+d;
     end;
    x:=d;
    for j:=1 to n do{исходное положение столбцов 2 массива}
     begin
      if not(j in mb) then rectangle(x,y2-round(b[j]*ms),x+w,y2);
      x:=x+w+d;
     end;
    x:=d;
    for k:=1 to i do{исходное положение столбцов 2 массива}
     begin
      rectangle(x,y3-round(c[k]*ms),x+w,y3);
      x:=x+w+d;
     end;
   end
  else
   begin
    c[i]:=a[ia];
    ma:=ma+[ia];
    ia:=ia+1;
    x:=d;
    for j:=1 to n do {новое положение столбцов 1 массива}
     begin
      if not(j in ma) then rectangle(x,y1-round(a[j]*ms),x+w,y1);
      x:=x+w+d;
     end;
    x:=d;
    for j:=1 to n do{исходное положение столбцов 2 массива}
     begin
      if not(j in mb) then rectangle(x,y2-round(b[j]*ms),x+w,y2);
      x:=x+w+d;
     end;
    x:=d;
    for k:=1 to i do{исходное положение столбцов 2 массива}
     begin
      rectangle(x,y3-round(c[k]*ms),x+w,y3);
      x:=x+w+d;
     end;
   end;
  delay(1000);
 end;
outtextXY(getmaxX div 2-150,20,'Sortirovka zakonchena, press any key');
readkey
end.