  We describe a method for calculating dynamical spin-spin correlation
functions in the finite isotropic and anisotropic antiferromagnetic Heisenberg
models. Our method is able to produce results with high accuracy over the full
parameter space.
