  An algorithm to compute polynomial conserved densities of polynomial
nonlinear lattices is presented. The algorithm is implemented in Mathematica
and can be used as an automated integrability test. With the code diffdens.m,
conserved densities are obtained for several well-known lattice equations. For
systems with parameters, the code allows one to determine the conditions on
these parameters so that a sequence of conservation laws exist.
