ToolPythonIntermediate Hückel theory with SymPy: exact orbital energies of butadiene and benzene print(sp.Poly(cubic, lam).nroots(n=4))
ToolPythonAdvanced Von Neumann stability analysis with SymPy: the leapfrog Courant limit return sp.collect(sp.Poly(expr, G).monic().as_expr(), G)