v unittestu ty vzorce mohu otestoval nějak tak
from sympy import solve, symbols, var, sqrt, simplify
a, b, c = symbols('a b c')
x = var('x')
def roots():
return [simplify((-b + sqrt(-4*a*c + b**2))/(2*a)),
simplify(-(b + sqrt(-4*a*c + b**2))/(2*a))]
def test():
assert roots() == solve(a*x**2 + b*x + c, x)
chci vidět, jak to stejné otestuješ typovým systémem