We would like to share the little Python notebook for solving simple Hückel molecular orbital (HMO) of π-conjugated hydrocarbons. The script can generate the Hückel matrix from a molecular SMILES string, and the eigenproblem is then solved numerically or symbolically. The former approach is recommended, as the symbolic solution in terms of α (Coulomb integrals) and β (resonance integrals) can look quite daunting. The script was prepared when analyzing the electronic effect of thionation (DOI: 10.1039/D2CP05186A); see also the blog post for more discussion.