Даны координаты вершин треугольника, найти центр вписанной окружности

uses crt;
function dlina(ax,ay,bx,by:real):real;//определение длин сторон
begin
dlina:=sqrt(sqr(ax-bx)+sqr(ay-by))
end;
var ax,ay,bx,by,cx,cy,a,b,c,k,x1,y1,x2,y2,xc,yc,k1,b1,k2,b2:real;
begin
repeat
writeln('Введите координаты вершин треугольника');
writeln('Вершина А');
readln(ax,ay);
writeln('Вершина B');
readln(bx,by);
writeln('Вершина C');
readln(cx,cy);
if (ax-cx)*(by-cy)-(bx-cx)*(ay-cy)=0 then
 writeln('Треугольник вырожденный, повторите ввод')
until (ax-cx)*(by-cy)-(bx-cx)*(ay-cy)<>0;
a:=dlina(bx,by,cx,cy);
b:=dlina(ax,ay,cx,cy);
c:=dlina(ax,ay,bx,by);
//найдем координаты концов 2х биссектрис
k:=c/a; //коэффициент пропорциональности
x1:=(cx+k*bx)/(1+k);
y1:=(cy+k*by)/(1+k);
k:=b/a;
x2:=(cx+k*ax)/(1+k);
y2:=(cy+k*ay)/(1+k);
//найдем координаты пересечения биссектрис=координаты центра окружности
if ax=x1 then  //если вдруг первая биссектриса вертикально
 begin
  xc:=x1;//х известен
  k2:=(y2-by)/(x2-bx);//находим коэффициенты второго уравнения
  b2:=(x2*by-bx*y2)/(x2-bx);
  yc:=k2*xc+b2;//координату У
 end
else if x2=x2 then //или вторая
 begin
  xc:=bx;
  k1:=(y1-ay)/(x1-ax);
  b1:=(x1*ay-ax*y1)/(x1-ax);
  yc:=k1*xc+b1;
 end
else //если обе не вертикальные
 begin
  k1:=(y1-ay)/(x1-ax); //коэффициенты 1 уравнения
  b1:=(x1*ay-ax*y1)/(x1-ax);
  k2:=(y2-by)/(x2-bx);//второго
  b2:=(x2*by-bx*y2)/(x2-bx);
  xc:=(b1-b2)/(k2-k1); //координаты центра
  yc:=(k2*b1-k1*b2)/(k2-k1);
 end;
writeln('Координаты центра вписанной окружности X=',xc:0:2,'  Y=',yc:0:2);
end.

uses crt;
var x1,y1,x2,y2,x3,y3:integer;
    a,b,c,k,x11,y11,x12,y12,k1,b1,k2,b2,xc,yc:real;
begin
{вводим координаты с проверкой}
writeln('Введите координаты 1 точки:');
readln(x1,y1);
writeln('Введите координаты 2 точки:');
readln(x2,y2);
writeln('Введите координаты 3 точки:');
readln(x3,y3);
//найдем точки пересечения 2х биссектрис со сторонами
//1. найдем длины  сторон
a:=sqrt(sqr(x1-x2)+sqr(y1-y2));
b:=sqrt(sqr(x2-x3)+sqr(y2-y3));
c:=sqrt(sqr(x1-x3)+sqr(y1-y3));
writeln('a=',a:0:1,' b=',b:0:1,' c=',c:0:1);
//2. найдем координаты концов 2х биссектрис
k:=c/a; //коэффициент пропорциональности
x11:=(x3+k*x2)/(1+k);
y11:=(y3+k*y2)/(1+k);
k:=b/a;
x12:=(x3+k*x1)/(1+k);
y12:=(y3+k*y1)/(1+k);
//найдем координаты пересечения биссектрис=координаты центра окружности
if x1=x11 then  //если вдруг первая биссектриса вертикально
 begin
  xc:=x1;//х известен
  k2:=(y12-y2)/(x12-x2);//находим коэффициенты второго уравнения
  b2:=(x12*y2-x2*y12)/(x12-x2);
  yc:=k2*xc+b2;//координату У
 end
else if x2=x12 then //или вторая
 begin
  xc:=x2;
  k1:=(y11-y1)/(x11-x1);
  b1:=(x11*y1-x1*y11)/(x11-x1);
  yc:=k1*xc+b1;
 end
else //если обе не вертикальные
 begin
  k1:=(y11-y1)/(x11-x1); //коэффициенты 1 уравнения
  b1:=(x11*y1-x1*y11)/(x11-x1);
  k2:=(y12-y2)/(x12-x2);//второго
  b2:=(x12*y2-x2*y12)/(x12-x2);
  xc:=(b1-b2)/(k2-k1); //координаты центра
  yc:=(k2*b1-k1*b2)/(k2-k1);
 end;
writeln('Координаты центра вписанной окружности X=',xc:0:2,'  Y=',yc:0:2)
end.

else if bx=x2 then //èëè âòîðàÿ