How to calculate correlation functions of Heisenberg chains

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.