for
a = 1 b = 0 c = -4 disc = b**2 - 4*a*c print(disc) x1 = (-b + disc**0.5)/(2*a) x2 = (-b - disc**0.5)/(2*a) print(x1, x2)