На плоскости находятся точки. Найти те точки, образующие замкнутый контур этого множества

uses crt,graphABC;
const nmax=15;
{угол направления от левой точки к правой
от оси Х}
type point=record
           x,y:integer
           end;
function Ugol(a,b:point):real;
begin
if a.x=b.x then
 begin
  if a.y<b.y then Ugol:=pi/2
  else Ugol:=-pi/2;
 end
else
 if a.y=b.y then
  begin
   if a.x<b.x then Ugol:=0
   else Ugol:=pi;
  end
else Ugol:=arctan((b.y-a.y)/(b.x-a.x));
end;
 
var a:array[1..nmax] of point;{исходный массив точек}
    p:array[1..nmax+1] of point;{замкнутая ломаная=полигон}
    b:array[1..nmax] of integer;{массив номеров точек исходного массива}
    c:array[1..nmax+1] of integer;{массив номеров в порядке обхода}
    n:integer;{количество точек в исходном}
    np:integer;{то же в полигоне}
    i,j,imx,k,y,xc,yc:integer;
    x:point;{буфер для обмена при сортировке}
    mx:real;{максимальный угол}
    s:string;
begin
randomize;
repeat
write('2<n<=',nmax,' n=');
readln(n);
until n in [3..nmax];
for i:=1 to n do{координаты точек в пределах экрана}
 begin
  a[i].x:=random(600)-300;
  a[i].y:=random(400)-200;
 end;
for i:=1 to n do{номера по порядку}
b[i]:=i;
writeln('Ishodnyj massiv:');
write(' ':2);
for i:=1 to n do
write(b[i]:5);
writeln;
write('x:');
for i:=1 to n do
write(a[i].x:5);
writeln;
write('y:');
for i:=1 to n do
write(a[i].y:5);
writeln;
{сортировка по возрастанию координаты Х}
for i:=1 to n-1 do
for j:=i+1 to n do
if a[i].x>a[j].x then
 begin
  x:=a[i];
  a[i]:=a[j];{обмениваем точки}
  a[j]:=x;
  y:=b[i];
  b[i]:=b[j];{и их номера}
  b[j]:=y;
 end;
writeln;
writeln('Sortirovannyj massiv:');
write(' ':2);
for i:=1 to n do
write(b[i]:5);
writeln;
write('x:');
for i:=1 to n do
write(a[i].x:5);
writeln;
write('y:');
for i:=1 to n do
write(a[i].y:5);
writeln;
{Построение замкнутой ломаной без пересечений=многоугольника}
np:=1;{первая левая точка}
p[1]:=a[1];
c[1]:=b[1];{ее номер}
i:=1;
while i<n do{до предпоследней}
 begin
  mx:=Ugol(a[i],a[i+1]);{пока макс угол первые 2 соседних}
  imx:=i+1;{номер точки для построения}
  for j:=i+1 to n do{проверяем все точки впереди}
  if Ugol(a[i],a[j])>mx then
   begin
    mx:=Ugol(a[i],a[j]);{выбираем точку с макс углом}
    imx:=j;
   end;
  np:=np+1;{добавляем в массив для полигона координаты и номер}
  p[np]:=a[imx];
  c[np]:=b[imx];
  {удаляем ее из исходного массива}
  if imx=n then n:=n-1{если последняя, обрезаем}
  else{иначе сдвигаем на нее хвост массива}
   begin
    for k:=imx to n-1 do
     begin
      a[k]:=a[k+1];
      b[k]:=b[k+1];
     end;
    n:=n-1;{обрезаем}
   end;
  i:=imx-1;{встаем перед удаленной}
 end;
for i:=n downto 1 do{пишем начиная с конца все оставшиеся точки}
 begin
  np:=np+1;
  p[np]:=a[i];
  c[np]:=b[i];
 end;
writeln('Mnogougolnik:');
write(' ':2);
for i:=1 to np-1 do {первую точку полигона вывожим 1 раз}
write(c[i]:5);
writeln;
write('x:');
for i:=1 to np-1 do
write(p[i].x:5);
writeln;
write('y:');
for i:=1 to np-1 do
write(p[i].y:5);
readln;
clearwindow;
i:=0;
j:=0;
xc:=windowwidth div 2;
yc:=windowheight div 2;
for i:=1 to np do{определяем экранные координаты точек}
 begin
  p[i].x:=xc+p[i].x;
  p[i].y:=yc-p[i].y;
 end;
setpencolor(clRed);
polygon(p,np);{рисуем полигон}
setpencolor(clGreen);
for i:=1 to np do{выводим точки и их номера}
 begin
   circle(p[i].x,p[i].y,2);
   //str(c[i],s);
   textout(p[i].x+10,p[i].y,inttostr(c[i]));
 end;
end.