How can simultaneous sinusoidal equations be solved? I have come across a set of simultaneous equations which I can't figure out how to solve. I have 3 equations and only two unknowns, but they are angular quantities and feature in the equations as sinusoidal functions of the angular quantities.
The system of equations is:
$$ \begin{Bmatrix}\cos\psi\sin\theta+\cos\theta\sin\phi\sin\psi \\ \sin\psi\sin\theta-\cos\psi\cos\theta\sin\phi \\ \cos\phi\cos\theta \end{Bmatrix} = \begin{Bmatrix} x \\ y \\ z \end{Bmatrix} $$
Where: x, y, z and psi are known and phi and theta are unknown.
Is it possible to rearrange these equations to solve for theta and phi using the other terms?

Note that
$$
\begin{array}{l}
 \left( {\begin{array}{*{20}c}
   x  \\
   y  \\
   z  \\
\end{array}} \right) = \left( {\begin{array}{*{20}c}
   {\cos \psi \sin \theta  + \cos \theta \sin \phi \sin \psi }  \\
   {\sin \psi \sin \theta  - \cos \psi \cos \theta \sin \phi }  \\
   {\cos \phi \cos \theta }  \\
\end{array}} \right) =  \\ 
  = \left( {\begin{array}{*{20}c}
   {\sin \theta \cos \psi  + \cos \theta \sin \phi \sin \psi }  \\
   {\sin \theta \sin \psi  - \cos \theta \sin \phi \cos \psi }  \\
   {\cos \theta \cos \phi }  \\
\end{array}} \right) =  \\ 
  = \sin \theta \left( {\begin{array}{*{20}c}
   {\cos \psi }  \\
   {\sin \psi }  \\
   0  \\
\end{array}} \right) + \cos \theta \left( {\begin{array}{*{20}c}
   {\sin \phi \sin \psi }  \\
   { - \sin \phi \cos \psi }  \\
   {\cos \phi }  \\
\end{array}} \right) =  \\ 
  = \sin \theta \left( {\begin{array}{*{20}c}
   {\cos \psi }  \\
   {\sin \psi }  \\
   0  \\
\end{array}} \right) + \cos \theta \left( {\begin{array}{*{20}c}
   {\sin \phi \cos \left( {\psi  + \pi /2} \right)}  \\
   {\sin \phi \sin \left( {\psi  + \pi /2} \right)}  \\
   {\cos \phi }  \\
\end{array}} \right) =  \\ 
  = \sin \theta \;{\bf u} + \cos \theta \;{\bf v} \\ 
 \end{array}
$$
where $\bf u$ and $\bf v$ are  given in spherical coordinates, and can be seen to be
unitary vectors and orthogonal to each other.
Their combination with $sin \theta$ and $cos \theta$ is therefore a circle in their plane.
For the system to have solutions, the vector $\bf p=(x,y,z)$ shall be on the unitary sphere.
That given, there will be multiple solutions for the three angles, corresponding to the 
infinite main circles passing through $\bf p$.
