Дан треугольник по заданным координатам своих вершин.Определить радиус окружности возле треугольника

function radius(x1,y1,x2,y2,x3,y3:real):real;
var x4,y4,x5,y5,a1,b1,c1,a2,b2,c2,xc,yc:real;
begin
//середины сторон
x4:=(x1+x2)/2;
y4:=(y1+y2)/2;
x5:=(x1+x3)/2;
y5:=(y1+y3)/2;
//координаты центра
a1:=x2-x1; //коэффициенты уравнения 1 прямой
b1:=y2-y1;
c1:=x4*(x2-x1)+y4*(y2-y1);
a2:=x3-x1; //коэффициенты уравнения 2 прямой
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);
radius:=sqrt(sqr(x1-xc)+sqr(y1-yc));//радиус
end;
var x1,y1,x2,y2,x3,y3:real;
begin
writeln('Введите координаты 1 вершины');
readln(x1,y1);
writeln('Введите координаты 2 вершины');
readln(x2,y2);
writeln('Введите координаты 3 вершины');
readln(x3,y3);
write('Радиус описанной окружности=',radius(x1,y1,x2,y2,x3,y3):0:2)
end.

var
  x1, y1, x2, y2, x3, y3, a, b, c, p, r,s: real;
function rad(a,b,c,p: real): real;
begin
  rad := (a+b+c/4*s);
end;
begin
  write('A='); readln(x1, y1);
  write('B='); readln(x2, y2);
  write('C='); readln(x3, y3);
  a := sqrt(sqr(x1 - x2) + sqr(y1 - y2)); //длины сторон
writeln('a=',a:5:2);
  b := sqrt(sqr(x2 - x3) + sqr(y2 - y3));
writeln('b=',b:5:2);
  c := sqrt(sqr(x1 - x3) + sqr(y1 - y3));
writeln('c=',c:5:2);
s:=sqrt(p*(p-a)*(p-b)-(p-c));
writeln('Ploshad=',s:5:2);
p := (a + b + c) / 2;
writeln('Perimetr=',p:5:2);
  r := rad(a,b,c,s);
  writeln('ответ : ', r:5:5);
end.

var  x1, y1, x2, y2, x3, y3, a, b, c, p, r,s: real;
function rad(a,b,c:real): real;
var p,s:real;
begin
p:=(a+b+c)/2;
s:=sqrt(p*(p-a)*(p-b)*(p-c));
rad:=a*b*c/(4*s);
end;
begin
writeln('Введите координаты вершины A=');
readln(x1, y1);
writeln('Введите координаты вершины B=');
readln(x2, y2);
writeln('Введите координаты вершины C=');
readln(x3, y3);
a := sqrt(sqr(x1 - x2) + sqr(y1 - y2)); //длины сторон
b := sqrt(sqr(x2 - x3) + sqr(y2 - y3));
c := sqrt(sqr(x1 - x3) + sqr(y1 - y3));
r := rad(a,b,c);
writeln('ответ : ', r:0:5);
end.

function dlina(x1,y1,x2,y2:real):real;
begin
dlina:=sqrt(sqr(x1-x2)+sqr(y1-y2));
end;
function plosh(x1,y1,x2,y2,x3,y3:real):real;
begin
plosh:=0.5*abs(((x1-x3)*(y2-y3)-(x2-x3)*(y1-y3)));
end;
procedure rad(a,b,c,s:real;var r:real);
begin
r:=a*b*c/(4*s);
write('Ответ:',r:0:5);
end;
var x1,y1,x2,y2,x3,y3,r:real;
begin
writeln('Введите координаты вершины A=');
readln(x1,y1);
writeln('Введите координаты вершины B=');
readln(x2,y2);
writeln('Введите координаты вершины C=');
readln(x3,y3);
rad(dlina(x1,y1,x2,y2),dlina(x2,y2,x3,y3),dlina(x1,y1,x3,y3),plosh(x1,y1,x2,y2,x3,y3),r)
end.