Точки на плоскости. Пересечение треугольников.

const nmax=20
{ Используемые типы }
type
   float = double; { Для простоты замены используемого типа в зависимости
                     от необходимой точности и ограничения на память. }
   tpoint    = record x, y : float; end;
   tsegment    = record { Отрезок, задается кому как нравится }
     case byte of
     0: (a: array [1..2] of tpoint;);
     1: (p1, p2: tpoint);
     2: (x1, y1, x2, y2: float);
   end;
   tpoly = array[0..maxn] of tpoint;
     { Первый способ задания полигона, количество вершин
       храниться в отдельной переменной, лучше использовать
       когда в задаче рассматривается 1-2 полигона }
   tpolygon = record
     n : integer;
     p : tpoly;
     { Второй способ задания полигона, используется когда
       полигонов много}
   end;
   ttriangle = array[1..3] of tpoint;
   tline = record a,b,c: float; end; { Для Ax + By + C = 0 }
 
 
{ _Проверка вхождения точки в треугольник._
 Опять же векторное произведение. Проверяем лежит ли точка по одну сторону от
 всех сторон треугольника.}
 
function intriangle(q: tpoint; t: ttriangle): boolean;
var v: float;
begin
 v:= (q.x-t[1].x)*(t[2].y-t[1].y)-(q.y-t[1].y)*(t[2].x-t[1].x);
 intriangle:= false;
 if((q.x-t[2].x)*(t[3].y-t[2].y)-(q.y-t[2].y)*(t[3].x-t[2].x))*v>-eps then
   if((q.x-t[3].x)*(t[1].y-t[3].y)-(q.y-t[3].y)*(t[1].x-t[3].x))*v>-eps then
   intriangle:= true;
end;

var
begin
.........
end.

Можно так.
uses wincrt,graph;
//удвоенная площадь треугольника
function Square(a,b,c:pointtype):longint;
begin
Square:=abs((a.x-c.x)*(b.y-c.y)-(b.x-c.x)*(a.y-c.y));
end;
//принадлежность точки треугольнику
function InTriangle(q,a,b,c:pointtype):boolean;
begin
InTriangle:=Square(q,a,b)+Square(q,a,c)+Square(q,b,c)=Square(a,b,c);
end;
 
var a:array[1..60] of pointtype;
    n,i,j,k,g,m:integer;
    x:pointtype;
    f:boolean;
begin
randomize;
repeat
write('Введите количество треугольников от 2 до 20 n=');
readln(n);
until n in [2..20];
write('Нажмите Enter');
readln;
g:=0;
initgraph(g,m,'');
for i:=1 to 2 do
 begin
  a[i].x:=1+random(getmaxX-2);
  a[i].y:=1+random(getmaxY-2);
  circle(a[i].x,a[i].y,3);
 end;
for i:=3 to 3*n do
 begin
  repeat
  a[i].x:=1+random(getmaxX-2);
  a[i].y:=1+random(getmaxY-2);
  f:=true;
  j:=1;
  while(j<i-1)and f do
   begin
    k:=j+1;
    while(k<i)and f do
    if Square(a[i],a[j],a[k])=0 then f:=false
    else inc(k);
    if f then inc(j);
   end;
  until f;
  circle(a[i].x,a[i].y,3);
 end;
for i:=1 to 3*n-1 do
for j:=i+1 to 3*n do
if a[i].x>a[j].x then
 begin
  x:=a[i];
  a[i]:=a[j];
  a[j]:=x;
 end;
for i:=1 to n do
 begin
  moveto(a[i*3].x,a[i*3].y);
  lineto(a[i*3-1].x,a[i*3-1].y);
  lineto(a[i*3-2].x,a[i*3-2].y);
  lineto(a[i*3].x,a[i*3].y);
 end;
readkey
end.