This is a comment that’s too long for the usual format.
LAcarguy’s method generalizes to any values $a,b,c$ instead
of $\frac14,\frac15,\frac16$.
If we put
$$
\phi(a,b,c)=ab+ac-bc, \ \psi(a,b,c)=ab+ac+bc \tag{1}
$$
and
$$
G(a,b,c)=\phi(a,b,c)\phi(b,c,a)\phi(c,a,b)\psi(a,b,c) \tag{2}
$$
then LAcarguy’s method leads to the following explicit formulas for the
solutions :
$$
x=\frac{\varepsilon}{\phi(a^2,b^2,c^2)}\sqrt{G(a,b,c)},\
y=\frac{\varepsilon}{\phi(b^2,c^2,a^2)}\sqrt{G(a,b,c)},\
z=\frac{\varepsilon}{\phi(c^2,a^2,b^2)}\sqrt{G(a,b,c)} \tag{3}
$$
where $\varepsilon$ is $+1$ or $-1$.
In your initial question where $a=\frac14,b=\frac15,c=\frac16$, one finds
$G=\frac{7}{921600}$,
$$
x=\varepsilon\frac{\sqrt{7}}{3}, \
y=\varepsilon\frac{5\sqrt{7}}{9}, \
z=3\varepsilon\sqrt{7} \tag{4}
$$
$\frac{\sin(\theta_1)}{4}=\frac{\sin(\theta_2)}{5}=\frac{\sin(\theta_3)}{6} \implies$
$15\sin(\theta_1)=12\sin(\theta_2)=10\sin(\theta_3)$ Not sure where to take it from there though. I haven't messed with trig functions for a while.
– ruler501 Mar 27 '14 at 06:34