Геометрическая задача. Объясните как работает (нужно описание алгоритма)

program convex_quadrilateral;
const MaxDotCount = 100;
type dot=record x,y:real; end;
var
   epsilon:real;
   space:array [0..MaxDotCount] of dot;
   dotCount:integer;
   a,b,c,d,count,i:integer;
 
function getY(const a,b,c:integer):real;
         begin
              getY:=(space[c].x * (space[b].y - space[a].y) + space[a].y * space[b].x -
              space[b].y * space[a].x) / (space[b].x - space[a].x);
         end;
 
function onOneLine(const a,b,c:integer):boolean;
         begin
              if (abs(space[a].x-space[b].x)<=epsilon) then 
                  begin
                       if (abs (space[b].x - space[c].x) <= epsilon) then
                          onOneLine:=true else onOneLine:=false; 
                       exit;
                  end                  
                  else if (abs(space[a].y-space[b].y)<=epsilon) then
                           begin
                                if (abs (space[b].y - space[c].y) <= epsilon) then
                                    onOneLine:=true else onOneLine:=false;
                                exit;
                           end                                
                           else if (abs (space[c].y - getY (a, b, c)) <= epsilon) then
                                    onOneLine:=true
                                    else onOneLine:=false;
         end;
function diagonalRule(const a1,a2,b1,b2:integer):boolean;
         var
            y1,y2:real;
         begin
              if abs(space[a1].x - space[a2].x)<= epsilon then 
                  begin
                        if ((space[a1].x - space[b1].x) * (space[a1].x - space[b2].x)) < 0 then
                            diagonalRule:=true else diagonalRule:=false; 
                        exit;                        
                  end                      
                 else if (abs (space[a1].y - space[a2].y) <= epsilon) then begin
                          y1:=space[a1].y;
                          y2:=space[a1].y;
                      end
                         else begin
                                   y1:= getY (a1, a2, b1);
                                   y2:= getY (a1, a2, b2);
                              end;
              if ((y1 - space[b1].y) * (y2 - space[b2].y)) < 0 then
                  diagonalRule:=true else diagonalRule:=false;
end;
begin
     write('Vedite tochnost vichislenii epsilon=');
     readln(epsilon);
     writeln('Vvedite kol-vo tochek');
     readln(dotCount);
     for i:=0 to dotCount-1 do begin
         write('Vvedite x[',i,'] y[',i,']');
         readln(space[i].x,space[i].y);
         end;
     writeln('Resultat');
         if (dotCount < 4) then
             writeln('Ne hvataet tochek,minimum 4')
             else begin
                       count:=0;
                       for a:=0 to dotCount - 4 do
                begin
                           for b:= a + 1 to dotCount - 3 do
                                begin
                                 for c:= b + 1 to dotCount - 2 do
                                        begin
                                         if onOneLine (a, b, c)= false then
                                                begin
                                                  for d:= c + 1 to dotCount-1 do
                                                        begin
                                                      if (diagonalRule (a, c, b, d)=true) and (diagonalRule (b, d, a, c)=true) then begin
                                                          count:=count+1;
                                                          writeln (count,':',a, b, c, d);                                                        
                                                          end
 
                                                          else if (diagonalRule (a, b, c, d)=true) and (diagonalRule (c, d, a, b)=true) then begin
                                                                   count:=count+1;
                                                                   writeln (count,':',a, c, b, d);
                                                                   end
                                                                   else if (diagonalRule (a, d, b, c)=true) and (diagonalRule (b, c, a, d)=true) then 
                                                                        begin
                                                                             count:=count+1;
                                                                             writeln(count,':',a, b, d, c);
                                                                        end; 
                                                        end;
                                                end;
                                        end;
                                end;                                                                         
                end;
                      if count = 0 then 
                         writeln ('Net Vipuklih');
     end;
     writeln ('Nazhmite lubuu klavishu');
     readln;
end.