DeterminingPDE - Maple Help

PDEtools

 DeterminingPDE
 computes the determining PDE system satisfied by the infinitesimals of the symmetry groups, or the integrating factors or the conserved currents, admitted by a given PDE system

 Calling Sequence DeterminingPDE(PDESYS, DepVars, S, options=value)

Parameters

 PDESYS - PDE or a set or list of PDEs; it can include ODEs and non-differential equations DepVars - (optional - may be required) function or list of functions indicating the dependent variables of the problem S - (optional) list with the functional form of the infinitesimals of a symmetry generator checkconsistency = ... - (optional) can be true or false (default); to check the consistency of PDESYS before proceeding closesystem = ... - (optional) can be true or false (default); to derive and then include in PDESYS, explicitly, all its integrability conditions, necessary to compute all of PDESYS symmetries conservedcurrents = ... - (optional) indicates the computation of the determining system for the conserved currents or their algebraic form dependency = ... - (optional) name or a set of names; to indicate (extend or restrict) the dependency of the infinitesimals, integrating factors and conserved currents integrabilityconditions = ... - (optional) can be true (default) or false; to reduce, or not reduce, the determining PDE output system taking into account its integrability conditions integratingfactors = ... - (optional) indicates the computation of the determining system for the integrating factors or their algebraic form jetnotation = ... - (optional) can be jetvariables (default), jetvariableswithbrackets, jetnumbers or jetODE; to respectively return or not using the different jet notations available moduloPDESYS = ... - (optional) can be true (default) or false; to reduce or not the determining PDE output system taking into account the given PDESYS simplifier = ... - (optional) indicates the simplifier to be used instead of the default simplify/size split = ... - (optional) can be true (default) or false; to split or not the returned determining PDE system by taking the coefficients of power of independent objects not entering the dependency of the infinitesimals typeofsymmetry = ... - (optional) can be pointlike (default), contact, evolutionary or general; indicates the type of symmetry for which a determining PDE system is to be returned xi_eta = ... - (optional) list of two names to be used as roots for the ${\mathrm{\xi }}_{n}$ and ${\mathrm{\eta }}_{m}$ components of the infinitesimals

Options

 • checkconsistency = true | false

The default value is false; when true, the consistency of PDESYS is checked before proceeding. Note that one could construct the determining system for the symmetries of an inconsistent PDESYS, in which case the solution of the determining system would be meaningless.

 • closesystem = true | false

By default DeterminingPDE processes PDESYS as given. If, in PDESYS, all the integrability conditions that can be derived from its equations are explicitly present, then the resulting determining system for the symmetries is assured to have for solution all the symmetries of PDESYS. Otherwise, the returned determining system may in rare situations also (correctly) include equations that however may restrict the possible values of the symmetry solutions, and hence some of the symmetries be missed. Note also that, for a system to be closed, it is not sufficient to have it in reduced involutive or equivalent form produced by the casesplit and rifsimp commands or the DifferentialAlgebra package. To assure that all the integrability conditions are explicitly present use the option closesystem - see related example in the Examples section. See also the Library routine CloseSystem. This option is ignored when computing integrating factors or conserved currents

 • conservedcurrents

Use this option, giving just the keyword conservedcurrents, to request the determining system satisfied by the conserved currents. Each component of the conserved current will be represented by a different function depending on all the jet space variables including derivatives up to the maximum differential order found in PDESYS minus one. This option can also be used as an equation, in which case the right-hand side can be a single function, and the components of the conserved currents will be taken to be different functions all of them having the same functionality as the one indicated in the right-hand side. Alternatively you can specify the right-hand side as a list directly indicating the functionality of each of the components of the conserved currents.

 • dependency = x | {x, t, ...}

Indicates the dependency of the infinitesimals, integrating factors or conserved currents. Default value is all the jet space variables associated to point transformations. To change this default use this option with a name variable or a set or list of name variables in the right-hand side.

 • integrabilityconditions = true | false

Default value is true, and so the output system is simplified with respect to its integrability conditions before being returned. This default frequently simplifies the determining PDE system significantly, but in same cases the mathematical computations to achieve the result are lengthy. To skip that simplification and have DeterminingPDE returning faster pass this option with false in the right-hand side.

 • integratingfactors

Use this option, giving just the keyword integratingfactors, to request the determining system satisfied by the integrating factors. Each component of the conserved current will be represented by a different function depending on all the jet space variables including derivatives up to the maximum differential order found in PDESYS minus one. This option can also be used as an equation, in which case the right-hand side can be a single function, and the components of the integrating factors will be taken to be different functions all of them having the same functionality as the one indicated in the right-hand side. Alternatively you can specify the right-hand side as a list directly indicating the functionality of each of the components of the integrating factors.

 • jetnotation = jetvariables | jetvariableswithbrackets | jetODE | jetnumbers

The default value is jetvariables so that the output has all the dependent variables and its derivatives expressed in jet notation. The value jetODE can be used provided that the input is a single ODE; the output is then expressed using the jetODE notation of DEtools,Lie, the package for ODE symmetries developed before the introduction of symmetry commands in PDEtools.

 • moduloPDESYS = true | false

By default the output system is computed modulo the given PDESYS; when this option is passed with false in the right-hand side, this computational step is skipped.

 • simplifier = anything

The default value is simplify/size. Use this option to specify the simplifier of your choice. To indicate that no simplification is preferred pass the keyword none in the right-hand side.

 • split = true | false

