a = 1 b = -9.5 c = 15 D = b**2 - 4*a*c x1 = (-b - D**0.5) / (2*a) x2 = (-b + D**0.5) / (2*a) print('x1 =', x1) print('x2 =', x2)