Hückel theory with SymPy: exact orbital energies of butadiene and benzene

Afterwards you can compute exact Hückel energies, delocalization energies, and bond orders with SymPy's eigenvals and eigenvects, and prove the 4n + 2 rule.