Определить, попадает ли точка в окружность, образованную тремя точками

program lab1zd6;
 
var
  x, y, x0, x1, x2, x3, y0, y1, y2, y3, ma, mb, r, d: real;
 
begin
  writeln('vvedite kordinati (x0,y0)');
  readln(x0, y0);
  writeln('vvedite kordinati (x1,y1)');
  readln(x1, y1);
  writeln('vvedite kordinati (x2,y2)');
  readln(x2, y2);
  writeln('vvedite kordinati (x3,y3)');
  readln(x3, y3);
  //угловые коэфиценты прямых соединяющих точки
  ma := (y2 - y1) / (x2 - x1);
  mb := (y3 - 21) / (x3 - x2);
  //проведём перпендикуляры к этим прямым они будут пересекаться в центре окружности
  x := (ma * mb * (y1 - y3) + mb * (x1 + x2) - ma * (x2 + x3)) / (2 * (mb - ma));
  //подставим x в одно из уравнений перпендикуляров и найдём y
  y := (-1 / ma) * (x - ((x1 + x2) / 2)) + ((y1 + y2) / 2);
  //найдём радиус
  r := sqrt(sqr(x - x1) + sqr(y - y1));
  //теперь найдём растояние от центра окружности до нашей точки
  d := sqrt(sqr(x - x0) + sqr(y - y0));
  if d < r then
    writeln('YES')
  else
    writeln('NO');
  readln;
end.

uses crt;
var  x0,y0,x1,y1,x2,y2,x3,y3,p,x4,y4,x5,y5,a1,b1,c1,a2,b2,c2,xc,yc,r:real;
begin
clrscr;
repeat
writeln('Введите координты 3х точек, лежащих на окружности');
readln(x1,y1,x2,y2,x3,y3);
p:=(x1-x3)*(y2-y3)-(x2-x3)*(y1-y3);
if p=0 then writeln('Точки на 1 прямой, повторите ввод');
until p<>0;
writeln('Введите координаты отдельной точки');
readln(x0,y0);
//середины сторон
x4:=(x1+x2)/2;
y4:=(y1+y2)/2;
x5:=(x1+x3)/2;
y5:=(y1+y3)/2;
//координаты центра и радиус
 //-коэффициенты прямых, перпендикулярных прямым, проходящим через точки 1-2,1-3
a1:=x2-x1;
b1:=y2-y1;
c1:=x4*(x2-x1)+y4*(y2-y1);
a2:=x3-x1;
b2:=y3-y1;
c2:=x5*(x3-x1)+y5*(y3-y1);
 //-координаты точки их перемечения=центр круга
xc:=(c1*b2-c2*b1)/(a1*b2-a2*b1);
yc:=(a1*c2-a2*c1)/(a1*b2-a2*b1);
 //-радиус
r:=sqrt(sqr(x1-xc)+sqr(y1-yc));
if sqr(xc-x0)+sqr(yc-y0)<=r*r then
writeln('Точка принадлежит кругу')
else writeln('Точка не принадлежит кругу')
end.