solve for sqrt((0.5 (a^2 - b^2) + a + 1)^2 + (a b + b)^2) = 1 for b