Derivation of the central equations for the square wave potential.
Maple solutionof the above for k=0, gamma point.
Maple solution for k= edge of Brillouin Zone.
Maple solution for half way between center and zone boundary.