Carnot:=proc() local m,n,T,O,O1,O2,O3,R,r:
T:=Te
(m,n): O:=Circumcenter(T):
R:=Circumradius(T):
r:=Inradius(T):
O1:=Ft
(O , Le
(T[2],T[3])):
O2:=
Ft
(O ,
Le
(T[1],T[3])):
O3:=
Ft
(O ,
Le(T[1],T[2])):
SumSqRtsIsZero
(
[
DeSq(O,O1),
DeSq(O,O2),DeSq
(O,O3),R^2,r^2]):
end: