La traducción de la ecuación a la forma
$$(xy)^2+(yz)^2+(zx)^2=15xyz,\qquad(1)$$
fácil ver que $xyz>0.$
Eso significa que todas las soluciones con diferentes signos pueden ser obtenidos a partir de soluciones positivas del signo colocación de dos incógnitas al mismo tiempo.
Así que podemos buscar soluciones positivas,
$$x,y,z>0.$$
Todas las raíces racionales de la ecuación con coeficientes enteros
$$\left(\dfrac{yz}x\right)^2-15\dfrac{yz}x+y^2+z^2=0$$
tiene que ser entero, por lo $x\ |\ yz$ y de manera similar a $y\ |\ zx,\ z\ |\ xy.$
Vamos
$$u=\dfrac{xy}z,\quad v=\dfrac{yz}x,\quad w=\dfrac{zx}y,$$
a continuación, obtenemos
$$u+v+w=15,\quad u,v,w,\sqrt{uv},\sqrt{vw},\sqrt{wu}\in\mathbb N$$
con la única solución
$$u=v=w=5,$$
$$x=\sqrt{uw}=5,\quad y=\sqrt{uv}=5,\quad z=\sqrt{vw}=5.$$
Por lo que el conjunto solución es
$$\boxed{(x,y,z)\in\{(5,5,5),(5,-5,-5),(-5,5,-5),(-5,-5,5)\}.}$$
Nota
Para la ecuación de $LHS=9$ hemos
$$u+v+w=9,\quad u,v,w,\sqrt{uv},\sqrt{vw},\sqrt{wu}\in\mathbb N$$
con las soluciones
$$(u,v,w)=\{(4,4,1), (4,1,4), (1,4,4), (3,3,3)\},$$
así que las soluciones positivas de $LHS=9$
$$(x,y,z)=\{(2,4,2),(4,2,2),(2,2,4),(3,3,3)\}$$
y con la cuenta a la par de signos de permutaciones tenemos 16 soluciones.