Use the code posted at www.elsevierdirect.com to evaluate equations (13.4.56) through (13.4.59) to determine the low-temperature series coefficients for the two-dimensional Ising model for an \(8 \times 8\) lattice. Plot the internal energy and the specific heat as a function of temperature. Repeat your calculation for \(16 \times 16\) and \(32 \times 32\) lattices.