I have the following system of equations:
$$ \begin{cases} x_0=\cos\theta_1\cos\theta_0 \\[4pt] x_1=\sin\theta_1\cos\theta_0 \\[4pt] x_2=\cos\theta_2\sin\theta_0 \\[4pt] x_3=\sin\theta_2\sin\theta_0 \end{cases} $$
The $x_i$ are known and $\sum_i x_i^2 =1$, while the $\theta_i$ are between -$\pi$ and $\pi$ but not known. How can they be found?
It is easy to get $\theta_1$ and $\theta_2$ by using $\arcsin$ and $\arccos$, but both depend on $\theta_0$, and I do not know how to continue. I have also tried using SymPy and Wolfram Alpha, but they were not successful.
Assuming $x_0x_2\neq 0$, dividing the second one by the first one and the fourth one by the third one we obtain
$$\tan \theta_1=\frac{x_1}{x_0} \implies \theta_1=\arctan\left(\frac{x_1}{x_0}\right)+k_1\pi$$
$$\tan \theta_2=\frac{x_3}{x_2} \implies \theta_2=\arctan\left(\frac{x_3}{x_2}\right)+k_2\pi$$
then we can plug these values in the equations to obtain the right value for $k_1\in\mathbb Z$, $k_2\in\mathbb Z$ and $\theta_0$.