b: x^(2) / 2 + y^(2) / 2 + z^(2) = 1
A = PointIn(b)
c = Plane(A, zAxis)
d = IntersectPath(c, b)
f = Tangent(A, d)
g = OrthogonalLine(A, f, space)
g_1 = OrthogonalLine(A, f, c)
α = Angle(g_1, Plane[xAxis, yAxis])
β = Angle(c, Plane[xAxis, zAxis])
O1 = Point({0.2, 0.2, 0.2})
X1 = Line(O1, xAxis)
Y1 = Line(O1, yAxis)
Z1 = Line(O1, zAxis)
Ecuador1 = Plane(X1, Y1)
Meridiano1 = Plane(A, Z1)
Datum1: (0.6 * (x - x(O1))^(2)) + (0.6 * (y - y(O1))^(2)) + (z - z(O1))^(2) = 1
CurvaMeridiano1 = IntersectPath(Meridiano1, Datum1)
lon1 = Angle(Meridiano1, Plane[X1, Z1])
A1 = ClosestPoint(CurvaMeridiano1, A)
T1 = Tangent(A1, CurvaMeridiano1)
Per1 = OrthogonalLine(A, T1, Meridiano1)
lat1 = Angle(Per1, Ecuador1)