The system of PDEs characterizing the symmetries of a system of PDEs will often include a number of equations which are particularly easy to integrate (for example, first order partials equated with 0). This command will automatically carry out many of these simple integrations.  If no such integrations can be found then the system is left alone.

When simplifications are completed the result is returned in a form of  eqns &where equations where the first operand is the set of remaining equations and the second operand is a set of algebraic equations indicating the relationship between the new functions appearing in the set of simplified equations and the original functions.

This routine is ordinarily loaded via with(liesymm) but can be used in the ``package style'' as liesymm[autosimp]()