Computation of Higher-order Symmetries for Nonlinear Evolution and   Lattice Equations
  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. 
Unal Goktas (1), Willy Hereman (1) ((1) Colorado School of Mines)