How to solve a system of 3 trigonometric equations How to solve  system of three trigonometric equations:
$(\sin x)^2 (\cos y)^2 = 4 \cos x \sin y\tag1$
$(\sin y)^2 (\cos z)^2  = 4 \cos y \sin z \tag2$
$1- \sqrt{\sin z}(1+\sqrt{\cos x})=\sqrt{\frac{1-\sin y}{1+\sin y}}\tag3$
and to verify that
$\sin x=\sqrt{2}(\sqrt{2}+1)\sqrt{\sqrt{10}-3}(\sqrt{5}-2)\\
\sin y=(\sqrt{2}-1)^2(\sqrt{10}-3)\\
\sin z=(2\sqrt{2}+\sqrt{5}-\sqrt{12+4\sqrt{10}})^2(\sqrt{5}+\sqrt{2}-\sqrt{6+2\sqrt{10}})^2$
not is the only solution, but that many others exist?
The simplest solution is:
$\sin x=2^{5/4}\big(\sqrt{2}-1\big)$
$\sin y=\frac{1}{\sqrt{2}}$
$\sin z=3-2\sqrt{2}$.
Three other solutions are:
1)
$\sin x=\sqrt{2\sqrt{2}-2}$
$\sin y=\sqrt{2}-1$
$\sin z=5+4\sqrt{2}-\sqrt{56+40\sqrt{2}}$,
2)
$\sin x=\frac{\sqrt{2}}{2}$
$\sin y=3-2\sqrt{2}$
$\sin z=\big(\sqrt{2}+1\big)^{2}\big(2^{1/4}-1\big)^{4}$
3)
$\sin x=\frac{2^{5/4}a^{2}c. d^{6}\varphi^{3}}{b^{2}}$
$\sin y=\frac{d^4}{\sqrt{2}\varphi^{4}}$
$\sin z=a^{4}b^{4}c^{2}\varphi^{6}$
where
$$a=\sqrt{2}-1$$
$$b=5^{1/4}-\sqrt{2}$$
$$c=\sqrt{10}-3$$
$$d=\frac{5^{1/4}-1}{\sqrt{2}}$$
$\varphi$ is the golden ratio.

This is based on modular equations of degree 2. If $k, l$ are elliptic moduli with $l$ of degree 2 over $k$ then we have $$l=\frac{1-k'}{1+k'}\tag{1}$$ where $k'=\sqrt{1-k^2}$.
One should observe that the relation between $\sin x, \sin y$ is such that $\sin y$ is of degree 2 over $\sin x$. This is easily verified by writing $k, l$ in place of $\sin x, \sin y$ in the equation given in question. Doing so we get $$k^2(1-l^2)=4k'l$$ This gives us a quadratic equation in $l$ as $$k^2l^2+4k'l-k^2=0$$ or $$l=\frac{-2k'\pm\sqrt{4k'^2+k^4}}{k^2}$$ Choosing $+$ sign as $l>0$ we get $$l=\frac{2-k^2-2k'}{k^2}=\frac{(1-k')^2}{1-k'^2}=\frac{1-k'}{1+k'}$$ as in equation $(1)$.
Similarly $\sin z$ of degree 2 over $\sin y$ and hence $\sin z$ is of degree 4 over $\sin x$. Let $l_1=\sin z$ so that $$l_1=\frac{1-l'}{1+l'}$$ Next we have $$1+\sqrt{\cos x} =1+\sqrt{k'}=1+\sqrt{\frac{1-l}{1+l}}$$ The equation to be proved  can now be rewritten as $$\sqrt{\sin z} =\dfrac{1-\sqrt{\dfrac{1-l}{1+l}}} {1+\sqrt{\dfrac{1-l} {1+l} } } $$ or $$l_1=\sin z=\frac{2-2\sqrt{1-l^2} } {2+2\sqrt{1-l^2}}=\frac{1-l'}{1+l'} $$ as expected.
The given numerical values for the moduli are taken from a table of singular moduli. In particular the value $$\sin y=(\sqrt{10}-3)(\sqrt{2}-1)^2$$ is $k_{10}$ where $$k_n=\frac{\vartheta_2^2(e^{-\pi\sqrt{n}})}{\vartheta_3^2(e^{-\pi\sqrt{n}})} $$ and $$\vartheta_2(q)=\sum_{r\in\mathbb{Z}}q^{(r+(1/2))^2},\,\vartheta_3(q)=\sum_{r\in\mathbb {Z}} q^{r^2} $$ And therefore the values of $\sin x, \sin z$ are $k_{5/2},k_{40}$ respectively. 
