a = a1b1 - a2b2

b = 2b1b2(a1 + a2)

c = b1b2(a1b2-a2b1)

ax^2 + bx +c = 0

x = -b1b2(a1 + a2) / (a1b1 - a2b2) ± b1b2√(a1^2 + a2^2 + 2a1a2 - a1b1b2 - a2b1b2) / (a1b1 - a2b2)