I am simply trying to find the values of the two variables X and Y using the following system of two equations, in Python:
z = arccos [(sin(a) * sin(X)) + (cos(a) * cos(X) * cos(Y - b))]
Az = arctan2[(cos(X) * sin(Y - b)) , (cos(a) * sin(X)) - (sin(a) * cos(X) * cos(Y-b))]
where the values for the variables z, Az, a, and b are already known. The only unknown variables are X and Y.
How can I do this, using Python?