  Thermally excited waves in a Plasma crystal were numerically simulated using
a Box_Tree code. The code is a Barnes_Hut tree code proven effective in
modeling systems composed of large numbers of particles. Interaction between
individual particles was assumed to conform to a Yukawa potential. Particle
charge, mass, density, Debye length and output data intervals are all
adjustable parameters in the code. Employing a Fourier transform on the output
data, dispersion relations for both longitudinal and transverse wave modes were
determined. These were compared with the dispersion relations obtained from
experiment as well as a theory based on a harmonic approximation to the
potential. They were found to agree over a range of 0.9<k<5, where k is the
shielding parameter, defined by the ratio between interparticle distance a and
dust Debye length lD. This is an improvement over experimental data as current
experiments can only verify the theory up to k = 1.5.
