Student Project by Christian Neubauer and Gernot Kraberger (2010)
The form below generates a table of where the first column is the angular frequency ω in rad/s and the second column is the density of states D(ω) in units of s/(rad m³).

The Matlab files used to calculate the data are base.m, kvals.m, disp_solver.m and analy.m.
The density of states can be used to calculate the temperature dependence of thermodynamic quantities.