Brianchon:=proc() local L,t,i,c,d,P:
for i from 0 to 5 do
L[i]:=
TangentToEllipse(c,d,t[i]):
od:
for i from 0 to 5 do
P[i]:=
Pt
(L[i],L[i+1 mod 6]):
od:
Concurrent
(
Le (P[0],P[3]) ,
Le
(P[1],P[4]) ,
Le (P[2],P[5])
):
end: