The main objective of this project is the advancement of methods for the exact and efficient calculation of correlation functions of integrable lattice models. In previous work we (and others) have shown that the correlation functions of the spin-1/2 Heisenberg chain and of related integrable models are characterized by a specific structure which we have called factorization: longer-range correlation functions are polynomials in one-point functions and neighbour-correlators, whose coefficients are determined by the underlying infinite-dimensional symmetry algebra. Here we plan to develop efficient computer algebra programs for the calculation of these coefficients. We want to show that the correlation functions of other integrable models factorize as well, and we want to extend the notion of factorization to the matrix elements occurring in the spectral representations of correlation functions.