CentralProjection.py

from pygeo import *

v=display(scale=70)

p1 = FreePoint(-33.,-33,-3,color=WHITE)
p2 = FreePoint(12,-15,4,color=WHITE)
p3 = FreePoint(40,-3,0,color=WHITE)
plane=Plane(p1,p2,p3,color=BLUE,level=2)
point = FreePoint(2,36,1,color=BLUE,pointsize=1.5)

p1a = FreePoint(-3.,-3,13,color=WHITE)
p2a = FreePoint(1,11,-11,color=WHITE)
p3a = FreePoint(17,-6,7,color=WHITE)

circle=CircumCircle(p1a,p2a,p3a,color=BLUE)
CentralProjection(plane,point,[circle],color=RED)

pa=PointArray(circle,color=WHITE,density=200)
LineArray(pa,point,color=WHITE)
v.pickloop()