M = (0, 0, 0)
kx = Sphere(M, 4)
B = PointIn(kx)
r = Line(M, B)
T = OrthogonalPlane(B, r)
c = Circle(xAxis, B)
d = Circle(yAxis, B)
rx = OrthogonalLine(B, xAxis, space)
Mx = Intersect(rx, xAxis)
ax = OrthogonalPlane(Mx, xAxis)
A = PointIn(ax)
α = Angle(Mx, B, A)
tx = Line(B, A)
ry = OrthogonalLine(B, yAxis, space)
My = Intersect(ry, yAxis)
ay = OrthogonalPlane(My, yAxis)
C = PointIn(ay)
β = Angle(My, B, C)
ty = Line(C, B)