Morley's Theorem

# W.l.o.g. A=[0,0] and B=[1,0]

Morley:=proc() local a,b,c,A,B,C,D,E,F,TriA1,TriA2,AC,BC, TriB1,TriB2,TriC1,TriC2:

c:=Pi/3-a-b: A:=[0,0]: B:=[1,0]: AC:=LePtAngle(A,3*a): BC:=LePtAngle(B,-3*b): C:=Pt(AC,BC):

TriA1:=LePtAngle(A,a): TriA2:=LePtAngle(A,2*a): TriB1:=LePtAngle(B,-b):
TriB2:=LePtAngle(B,-2*b): TriC1:=LePtAngle(C,3*a+c): TriC2:=LePtAngle(C,3*a+2*c):

D:= Pt (TriA1,TriB1): E:= Pt (TriA2,TriC1): F:= Pt(TriB2,TriC2):

end: