На плоскости даны 1000 точек...

uses graph;
const t=0.001;
      n=500;{у меня при 1000 зависае, даже 500 долго считает}
{oпределение площади}
function Plosh(a,b,c:pointtype):real;
begin
Plosh:=abs(a.x*(b.y-c.y)+b.x*(c.y-a.y)+c.x*(a.y-b.y))/2;
end;
{проверка внутри или нет}
function Prin(a,b,c,d:pointtype):boolean;
var s,s1,s2,s3:real;
begin
s:=Plosh(a,b,c);{площадь данного треугольника}
s1:=Plosh(a,b,d);{1 маленький}
s2:=Plosh(a,c,d);{2}
s3:=Plosh(b,c,d);{3}
if abs((s1+s2+s3)-s)<t then Prin:=true
{если площади равны с заданной точностью}
else Prin:=false;{не равны}
end;
var m:array[1..n] of pointtype;
    xc,yc,ms,i,j,k,imx,jmx,kmx:integer;
    mx:real;
begin
randomize;
i:=0;
initgraph(i,j,'');
{центр экрана}
xc:=getmaxX div 2;
yc:=getmaxY div 2;
{создаем массив точек, чтобы не считать большие числа, возьмем
координаты по х=[-32..32] по y=[-24..24], потом умножим на 10 при рисовании}
for i:=1 to n do
 begin
  m[i].x:=random(xc div 5)-(xc div 10);
  m[i].y:=random(yc div 5)-(yc div 10);
 end;
{найдем вершины максимального по площади}
mx:=Plosh(m[1],m[2],m[3]);
imx:=1;
for i:=1 to n-2 do
for j:=i+1 to n-1 do
for k:=j+1 to n do
if Plosh(m[i],m[j],m[k])>mx then
 begin
  mx:=Plosh(m[i],m[j],m[k]);
  imx:=i;
  jmx:=j;
  kmx:=k;
 end;
ms:=10;{масштаб}
{нарисуем треугольник}
setcolor(10);
circle(xc+m[imx].x*10,yc-m[imx].y*10,3);
setfillstyle(1,10);
floodfill(xc+m[imx].x*10,yc-m[imx].y*10,10);
circle(xc+m[jmx].x*10,yc-m[jmx].y*10,3);
floodfill(xc+m[jmx].x*10,yc-m[jmx].y*10,10);
circle(xc+m[kmx].x*10,yc-m[kmx].y*10,3);
floodfill(xc+m[kmx].x*10,yc-m[kmx].y*10,10);
line(xc+m[imx].x*10,yc-m[imx].y*10,xc+m[jmx].x*10,yc-m[jmx].y*10);
line(xc+m[jmx].x*10,yc-m[jmx].y*10,xc+m[kmx].x*10,yc-m[kmx].y*10);
line(xc+m[imx].x*10,yc-m[imx].y*10,xc+m[kmx].x*10,yc-m[kmx].y*10);
{нарисуем точки}
for i:=1 to n do
 begin
  if Prin(m[imx],m[jmx],m[kmx],m[i]) then setcolor(9){внутри}
  else setcolor(12);{снаружи}
  circle(xc+m[i].x*10,yc-m[i].y*10,1);
 end;
readln
end.