(a - x)·x + a·y = a + b
Wir können einfach nach y auflösen
y = x^2/a - x + b/a + 1
Das setzt man für y in die 2. Gleichung ein
a·x + (a + b)·y = a + 2·b
a·x + (a + b)·(x^2/a - x + b/a + 1) = a + 2·b
b·x^2/a + x^2 - b·x + a + b^2/a + 2·b = a + 2·b
b·x^2/a + x^2 - b·x + b^2/a = 0
(x^2·(a + b) - a·b·x + b^2)/a = 0
x^2·(a + b) - a·b·x + b^2 = 0
Das löst man mit der quadratischen Lösungsformel
x = b·(a - √(a^2 - 4·a - 4·b))/(2·(a + b))
x = b·(√(a^2 - 4·a - 4·b) + a)/(2·(a + b))