The command ConservedCurrents computes conserved currents for a given differential equation system PDESYS. Given a system  (PDESYS) consisting of  equations , , involving  independent variables , and where the dependent variables are , with  denoting the set of partial derivatives of , conserved currents are expressions  each of which satisfy:  = 0, where  and  represents the total derivative with respect to . The conserved currents coincide with the traditional first integrals when there is only one independent variable, so that  is a system of ODEs.

The command ConservedCurrentTest verifies whether a given expression J is a conserved current of the system PDESYS. Note: Herein,  is assumed to mean that the left hand side is taken modulo the given PDE system; i.e.  may result in a function of  that becomes zero whenever .

The conserved currents  are computed by first constructing the PDE system  modulo PDESYS that they satisfy (see DeterminingPDE), then solving this system using pdsolve.

By default, the conserved currents are searched as functions depending on the derivatives of the unknowns of the system (they are automatically detected or can be specified as DepVars) up to the order d-1, where d is the highest order of derivatives entering PDESYS. This default can be changed by optionally passing the argument order = m, where  is a nonnegative integer.

By default, the conserved currents are searched as functions with no pre-especified form, just with the depencency explained in the previous paragraph. This default can be changed with the option typeofconservedcurrent = ... where the right-hand-side can be polynomial or functionfield, respectively indicating a conserved current of polynomial type, or of a functionfield type with the meaning explained in FunctionField.

By default, the functionality of , entering the left-hand-sides of each element in the returned lists, is displayed, the output is presented in functional notation instead of jet notation and is simplified with respect to its size. The PDE system solved to compute the  is split, when that is possible, before being tackled. All these defaults can be changed by passing the optional arguments displayfunctionality = ..., jetnotation = ..., simplifier = ..., split = false.

It is also possible to directly specify the functionality expected for the  using _J = .... See the examples for a demonstration of the use of this parameter.

To avoid having to remember the optional keywords, if you type the keyword misspelled, or just a portion of it, a matching against the correct keywords is performed, and when there is only one match, the input is automatically corrected.

[1] Olver, P.J. Applications of Lie Groups to Differential Equations. Graduate Texts in Mathematics. Springer-Verlag, 1993.