program tri;
uses crt;

function len(a1,b1,a2,b2:byte):real;{длинна отрезка с коорд (a1;b1) и (a2;b2)}
begin
len:=sqrt(sqr(a1-a2)+sqr(b1-b2));
end;

function S(l1,l2,l3:real):real;{находим площадь фигуры(треугольника) по формуле Герона}
var p:real;
begin
p:=(l1+l2+l3)/2;
S:=sqrt(p*(p-l1)*(p-l2)*(p-l3));
end;

var
x1,x2,x3,x4,y1,y2,y3,y4:byte;
flag:boolean;
begin
clrscr;

write('x1,y1> ');
readln(x1,y1);

write('x2,y2> ');
readln(x2,y2);

write('x3,y3> ');
readln(x3,y3);

write('x4,y4> ');
readln(x4,y4);

flag:=
(
S(len(x1,y1,x2,y2),len(x2,y2,x3,y3),len(x3,y3,x1,y1))=
S(len(x1,y1,x2,y2),len(x2,y2,x4,y4),len(x4,y4,x1,y1))+
S(len(x2,y2,x3,y3),len(x3,y3,x4,y4),len(x4,y4,x2,y2))+
S(len(x1,y1,x4,y4),len(x4,y4,x3,y3),len(x3,y3,x1,y1))
);
writeln('[',x4,';',y4,']',flag);
end.

207 Invalid floating point operation(Недопустимая операция с плавающей запятой) .

Возможные причины сообщения:

аргумент функций TRUNC или ROUND не может быть преобразован в целое число, находящееся внутри диапазона типа LONGINT (от -2147483648 до +2147483647);
отрицательный аргумент функции SQRT (извлечение квадратного корня);
аргумент функции LN (логарифм) равен нулю или имеет отрицательное значение;
произошло переполнение стека сопроцессора.

function S(l1,l2,l3:real):real;{находим площадь фигуры(треугольника) по формуле Герона}
var p:real;
begin
p:=(l1+l2+l3)/2; { по формуле Герона p - полупериметр }
S:=sqrt(abs((p-l1)*(p-l2)*(p-l3)/p)); {abs - избегаем взятия корня из отриц. числа}
end;

Type
 TPoint =
   Record
     X, Y: Integer;
   End;

Function min(a, b: Integer): Integer;
 Begin
   min := a;
   if b < a Then min := b
 End;

Function max(a, b: Integer): Integer;
 Begin
   max := a;
   if b > a Then max := b
 End;

Function isInside(Const a: Array Of TPoint;
        b: TPoint; n: Integer): Boolean;
 Var
   i, j, Count: Word;
   T: Double;
 Begin
   Count := 0;
   For i := 0 To Pred(n) Do
     Begin
       j := (i+1) mod n;
       If a[i].Y = a[j].Y Then Continue;

       If (a[i].Y > b.Y) and (a[j].Y > b.Y) Then Continue;
       If (a[i].Y < b.Y) and (a[j].Y < b.Y) Then Continue;
       If (max(a[i].Y, a[j].Y) = b.Y)
         Then Inc(Count)
         Else
           Begin
             If min(a[i].Y, a[j].Y) = b.y Then Continue
             Else
               Begin
                 T := (b.Y - a[i].Y) / (a[j].Y - a[i].Y);
                 If a[i].X + T * (a[j].X - a[i].X) >= b.X Then Inc(Count)
               End
           End;
     End;
   IsInside := ((Count and $0001) <> 0)
 End;

Const
 arr: Array[0 .. 3] Of TPoint =
   (
     { Точки - в порядке соединения в многоугольник !!!
        Если многоугольник = T1-T3-T4-T2-T1, то так: }
     (X:2;Y:2), {T1}
     (X:4;Y:6), {T3}
     (X:6;Y:4), {T4}
     (X:4;Y:3)  {T2}
   );

Const
 Check1: TPoint = (X:7; Y:2);
 Check2: TPoint = (X:1; Y:7);

Var
 b: Boolean;
Begin
 b := IsInside(arr, Check1, 4);
 b := IsInside(arr, Check2, 4);
