Найти точки образующие замкнутый контур плоскости

var x0,x1,x2,y0,y1,y2,i,j,u,n,k1,k2: integer; 
x: array [1..20] of integer;
y: array [1..20] of integer;
xx: array [1..20] of integer;
yy: array [1..20] of integer;
 
procedure lin(x0,x1,x2,y0,y1,y2,n,k1,k2: integer);
var f,i:integer;
begin
  f:=(y1-y1)*x0+(x2-x1)*y0+(x1*y2-x2*y1);
  if f < 0 then inc (k1);
  if f > 0 then inc (k2);
end;
 
begin
  writeln('Введите количество точок на плоскости');
  readln(n);
  for i:= 1 to n do
  begin
    x[i]:= random(10);
    y[i]:= random(10);
  end;
  for j:= 1 to n do
  writeln('(',x[j],';',y[j],')');
  writeln('Точки, что образуют замкнутый контур плоскости');
  for j:=i to n do
  begin
    x1:=x[j];
    y1:=y[j];
    for u:= 1 to n do
      x2:=x[u];
      y2:=x[u];
      for i:= 1 to n do
      begin
        x0:=x[i];
        y0:=x[i];
        lin(x0,x1,x2,y0,y1,y2,n,k1,k2);
      end;
      if (k1 = 0) or (k2 = 0) then
      begin
        xx[j]:= x[j];
        yy[j]:= y[j];
      end;
  end;
  for j:= 1 to n do
  writeln('(',xx[j],';',yy[j],')');
end.

const armax=100;
var x,y,x1,y1:array[1..armax] of integer; i,j,n,k,gd,gm:integer; max_it:real;
function it(x1,y1,x2,y2,x3,y3:integer):real;
var l1,l2,l3,cosx:real;
begin
  l1:=sqrt(sqr(x2-x1)+sqr(y2-y1));
  l2:=sqrt(sqr(x3-x2)+sqr(y3-y2));
  l3:=sqrt(sqr(x3-x1)+sqr(y3-y1));
  cosx:=(l1*l1+l2*l2-l3*l3)/(2*l1*l2);
  it:=arccos(cosx);
end;
begin
 writeln('Количество точек:');
 read(n);
 j:=2;
 max_it:=0;
 for i:=1 to n do
    begin
      writeln('x[',i,']');
      read(x[i]);
      writeln('y[',i,']');
      read(y[i]);
    end;
 y1[1]:=y[1];
 x1[1]:=x[1];
  for i:=1 to n do
        begin
         if y[i]<y1[1] then begin y1[1]:=y[i]; x1[1]:=x[i]; end
         else if (y[i]=y1[1]) and (x[i]>x1[1]) then
                begin
                  y1[1]:=y[i]; x1[1]:=x[i];
                end;
        end;
  for i:=1 to n do
        begin
         if (x1[1]<>x[i]) and (y1[1]<>y[i]) then if (it((x1[1]-3),y1[1],x1[1],y1[1],x[i],y[i])>max_it) then
           begin
              max_it:=it((x1[1]-3),y1[1],x1[1],y1[1],x[i],y[i]);
              x1[2]:=x[i];
              y1[2]:=y[i];
           end;
        end;
  repeat
    inc(j);
    max_it:=0;
    for i:=1 to n do
        begin
          if ((x1[j-2]=x[i]) and (y1[j-2]=y[i])) or
             ((x1[j-1]=x[i]) and (y1[j-1]=y[i])) then continue
          else
           begin
            if (it(x1[j-2],y1[j-2],x1[j-1],y1[j-1],x[i],y[i])>max_it) then
               begin
                 max_it:=it(x1[j-2],y1[j-2],x1[j-1],y1[j-1],x[i],y[i]);
                 x1[j]:=x[i];
                 y1[j]:=y[i];
               end;
           end;
        end;
  until (x1[j]=x1[1]) and (y1[j]=y1[1]);
  writeln('Оболочка идет через точки:');
  for i:=1 to j-1 do writeln('x[',i,']=',x1[i],'   ','y[',i,']=',y1[i]);
end.