(a - x)·x + a·y = a + b
Wir können einfach nach y auflösen
y = x2/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)·(x2/a - x + b/a + 1) = a + 2·b
b·x2/a + x2 - b·x + a + b2/a + 2·b = a + 2·b
b·x2/a + x2 - b·x + b2/a = 0
(x2·(a + b) - a·b·x + b2)/a = 0
x2·(a + b) - a·b·x + b2 = 0
Das löst man mit der quadratischen Lösungsformel
x = b·(a - √(a2 - 4·a - 4·b))/(2·(a + b))
x = b·(√(a2 - 4·a - 4·b) + a)/(2·(a + b))