Задача на построение выпуклой оболочки

m:=(Image1.width-40)/max;//40 для отступа от краев

unit Unit1;
 
interface
 
uses
  Windows, Messages, SysUtils, Variants, Classes, Graphics, Controls, Forms,
  Dialogs, Grids, StdCtrls, Spin, ExtCtrls;
 
type
  TForm1 = class(TForm)
    Button1: TButton;
    Button2: TButton;
    SpinEdit1: TSpinEdit;
    Label1: TLabel;
    Sg1: TStringGrid;
    Label2: TLabel;
    Label3: TLabel;
    Sg2: TStringGrid;
    Image1: TImage;
    procedure Button1Click(Sender: TObject);
    procedure Button2Click(Sender: TObject);
  private
    { Private declarations }
  public
    { Public declarations }
  end;
  Tvyp=record  //точки оболочки с номерами
       x,y:integer;
       k:integer;{номер точки}
       end;
var
  Form1: TForm1;
  n:integer;
  p:array[1..100] of TPoint; //массив всех точек
  p1:array[1..100] of Tvyp; //выпуклая оболочка
implementation
 
{$R *.dfm}
{находится ли точка левее луча}
function IsLeft(p,p1,p2:TPoint):boolean;
begin
IsLeft:=(p1.x-p.x)*(p2.y-p.y)-(p1.y-p.y)*(p2.x-p.x)>0;
end;
 
procedure TForm1.Button1Click(Sender: TObject);
var i:integer;
begin
n:=SpinEdit1.Value;
if n=0 then
 begin
  Showmessage('Введите количество точек');
  exit;
 end;
with sg1 do
 begin
  ColCount:=3;
  RowCount:=n+1;
  FixedCols:=1;
  FixedRows:=1;
  defaultcolwidth:=40;
  for i:=1 to n do
  Cells[0,i]:=IntToStr(i);
  Cells[0,0]:='№';
  Cells[1,0]:='X';
  Cells[2,0]:='Y';
 end;
with sg2 do
 begin
  ColCount:=3;
  RowCount:=2;
  FixedCols:=1;
  FixedRows:=1;
  defaultcolwidth:=40;
  Cells[0,0]:='№';
  Cells[1,0]:='X';
  Cells[2,0]:='Y';
 end;
Showmessage('Введите координаты в таблицу')  
end;
 
procedure TForm1.Button2Click(Sender: TObject);
var i,j,k,t,t1,m,newt:integer;
    xmx,xmn,ymx,ymn,x0,y0:integer;
    ms:real;
    s:string;
begin
k:=0;
for i:=1 to n do
for j:=1 to 2 do
if StrToIntDef(Sg1.Cells[j,i],-1)=-1 then k:=1;
if k=1 then
 begin
  Showmessage('Не все координаты введены или неправильный формат');
  exit
 end;
for i:=1 to n do
 begin
  p[i].x:=StrToInt(Sg1.Cells[1,i]);
  p[i].y:=StrToInt(Sg1.Cells[2,i]);
 end;
//строим оболочку
t1:=1;{выбираем первую точку}
for i:=1 to n do
if (p[i].y<p[t1].y)or((p[i].y=p[t1].y) and (p[i].x<p[t1].x)) then t1:=i;
t:=t1;
m:=0;
repeat
m:=m+1;//запоминаем координаты и номер выбранной точки
p1[m].x:=p[t].x;p1[m].y:=p[t].y;p1[m].k:=t;
newt:=0;
{ Поиск следующей вершины }
for i:=1 to n do
if i<>t then
if (newt=0) or (IsLeft(p[t], p[i], p[newt]))then newt:=i;
t:=newt;
until t=t1;//дошли обратно до первой выбранной точки
for i:=1 to m do
with sg2 do
 begin
  Cells[0,i]:=IntToStr(i);
  Cells[1,i]:=IntToStr(p1[i].x);
  Cells[2,i]:=IntToStr(p1[i].y);
  RowCount:=RowCount+1;
 end;
m:=m+1;
p1[m].x:=p1[1].x;p1[m].y:=p1[1].y;{замкнем оболочку}
with sg2 do
 begin
  Cells[0,m]:='1';
  Cells[1,m]:=IntToStr(p1[m].x);
  Cells[2,m]:=IntToStr(p1[m].y);
 end;
//найдем мин и макс по Х и У  
xmn:=p[1].x;ymn:=p[1].x;
ymn:=p[1].Y;ymx:=p[1].Y;
for i:=1 to n do
 begin
  if p[i].x<xmn then xmn:=p[i].X;
  if p[i].x>xmx then xmx:=p[i].X;
  if p[i].y<ymn then ymn:=p[i].y;
  if p[i].y>ymx then ymx:=p[i].y;
 end;
Image1.Height:=Image1.Width; //квадратный Image
//определим масштаб
if xmx-xmn>ymx-ymn then ms:=(Image1.Width-40)/(xmx-xmn)
else ms:=(Image1.Width-40)/(ymx-ymn);
x0:=20; //начало отсчета по Х и Y
y0:=20;
with Image1.Canvas do
 begin
   //нарисуем выпуклую оболочку
  pen.Color:=clBlue;
  moveto(x0+round((p1[1].x-xmn)*ms),y0+round((p1[1].y-ymn)*ms));
  for i:=1 to m do
  lineto(x0+round((p1[i].x-xmn)*ms),y0+round((p1[i].y-ymn)*ms));
  //нарисуем все точки
  pen.Color:=clRed;
  for i:=1 to n do
   begin
    ellipse(x0+round((p[i].x-xmn)*ms)-2,y0+round((p[i].y-ymn)*ms)-2,
            x0+round((p[i].x-xmn)*ms)+2,y0+round((p[i].y-ymn)*ms)+2);
    str(i,s);
    textout(x0+round((p[i].x-xmn)*ms)+5,y0+round((p[i].y-ymn)*ms),s);
   end;
 end;
end; 
end.