Solución de la ecuación cuadratica (solo raíces reales)

Cuadratica.py

# Módulo para funciones matematicas en los reales
# https://docs.python.org/3/library/math.html
import math

a = float(input("Ingrese a: "))
b = float(input("Ingrese b: "))
c = float(input("Ingrese c: "))

# Comprobar la precondición
if b**2-4*a*c>=0:
    x1 = (-b+math.sqrt(b**2-4*a*c))/(2*a)
    x2 = (-b-math.sqrt(b**2-4*a*c))/(2*a)
    print('Primera raiz ', x1)
    print('Segunda raiz ', x2)

    # Comprobar las postcondiciones
    assert math.fabs(a*x1**2 + b*x1 + c) < 0.000000001, "La primera raiz no cumple"
    assert math.fabs(a*x2**2 + b*x2 + c) < 0.000000001, "La segunda raiz no cumple"

else:
    print('Ecuación sin raices reales')