By default the output system is split, that is, conformed by the coefficients of different powers of independent objects (e.g.: when computing point symmetries, the system's equations are the coefficients of different powers of the derivatives of the unknowns). That frequently makes the output to be an overdetermined system. When this option is passed with false in the right-hand side, the output system consists of a single partial differential equation.

 • typeofsymmetry = pointlike | contact | evolutionary | general

By default DeterminingPDE returns a system whose solutions are the point symmetries of the given PDESYS. You can also request the determining system for symmetries of evolutionary, contact or general types, where in this context general means infinitesimals depending on derivatives of each unknown of the system up to order $n-1$, where $n$ is the differential order of the unknown as found in the system, and evolutionary is like general but with the $\mathrm{\xi }$ infinitesimals (related to transformations of the independent variables) all equal to zero. By differential order of the unknown as found in the system it is meant: take the maximum differential order of an unknown within each equation of the PDE system, then take the minimum of those numbers, so that it is possible to return a determining system moduloPDESYS.

 • xi_eta = true | false

This option is used to indicate the root names to be used to represent the infinitesimals instead of the default ${\mathrm{_ξ}}_{n}$ and ${\mathrm{_η}}_{m}$.

Description

 • Given a PDE problem (PDESYS), as an equation, or a set or list of them, the DeterminingPDE command computes the determining PDE system for the symmetries, integrating factors or conserved currents of PDESYS. The output is thus a system of partial differential equations whose solution gives, respectively, the infinitesimals of the symmetry generators leaving invariant PDESYS, the integrating factors or the conserved currents of PDESYS. DeterminingPDE also works with anticommutative variables set using the Physics package.
 • By default, DeterminingPDE computes the system satisfied by the symmetry infinitesimals of PDESYS and processes PDESYS as given. To request the determining system for the integrating factors or the conserved currents pass the keywords integratingfactors or conservedcurrents. Note that in the case of symmetries, if PDESYS is not closed in the sense of including explicitly all its integrability conditions, the resulting determining system may not include for solutions all the symmetries of PDESYS. To close the system in order to compute all the symmetries use the option closesystem - see details in the Optional arguments section below - see also the related PDEtools Library routine CloseSystem.
 • Within the framework of computing the so-called group-invariant (GI) solutions of PDESYS, the first step is to compute the determining PDE system, next you attempt solving this system to obtain the infinitesimals, with them you compute an invariant transformation to reduce the number of independent variables of PDESYS, finally you attempt solving the reduced system to change variables back resulting in GI solutions for PDESYS. You can compute any of these steps directly departing from PDESYS, respectively using the commands: DeterminingPDE, Infinitesimals, InvariantTransformation and InvariantSolutions.
 • If DepVars is not given, DeterminingPDE will consider all the differentiated unknown functions in PDESYS as unknowns of the problem. Specifying DepVars however permits not only restricting the unknowns in different ways but also specifying unknowns of the problems which do not appear differentiated in PDESYS.
 • The argument S is optionally used to indicate the functional form of the symmetries - see the examples - if not given, DeterminingPDE will return a determining system whose solution are the infinitesimals of the point symmetries of PDESYS. To override that behavior when not passing S, use the option typeofsymmetry = ... where the right-hand side can be any of pointlike (default), contact, evolutionary or general (see definitions in the Options section).
 • Among the most commonly used options explained in the Options section, the PDE system returned by DeterminingPDE is automatically reduced taking into account its integrability conditions, a frequently desired but in some cases expensive mathematical computation, unless explicitly requested otherwise by passing the optional argument integrabilityconditions = false, in which case DeterminingPDE will return faster. Also, the output is by default in jetnotation = jetvariables, modulo the given PDESYS, split with respect to powers of independent objects not entering the dependency of the infinitesimals, and simplified in size. That is frequently the most compact and useful form of the determining system. You can override any of these defaults by passing any of the optional arguments integrabilityconditions = false, jetnotation = ..., moduloPDESYS = false, simplifier = none (or the one you prefer).
 • 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.

Examples

 > $\mathrm{with}\left(\mathrm{PDEtools},\mathrm{DeterminingPDE},\mathrm{declare},\mathrm{diff_table},\mathrm{casesplit},\mathrm{InfinitesimalGenerator},\mathrm{Infinitesimals},\mathrm{SymmetryTest}\right):$

Consider first a simple case; use diff_table and declare to avoid redundancies in the input (you can enter things using jetvariableswithbrackets jet notation) and in the display of the output (expressed using jetvariables jet notation)

 > $\mathrm{declare}\left(u\left(x,t\right)\right)$
 ${u}{}\left({x}{,}{t}\right){}{\mathrm{will now be displayed as}}{}{u}$ (1)
 > $U≔\mathrm{diff_table}\left(u\left(x,t\right)\right):$
 > $\mathrm{PDE1}≔{U}_{x,x}-{U}_{t}=0$
 ${\mathrm{diff}}{}\left({\mathrm{diff}}{}\left({u}{}\left({x}{,}{t}\right){,}{x}\right){,}{x}\right){-}\left({\mathrm{diff}}{}\left({u}{}\left({x}{,}{t}\right){,}{t}\right)\right){=}{0}$ (2)

The determining PDE system whose solutions are the infinitesimals

 > $\mathrm{DetSys}≔\mathrm{DeterminingPDE}\left(\mathrm{PDE1}\right)$
 $\left\{{\mathrm{diff}}{}\left({\mathrm{diff}}{}\left({\mathrm{diff}}{}\left({{\mathrm{_ξ}}}_{{t}}{}\left({x}{,}{t}{,}{u}\right){,}{t}\right){,}{t}\right){,}{t}\right){=}{0}{,}{\mathrm{diff}}{}\left({\mathrm{diff}}{}\left({{\mathrm{_η}}}_{{u}}{}\left({x}{,}{t}{,}{u}\right){,}{t}\right){,}{u}\right){=}{-}\frac{{1}}{{4}}{}{\mathrm{diff}}{}\left({\mathrm{diff}}{}\left({{\mathrm{_ξ}}}_{{t}}{}\left({x}{,}{t}{,}{u}\right){,}{t}\right){,}{t}\right){,}{\mathrm{diff}}{}\left({\mathrm{diff}}{}\left({{\mathrm{_η}}}_{{u}}{}\left({x}{,}{t}{,}{u}\right){,}{u}\right){,}{u}\right){=}{0}{,}{\mathrm{diff}}{}\left({\mathrm{diff}}{}\left({{\mathrm{_η}}}_{{u}}{}\left({x}{,}{t}{,}{u}\right){,}{x}\right){,}{x}\right){=}{\mathrm{diff}}{}\left({{\mathrm{_η}}}_{{u}}{}\left({x}{,}{t}{,}{u}\right){,}{t}\right){,}{\mathrm{diff}}{}\left({{\mathrm{_ξ}}}_{{t}}{}\left({x}{,}{t}{,}{u}\right){,}{u}\right){=}{0}{,}{\mathrm{diff}}{}\left({{\mathrm{_ξ}}}_{{t}}{}\left({x}{,}{t}{,}{u}\right){,}{x}\right){=}{0}{,}{\mathrm{diff}}{}\left({{\mathrm{_ξ}}}_{{x}}{}\left({x}{,}{t}{,}{u}\right){,}{t}\right){=}{-}{2}{}\left({\mathrm{diff}}{}\left({\mathrm{diff}}{}\left({{\mathrm{_η}}}_{{u}}{}\left({x}{,}{t}{,}{u}\right){,}{u}\right){,}{x}\right)\right){,}{\mathrm{diff}}{}\left({{\mathrm{_ξ}}}_{{x}}{}\left({x}{,}{t}{,}{u}\right){,}{u}\right){=}{0}{,}{\mathrm{diff}}{}\left({{\mathrm{_ξ}}}_{{x}}{}\left({x}{,}{t}{,}{u}\right){,}{x}\right){=}\frac{{1}}{{2}}{}{\mathrm{diff}}{}\left({{\mathrm{_ξ}}}_{{t}}{}\left({x}{,}{t}{,}{u}\right){,}{t}\right)\right\}$ (3)

The system returned has the leading derivatives isolated on the left-hand sides, and the jet notation used by default is jetvariables, so $u$ is represented by $u$, $\frac{\partial u}{\partial x}$ by ${u}_{x}$, and so on. The derivatives of the $\mathrm{_ξ}$ and $\mathrm{_η}$ infinitesimals are displayed indexed, in jetvariableswithbrackets notation, because of the use of declare.

Note the labels for the infinitesimals: they are of the form ${\mathrm{\xi }}_{{x}_{n}}$ and ${\mathrm{\eta }}_{{u}_{m}}$, where ${x}_{n}$ and ${u}_{m}$ are respectively the names of independent and dependent variables; these are the infinitesimals associated to symmetry transformations of the corresponding (independent or dependent) variable.

You can also specify that you want a different jet notation, for instance jetnumbers (a less mathematical notation, sometimes preferred for programming purposes), in which case the labels identifying each infinitesimal will also appear in jetnumbers notation, as in ${\mathrm{\xi }}_{n}$ and ${\mathrm{\eta }}_{m}$

 > $\mathrm{DeterminingPDE}\left(\mathrm{PDE1},\mathrm{jetnotation}=\mathrm{jetnumbers}\right)$
 $\left\{{\mathrm{diff}}{}\left({\mathrm{diff}}{}\left({\mathrm{diff}}{}\left({{\mathrm{_ξ}}}_{{2}}{}\left({x}{,}{t}{,}{{u}}_{{[}{]}}\right){,}{t}\right){,}{t}\right){,}{t}\right){=}{0}{,}{\mathrm{diff}}{}\left({\mathrm{diff}}{}\left({{\mathrm{_η}}}_{{1}}{}\left({x}{,}{t}{,}{{u}}_{{[}{]}}\right){,}{t}\right){,}{{u}}_{{[}{]}}\right){=}{-}\frac{{1}}{{4}}{}{\mathrm{diff}}{}\left({\mathrm{diff}}{}\left({{\mathrm{_ξ}}}_{{2}}{}\left({x}{,}{t}{,}{{u}}_{{[}{]}}\right){,}{t}\right){,}{t}\right){,}{\mathrm{diff}}{}\left({\mathrm{diff}}{}\left({{\mathrm{_η}}}_{{1}}{}\left({x}{,}{t}{,}{{u}}_{{[}{]}}\right){,}{x}\right){,}{x}\right){=}{\mathrm{diff}}{}\left({{\mathrm{_η}}}_{{1}}{}\left({x}{,}{t}{,}{{u}}_{{[}{]}}\right){,}{t}\right){,}{\mathrm{diff}}{}\left({\mathrm{diff}}{}\left({{\mathrm{_η}}}_{{1}}{}\left({x}{,}{t}{,}{{u}}_{{[}{]}}\right){,}{{u}}_{{[}{]}}\right){,}{{u}}_{{[}{]}}\right){=}{0}{,}{\mathrm{diff}}{}\left({{\mathrm{_ξ}}}_{{1}}{}\left({x}{,}{t}{,}{{u}}_{{[}{]}}\right){,}{t}\right){=}{-}{2}{}\left({\mathrm{diff}}{}\left({\mathrm{diff}}{}\left({{\mathrm{_η}}}_{{1}}{}\left({x}{,}{t}{,}{{u}}_{{[}{]}}\right){,}{x}\right){,}{{u}}_{{[}{]}}\right)\right){,}{\mathrm{diff}}{}\left({{\mathrm{_ξ}}}_{{1}}{}\left({x}{,}{t}{,}{{u}}_{{[}{]}}\right){,}{x}\right){=}\frac{{1}}{{2}}{}{\mathrm{diff}}{}\left({{\mathrm{_ξ}}}_{{2}}{}\left({x}{,}{t}{,}{{u}}_{{[}{]}}\right){,}{t}\right){,}{\mathrm{diff}}{}\left({{\mathrm{_ξ}}}_{{1}}{}\left({x}{,}{t}{,}{{u}}_{{[}{]}}\right){,}{{u}}_{{[}{]}}\right){=}{0}{,}{\mathrm{diff}}{}\left({{\mathrm{_ξ}}}_{{2}}{}\left({x}{,}{t}{,}{{u}}_{{[}{]}}\right){,}{x}\right){=}{0}{,}{\mathrm{diff}}{}\left({{\mathrm{_ξ}}}_{{2}}{}\left({x}{,}{t}{,}{{u}}_{{[}{]}}\right){,}{{u}}_{{[}{]}}\right){=}{0}\right\}$ (4)

The system returned by DeterminingPDE is automatically simplified with respect to its integrability conditions. This simplification is normally desirable, but in some cases it can be expensive in time or depending on the problem even a bottleneck in the computation of symmetries. To skip this simplification use the option integrabilityconditions = false, in which case the leading derivatives are not isolated; compare for instance $\mathrm{DetSys}$ with this other equivalent determining system

 > $\mathrm{DeterminingPDE}\left(\mathrm{PDE1},\mathrm{integrabilityconditions}=\mathrm{false}\right)$
 $\left\{{\mathrm{diff}}{}\left({\mathrm{diff}}{}\left({{\mathrm{_η}}}_{{u}}{}\left({x}{,}{t}{,}{u}\right){,}{x}\right){,}{x}\right){-}\left({\mathrm{diff}}{}\left({{\mathrm{_η}}}_{{u}}{}\left({x}{,}{t}{,}{u}\right){,}{t}\right)\right){,}{-}{2}{}\left({\mathrm{diff}}{}\left({\mathrm{diff}}{}\left({{\mathrm{_ξ}}}_{{t}}{}\left({x}{,}{t}{,}{u}\right){,}{u}\right){,}{x}\right)\right){-}{2}{}\left({\mathrm{diff}}{}\left({{\mathrm{_ξ}}}_{{x}}{}\left({x}{,}{t}{,}{u}\right){,}{u}\right)\right){,}{-}{2}{}\left({\mathrm{diff}}{}\left({\mathrm{diff}}{}\left({{\mathrm{_ξ}}}_{{x}}{}\left({x}{,}{t}{,}{u}\right){,}{u}\right){,}{x}\right)\right){+}{\mathrm{diff}}{}\left({\mathrm{diff}}{}\left({{\mathrm{_η}}}_{{u}}{}\left({x}{,}{t}{,}{u}\right){,}{u}\right){,}{u}\right){,}{2}{}\left({\mathrm{diff}}{}\left({\mathrm{diff}}{}\left({{\mathrm{_η}}}_{{u}}{}\left({x}{,}{t}{,}{u}\right){,}{u}\right){,}{x}\right)\right){-}\left({\mathrm{diff}}{}\left({\mathrm{diff}}{}\left({{\mathrm{_ξ}}}_{{x}}{}\left({x}{,}{t}{,}{u}\right){,}{x}\right){,}{x}\right)\right){+}{\mathrm{diff}}{}\left({{\mathrm{_ξ}}}_{{x}}{}\left({x}{,}{t}{,}{u}\right){,}{t}\right){,}{-}\left({\mathrm{diff}}{}\left({\mathrm{diff}}{}\left({{\mathrm{_ξ}}}_{{t}}{}\left({x}{,}{t}{,}{u}\right){,}{x}\right){,}{x}\right)\right){-}{2}{}\left({\mathrm{diff}}{}\left({{\mathrm{_ξ}}}_{{x}}{}\left({x}{,}{t}{,}{u}\right){,}{x}\right)\right){+}{\mathrm{diff}}{}\left({{\mathrm{_ξ}}}_{{t}}{}\left({x}{,}{t}{,}{u}\right){,}{t}\right){,}{\mathrm{diff}}{}\left({\mathrm{diff}}{}\left({{\mathrm{_ξ}}}_{{t}}{}\left({x}{,}{t}{,}{u}\right){,}{u}\right){,}{u}\right){,}{\mathrm{diff}}{}\left({\mathrm{diff}}{}\left({{\mathrm{_ξ}}}_{{x}}{}\left({x}{,}{t}{,}{u}\right){,}{u}\right){,}{u}\right){,}{\mathrm{diff}}{}\left({{\mathrm{_ξ}}}_{{t}}{}\left({x}{,}{t}{,}{u}\right){,}{u}\right){,}{\mathrm{diff}}{}\left({{\mathrm{_ξ}}}_{{t}}{}\left({x}{,}{t}{,}{u}\right){,}{x}\right)\right\}$ (5)

The fastest (and less simplified) way to compute the determining PDE for the symmetries includes skipping the automatic splitting into coefficients of independent objects, in which case the output consists of a set with a single determining PDE, not a system

 > $\mathrm{DeterminingPDE}\left(\mathrm{PDE1},\mathrm{integrabilityconditions}=\mathrm{false},\mathrm{split}=\mathrm{false}\right)$
 $\left\{{-}\left({\mathrm{diff}}{}\left({\mathrm{diff}}{}\left({{\mathrm{_ξ}}}_{{x}}{}\left({x}{,}{t}{,}{u}\right){,}{u}\right){,}{u}\right)\right){}{{u}}_{{x}}^{{3}}{-}\left({\mathrm{diff}}{}\left({\mathrm{diff}}{}\left({{\mathrm{_ξ}}}_{{t}}{}\left({x}{,}{t}{,}{u}\right){,}{u}\right){,}{u}\right)\right){}{{u}}_{{x}}^{{2}}{}{{u}}_{{t}}{-}{2}{}\left({\mathrm{diff}}{}\left({\mathrm{diff}}{}\left({{\mathrm{_ξ}}}_{{x}}{}\left({x}{,}{t}{,}{u}\right){,}{u}\right){,}{x}\right)\right){}{{u}}_{{x}}^{{2}}{-}{2}{}\left({\mathrm{diff}}{}\left({\mathrm{diff}}{}\left({{\mathrm{_ξ}}}_{{t}}{}\left({x}{,}{t}{,}{u}\right){,}{u}\right){,}{x}\right)\right){}{{u}}_{{x}}{}{{u}}_{{t}}{+}\left({\mathrm{diff}}{}\left({\mathrm{diff}}{}\left({{\mathrm{_η}}}_{{u}}{}\left({x}{,}{t}{,}{u}\right){,}{u}\right){,}{u}\right)\right){}{{u}}_{{x}}^{{2}}{-}{2}{}\left({\mathrm{diff}}{}\left({{\mathrm{_ξ}}}_{{x}}{}\left({x}{,}{t}{,}{u}\right){,}{u}\right)\right){}{{u}}_{{x}}{}{{u}}_{{t}}{-}{2}{}\left({\mathrm{diff}}{}\left({{\mathrm{_ξ}}}_{{t}}{}\left({x}{,}{t}{,}{u}\right){,}{u}\right)\right){}{{u}}_{{x}}{}{{u}}_{{x}{,}{t}}{+}{2}{}{{u}}_{{x}}{}\left({\mathrm{diff}}{}\left({\mathrm{diff}}{}\left({{\mathrm{_η}}}_{{u}}{}\left({x}{,}{t}{,}{u}\right){,}{u}\right){,}{x}\right)\right){-}\left({\mathrm{diff}}{}\left({\mathrm{diff}}{}\left({{\mathrm{_ξ}}}_{{x}}{}\left({x}{,}{t}{,}{u}\right){,}{x}\right){,}{x}\right)\right){}{{u}}_{{x}}{-}\left({\mathrm{diff}}{}\left({\mathrm{diff}}{}\left({{\mathrm{_ξ}}}_{{t}}{}\left({x}{,}{t}{,}{u}\right){,}{x}\right){,}{x}\right)\right){}{{u}}_{{t}}{-}{2}{}\left({\mathrm{diff}}{}\left({{\mathrm{_ξ}}}_{{x}}{}\left({x}{,}{t}{,}{u}\right){,}{x}\right)\right){}{{u}}_{{t}}{-}{2}{}\left({\mathrm{diff}}{}\left({{\mathrm{_ξ}}}_{{t}}{}\left({x}{,}{t}{,}{u}\right){,}{x}\right)\right){}{{u}}_{{x}{,}{t}}{+}{\mathrm{diff}}{}\left({\mathrm{diff}}{}\left({{\mathrm{_η}}}_{{u}}{}\left({x}{,}{t}{,}{u}\right){,}{x}\right){,}{x}\right){+}\left({\mathrm{diff}}{}\left({{\mathrm{_ξ}}}_{{x}}{}\left({x}{,}{t}{,}{u}\right){,}{t}\right)\right){}{{u}}_{{x}}{+}\left({\mathrm{diff}}{}\left({{\mathrm{_ξ}}}_{{t}}{}\left({x}{,}{t}{,}{u}\right){,}{t}\right)\right){}{{u}}_{{t}}{-}\left({\mathrm{diff}}{}\left({{\mathrm{_η}}}_{{u}}{}\left({x}{,}{t}{,}{u}\right){,}{t}\right)\right)\right\}$ (6)

Note in the output above the presence of ${u}_{x}$ and ${u}_{t}$ representing derivatives of $u$ written in jetnotation. To see the last computed output with all its functionality and structure explicit, not using the enhanced compact mathematical display due to declare, use show

 > $\mathrm{show}$
 $\left\{{-}\left(\frac{{{\partial }}^{{2}}}{{\partial }{{u}}^{{2}}}\phantom{\rule[-0.0ex]{0.4em}{0.0ex}}{{\mathrm{_ξ}}}_{{x}}{}\left({x}{,}{t}{,}{u}\right)\right){}{{u}}_{{x}}^{{3}}{-}\left(\frac{{{\partial }}^{{2}}}{{\partial }{{u}}^{{2}}}\phantom{\rule[-0.0ex]{0.4em}{0.0ex}}{{\mathrm{_ξ}}}_{{t}}{}\left({x}{,}{t}{,}{u}\right)\right){}{{u}}_{{x}}^{{2}}{}{{u}}_{{t}}{-}{2}{}\left(\frac{{{\partial }}^{{2}}}{{\partial }{u}{\partial }{x}}\phantom{\rule[-0.0ex]{0.4em}{0.0ex}}{{\mathrm{_ξ}}}_{{x}}{}\left({x}{,}{t}{,}{u}\right)\right){}{{u}}_{{x}}^{{2}}{-}{2}{}\left(\frac{{{\partial }}^{{2}}}{{\partial }{u}{\partial }{x}}\phantom{\rule[-0.0ex]{0.4em}{0.0ex}}{{\mathrm{_ξ}}}_{{t}}{}\left({x}{,}{t}{,}{u}\right)\right){}{{u}}_{{x}}{}{{u}}_{{t}}{+}\left(\frac{{{\partial }}^{{2}}}{{\partial }{{u}}^{{2}}}\phantom{\rule[-0.0ex]{0.4em}{0.0ex}}{{\mathrm{_η}}}_{{u}}{}\left({x}{,}{t}{,}{u}\right)\right){}{{u}}_{{x}}^{{2}}{-}{2}{}\left(\frac{{\partial }}{{\partial }{u}}\phantom{\rule[-0.0ex]{0.4em}{0.0ex}}{{\mathrm{_ξ}}}_{{x}}{}\left({x}{,}{t}{,}{u}\right)\right){}{{u}}_{{x}}{}{{u}}_{{t}}{-}{2}{}\left(\frac{{\partial }}{{\partial }{u}}\phantom{\rule[-0.0ex]{0.4em}{0.0ex}}{{\mathrm{_ξ}}}_{{t}}{}\left({x}{,}{t}{,}{u}\right)\right){}{{u}}_{{x}}{}{{u}}_{{x}{,}{t}}{+}{2}{}{{u}}_{{x}}{}\left(\frac{{{\partial }}^{{2}}}{{\partial }{u}{\partial }{x}}\phantom{\rule[-0.0ex]{0.4em}{0.0ex}}{{\mathrm{_η}}}_{{u}}{}\left({x}{,}{t}{,}{u}\right)\right){-}\left(\frac{{{\partial }}^{{2}}}{{\partial }{{x}}^{{2}}}\phantom{\rule[-0.0ex]{0.4em}{0.0ex}}{{\mathrm{_ξ}}}_{{x}}{}\left({x}{,}{t}{,}{u}\right)\right){}{{u}}_{{x}}{-}\left(\frac{{{\partial }}^{{2}}}{{\partial }{{x}}^{{2}}}\phantom{\rule[-0.0ex]{0.4em}{0.0ex}}{{\mathrm{_ξ}}}_{{t}}{}\left({x}{,}{t}{,}{u}\right)\right){}{{u}}_{{t}}{-}{2}{}\left(\frac{{\partial }}{{\partial }{x}}\phantom{\rule[-0.0ex]{0.4em}{0.0ex}}{{\mathrm{_ξ}}}_{{x}}{}\left({x}{,}{t}{,}{u}\right)\right){}{{u}}_{{t}}{-}{2}{}\left(\frac{{\partial }}{{\partial }{x}}\phantom{\rule[-0.0ex]{0.4em}{0.0ex}}{{\mathrm{_ξ}}}_{{t}}{}\left({x}{,}{t}{,}{u}\right)\right){}{{u}}_{{x}{,}{t}}{+}\frac{{{\partial }}^{{2}}}{{\partial }{{x}}^{{2}}}\phantom{\rule[-0.0ex]{0.4em}{0.0ex}}{{\mathrm{_η}}}_{{u}}{}\left({x}{,}{t}{,}{u}\right){+}\left(\frac{{\partial }}{{\partial }{t}}\phantom{\rule[-0.0ex]{0.4em}{0.0ex}}{{\mathrm{_ξ}}}_{{x}}{}\left({x}{,}{t}{,}{u}\right)\right){}{{u}}_{{x}}{+}\left(\frac{{\partial }}{{\partial }{t}}\phantom{\rule[-0.0ex]{0.4em}{0.0ex}}{{\mathrm{_ξ}}}_{{t}}{}\left({x}{,}{t}{,}{u}\right)\right){}{{u}}_{{t}}{-}\frac{{\partial }}{{\partial }{t}}\phantom{\rule[-0.0ex]{0.4em}{0.0ex}}{{\mathrm{_η}}}_{{u}}{}\left({x}{,}{t}{,}{u}\right)\right\}$ (7)

Note in the above that although the infinitesimals and its derivatives are not displayed in jet notation anymore, the unknown $u$ and its derivatives ${u}_{x}$, ${u}_{t}$ are still in jet notation - they enter the expression actually as names, not functions, but representing the functions. As usual, to convert into function notation use FromJet.

The infinitesimals are obtained solving this determining system - we take the most simplified version of it, $\mathrm{DetSys}$, computed lines above

 > $\mathrm{pdsolve}\left(\mathrm{DetSys}\right)$
 $\left\{{{\mathrm{_η}}}_{{u}}{}\left({x}{,}{t}{,}{u}\right){=}\left(\frac{\left({-}{{x}}^{{2}}{-}{2}{}{t}\right){}\mathrm{c__1}}{{8}}{+}\mathrm{c__4}{}{x}{+}\mathrm{c__5}\right){}{u}{+}\mathrm{c__8}{}{{ⅇ}}^{{{\mathrm{_c}}}_{{1}}{}{t}}{}\mathrm{c__6}{}{{ⅇ}}^{\sqrt{{{\mathrm{_c}}}_{{1}}}{}{x}}{+}\frac{\mathrm{c__8}{}{{ⅇ}}^{{{\mathrm{_c}}}_{{1}}{}{t}}{}\mathrm{c__7}}{{{ⅇ}}^{\sqrt{{{\mathrm{_c}}}_{{1}}}{}{x}}}{,}{{\mathrm{_ξ}}}_{{t}}{}\left({x}{,}{t}{,}{u}\right){=}\frac{{1}}{{2}}{}\mathrm{c__1}{}{{t}}^{{2}}{+}\mathrm{c__2}{}{t}{+}\mathrm{c__3}{,}{{\mathrm{_ξ}}}_{{x}}{}\left({x}{,}{t}{,}{u}\right){=}\frac{\left({x}{}\mathrm{c__1}{-}{4}{}\mathrm{c__4}\right){}{t}}{{2}}{+}\frac{\mathrm{c__2}{}{x}}{{2}}{+}\mathrm{c__9}\right\}$ (8)

You can compute these infinitesimals directly from PDE1 as well. By default, whenever, it is possible Infinitesimals will present the symmetry groups with all its independent members explicitly, that is, with the integration constants specialized, as in

 > $\mathrm{Infinitesimals}\left(\mathrm{PDE1}\right)$
 $\left[{{\mathrm{_ξ}}}_{{x}}{}\left({x}{,}{t}{,}{u}\right){=}{0}{,}{{\mathrm{_ξ}}}_{{t}}{}\left({x}{,}{t}{,}{u}\right){=}{1}{,}{{\mathrm{_η}}}_{{u}}{}\left({x}{,}{t}{,}{u}\right){=}{0}\right]{,}\left[{{\mathrm{_ξ}}}_{{x}}{}\left({x}{,}{t}{,}{u}\right){=}{1}{,}{{\mathrm{_ξ}}}_{{t}}{}\left({x}{,}{t}{,}{u}\right){=}{0}{,}{{\mathrm{_η}}}_{{u}}{}\left({x}{,}{t}{,}{u}\right){=}{0}\right]{,}\left[{{\mathrm{_ξ}}}_{{x}}{}\left({x}{,}{t}{,}{u}\right){=}{0}{,}{{\mathrm{_ξ}}}_{{t}}{}\left({x}{,}{t}{,}{u}\right){=}{0}{,}{{\mathrm{_η}}}_{{u}}{}\left({x}{,}{t}{,}{u}\right){=}{u}\right]{,}\left[{{\mathrm{_ξ}}}_{{x}}{}\left({x}{,}{t}{,}{u}\right){=}\frac{{x}}{{2}}{,}{{\mathrm{_ξ}}}_{{t}}{}\left({x}{,}{t}{,}{u}\right){=}{t}{,}{{\mathrm{_η}}}_{{u}}{}\left({x}{,}{t}{,}{u}\right){=}{0}\right]{,}\left[{{\mathrm{_ξ}}}_{{x}}{}\left({x}{,}{t}{,}{u}\right){=}{-}{2}{}{t}{,}{{\mathrm{_ξ}}}_{{t}}{}\left({x}{,}{t}{,}{u}\right){=}{0}{,}{{\mathrm{_η}}}_{{u}}{}\left({x}{,}{t}{,}{u}\right){=}{x}{}{u}\right]{,}\left[{{\mathrm{_ξ}}}_{{x}}{}\left({x}{,}{t}{,}{u}\right){=}{0}{,}{{\mathrm{_ξ}}}_{{t}}{}\left({x}{,}{t}{,}{u}\right){=}{0}{,}{{\mathrm{_η}}}_{{u}}{}\left({x}{,}{t}{,}{u}\right){=}{{ⅇ}}^{{{\mathrm{_c}}}_{{1}}{}{t}}{}{{ⅇ}}^{\sqrt{{{\mathrm{_c}}}_{{1}}}{}{x}}\right]{,}\left[{{\mathrm{_ξ}}}_{{x}}{}\left({x}{,}{t}{,}{u}\right){=}{0}{,}{{\mathrm{_ξ}}}_{{t}}{}\left({x}{,}{t}{,}{u}\right){=}{0}{,}{{\mathrm{_η}}}_{{u}}{}\left({x}{,}{t}{,}{u}\right){=}\frac{{{ⅇ}}^{{{\mathrm{_c}}}_{{1}}{}{t}}}{{{ⅇ}}^{\sqrt{{{\mathrm{_c}}}_{{1}}}{}{x}}}\right]{,}\left[{{\mathrm{_ξ}}}_{{x}}{}\left({x}{,}{t}{,}{u}\right){=}\frac{{x}{}{t}}{{2}}{,}{{\mathrm{_ξ}}}_{{t}}{}\left({x}{,}{t}{,}{u}\right){=}\frac{{{t}}^{{2}}}{{2}}{,}{{\mathrm{_η}}}_{{u}}{}\left({x}{,}{t}{,}{u}\right){=}{-}\frac{\left({{x}}^{{2}}{+}{2}{}{t}\right){}{u}}{{8}}\right]$ (9)

The determining PDE system for the integrating factors of PDE1

 > $\mathrm{DeterminingPDE}\left(\mathrm{PDE1},\mathrm{integratingfactors}\right)$
 $\left\{{\mathrm{diff}}{}\left({\mathrm{diff}}{}\left({{\mathrm{_μ}}}_{{1}}{}\left({x}{,}{t}{,}{u}{,}{{u}}_{{x}}{,}{{u}}_{{t}}\right){,}{x}\right){,}{x}\right){=}{-}\left({\mathrm{diff}}{}\left({{\mathrm{_μ}}}_{{1}}{}\left({x}{,}{t}{,}{u}{,}{{u}}_{{x}}{,}{{u}}_{{t}}\right){,}{t}\right)\right){,}{\mathrm{diff}}{}\left({{\mathrm{_μ}}}_{{1}}{}\left({x}{,}{t}{,}{u}{,}{{u}}_{{x}}{,}{{u}}_{{t}}\right){,}{u}\right){=}{0}{,}{\mathrm{diff}}{}\left({{\mathrm{_μ}}}_{{1}}{}\left({x}{,}{t}{,}{u}{,}{{u}}_{{x}}{,}{{u}}_{{t}}\right){,}{{u}}_{{t}}\right){=}{0}{,}{\mathrm{diff}}{}\left({{\mathrm{_μ}}}_{{1}}{}\left({x}{,}{t}{,}{u}{,}{{u}}_{{x}}{,}{{u}}_{{t}}\right){,}{{u}}_{{x}}\right){=}{0}\right\}$ (10)

The determining PDE system for the conserved currents of PDE1 that depend only on $x$ and $t$

 > $\mathrm{DeterminingPDE}\left(\mathrm{PDE1},\mathrm{conservedcurrents},\mathrm{dependency}=\left[x,t\right]\right)$
 $\left\{{\mathrm{diff}}{}\left({{\mathrm{_J}}}_{{x}}{}\left({x}{,}{t}\right){,}{x}\right){+}{\mathrm{diff}}{}\left({{\mathrm{_J}}}_{{t}}{}\left({x}{,}{t}\right){,}{t}\right)\right\}$ (11)

To compute the actual integrating factors or conserved currents - so not just their determining equations- see IntegratingFactors and ConservedCurrents.

The Laplace equation  in spherical coordinates with a source term ${u}_{t,t}$ symmetries

 > $\mathrm{declare}\left(u\left(r,\mathrm{φ},\mathrm{θ},t\right),\left(\mathrm{_ξ},\mathrm{_η}\right)\left(r,\mathrm{φ},\mathrm{θ},t,u\right)\right)$
 ${u}{}\left({r}{,}{\mathrm{\phi }}{,}{\mathrm{\theta }}{,}{t}\right){}{\mathrm{will now be displayed as}}{}{u}$
 ${\mathrm{_ξ}}{}\left({r}{,}{\mathrm{\phi }}{,}{\mathrm{\theta }}{,}{t}{,}{u}\right){}{\mathrm{will now be displayed as}}{}{\mathrm{_ξ}}$
 ${\mathrm{_η}}{}\left({r}{,}{\mathrm{\phi }}{,}{\mathrm{\theta }}{,}{t}{,}{u}\right){}{\mathrm{will now be displayed as}}{}{\mathrm{_η}}$ (12)
 > $U≔\mathrm{diff_table}\left(u\left(r,\mathrm{φ},\mathrm{θ},t\right)\right):$
 > $\mathrm{PDE2}≔\frac{1\left(2r{U}_{r}+{r}^{2}{U}_{r,r}\right)}{{r}^{2}}+\frac{1\left(\mathrm{cos}\left(\mathrm{θ}\right){U}_{\mathrm{θ}}+\mathrm{sin}\left(\mathrm{θ}\right){U}_{\mathrm{θ},\mathrm{θ}}\right)}{{r}^{2}\mathrm{sin}\left(\mathrm{θ}\right)}+\frac{1{U}_{\mathrm{φ},\mathrm{φ}}}{{r}^{2}{\mathrm{sin}\left(\mathrm{θ}\right)}^{2}}={U}_{t,t}$
 $\frac{{{r}}^{{2}}{}\left({\mathrm{diff}}{}\left({\mathrm{diff}}{}\left({u}{}\left({r}{,}{\mathrm{φ}}{,}{\mathrm{θ}}{,}{t}\right){,}{r}\right){,}{r}\right)\right){+}{2}{}{r}{}\left({\mathrm{diff}}{}\left({u}{}\left({r}{,}{\mathrm{φ}}{,}{\mathrm{θ}}{,}{t}\right){,}{r}\right)\right)}{{{r}}^{{2}}}{+}\frac{{\mathrm{cos}}{}\left({\mathrm{θ}}\right){}\left({\mathrm{diff}}{}\left({u}{}\left({r}{,}{\mathrm{φ}}{,}{\mathrm{θ}}{,}{t}\right){,}{\mathrm{θ}}\right)\right){+}{\mathrm{sin}}{}\left({\mathrm{θ}}\right){}\left({\mathrm{diff}}{}\left({\mathrm{diff}}{}\left({u}{}\left({r}{,}{\mathrm{φ}}{,}{\mathrm{θ}}{,}{t}\right){,}{\mathrm{θ}}\right){,}{\mathrm{θ}}\right)\right)}{{{r}}^{{2}}{}{\mathrm{sin}}{}\left({\mathrm{θ}}\right)}{+}\frac{{\mathrm{diff}}{}\left({\mathrm{diff}}{}\left({u}{}\left({r}{,}{\mathrm{φ}}{,}{\mathrm{θ}}{,}{t}\right){,}{\mathrm{φ}}\right){,}{\mathrm{φ}}\right)}{{{r}}^{{2}}{}{{\mathrm{sin}}{}\left({\mathrm{θ}}\right)}^{{2}}}{=}{\mathrm{diff}}{}\left({\mathrm{diff}}{}\left({u}{}\left({r}{,}{\mathrm{φ}}{,}{\mathrm{θ}}{,}{t}\right){,}{t}\right){,}{t}\right)$ (13)

The infinitesimals are computed either by Infinitesimals by calling DeterminingPDE to obtain the determining system then solving it calling pdsolve, finally splitting the result according to the integration constants found in pdsolve's solution

 > $\mathrm{Infinitesimals}\left(\mathrm{PDE2}\right)$
 $\left[{{\mathrm{_ξ}}}_{{r}}{}\left({r}{,}{\mathrm{φ}}{,}{\mathrm{θ}}{,}{t}{,}{u}\right){=}{0}{,}{{\mathrm{_ξ}}}_{{\mathrm{φ}}}{}\left({r}{,}{\mathrm{φ}}{,}{\mathrm{θ}}{,}{t}{,}{u}\right){=}{0}{,}{{\mathrm{_ξ}}}_{{\mathrm{θ}}}{}\left({r}{,}{\mathrm{φ}}{,}{\mathrm{θ}}{,}{t}{,}{u}\right){=}{0}{,}{{\mathrm{_ξ}}}_{{t}}{}\left({r}{,}{\mathrm{φ}}{,}{\mathrm{θ}}{,}{t}{,}{u}\right){=}{1}{,}{{\mathrm{_η}}}_{{u}}{}\left({r}{,}{\mathrm{φ}}{,}{\mathrm{θ}}{,}{t}{,}{u}\right){=}{0}\right]{,}\left[{{\mathrm{_ξ}}}_{{r}}{}\left({r}{,}{\mathrm{φ}}{,}{\mathrm{θ}}{,}{t}{,}{u}\right){=}{0}{,}{{\mathrm{_ξ}}}_{{\mathrm{φ}}}{}\left({r}{,}{\mathrm{φ}}{,}{\mathrm{θ}}{,}{t}{,}{u}\right){=}{1}{,}{{\mathrm{_ξ}}}_{{\mathrm{θ}}}{}\left({r}{,}{\mathrm{φ}}{,}{\mathrm{θ}}{,}{t}{,}{u}\right){=}{0}{,}{{\mathrm{_ξ}}}_{{t}}{}\left({r}{,}{\mathrm{φ}}{,}{\mathrm{θ}}{,}{t}{,}{u}\right){=}{0}{,}{{\mathrm{_η}}}_{{u}}{}\left({r}{,}{\mathrm{φ}}{,}{\mathrm{θ}}{,}{t}{,}{u}\right){=}{0}\right]{,}\left[{{\mathrm{_ξ}}}_{{r}}{}\left({r}{,}{\mathrm{φ}}{,}{\mathrm{θ}}{,}{t}{,}{u}\right){=}{r}{,}{{\mathrm{_ξ}}}_{{\mathrm{φ}}}{}\left({r}{,}{\mathrm{φ}}{,}{\mathrm{θ}}{,}{t}{,}{u}\right){=}{0}{,}{{\mathrm{_ξ}}}_{{\mathrm{θ}}}{}\left({r}{,}{\mathrm{φ}}{,}{\mathrm{θ}}{,}{t}{,}{u}\right){=}{0}{,}{{\mathrm{_ξ}}}_{{t}}{}\left({r}{,}{\mathrm{φ}}{,}{\mathrm{θ}}{,}{t}{,}{u}\right){=}{t}{,}{{\mathrm{_η}}}_{{u}}{}\left({r}{,}{\mathrm{φ}}{,}{\mathrm{θ}}{,}{t}{,}{u}\right){=}{0}\right]{,}\left[{{\mathrm{_ξ}}}_{{r}}{}\left({r}{,}{\mathrm{φ}}{,}{\mathrm{θ}}{,}{t}{,}{u}\right){=}{0}{,}{{\mathrm{_ξ}}}_{{\mathrm{φ}}}{}\left({r}{,}{\mathrm{φ}}{,}{\mathrm{θ}}{,}{t}{,}{u}\right){=}{0}{,}{{\mathrm{_ξ}}}_{{\mathrm{θ}}}{}\left({r}{,}{\mathrm{φ}}{,}{\mathrm{θ}}{,}{t}{,}{u}\right){=}{0}{,}{{\mathrm{_ξ}}}_{{t}}{}\left({r}{,}{\mathrm{φ}}{,}{\mathrm{θ}}{,}{t}{,}{u}\right){=}{0}{,}{{\mathrm{_η}}}_{{u}}{}\left({r}{,}{\mathrm{φ}}{,}{\mathrm{θ}}{,}{t}{,}{u}\right){=}{u}\right]{,}\left[{{\mathrm{_ξ}}}_{{r}}{}\left({r}{,}{\mathrm{φ}}{,}{\mathrm{θ}}{,}{t}{,}{u}\right){=}{-}{\mathrm{cos}}{}\left({\mathrm{\theta }}\right){,}{{\mathrm{_ξ}}}_{{\mathrm{φ}}}{}\left({r}{,}{\mathrm{φ}}{,}{\mathrm{θ}}{,}{t}{,}{u}\right){=}{0}{,}{{\mathrm{_ξ}}}_{{\mathrm{θ}}}{}\left({r}{,}{\mathrm{φ}}{,}{\mathrm{θ}}{,}{t}{,}{u}\right){=}\frac{{\mathrm{sin}}{}\left({\mathrm{\theta }}\right)}{{r}}{,}{{\mathrm{_ξ}}}_{{t}}{}\left({r}{,}{\mathrm{φ}}{,}{\mathrm{θ}}{,}{t}{,}{u}\right){=}{0}{,}{{\mathrm{_η}}}_{{u}}{}\left({r}{,}{\mathrm{φ}}{,}{\mathrm{θ}}{,}{t}{,}{u}\right){=}{0}\right]{,}\left[{{\mathrm{_ξ}}}_{{r}}{}\left({r}{,}{\mathrm{φ}}{,}{\mathrm{θ}}{,}{t}{,}{u}\right){=}{r}{}{t}{,}{{\mathrm{_ξ}}}_{{\mathrm{φ}}}{}\left({r}{,}{\mathrm{φ}}{,}{\mathrm{θ}}{,}{t}{,}{u}\right){=}{0}{,}{{\mathrm{_ξ}}}_{{\mathrm{θ}}}{}\left({r}{,}{\mathrm{φ}}{,}{\mathrm{θ}}{,}{t}{,}{u}\right){=}{0}{,}{{\mathrm{_ξ}}}_{{t}}{}\left({r}{,}{\mathrm{φ}}{,}{\mathrm{θ}}{,}{t}{,}{u}\right){=}\frac{{{r}}^{{2}}}{{2}}{+}\frac{{{t}}^{{2}}}{{2}}{,}{{\mathrm{_η}}}_{{u}}{}\left({r}{,}{\mathrm{φ}}{,}{\mathrm{θ}}{,}{t}{,}{u}\right){=}{-}{t}{}{u}\right]{,}\left[{{\mathrm{_ξ}}}_{{r}}{}\left({r}{,}{\mathrm{φ}}{,}{\mathrm{θ}}{,}{t}{,}{u}\right){=}{0}{,}{{\mathrm{_ξ}}}_{{\mathrm{φ}}}{}\left({r}{,}{\mathrm{φ}}{,}{\mathrm{θ}}{,}{t}{,}{u}\right){=}{\mathrm{cos}}{}\left({\mathrm{\phi }}\right){}{\mathrm{cot}}{}\left({\mathrm{\theta }}\right){,}{{\mathrm{_ξ}}}_{{\mathrm{θ}}}{}\left({r}{,}{\mathrm{φ}}{,}{\mathrm{θ}}{,}{t}{,}{u}\right){=}{\mathrm{sin}}{}\left({\mathrm{\phi }}\right){,}{{\mathrm{_ξ}}}_{{t}}{}\left({r}{,}{\mathrm{φ}}{,}{\mathrm{θ}}{,}{t}{,}{u}\right){=}{0}{,}{{\mathrm{_η}}}_{{u}}{}\left({r}{,}{\mathrm{φ}}{,}{\mathrm{θ}}{,}{t}{,}{u}\right){=}{0}\right]{,}\left[{{\mathrm{_ξ}}}_{{r}}{}\left({r}{,}{\mathrm{φ}}{,}{\mathrm{θ}}{,}{t}{,}{u}\right){=}{0}{,}{{\mathrm{_ξ}}}_{{\mathrm{φ}}}{}\left({r}{,}{\mathrm{φ}}{,}{\mathrm{θ}}{,}{t}{,}{u}\right){=}{\mathrm{sin}}{}\left({\mathrm{\phi }}\right){}{\mathrm{cot}}{}\left({\mathrm{\theta }}\right){,}{{\mathrm{_ξ}}}_{{\mathrm{θ}}}{}\left({r}{,}{\mathrm{φ}}{,}{\mathrm{θ}}{,}{t}{,}{u}\right){=}{-}{\mathrm{cos}}{}\left({\mathrm{\phi }}\right){,}{{\mathrm{_ξ}}}_{{t}}{}\left({r}{,}{\mathrm{φ}}{,}{\mathrm{θ}}{,}{t}{,}{u}\right){=}{0}{,}{{\mathrm{_η}}}_{{u}}{}\left({r}{,}{\mathrm{φ}}{,}{\mathrm{θ}}{,}{t}{,}{u}\right){=}{0}\right]{,}\left[{{\mathrm{_ξ}}}_{{r}}{}\left({r}{,}{\mathrm{φ}}{,}{\mathrm{θ}}{,}{t}{,}{u}\right){=}{-}{t}{}{\mathrm{cos}}{}\left({\mathrm{\theta }}\right){,}{{\mathrm{_ξ}}}_{{\mathrm{φ}}}{}\left({r}{,}{\mathrm{φ}}{,}{\mathrm{θ}}{,}{t}{,}{u}\right){=}{0}{,}{{\mathrm{_ξ}}}_{{\mathrm{θ}}}{}\left({r}{,}{\mathrm{φ}}{,}{\mathrm{θ}}{,}{t}{,}{u}\right){=}\frac{{t}{}{\mathrm{sin}}{}\left({\mathrm{\theta }}\right)}{{r}}{,}{{\mathrm{_ξ}}}_{{t}}{}\left({r}{,}{\mathrm{φ}}{,}{\mathrm{θ}}{,}{t}{,}{u}\right){=}{-}{r}{}{\mathrm{cos}}{}\left({\mathrm{\theta }}\right){,}{{\mathrm{_η}}}_{{u}}{}\left({r}{,}{\mathrm{φ}}{,}{\mathrm{θ}}{,}{t}{,}{u}\right){=}{0}\right]{,}\left[{{\mathrm{_ξ}}}_{{r}}{}\left({r}{,}{\mathrm{φ}}{,}{\mathrm{θ}}{,}{t}{,}{u}\right){=}{\mathrm{sin}}{}\left({\mathrm{\theta }}\right){}{\mathrm{sin}}{}\left({\mathrm{\phi }}\right){,}{{\mathrm{_ξ}}}_{{\mathrm{φ}}}{}\left({r}{,}{\mathrm{φ}}{,}{\mathrm{θ}}{,}{t}{,}{u}\right){=}\frac{{\mathrm{cos}}{}\left({\mathrm{\phi }}\right){}{\mathrm{csc}}{}\left({\mathrm{\theta }}\right)}{{r}}{,}{{\mathrm{_ξ}}}_{{\mathrm{θ}}}{}\left({r}{,}{\mathrm{φ}}{,}{\mathrm{θ}}{,}{t}{,}{u}\right){=}\frac{{\mathrm{sin}}{}\left({\mathrm{\phi }}\right){}{\mathrm{cos}}{}\left({\mathrm{\theta }}\right)}{{r}}{,}{{\mathrm{_ξ}}}_{{t}}{}\left({r}{,}{\mathrm{φ}}{,}{\mathrm{θ}}{,}{t}{,}{u}\right){=}{0}{,}{{\mathrm{_η}}}_{{u}}{}\left({r}{,}{\mathrm{φ}}{,}{\mathrm{θ}}{,}{t}{,}{u}\right){=}{0}\right]{,}\left[{{\mathrm{_ξ}}}_{{r}}{}\left({r}{,}{\mathrm{φ}}{,}{\mathrm{θ}}{,}{t}{,}{u}\right){=}{-}{\mathrm{cos}}{}\left({\mathrm{\phi }}\right){}{\mathrm{sin}}{}\left({\mathrm{\theta }}\right){,}{{\mathrm{_ξ}}}_{{\mathrm{φ}}}{}\left({r}{,}{\mathrm{φ}}{,}{\mathrm{θ}}{,}{t}{,}{u}\right){=}\frac{{\mathrm{sin}}{}\left({\mathrm{\phi }}\right){}{\mathrm{csc}}{}\left({\mathrm{\theta }}\right)}{{r}}{,}{{\mathrm{_ξ}}}_{{\mathrm{θ}}}{}\left({r}{,}{\mathrm{φ}}{,}{\mathrm{θ}}{,}{t}{,}{u}\right){=}{-}\frac{{\mathrm{cos}}{}\left({\mathrm{\phi }}\right){}{\mathrm{cos}}{}\left({\mathrm{\theta }}\right)}{{r}}{,}{{\mathrm{_ξ}}}_{{t}}{}\left({r}{,}{\mathrm{φ}}{,}{\mathrm{θ}}{,}{t}{,}{u}\right){=}{0}{,}{{\mathrm{_η}}}_{{u}}{}\left({r}{,}{\mathrm{φ}}{,}{\mathrm{θ}}{,}{t}{,}{u}\right){=}{0}\right]{,}\left[{{\mathrm{_ξ}}}_{{r}}{}\left({r}{,}{\mathrm{φ}}{,}{\mathrm{θ}}{,}{t}{,}{u}\right){=}\left({-}\frac{{{r}}^{{2}}}{{2}}{-}\frac{{{t}}^{{2}}}{{2}}\right){}{\mathrm{cos}}{}\left({\mathrm{\theta }}\right){,}{{\mathrm{_ξ}}}_{{\mathrm{φ}}}{}\left({r}{,}{\mathrm{φ}}{,}{\mathrm{θ}}{,}{t}{,}{u}\right){=}{0}{,}{{\mathrm{_ξ}}}_{{\mathrm{θ}}}{}\left({r}{,}{\mathrm{φ}}{,}{\mathrm{θ}}{,}{t}{,}{u}\right){=}{-}\frac{\left({{r}}^{{2}}{-}{{t}}^{{2}}\right){}{\mathrm{sin}}{}\left({\mathrm{\theta }}\right)}{{2}{}{r}}{,}{{\mathrm{_ξ}}}_{{t}}{}\left({r}{,}{\mathrm{φ}}{,}{\mathrm{θ}}{,}{t}{,}{u}\right){=}{-}{r}{}{\mathrm{cos}}{}\left({\mathrm{\theta }}\right){}{t}{,}{{\mathrm{_η}}}_{{u}}{}\left({r}{,}{\mathrm{φ}}{,}{\mathrm{θ}}{,}{t}{,}{u}\right){=}{r}{}{\mathrm{cos}}{}\left({\mathrm{\theta }}\right){}{u}\right]{,}\left[{{\mathrm{_ξ}}}_{{r}}{}\left({r}{,}{\mathrm{φ}}{,}{\mathrm{θ}}{,}{t}{,}{u}\right){=}{\mathrm{sin}}{}\left({\mathrm{\phi }}\right){}{t}{}{\mathrm{sin}}{}\left({\mathrm{\theta }}\right){,}{{\mathrm{_ξ}}}_{{\mathrm{φ}}}{}\left({r}{,}{\mathrm{φ}}{,}{\mathrm{θ}}{,}{t}{,}{u}\right){=}\frac{{t}{}{\mathrm{cos}}{}\left({\mathrm{\phi }}\right){}{\mathrm{csc}}{}\left({\mathrm{\theta }}\right)}{{r}}{,}{{\mathrm{_ξ}}}_{{\mathrm{θ}}}{}\left({r}{,}{\mathrm{φ}}{,}{\mathrm{θ}}{,}{t}{,}{u}\right){=}\frac{{t}{}{\mathrm{sin}}{}\left({\mathrm{\phi }}\right){}{\mathrm{cos}}{}\left({\mathrm{\theta }}\right)}{{r}}{,}{{\mathrm{_ξ}}}_{{t}}{}\left({r}{,}{\mathrm{φ}}{,}{\mathrm{θ}}{,}{t}{,}{u}\right){=}{\mathrm{sin}}{}\left({\mathrm{\theta }}\right){}{r}{}{\mathrm{sin}}{}\left({\mathrm{\phi }}\right){,}{{\mathrm{_η}}}_{{u}}{}\left({r}{,}{\mathrm{φ}}{,}{\mathrm{θ}}{,}{t}{,}{u}\right){=}{0}\right]{,}\left[{{\mathrm{_ξ}}}_{{r}}{}\left({r}{,}{\mathrm{φ}}{,}{\mathrm{θ}}{,}{t}{,}{u}\right){=}{-}{t}{}{\mathrm{cos}}{}\left({\mathrm{\phi }}\right){}{\mathrm{sin}}{}\left({\mathrm{\theta }}\right){,}{{\mathrm{_ξ}}}_{{\mathrm{φ}}}{}\left({r}{,}{\mathrm{φ}}{,}{\mathrm{θ}}{,}{t}{,}{u}\right){=}\frac{{t}{}{\mathrm{sin}}{}\left({\mathrm{\phi }}\right){}{\mathrm{csc}}{}\left({\mathrm{\theta }}\right)}{{r}}{,}{{\mathrm{_ξ}}}_{{\mathrm{θ}}}{}\left({r}{,}{\mathrm{φ}}{,}{\mathrm{θ}}{,}{t}{,}{u}\right){=}{-}\frac{{t}{}{\mathrm{cos}}{}\left({\mathrm{\phi }}\right){}{\mathrm{cos}}{}\left({\mathrm{\theta }}\right)}{{r}}{,}{{\mathrm{_ξ}}}_{{t}}{}\left({r}{,}{\mathrm{φ}}{,}{\mathrm{θ}}{,}{t}{,}{u}\right){=}{-}{\mathrm{sin}}{}\left({\mathrm{\theta }}\right){}{r}{}{\mathrm{cos}}{}\left({\mathrm{\phi }}\right){,}{{\mathrm{_η}}}_{{u}}{}\left({r}{,}{\mathrm{φ}}{,}{\mathrm{θ}}{,}{t}{,}{u}\right){=}{0}\right]{,}\left[{{\mathrm{_ξ}}}_{{r}}{}\left({r}{,}{\mathrm{φ}}{,}{\mathrm{θ}}{,}{t}{,}{u}\right){=}\left({-}\frac{{{r}}^{{2}}}{{2}}{-}\frac{{{t}}^{{2}}}{{2}}\right){}{\mathrm{cos}}{}\left({\mathrm{\phi }}\right){}{\mathrm{sin}}{}\left({\mathrm{\theta }}\right){,}{{\mathrm{_ξ}}}_{{\mathrm{φ}}}{}\left({r}{,}{\mathrm{φ}}{,}{\mathrm{θ}}{,}{t}{,}{u}\right){=}{-}\frac{\left({{r}}^{{2}}{-}{{t}}^{{2}}\right){}{\mathrm{sin}}{}\left({\mathrm{\phi }}\right){}{\mathrm{csc}}{}\left({\mathrm{\theta }}\right)}{{2}{}{r}}{,}{{\mathrm{_ξ}}}_{{\mathrm{θ}}}{}\left({r}{,}{\mathrm{φ}}{,}{\mathrm{θ}}{,}{t}{,}{u}\right){=}\frac{\left({{r}}^{{2}}{-}{{t}}^{{2}}\right){}{\mathrm{cos}}{}\left({\mathrm{\phi }}\right){}{\mathrm{cos}}{}\left({\mathrm{\theta }}\right)}{{2}{}{r}}{,}{{\mathrm{_ξ}}}_{{t}}{}\left({r}{,}{\mathrm{φ}}{,}{\mathrm{θ}}{,}{t}{,}{u}\right){=}{-}{t}{}{\mathrm{cos}}{}\left({\mathrm{\phi }}\right){}{r}{}{\mathrm{sin}}{}\left({\mathrm{\theta }}\right){,}{{\mathrm{_η}}}_{{u}}{}\left({r}{,}{\mathrm{φ}}{,}{\mathrm{θ}}{,}{t}{,}{u}\right){=}{r}{}{\mathrm{sin}}{}\left({\mathrm{\theta }}\right){}{\mathrm{cos}}{}\left({\mathrm{\phi }}\right){}{u}\right]{,}\left[{{\mathrm{_ξ}}}_{{r}}{}\left({r}{,}{\mathrm{φ}}{,}{\mathrm{θ}}{,}{t}{,}{u}\right){=}\left(\frac{{{r}}^{{2}}}{{2}}{+}\frac{{{t}}^{{2}}}{{2}}\right){}{\mathrm{sin}}{}\left({\mathrm{\theta }}\right){}{\mathrm{sin}}{}\left({\mathrm{\phi }}\right){,}{{\mathrm{_ξ}}}_{{\mathrm{φ}}}{}\left({r}{,}{\mathrm{φ}}{,}{\mathrm{θ}}{,}{t}{,}{u}\right){=}{-}\frac{\left({{r}}^{{2}}{-}{{t}}^{{2}}\right){}{\mathrm{cos}}{}\left({\mathrm{\phi }}\right){}{\mathrm{csc}}{}\left({\mathrm{\theta }}\right)}{{2}{}{r}}{,}{{\mathrm{_ξ}}}_{{\mathrm{θ}}}{}\left({r}{,}{\mathrm{φ}}{,}{\mathrm{θ}}{,}{t}{,}{u}\right){=}{-}\frac{\left({{r}}^{{2}}{-}{{t}}^{{2}}\right){}{\mathrm{sin}}{}\left({\mathrm{\phi }}\right){}{\mathrm{cos}}{}\left({\mathrm{\theta }}\right)}{{2}{}{r}}{,}{{\mathrm{_ξ}}}_{{t}}{}\left({r}{,}{\mathrm{φ}}{,}{\mathrm{θ}}{,}{t}{,}{u}\right){=}{\mathrm{sin}}{}\left({\mathrm{\phi }}\right){}{t}{}{r}{}{\mathrm{sin}}{}\left({\mathrm{\theta }}\right){,}{{\mathrm{_η}}}_{{u}}{}\left({r}{,}{\mathrm{φ}}{,}{\mathrm{θ}}{,}{t}{,}{u}\right){=}{-}{r}{}{\mathrm{sin}}{}\left({\mathrm{\theta }}\right){}{\mathrm{sin}}{}\left({\mathrm{\phi }}\right){}{u}\right]$ (14)

An ODE example (see more commands for handling ODE symmetry problems also in DEtools,Lie) and a dynamical symmetry represented by infinitesimals depending on the first derivative y'

 > $\mathrm{declare}\left(y\left(x\right),\mathrm{prime}=x\right)$
 ${y}{}\left({x}\right){}{\mathrm{will now be displayed as}}{}{y}$
 ${\mathrm{derivatives with respect to}}{}{x}{}{\mathrm{of functions of one variable will now be displayed with \text{'}}}$ (15)
 > $\mathrm{ODE2}≔\frac{{ⅆ}^{2}}{ⅆ{x}^{2}}y\left(x\right)=\mathrm{Φ}\left(y\left(x\right)\right)$
 ${\mathrm{diff}}{}\left({\mathrm{diff}}{}\left({y}{}\left({x}\right){,}{x}\right){,}{x}\right){=}{\mathrm{Φ}}{}\left({y}{}\left({x}\right)\right)$ (16)
 > $\mathrm{DS}≔\left[\mathrm{_ξ}=f\left(x,y,\mathrm{_y1}\right),\mathrm{_η}=0\right]$
 ${\mathrm{DS}}{≔}\left[{\mathrm{_ξ}}{=}{f}{}\left({x}{,}{y}{,}{\mathrm{_y1}}\right){,}{\mathrm{_η}}{=}{0}\right]$ (17)

The corresponding determining system directly in jetODE notation, that is the notation used in DEtools,Lie, the package for symmetries and integrating factors for ODEs developed before the symmetry commands of PDEtools

 > $\mathrm{DeterminingPDE}\left(\mathrm{ODE2},y\left(x\right),\mathrm{DS},'\mathrm{jetnotation}'='\mathrm{jetODE}'\right)$
 $\left\{{\mathrm{diff}}{}\left({\mathrm{diff}}{}\left({f}{}\left({x}{,}{y}{,}{\mathrm{_y1}}\right){,}{x}\right){,}{x}\right){=}\frac{{-}\left({\mathrm{diff}}{}\left({\mathrm{diff}}{}\left({f}{}\left({x}{,}{y}{,}{\mathrm{_y1}}\right){,}{\mathrm{_y1}}\right){,}{\mathrm{_y1}}\right)\right){}{\mathrm{_y1}}{}{{\mathrm{Φ}}{}\left({y}\right)}^{{2}}{-}{2}{}\left({\mathrm{diff}}{}\left({\mathrm{diff}}{}\left({f}{}\left({x}{,}{y}{,}{\mathrm{_y1}}\right){,}{\mathrm{_y1}}\right){,}{y}\right)\right){}{{\mathrm{_y1}}}^{{2}}{}{\mathrm{Φ}}{}\left({y}\right){-}\left({\mathrm{diff}}{}\left({\mathrm{Φ}}{}\left({y}\right){,}{y}\right)\right){}{{\mathrm{_y1}}}^{{2}}{}\left({\mathrm{diff}}{}\left({f}{}\left({x}{,}{y}{,}{\mathrm{_y1}}\right){,}{\mathrm{_y1}}\right)\right){-}\left({\mathrm{diff}}{}\left({\mathrm{diff}}{}\left({f}{}\left({x}{,}{y}{,}{\mathrm{_y1}}\right){,}{y}\right){,}{y}\right)\right){}{{\mathrm{_y1}}}^{{3}}{-}{2}{}\left({\mathrm{diff}}{}\left({f}{}\left({x}{,}{y}{,}{\mathrm{_y1}}\right){,}{\mathrm{_y1}}\right)\right){}{{\mathrm{Φ}}{}\left({y}\right)}^{{2}}{-}{3}{}\left({\mathrm{diff}}{}\left({f}{}\left({x}{,}{y}{,}{\mathrm{_y1}}\right){,}{y}\right)\right){}{\mathrm{_y1}}{}{\mathrm{Φ}}{}\left({y}\right){-}{2}{}\left({\mathrm{diff}}{}\left({\mathrm{diff}}{}\left({f}{}\left({x}{,}{y}{,}{\mathrm{_y1}}\right){,}{\mathrm{_y1}}\right){,}{x}\right)\right){}{\mathrm{_y1}}{}{\mathrm{Φ}}{}\left({y}\right){-}{2}{}\left({\mathrm{diff}}{}\left({\mathrm{diff}}{}\left({f}{}\left({x}{,}{y}{,}{\mathrm{_y1}}\right){,}{x}\right){,}{y}\right)\right){}{{\mathrm{_y1}}}^{{2}}{-}{2}{}\left({\mathrm{diff}}{}\left({f}{}\left({x}{,}{y}{,}{\mathrm{_y1}}\right){,}{x}\right)\right){}{\mathrm{Φ}}{}\left({y}\right)}{{\mathrm{_y1}}}\right\}$ (18)

An ODE in disguise, appearing as a PDE, and the corresponding determining PDE system for a very particular form $\mathrm{S2}$ of its infinitesimals, depending on two arbitrary functions $\mathrm{f__1}$ and $\mathrm{f__2}$

 > $\mathrm{ODE2}≔\frac{\partial }{\partial x}z\left(x,y\right)=\frac{yz\left(x,y\right)}{x}$
 ${\mathrm{diff}}{}\left({z}{}\left({x}{,}{y}\right){,}{x}\right){=}\frac{{y}{}{z}{}\left({x}{,}{y}\right)}{{x}}$ (19)
 > $\mathrm{S2}≔\left[0,\mathrm{_F1}\left(y,z{x}^{-y}\right),{x}^{y}z{x}^{-y}\mathrm{_F1}\left(y,z{x}^{-y}\right)\mathrm{ln}\left(x\right)+{x}^{y}\mathrm{_F2}\left(y,z{x}^{-y}\right)\right]$
 ${\mathrm{S2}}{≔}\left[{0}{,}\mathrm{f__1}{}\left({y}{,}{z}{}{{x}}^{{-}{y}}\right){,}{{x}}^{{y}}{}{z}{}{{x}}^{{-}{y}}{}\mathrm{f__1}{}\left({y}{,}{z}{}{{x}}^{{-}{y}}\right){}{\mathrm{ln}}{}\left({x}\right){+}{{x}}^{{y}}{}\mathrm{f__2}{}\left({y}{,}{z}{}{{x}}^{{-}{y}}\right)\right]$ (20)
 > $\mathrm{DeterminingPDE}\left(\mathrm{ODE2},\mathrm{S2}\right)$
 $\left\{{0}\right\}$ (21)

The output above is zero because S2 is already a symmetry of ODE2, thus canceling completely the determining system. To test symmetries given as lists of infinitesimals use the SymmetryTest command

 > $\mathrm{SymmetryTest}\left(\mathrm{S2},\mathrm{ODE2}\right)$
 $\left\{{0}\right\}$ (22)

The determining system for symmetries in evolutionary form for an ODE, directly in jetODE notation

 > $\mathrm{ODE3}≔\frac{{ⅆ}^{2}}{ⅆ{x}^{2}}y\left(x\right)=xy\left(x\right)\left(\frac{ⅆ}{ⅆx}y\left(x\right)\right)$
 ${\mathrm{diff}}{}\left({\mathrm{diff}}{}\left({y}{}\left({x}\right){,}{x}\right){,}{x}\right){=}{x}{}{y}{}\left({x}\right){}\left({\mathrm{diff}}{}\left({y}{}\left({x}\right){,}{x}\right)\right)$ (23)
 > $\mathrm{DeterminingPDE}\left(\mathrm{ODE3},\mathrm{typeofsymmetry}=\mathrm{evolutionary},\mathrm{notation}=\mathrm{jetODE}\right)$
 $\left\{{\mathrm{diff}}{}\left({\mathrm{diff}}{}\left({{\mathrm{_η}}}_{{y}}{}\left({x}{,}{y}{,}{\mathrm{_y1}}\right){,}{x}\right){,}{x}\right){=}{-}\left({\mathrm{diff}}{}\left({\mathrm{diff}}{}\left({{\mathrm{_η}}}_{{y}}{}\left({x}{,}{y}{,}{\mathrm{_y1}}\right){,}{\mathrm{_y1}}\right){,}{\mathrm{_y1}}\right)\right){}{{x}}^{{2}}{}{{y}}^{{2}}{}{{\mathrm{_y1}}}^{{2}}{-}{2}{}\left({\mathrm{diff}}{}\left({\mathrm{diff}}{}\left({{\mathrm{_η}}}_{{y}}{}\left({x}{,}{y}{,}{\mathrm{_y1}}\right){,}{\mathrm{_y1}}\right){,}{y}\right)\right){}{{\mathrm{_y1}}}^{{2}}{}{x}{}{y}{-}\left({\mathrm{diff}}{}\left({{\mathrm{_η}}}_{{y}}{}\left({x}{,}{y}{,}{\mathrm{_y1}}\right){,}{\mathrm{_y1}}\right)\right){}{{\mathrm{_y1}}}^{{2}}{}{x}{-}{2}{}{x}{}{y}{}{\mathrm{_y1}}{}\left({\mathrm{diff}}{}\left({\mathrm{diff}}{}\left({{\mathrm{_η}}}_{{y}}{}\left({x}{,}{y}{,}{\mathrm{_y1}}\right){,}{\mathrm{_y1}}\right){,}{x}\right)\right){-}\left({\mathrm{diff}}{}\left({\mathrm{diff}}{}\left({{\mathrm{_η}}}_{{y}}{}\left({x}{,}{y}{,}{\mathrm{_y1}}\right){,}{y}\right){,}{y}\right)\right){}{{\mathrm{_y1}}}^{{2}}{-}{y}{}\left({\mathrm{diff}}{}\left({{\mathrm{_η}}}_{{y}}{}\left({x}{,}{y}{,}{\mathrm{_y1}}\right){,}{\mathrm{_y1}}\right)\right){}{\mathrm{_y1}}{+}\left({\mathrm{diff}}{}\left({{\mathrm{_η}}}_{{y}}{}\left({x}{,}{y}{,}{\mathrm{_y1}}\right){,}{x}\right)\right){}{x}{}{y}{+}{{\mathrm{_η}}}_{{y}}{}\left({x}{,}{y}{,}{\mathrm{_y1}}\right){}{x}{}{\mathrm{_y1}}{-}{2}{}{\mathrm{_y1}}{}\left({\mathrm{diff}}{}\left({\mathrm{diff}}{}\left({{\mathrm{_η}}}_{{y}}{}\left({x}{,}{y}{,}{\mathrm{_y1}}\right){,}{x}\right){,}{y}\right)\right)\right\}$ (24)

By default DeterminingPDE processes DE systems as given. If all the integrability conditions that can be derived from the system's equations are explicitly present, then the resulting determining system for the symmetry infinitesimals has for solution all the symmetries. This is not always true otherwise. Consider for instance:

 > $\mathrm{declare}\left(u\left(x,y\right)\right)$
 ${u}{}\left({x}{,}{y}\right){}{\mathrm{will now be displayed as}}{}{u}$ (25)
 > $U≔\mathrm{diff_table}\left(u\left(x,y\right)\right):$
 > $\mathrm{DE}≔\left[{U}_{x}-{U}_{y}=0,{U}_{y,y}=0\right]$
 $\left[{\mathrm{diff}}{}\left({u}{}\left({x}{,}{y}\right){,}{x}\right){-}\left({\mathrm{diff}}{}\left({u}{}\left({x}{,}{y}\right){,}{y}\right)\right){=}{0}{,}{\mathrm{diff}}{}\left({\mathrm{diff}}{}\left({u}{}\left({x}{,}{y}\right){,}{y}\right){,}{y}\right){=}{0}\right]$ (26)

The determining system for the point symmetries of DE and its corresponding solution are

 > $\mathrm{detsys}≔\mathrm{DeterminingPDE}\left(\mathrm{DE}\right)$
 $\left\{{\mathrm{diff}}{}\left({\mathrm{diff}}{}\left({\mathrm{diff}}{}\left({{\mathrm{_η}}}_{{u}}{}\left({x}{,}{y}{,}{u}\right){,}{u}\right){,}{u}\right){,}{u}\right){=}{0}{,}{\mathrm{diff}}{}\left({\mathrm{diff}}{}\left({\mathrm{diff}}{}\left({{\mathrm{_η}}}_{{u}}{}\left({x}{,}{y}{,}{u}\right){,}{u}\right){,}{u}\right){,}{y}\right){=}{0}{,}{\mathrm{diff}}{}\left({\mathrm{diff}}{}\left({{\mathrm{_η}}}_{{u}}{}\left({x}{,}{y}{,}{u}\right){,}{y}\right){,}{y}\right){=}{0}{,}{\mathrm{diff}}{}\left({\mathrm{diff}}{}\left({{\mathrm{_ξ}}}_{{y}}{}\left({x}{,}{y}{,}{u}\right){,}{u}\right){,}{u}\right){=}{0}{,}{\mathrm{diff}}{}\left({\mathrm{diff}}{}\left({{\mathrm{_ξ}}}_{{y}}{}\left({x}{,}{y}{,}{u}\right){,}{u}\right){,}{y}\right){=}\frac{{1}}{{2}}{}{\mathrm{diff}}{}\left({\mathrm{diff}}{}\left({{\mathrm{_η}}}_{{u}}{}\left({x}{,}{y}{,}{u}\right){,}{u}\right){,}{u}\right){,}{\mathrm{diff}}{}\left({\mathrm{diff}}{}\left({{\mathrm{_ξ}}}_{{y}}{}\left({x}{,}{y}{,}{u}\right){,}{y}\right){,}{y}\right){=}{2}{}\left({\mathrm{diff}}{}\left({\mathrm{diff}}{}\left({{\mathrm{_η}}}_{{u}}{}\left({x}{,}{y}{,}{u}\right){,}{u}\right){,}{y}\right)\right){,}{\mathrm{diff}}{}\left({{\mathrm{_η}}}_{{u}}{}\left({x}{,}{y}{,}{u}\right){,}{x}\right){=}{\mathrm{diff}}{}\left({{\mathrm{_η}}}_{{u}}{}\left({x}{,}{y}{,}{u}\right){,}{y}\right){,}{\mathrm{diff}}{}\left({{\mathrm{_ξ}}}_{{x}}{}\left({x}{,}{y}{,}{u}\right){,}{u}\right){=}{0}{,}{\mathrm{diff}}{}\left({{\mathrm{_ξ}}}_{{x}}{}\left({x}{,}{y}{,}{u}\right){,}{y}\right){=}{0}{,}{\mathrm{diff}}{}\left({{\mathrm{_ξ}}}_{{y}}{}\left({x}{,}{y}{,}{u}\right){,}{x}\right){=}{\mathrm{diff}}{}\left({{\mathrm{_ξ}}}_{{y}}{}\left({x}{,}{y}{,}{u}\right){,}{y}\right){-}\left({\mathrm{diff}}{}\left({{\mathrm{_ξ}}}_{{x}}{}\left({x}{,}{y}{,}{u}\right){,}{x}\right)\right)\right\}$ (27)
 > $\mathrm{pdsolve}\left(\mathrm{detsys}\right)$
 $\left\{{{\mathrm{_η}}}_{{u}}{}\left({x}{,}{y}{,}{u}\right){=}\frac{\mathrm{c__3}{}{{u}}^{{2}}}{{2}}{+}\left(\left({x}{+}{y}\right){}\mathrm{c__1}{+}\mathrm{c__4}\right){}{u}{+}\left({x}{+}{y}\right){}\mathrm{c__2}{+}\mathrm{c__5}{,}{{\mathrm{_ξ}}}_{{x}}{}\left({x}{,}{y}{,}{u}\right){=}\mathrm{f__1}{}\left({x}\right){,}{{\mathrm{_ξ}}}_{{y}}{}\left({x}{,}{y}{,}{u}\right){=}{-}\mathrm{f__1}{}\left({x}\right){+}\mathrm{c__1}{}{{x}}^{{2}}{+}\frac{\left({4}{}\mathrm{c__1}{}{y}{+}\mathrm{c__3}{}{u}{+}{2}{}\mathrm{c__7}\right){}{x}}{{2}}{+}\mathrm{c__1}{}{{y}}^{{2}}{+}\frac{\left(\mathrm{c__3}{}{u}{+}{2}{}\mathrm{c__7}\right){}{y}}{{2}}{+}\mathrm{c__6}{}{u}{+}\mathrm{c__8}\right\}$ (28)

The system DE, however, does not include the integrability condition ${u}_{x,x}=0$ that can be derived from its equations. To compute a version of DE that is closed, i.e., that it includes all its integrability conditions explicitly you can use the Library routine CloseSystem or directly pass the closesystem option to DeterminingPDE

 > $\mathrm{closed_DE}≔\mathrm{PDEtools}:-\mathrm{Library}:-\mathrm{CloseSystem}\left(\mathrm{DE},\left[u\left(x,y\right)\right]\right)$
 $\left\{{\mathrm{diff}}{}\left({u}{}\left({x}{,}{y}\right){,}{x}\right){-}\left({\mathrm{diff}}{}\left({u}{}\left({x}{,}{y}\right){,}{y}\right)\right){,}{\mathrm{diff}}{}\left({\mathrm{diff}}{}\left({u}{}\left({x}{,}{y}\right){,}{x}\right){,}{x}\right){,}{\mathrm{diff}}{}\left({\mathrm{diff}}{}\left({u}{}\left({x}{,}{y}\right){,}{x}\right){,}{y}\right){,}{\mathrm{diff}}{}\left({\mathrm{diff}}{}\left({u}{}\left({x}{,}{y}\right){,}{y}\right){,}{y}\right)\right\}$ (29)
 > $\mathrm{DeterminingPDE}\left(\mathrm{closed_DE}\right)$
 $\left\{{\mathrm{diff}}{}\left({\mathrm{diff}}{}\left({\mathrm{diff}}{}\left({{\mathrm{_η}}}_{{u}}{}\left({x}{,}{y}{,}{u}\right){,}{u}\right){,}{u}\right){,}{u}\right){=}{0}{,}{\mathrm{diff}}{}\left({\mathrm{diff}}{}\left({\mathrm{diff}}{}\left({{\mathrm{_η}}}_{{u}}{}\left({x}{,}{y}{,}{u}\right){,}{u}\right){,}{u}\right){,}{y}\right){=}{0}{,}{\mathrm{diff}}{}\left({\mathrm{diff}}{}\left({{\mathrm{_η}}}_{{u}}{}\left({x}{,}{y}{,}{u}\right){,}{y}\right){,}{y}\right){=}{0}{,}{\mathrm{diff}}{}\left({\mathrm{diff}}{}\left({{\mathrm{_ξ}}}_{{y}}{}\left({x}{,}{y}{,}{u}\right){,}{u}\right){,}{u}\right){=}{-}\left({\mathrm{diff}}{}\left({\mathrm{diff}}{}\left({{\mathrm{_ξ}}}_{{x}}{}\left({x}{,}{y}{,}{u}\right){,}{u}\right){,}{u}\right)\right){,}{\mathrm{diff}}{}\left({\mathrm{diff}}{}\left({{\mathrm{_ξ}}}_{{y}}{}\left({x}{,}{y}{,}{u}\right){,}{u}\right){,}{y}\right){=}{-}\left({\mathrm{diff}}{}\left({\mathrm{diff}}{}\left({{\mathrm{_ξ}}}_{{x}}{}\left({x}{,}{y}{,}{u}\right){,}{u}\right){,}{y}\right)\right){+}\frac{{1}}{{2}}{}{\mathrm{diff}}{}\left({\mathrm{diff}}{}\left({{\mathrm{_η}}}_{{u}}{}\left({x}{,}{y}{,}{u}\right){,}{u}\right){,}{u}\right){,}{\mathrm{diff}}{}\left({\mathrm{diff}}{}\left({{\mathrm{_ξ}}}_{{y}}{}\left({x}{,}{y}{,}{u}\right){,}{y}\right){,}{y}\right){=}{2}{}\left({\mathrm{diff}}{}\left({\mathrm{diff}}{}\left({{\mathrm{_η}}}_{{u}}{}\left({x}{,}{y}{,}{u}\right){,}{u}\right){,}{y}\right)\right){-}\left({\mathrm{diff}}{}\left({\mathrm{diff}}{}\left({{\mathrm{_ξ}}}_{{x}}{}\left({x}{,}{y}{,}{u}\right){,}{y}\right){,}{y}\right)\right){,}{\mathrm{diff}}{}\left({{\mathrm{_η}}}_{{u}}{}\left({x}{,}{y}{,}{u}\right){,}{x}\right){=}{\mathrm{diff}}{}\left({{\mathrm{_η}}}_{{u}}{}\left({x}{,}{y}{,}{u}\right){,}{y}\right){,}{\mathrm{diff}}{}\left({{\mathrm{_ξ}}}_{{y}}{}\left({x}{,}{y}{,}{u}\right){,}{x}\right){=}{\mathrm{diff}}{}\left({{\mathrm{_ξ}}}_{{x}}{}\left({x}{,}{y}{,}{u}\right){,}{y}\right){+}{\mathrm{diff}}{}\left({{\mathrm{_ξ}}}_{{y}}{}\left({x}{,}{y}{,}{u}\right){,}{y}\right){-}\left({\mathrm{diff}}{}\left({{\mathrm{_ξ}}}_{{x}}{}\left({x}{,}{y}{,}{u}\right){,}{x}\right)\right)\right\}$ (30)
 > $\mathrm{DeterminingPDE}\left(\mathrm{DE},\mathrm{closesystem}\right)$
 $\left\{{\mathrm{diff}}{}\left({\mathrm{diff}}{}\left({\mathrm{diff}}{}\left({{\mathrm{_η}}}_{{u}}{}\left({x}{,}{y}{,}{u}\right){,}{u}\right){,}{u}\right){,}{u}\right){=}{0}{,}{\mathrm{diff}}{}\left({\mathrm{diff}}{}\left({\mathrm{diff}}{}\left({{\mathrm{_η}}}_{{u}}{}\left({x}{,}{y}{,}{u}\right){,}{u}\right){,}{u}\right){,}{y}\right){=}{0}{,}{\mathrm{diff}}{}\left({\mathrm{diff}}{}\left({{\mathrm{_η}}}_{{u}}{}\left({x}{,}{y}{,}{u}\right){,}{y}\right){,}{y}\right){=}{0}{,}{\mathrm{diff}}{}\left({\mathrm{diff}}{}\left({{\mathrm{_ξ}}}_{{y}}{}\left({x}{,}{y}{,}{u}\right){,}{u}\right){,}{u}\right){=}{-}\left({\mathrm{diff}}{}\left({\mathrm{diff}}{}\left({{\mathrm{_ξ}}}_{{x}}{}\left({x}{,}{y}{,}{u}\right){,}{u}\right){,}{u}\right)\right){,}{\mathrm{diff}}{}\left({\mathrm{diff}}{}\left({{\mathrm{_ξ}}}_{{y}}{}\left({x}{,}{y}{,}{u}\right){,}{u}\right){,}{y}\right){=}{-}\left({\mathrm{diff}}{}\left({\mathrm{diff}}{}\left({{\mathrm{_ξ}}}_{{x}}{}\left({x}{,}{y}{,}{u}\right){,}{u}\right){,}{y}\right)\right){+}\frac{{1}}{{2}}{}{\mathrm{diff}}{}\left({\mathrm{diff}}{}\left({{\mathrm{_η}}}_{{u}}{}\left({x}{,}{y}{,}{u}\right){,}{u}\right){,}{u}\right){,}{\mathrm{diff}}{}\left({\mathrm{diff}}{}\left({{\mathrm{_ξ}}}_{{y}}{}\left({x}{,}{y}{,}{u}\right){,}{y}\right){,}{y}\right){=}{2}{}\left({\mathrm{diff}}{}\left({\mathrm{diff}}{}\left({{\mathrm{_η}}}_{{u}}{}\left({x}{,}{y}{,}{u}\right){,}{u}\right){,}{y}\right)\right){-}\left({\mathrm{diff}}{}\left({\mathrm{diff}}{}\left({{\mathrm{_ξ}}}_{{x}}{}\left({x}{,}{y}{,}{u}\right){,}{y}\right){,}{y}\right)\right){,}{\mathrm{diff}}{}\left({{\mathrm{_η}}}_{{u}}{}\left({x}{,}{y}{,}{u}\right){,}{x}\right){=}{\mathrm{diff}}{}\left({{\mathrm{_η}}}_{{u}}{}\left({x}{,}{y}{,}{u}\right){,}{y}\right){,}{\mathrm{diff}}{}\left({{\mathrm{_ξ}}}_{{y}}{}\left({x}{,}{y}{,}{u}\right){,}{x}\right){=}{\mathrm{diff}}{}\left({{\mathrm{_ξ}}}_{{x}}{}\left({x}{,}{y}{,}{u}\right){,}{y}\right){+}{\mathrm{diff}}{}\left({{\mathrm{_ξ}}}_{{y}}{}\left({x}{,}{y}{,}{u}\right){,}{y}\right){-}\left({\mathrm{diff}}{}\left({{\mathrm{_ξ}}}_{{x}}{}\left({x}{,}{y}{,}{u}\right){,}{x}\right)\right)\right\}$ (31)

The resulting determining system for the closed version of DE admits solutions more general than (6.28)

 > $\mathrm{pdsolve}\left(\right)$
 $\left\{{{\mathrm{_η}}}_{{u}}{}\left({x}{,}{y}{,}{u}\right){=}\frac{\mathrm{c__3}{}{{u}}^{{2}}}{{2}}{+}\left(\left({x}{+}{y}\right){}\mathrm{c__1}{+}\mathrm{c__4}\right){}{u}{+}\left({x}{+}{y}\right){}\mathrm{c__2}{+}\mathrm{c__5}{,}{{\mathrm{_ξ}}}_{{x}}{}\left({x}{,}{y}{,}{u}\right){=}{{\mathrm{_ξ}}}_{{x}}{}\left({x}{,}{y}{,}{u}\right){,}{{\mathrm{_ξ}}}_{{y}}{}\left({x}{,}{y}{,}{u}\right){=}{\left({x}{+}{y}\right)}^{{2}}{}\mathrm{c__1}{-}{{\mathrm{_ξ}}}_{{x}}{}\left({x}{,}{y}{,}{u}\right){+}\frac{\left(\left({x}{+}{y}\right){}\mathrm{c__3}{+}{2}{}\mathrm{c__6}\right){}{u}}{{2}}{+}\left({x}{+}{y}\right){}\mathrm{c__7}{+}\mathrm{c__8}\right\}$ (32)

Comparing, in (6.32) ${\mathrm{_ξ}}_{x}\left(x,y,u\right)$ is entirely arbitrary while in the infinitesimals (6.28) it is not. The closesystem option is also available when calling Infinitesimals, InvariantSolutions and SimilaritySolutions - all of them will forward the option to DeterminingPDE. For example: the same most general form of the infinitesimals (6.32) can be computed directly via

 > $\mathrm{Infinitesimals}\left(\mathrm{DE},\mathrm{closesystem},\mathrm{dis}=\mathrm{false}\right)$
 ${\mathrm{* Partial match of \text{'}dis\text{'} against keyword \text{'}displayfunctionality\text{'}}}$
 $\left[{{\mathrm{_ξ}}}_{{x}}{=}\mathrm{f__1}{}\left({x}{,}{y}{,}{u}\right){,}{{\mathrm{_ξ}}}_{{y}}{=}{-}\mathrm{f__1}{}\left({x}{,}{y}{,}{u}\right){,}{{\mathrm{_η}}}_{{u}}{=}{1}\right]{,}\left[{{\mathrm{_ξ}}}_{{x}}{=}\mathrm{f__2}{}\left({x}{,}{y}{,}{u}\right){,}{{\mathrm{_ξ}}}_{{y}}{=}{-}\mathrm{f__2}{}\left({x}{,}{y}{,}{u}\right){+}{1}{,}{{\mathrm{_η}}}_{{u}}{=}{0}\right]{,}\left[{{\mathrm{_ξ}}}_{{x}}{=}\mathrm{f__3}{}\left({x}{,}{y}{,}{u}\right){,}{{\mathrm{_ξ}}}_{{y}}{=}{-}\mathrm{f__3}{}\left({x}{,}{y}{,}{u}\right){,}{{\mathrm{_η}}}_{{u}}{=}{u}\right]{,}\left[{{\mathrm{_ξ}}}_{{x}}{=}\mathrm{f__4}{}\left({x}{,}{y}{,}{u}\right){,}{{\mathrm{_ξ}}}_{{y}}{=}{-}\mathrm{f__4}{}\left({x}{,}{y}{,}{u}\right){+}{x}{+}{y}{,}{{\mathrm{_η}}}_{{u}}{=}{0}\right]{,}\left[{{\mathrm{_ξ}}}_{{x}}{=}\mathrm{f__5}{}\left({x}{,}{y}{,}{u}\right){,}{{\mathrm{_ξ}}}_{{y}}{=}{-}\mathrm{f__5}{}\left({x}{,}{y}{,}{u}\right){+}{u}{,}{{\mathrm{_η}}}_{{u}}{=}{0}\right]{,}\left[{{\mathrm{_ξ}}}_{{x}}{=}\mathrm{f__6}{}\left({x}{,}{y}{,}{u}\right){,}{{\mathrm{_ξ}}}_{{y}}{=}{-}\mathrm{f__6}{}\left({x}{,}{y}{,}{u}\right){,}{{\mathrm{_η}}}_{{u}}{=}{x}{+}{y}\right]{,}\left[{{\mathrm{_ξ}}}_{{x}}{=}\mathrm{f__7}{}\left({x}{,}{y}{,}{u}\right){,}{{\mathrm{_ξ}}}_{{y}}{=}{-}\mathrm{f__7}{}\left({x}{,}{y}{,}{u}\right){+}\frac{\left({x}{+}{y}\right){}{u}}{{2}}{,}{{\mathrm{_η}}}_{{u}}{=}\frac{{{u}}^{{2}}}{{2}}\right]{,}\left[{{\mathrm{_ξ}}}_{{x}}{=}{\mathrm{_F8}}{}\left({x}{,}{y}{,}{u}\right){,}{{\mathrm{_ξ}}}_{{y}}{=}{-}{\mathrm{_F8}}{}\left({x}{,}{y}{,}{u}\right){+}{{x}}^{{2}}{+}{2}{}{x}{}{y}{+}{{y}}^{{2}}{,}{{\mathrm{_η}}}_{{u}}{=}\left({x}{+}{y}\right){}{u}\right]$ (33)

The DeterminingPDE command also works with anticommutative variables, natively, without using the approach explained in PerformOnAnticommutativeSystem.

 > $\mathrm{with}\left(\mathrm{Physics}\right)$
 $\left[{\mathrm{*}}{,}{\mathrm{.}}{,}{\mathrm{Annihilation}}{,}{\mathrm{AntiCommutator}}{,}{\mathrm{Antisymmetrize}}{,}{\mathrm{Assume}}{,}{\mathrm{Bra}}{,}{\mathrm{Bracket}}{,}{\mathrm{Check}}{,}{\mathrm{Christoffel}}{,}{\mathrm{Coefficients}}{,}{\mathrm{Commutator}}{,}{\mathrm{CompactDisplay}}{,}{\mathrm{Coordinates}}{,}{\mathrm{Creation}}{,}{\mathrm{D_}}{,}{\mathrm{Dagger}}{,}{\mathrm{Decompose}}{,}{\mathrm{Define}}{,}{\mathrm{Dγ}}{,}{\mathrm{Einstein}}{,}{\mathrm{EnergyMomentum}}{,}{\mathrm{Expand}}{,}{\mathrm{ExteriorDerivative}}{,}{\mathrm{Factor}}{,}{\mathrm{FeynmanDiagrams}}{,}{\mathrm{FeynmanIntegral}}{,}{\mathrm{Fundiff}}{,}{\mathrm{Geodesics}}{,}{\mathrm{GrassmannParity}}{,}{\mathrm{Gtaylor}}{,}{\mathrm{Intc}}{,}{\mathrm{Inverse}}{,}{\mathrm{Ket}}{,}{\mathrm{KillingVectors}}{,}{\mathrm{KroneckerDelta}}{,}{\mathrm{LagrangeEquations}}{,}{\mathrm{LeviCivita}}{,}{\mathrm{Library}}{,}{\mathrm{LieBracket}}{,}{\mathrm{LieDerivative}}{,}{\mathrm{Normal}}{,}{\mathrm{NumericalRelativity}}{,}{\mathrm{Parameters}}{,}{\mathrm{PerformOnAnticommutativeSystem}}{,}{\mathrm{Projector}}{,}{\mathrm{Psigma}}{,}{\mathrm{Redefine}}{,}{\mathrm{Ricci}}{,}{\mathrm{Riemann}}{,}{\mathrm{Setup}}{,}{\mathrm{Simplify}}{,}{\mathrm{SortProducts}}{,}{\mathrm{SpaceTimeVector}}{,}{\mathrm{StandardModel}}{,}{\mathrm{Substitute}}{,}{\mathrm{SubstituteTensor}}{,}{\mathrm{SubstituteTensorIndices}}{,}{\mathrm{SumOverRepeatedIndices}}{,}{\mathrm{Symmetrize}}{,}{\mathrm{TensorArray}}{,}{\mathrm{Tetrads}}{,}{\mathrm{ThreePlusOne}}{,}{\mathrm{ToContravariant}}{,}{\mathrm{ToCovariant}}{,}{\mathrm{ToFieldComponents}}{,}{\mathrm{ToSuperfields}}{,}{\mathrm{Trace}}{,}{\mathrm{TransformCoordinates}}{,}{\mathrm{Vectors}}{,}{\mathrm{Weyl}}{,}{\mathrm{^}}{,}{\mathrm{dAlembertian}}{,}{\mathrm{d_}}{,}{\mathrm{diff}}{,}{\mathrm{g_}}{,}{\mathrm{gamma_}}\right]$ (34)

Set first $\mathrm{\theta }$ and $Q$ as suffixes for variables of type/anticommutative (see Setup)

 > $\mathrm{Setup}\left(\mathrm{anticommutativepre}=\left\{\mathrm{θ},Q\right\}\right)$
 $\mathrm{* Partial match of \text{'}}\mathrm{anticommutativepre}\mathrm{\text{'} against keyword \text{'}}\mathrm{anticommutativeprefix}\text{'}$
 $\mathrm{_______________________________________________________}$
 $\left[{\mathrm{anticommutativeprefix}}{=}\left\{{Q}{,}{\mathrm{\theta }}\right\}\right]$ (35)

A PDE system example with two unknown anticommutative functions of four variables, two commutative and two anticommutative; to avoid redundant typing in the input that follows and redundant display of information on the screen let's use PDEtools:-diff_table and PDEtools:-declare.

 > $\mathrm{PDEtools}:-\mathrm{declare}\left(Q\left(x,y,{\mathrm{θ}}_{1},{\mathrm{θ}}_{2}\right)\right)$
 ${Q}{}\left({x}{,}{y}{,}{{\mathrm{\theta }}}_{{1}}{,}{{\mathrm{\theta }}}_{{2}}\right){}$