Vamos a mostrar que $(a,b)$ pertenece a una de las dos líneas con las ecuaciones de $b=\sqrt{a}$ e $b=\frac{1}{\sqrt{a}}$ como se muestra en la siguiente figura. Se le dará la respuesta, debido a la simetría del sistema de ecuaciones con respecto al grupo de variables $(a,b)$ vs $(c,d)$. Por otra parte, vamos a establecer (ver (*) en la parte inferior) que la última ecuación es superfluo.
Aquí está la explicación :
Vamos a :
$$S_1:=a+b, \ \ S_2:=c+d, \ \ P_1:=ab, \ \ P_2:=cd$$
El sistema constituido por las tres primeras ecuaciones pueden ser escritas, con estas variables, utilizando clásico transformaciones :
$$\begin{cases}
(A) \ &S_1S_2&=&1& \ &\\
(B) \ &(S_1^2-2P_1)(S_2^2-2P_2)&=&9 & \ \implies \ & (C) \ 1-2(P_1S_2^2+P_2S_1^2)+4(P_1P_2)=9\\
(D) \ &(S_1^3-2P_1S_1)(S_2^3-2P_2S_2)&=&7 & \ \implies \ & (E) \ 1-3S_1S_2(P_1S_2^2+P_2S_1^2)+9(P_1P_2)=7.
\end{casos}$$
(las ecuaciones (C) y (E) son obtenidos por la expansión de (B) y (D) resp., el uso de la relación (A)).
Configuración
$$\alpha := P_1P_2 \ \text{and} \ \beta := P_1S_2^2+P_2S_1^2,$$
las ecuaciones (C) y (E) se convierten en :
$$\begin{cases}
(C) & \ 2\alpha-\beta&=&4\\
(E) & \ 3\alpha-\beta&=&2
\end{casos} \ \ \implica \ \ \alpha=-2 \ \text{y} \ \beta=-8.$$
Utilizando el hecho de que $S_1S_2=1$ e $\alpha=P_1P_2=-2$, la ecuación de $\beta=-8$ se convierte en :
$$P_1 \frac{1}{S_1^2} - \frac{2}{P_1}S_1^2 = -8$$
es decir,
$$(F) \ \ \ \ P_1^2 + 8 P_1S_1^2 - 2 S_1^4 =0,$$
que puede ser considerado como una ecuación cuadrática en la variable $P_1$ da dos soluciones. Debido a la condición de clásico
$$(a+b)^2 \geq 2ab \ \iff \ S_1^2 \geq 2P_1,$$
sólo una de estas soluciones es elegible :
$$P_1=(-4+3\sqrt{2})S_1^2 \ \ \ \iff \ \ \ ab=(-4+3\sqrt{2})(a+b)^2 \ \ \ \iff \ \ \ (b-\sqrt{2}a)(b-\frac{\sqrt{2}}{2}a)=0$$
donde el resultado correspondiente a la figura.
Las ecuaciones paramétricas de las dos líneas son
$$(a,b)=(p,p \sqrt{2}) \ \ \text{and} \ \ (a,b)=(p,p \frac{\sqrt{2}}{2}), \ \ \text{for any} \ \ p \neq 0$$
Debido a la simetría de las ecuaciones, tenemos, también, para cualquier $q \neq 0$ :
$$(c,d)=(q,q \sqrt{2}) \ \ \text{and} \ \ (c,d)=(q,q \frac{\sqrt{2}}{2}).$$
Un rápido vistazo a cualquiera de las cuatro ecuaciones muestran que, necesariamente, $q=\frac{1}{p}$. Nos encontramos de nuevo en este modo, todas las soluciones dadas por @Claude Leibovici y @A. Pongrácz .
(*) De hecho, la cuarta ecuación es una consecuencia de los tres primeros. Aquí es por qué :
Primero de todo, la relación (F) es equivalente a :
$$(G) \ \ \ \ S_1^4=\frac12P_1^2+4P_1S_1^2.$$
Como la cuarta ecuación puede ser escrita :
$$(H) \ \ \ \ (S_1^4+2P_1^2-4P_1S_1^2)(S_2^4+2P_2^2-4P_2S_2^2)=25,$$
el uso de (G) (H), obtenemos :
$$\frac52P_1^2 \frac52P_2^2=25,$$
que es una tautología debido al hecho de que $\alpha=P_1P_2=-2.$