We calculate thermodynamic properties of a disordered model insulator, starting from the ideal simple-cubic lattice (g = 0) and increasing the disorder parameter g to » 1/2. As in the earlier Einstein and Debye approximations, ...