CirclePole.py

from pygeo import *

v=display(scale=70)

p1 = FreePoint(-33.,-33,-3,color=WHITE)
p2 = FreePoint(-23,33,4,color=WHITE)
p3 = FreePoint(12,9,-14,color=WHITE)

circle=CircumCircle(p1,p2,p3,color=BLUE)
plane=Plane(p1,p2,p3,level=3)
p4 = Slider(plane,1,9,-1,color=WHITE)
p5 = Slider(plane,1,7,-1,color=WHITE)
line=Line(p4,p5,color=BLUE)

cp=CirclePole(line,circle,color=RED,pointsize=1.5)


#check against CirclePolar algorithm
CirclePole(CirclePolar(circle,cp,level=2),circle,level=2,color=GREEN,pointsize=1.5)

v.pickloop()