zRotation.py

from pygeo import *

v=display(scale=5, camera_vector=(0,1,-1))

zpoint1=zFreePoint(1,3,color=BLUE)
zpoint2=zFreePoint(-2,-2,color=BLUE)

zrotation1=zRotation(zpoint1,zpoint2,angle=PI/3,color=RED)
zrotation2=zRotation(zpoint1,zpoint2,angle=PI*2/3,color=RED)
#reference objects
usphere=uSphere()

axis1=z_to_uPoint(zpoint1,color=WHITE)
axis2=uAntiPodal(axis1,color=WHITE)

up1=z_to_uPoint(zpoint2,color=WHITE)
up2=z_to_uPoint(zrotation1,color=WHITE)
up3=z_to_uPoint(zrotation2,color=WHITE)

rLine(axis1,axis2,seg=True,color=WHITE)
rLine(usphere.N,axis1,color=WHITE)
rLine(usphere.N,up1,color=WHITE)
rLine(usphere.N,up2,color=WHITE)
rLine(usphere.N,up3,color=WHITE)

uCircle(up1,up2,up3,color=WHITE)

zPlane(color=LIGHTGRAY)
v.pickloop()