Band structure with plane waves in NumPy: an electron in a 1D crystal

Afterwards you can compute the energy bands of a 1D periodic potential in a plane-wave basis with eigh and check the gaps against nearly-free-electron theory.