A = (4.157210710474489, 3.1445713919670464, 4)
B = Intersect(xAxis, zAxis)
a = Sphere(B, A)
C = PointIn(a)
X = Ray(B, A)
Y = Ray(B, C)
α = Angle(A, B, C)
p = Plane(C, B, A)
N = IntersectPath(p, xOyPlane)
c = IntersectPath(xOyPlane, a)
D = Point(c)
E = Intersect(N, c, 2)
Z = Ray(B, D)
β = Angle(D, B, E)
F = Intersect(Z, yAxis)
G = Point(xAxis)
f = Line(F, G)
H = Point(N)
ϕ = Angle(G, F, H)
θ = Angle(zAxis, Z)
I = Point(X)
J = Point(N)
ψ = Angle(I, F, J)