Вызывать функцию в программе несколько раз

uses crt;
const t=0.001;
function IsPixelInRegion(i,j:real; x1,x2,x3,x4,y1,y2,y3,y4:real):boolean;
var s1,s2,s3,s4,s41,s42,p:real;
    a,b,c,d,a1,b1,c1,d1,diag1,diag2:real;
begin
a:=sqrt(sqr(x2-x1)+sqr(y2-y1));
b:=sqrt(sqr(x3-x2)+sqr(y3-y2));
c:=sqrt(sqr(x4-x3)+sqr(y4-y3));
d:=sqrt(sqr(x4-x1)+sqr(y4-y1));
 
a1:=sqrt(sqr(x1-i)+sqr(y1-j));
b1:=sqrt(sqr(x2-i)+sqr(y2-j));
c1:=sqrt(sqr(x3-i)+sqr(y3-j));
d1:=sqrt(sqr(x4-i)+sqr(y4-j));
 
diag1:=sqrt(sqr(x3-x1)+sqr(y3-y1));
diag2:=sqrt(sqr(x4-x2)+sqr(y4-y2));
 
{1}p:=(a+a1+b1)/2;
s1:=sqrt(p*(p-a)*(p-a1)*(p-b1));
{2}p:=(b+b1+c1)/2;
s2:=sqrt(p*(p-b)*(p-b1)*(p-c1));
{3}p:=(c+c1+d1)/2;
s3:=sqrt(p*(p-c)*(p-c1)*(p-d1));
{4}p:=(d+d1+a1)/2;
s4:=sqrt(p*(p-d)*(p-d1)*(p-a1));
 
p:=(diag1+a+b)/2;
s41:=sqrt(p*(p-a)*(p-b)*(p-diag1));
p:=(diag2+c+d)/2;
s42:=sqrt(p*(p-c)*(p-d)*(p-diag2));
IsPixelInRegion:=abs(s41+s42-s1-s2-s3-s4)<t
end;
var x1,y1,x2,y2,x3,y3,x4,y4,tx,ty:real;
begin
writeln('Введите координаты вершин 4х-угольника в порядке обхода');
read(x1,y1,x2,y2,x3,y3,x4,y4);
writeln('Введите координаты точки');
read(tx,ty);
if IsPixelInRegion(tx,ty,x1,x2,x3,x4,y1,y2,y3,y4) then
write('Точка в 4х-угольнике') else write('Точка не в 4х-угольнике')
end.

function kor( x1,x2,x3,x4:real):real;
    var a:real;
  begin
    a:=sqrt(sqr(x2-x1)+sqr(x3-x4));
  end;
  
  function led( a,a1,b1:real):real;
    var p:real;
  begin
    p:=(a+a1+b1)/2;
  end;
  
  function bad( p,a,a1,b1:real):real;
    var s1:real;
  begin
    s1:=sqrt(p*(p-a)*(p-a1)*(p-b1));
  end;
 
var x1,y1,x2,y2,x3,y3,x4,y4,tx,ty:real;
                              i,j:real; 
            s1,s2,s3,s4,s41,s42,p:real;
  a,b,c,d,a1,b1,c1,d1,diag1,diag2:real;
  
begin
  writeln('Введите координаты вершин 4х-угольника в порядке обхода');
  read(x1,y1,x2,y2,x3,y3,x4,y4);
  writeln('Введите координаты точки');
  read(tx,ty);
a:=kor(x2,x1,y2,y1);
b:=kor(x3,x2,y3,y2);
c:=kor(x4,x3,y4,y3);
d:=kor(x4,x1,y4,y1);
 
a1:=kor(x1,i,y1,j);
b1:=kor(x2,i,y2,j);
c1:=kor(x3,i,y3,j);
d1:=kor(x4,i,y4,j);
 
diag1:=kor(x3,x1,y3,y1);
diag2:=kor(x4,x2,y4,y2);
 
{1}p:=led(a,a1,b1);
  s1:=bad(p,p,a1,b1);
    {2}p:=led(b,b1,c1);
        s2:=bad(p,p,b1,c1);
         {3}p:=led(c,c1,d1);
            s3:=bad(p,c,c1,d1);
              {4}p:=led(d,d1,a1);
                  s4:=bad(p,d,d1,a1);
 
p:=led(diag1,a,b);
  s41:=bad(p,diag1,a,b);
     p:=led(diag2,c,d);
      s42:=bad(p,diag2,c,d);
       
 
if kor(tx,ty,x1,x2,x3,x4,y1,y2,y3,y4) then
write('Точка в 4х-угольнике') else write('Точка не в 4х-угольнике')
end.