Eine etwas andere Idee:
Man kann die Gleichung durch zwei Gleichungen ersetzen:
\((1+a)^2+(1+b)^2=9\) und \(ab=1\), also
\((\frac{1+a}{3})^2+(\frac{1+b}{3})^2=1\) und \(ab=1\).
Wegen "Pythagoras" kann ich dies so substituieren:
\(\sin(\alpha)=\pm \frac{1+a}{3}\) und \(\cos(\alpha)=\pm\frac{1+b}{3}\) und \(ab=1\).
Dies ergibt für die 4 Lösungen die jeweiligen Gleichungen
\(3\sin(\alpha)\cos(\alpha)=+\sin(\alpha)+\cos(\alpha)\)
\(3\sin(\alpha)\cos(\alpha)=+\sin(\alpha)-\cos(\alpha)\)
\(3\sin(\alpha)\cos(\alpha)=-\sin(\alpha)+\cos(\alpha)\)
\(3\sin(\alpha)\cos(\alpha)=-\sin(\alpha)-\cos(\alpha)\)
Hiermit wollte ich keine exakte Lösung liefern, sondern nur auf
eine innere Struktur der Lösungen hinweisen ...