  A straightforward algorithm for the symbolic computation of higher-order
symmetries of nonlinear evolution equations and lattice equations is presented.
The scaling properties of the evolution or lattice equations are used to
determine the polynomial form of the higher-order symmetries. The coefficients
of the symmetry can be found by solving a linear system. The method applies to
polynomial systems of PDEs of first-order in time and arbitrary order in one
space variable. Likewise, lattices must be of first order in time but may
involve arbitrary shifts in the discretized space variable.
  The algorithm is implemented in Mathematica and can be used to test the
integrability of both nonlinear evolution equations and semi-discrete lattice
equations. With our Integrability Package, higher-order symmetries are obtained
for several well-known systems of evolution and lattice equations. For PDEs and
lattices with parameters, the code allows one to determine the conditions on
these parameters so that a sequence of higher-order symmetries exist. The
existence of a sequence of such symmetries is a predictor for integrability.