End.

flag:=
(
S(len(x1,y1,x2,y2),len(x2,y2,x3,y3),len(x3,y3,x1,y1))=
S(len(x1,y1,x2,y2),len(x2,y2,x4,y4),len(x4,y4,x1,y1))+
S(len(x2,y2,x3,y3),len(x3,y3,x4,y4),len(x4,y4,x2,y2))+
S(len(x1,y1,x4,y4),len(x4,y4,x3,y3),len(x3,y3,x1,y1))
);

...
flag:=
(
1.000001*S(len(x1,y1,x2,y2),len(x2,y2,x3,y3),len(x3,y3,x1,y1))<
S(len(x1,y1,x2,y2),len(x2,y2,x4,y4),len(x4,y4,x1,y1))+
S(len(x2,y2,x3,y3),len(x3,y3,x4,y4),len(x4,y4,x2,y2))+
S(len(x1,y1,x4,y4),len(x4,y4,x3,y3),len(x3,y3,x1,y1))
);
writeln('[',x4,';',y4,'] - ',flag);
...

program tri;
uses crt;

function len(a1,b1,a2,b2:byte):real;
begin
len:=sqrt(sqr(a1-a2)+sqr(b1-b2));
end;
function S(l1,l2,l3:real):real;
var p:real;
begin
p:=(l1+l2+l3)/2;
S:=sqrt(abs(p*((p-l1)*(p-l2)*(p-l3)));
end;

var
x1,x2,x3,x4,y1,y2,y3,y4:byte;
flag:boolean;
begin
clrscr;

write('x1,y1> ');
readln(x1,y1);

write('x2,y2> ');
readln(x2,y2);

write('x3,y3> ');
readln(x3,y3);

write('x4,y4> ');
readln(x4,y4);

flag:=
(
1.000001*S(len(x1,y1,x2,y2),len(x2,y2,x3,y3),len(x3,y3,x1,y1))<
S(len(x1,y1,x2,y2),len(x2,y2,x4,y4),len(x4,y4,x1,y1))+
S(len(x2,y2,x3,y3),len(x3,y3,x4,y4),len(x4,y4,x2,y2))+
S(len(x1,y1,x4,y4),len(x4,y4,x3,y3),len(x3,y3,x1,y1))
);
writeln('[',x4,';',y4,'] - ',flag);
end.

Function PointInSidebar(Const SideBar: Array Of TPoint;
        p: TPoint; Count: Integer): Boolean;
 Var
   x1,y1,x2,y2: Double;
   x0, y0, sum, divisor: Double;
   i: integer;
 Begin
   sum := 0;
   y0 := p.y;
   x0 := p.x;
   For i := 0 to Pred(count-1) do
     Begin
       x1:=sidebar[i].x;
       y1:=sidebar[i].y;
       x2:=sidebar[i+1].x;
       y2:=sidebar[i+1].y;

       divisor:=y1*x2-x1*y2-y0*x2+y0*x1+x0*y2-x0*y1;
       divisor := divisor + 1.0E-10;

       Sum := sum + (
       -arctan( (-x2*x1+x1*x0+y0*y1-x2*x0+x2*x2-y2*y0+y2*y2-y1*y2 )/divisor)+
       arctan( (-x1*x1+x2*x1+x1*x0+y0*y1+y1*y2-x2*x0-y1*y1-y2*y0 )/divisor)
       );
     End;

   x1:=sidebar[count-1].x;
   y1:=sidebar[count-1].y;
   x2:=sidebar[0].x;
   y2:=sidebar[0].y;

   divisor:=y1*x2-x1*y2-y0*x2+y0*x1+x0*y2-x0*y1;
   divisor := divisor + 1.0E-10;
   Sum := sum + (
   -arctan((-x2*x1+x1*x0+y0*y1-x2*x0+x2*x2-y2*y0+y2*y2-y1*y2)/divisor)+
   arctan((-x1*x1+x2*x1+x1*x0+y0*y1+y1*y2-x2*x0-y1*y1-y2*y0)/divisor)
   );

   pointinsidebar := not (abs(Sum) < 0.0001)
 End;