Make a plot of E vs. q for the band structure in Section 5.3.2. Use α = 1 (in units where m=\hbar=a=1 ). Hint: In Mathematica, ContourPlot will graph E(q) as defined implicitly by Equation 5.71. On other platforms the plot can be obtained as follows:
Choose a large number (say 30,000) of equally-spaced values for the energy in the range E = 0 and E = 30.
For each value of E, compute the right-hand side of Equation 5.71. If the result is between -1 and 1 , solve for q from Equation 5.71 and record the pair of values \{q, E\} \text { and }\{-q, E\} (there are two solutions for each energy).
You will then have a list of pairs \{\{q 1, E 1\},\{q 2, E 2\}, \ldots\} which you can plot.
\cos (q a)=\cos (k a)+\frac{m \alpha}{\hbar^{2} k} \sin (k a) (5.71).