Euler's Triangle Formula

# d^2=R^2-2*r*R

EulerTriangleFormula:=proc() local T,m,n,A,B,C,d,R,r,O,I:

T:= Te (m,n): A:=T[1]: B:=T[2]: C:=T[3]: I:= Incenter (m,n): O:=Circumcenter (A,B,C):
r:=Inradius(m,n): R:= Circumradius(A,B,C): d:=sqrt(DeSq (O,I)):

ItIsZero ((R^2-d^2)^2-(2*r*R)^2):

end: