Por qué no utilizar la separación de variables desde el principio. Introduzca $u=X(x)Y(y)$ para conseguir
$$x^2X'^2Y^2+y^2X^2Y'^2=X^2Y^2 \implies \frac{x^2X'^2}{X^2}+\frac{y^2Y'^2}{Y^2}=1$$
Está claro que la ecuación sólo puede ser cierta si $\frac{x^2X'^2}{X^2}$ y $\frac{y^2Y'^2}{Y^2}$ son constantes tales que su suma es igual a 1. Así obtenemos dos ecuaciones diferenciales ordinarias:
$$X'(x)^2=k^2\frac{X^2(x)}{x^2}$$ $$Y'(y)^2=(1-k^2)\frac{Y(y)^2}{y^2}$$
Ahora tomamos la raíz cuadrada de ambas expresiones
$$X'(x)=\pm k\frac{X(x)}{x} \implies X(x)=c_1x^{\pm k}$$ $$Y'(y)=\pm \sqrt{1-k^2}\frac{Y(y)}{y} \implies Y(y)=c_2y^{\pm \sqrt{1-k^2}}.$$
Ahora, la solución por separación de variables viene dada por $$u(x,y)=c_1c_2x^{\pm k}y^{\pm \sqrt{1-k^2}}.$$
Tenga en cuenta que no puede superponer diferentes soluciones para $k$ ya que la ecuación diferencial no es lineal.
EDITAR : Si quieres resolver
$$f'(x)^2=\frac{\lambda^2}{x^{2}}\implies f'(x)=\pm\frac{\lambda}{x} \implies f(x)=\pm \lambda\ln x +c_1$$
$$g'(y)^2=\frac{1-\lambda^2}{y^2} \implies g'(y)=\pm\frac{\sqrt{1-\lambda^2}}{y} \implies g(y)=\pm\sqrt{1-\lambda^2}\ln y +c_2$$
Verá que $|\lambda|\leq 1$ es una restricción lateral tal que las soluciones siguen siendo reales.