Треугольник, биссектрисы, вписанная окружность

program triangle;
uses Graph;
const
h=1e-4;
x=13;
y=8;
type
tri=record
x,y:real;
end;
long=record
a,b,c:real;
end;
straight=array[1..3] of real;
mas=array[1..3] of tri;
lo=array[1..3] of long;
 
 
procedure SetScreen (var mx,my:real; var x0,y0:integer);
begin
mx:=600/(x*2);
my:=420/(y*2);
x0:=GetMaxX div 2;
y0:=GetMaxY div 2;
end;
 
procedure countall (p:mas; var R:real; var p1:mas);
 
function D(x1,y1,x2,y2:real):real;
begin
D:=sqrt(sqr(x2-x1)+sqr(y2-y1));
end;
 
function Radius (a,b,c,S:real):real;
begin
Radius:=S/((a+b+c)/2);
end;
 
function Geron (a,b,c:real):real;
var
p:real;
begin
p:=(a+b+c)/2;
Geron:=sqrt(p*(p-a)*(p-b)*(p-c));
end;
var
S,k:real;
i,j:integer;
l:straight;
begin
l[1]:=d(p[1].x,p[1].y,p[2].x,p[2].y);
l[2]:=d(p[2].x,p[2].y,p[3].x,p[3].y);
l[3]:=d(p[3].x,p[3].y,p[1].x,p[1].y);
k:=l[3]/l[1];
p1[1].x:=(p[3].x+(k*p[2].x))/(1+k);
p1[1].y:=(p[3].y+(k*p[2].y))/(1+k);
k:=l[2]/l[1];
p1[2].x:=(p[3].x+(k*p[1].x))/(1+k);
p1[2].y:=(p[3].y+(k*p[1].y))/(1+k);
k:=l[3]/l[2];
p1[3].x:=(p[1].x+(k*p[2].x))/(1+k);
p1[3].y:=(p[1].y+(k*p[2].y))/(1+k);
S:=Geron(l[1],l[2],l[3]);
R:=Radius(l[1],l[2],l[3],S);
end;
 
procedure Xcircle (x1,y1,r:real);
var
mx,my,a,b,i:real;
x0,y0,px,py:integer;
begin
SetScreen(mx,my,x0,y0);
i:=0;
while i<=360 do begin
a:=r*cos(i)+x1;
b:=r*sin(i)+y1;
px:=x0+round(a*mx);
py:=y0-round(b*my);
putpixel(px,py,3);
i:=i+1e-3;
end;
px:=x0+round(x1*mx);
py:=y0-round(y1*my);
putpixel(px,py,2);
end;
 
var
gm,gd:integer;
i:integer;
x0,y0:integer;
px,py,px1,py1:integer;
x1,x2,x_1,y_1:real;
mx,my:real;
R:real;
p,p1:mas;
s:string;
n:lo;
BEGIN
writeln('Coordinates :');
for i:=1 to 3 do
readln(p[i].x,p[i].y);
gd:=detect;
initgraph(gd,gm,'');
setbkcolor(7);
SetScreen(mx,my,x0,y0);
setcolor(blue);
setlinestyle(0,1,1);
line(0,y0,GetMaxX,y0);
line(x0,0,x0,GetMaxY);
setlinestyle(1,0,0);
for i:=1 to x*2 do begin
setcolor(9);
line(0,y0-round(i*my),GetMaxX,y0-round(i*my));
line(0,y0+round(i*my),GetMaxX,y0+round(i*my));
line(x0-round(i*mx),0,x0-round(i*mx),GetMaxY);
line(x0+round(i*mx),0,x0+round(i*mx),GetMaxY);
setcolor(blue);
if i<=y then begin
str(i,s);
OutTextXY(x0-10,y0-round(i*my),s);
OutTextXY(x0-20,y0+round(i*my),'-'+s);
end;
str(i,s);
OutTextXY(x0+round(i*mx),y0+5,s);
OutTextXY(x0-round(i*mx)-10,y0+5,'-'+s);
end;
setcolor(Yellow);
OutTextXY(x0-10,y0+5,'0');
OutTextXY(x0-10,5,'Y');
OutTextXY(GetMaxX-12,y0-10,'X');
setcolor(red);
setlinestyle(0,0,1);
for i:=1 to 3 do begin
px:=x0+round(p[i].x*mx);
py:=y0-round(p[i].y*my);
px1:=x0+round(p[i mod 3+1].x*mx);
py1:=y0-round(p[i mod 3+1].y*my);
line(px,py,px1,py1);
end;
CountAll(p,R,p1);
for i:=1 to 3 do begin
px:=x0+round(p[i].x*mx);
py:=y0-round(p[i].y*my);
px1:=x0+round(p1[i].x*mx);
py1:=y0-round(p1[i].y*my);
line(px,py,px1,py1);
end;
n[1].a:=p[1].y-p1[1].y;
n[1].b:=p1[1].x-p[1].x;
n[1].c:=-(n[1].a*p[1].x+n[1].b*p[1].y);
n[2].a:=p[2].y-p1[2].y;
n[2].b:=p1[2].x-p[2].x;
n[2].c:=-(n[2].a*p[2].x+n[2].b*p[2].y);
x_1:=-((n[1].c*n[2].b)-(n[2].c*n[1].b))/((n[1].a*n[2].b)-(n[2].a*n[1].b));
y_1:=-((n[1].a*n[2].c)-(n[2].a*n[1].c))/((n[1].a*n[2].b)-(n[2].a*n[1].b));
Xcircle(x_1,y_1,R);
readln
END.