Viviani's Theorem

Viviani:=proc() local A,B,C,A1,B1,C1,C2,T,P,f,X,Y:

T:=Te(1/sqrt(3),1/sqrt(3)):   A:=T[1]:   B:=T[2]:   C:=T[3]:

A1:=Ft(P,Le(A,B)):   B1:=Ft(P,Le(A,C)):   C1:=Ft(P,Le(B,C)):   C2:=Ft(C,Le(A,B)):

SumSqRtsIsZero([ DeSq(P,A1), DeSq(P,B1), DeSq(P,C1),DeSq(C,C2)]):

end: