This routine is part of the liesymm package and is loaded via with(liesymm).

Given one or more partial differential equations this routine constructs a set of differential forms which after closure is equivalent to the original system of equations in the sense of Cartan. The forms are obtained by first reducing the system of PDEs to a system of first order equations.  It will handle systems which are quasi-linear.

If more than one equation is given, they must be given as a list or a set.

The second argument specifies the dependent and independent variables as in .  If there is more than one dependent variable, these must be specified as a set or list.  The Maple expression  can be used to construct the required list when more than one dependent variable is involved.

The third argument is used to construct the names of the extended variables. If it is the name k then the names k1, k2, k3, ..., kn are constructed as needed.  If a list of names is provided then these names are used in place of k1, k2, ..., kn.

The dependent variables (as defined by ) are processed in the order given.  For each dependent variable we introduce extended variables corresponding to the partial derivatives with respect to  t and x in that order.  These extended variables are in turn treated as dependent variables until the reduction to first order is achieved. The resulting coordinate list  is  with k1 and k2 being the partials of h with respect to t and x and w3 and w4 being the partials of u with respect to t and x.