Find a Non-negative Scaling Invariant Solution Having the Form

PDEtools

InvariantSolutions

computes the group invariant (symmetry) solutions for a given PDE system

Calling Sequence Parameters Description Examples Compatibility

Calling Sequence InvariantSolutions(PDESYS, DepVars, S, 'options'='value') Parameters PDESYS - a PDE or a set or list of them; it can include ODEs and non-differential equations DepVars - optional - may be required; a function or a list of them indicating the dependent variables of the problem S - optional - a list with the functional form of the infinitesimals of a symmetry generator, or a list of lists of them 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 before computing the determining system for its symmetries degreeofinfinitesimals = ... - optional - related to the option typeofsymmetry = polynomial , indicates the upper bound degree of the polynomial infinitesimals with respect to all or each of the dependent variables dependency = ... - optional - indicates the dependency of the invariant solutions, as a number (of jet variables), a range of them, a name, or a set of dependencies each of which can be a set or list of variables or numbers dependencyofinfinitesimals = ... - optional - same as dependency  but for the infinitesimals with which the invariant solutions are constructed discardinfinitesimalswithintegrals = ... - optional - can be true  or false  (default), to discard infinitesimals that contain uncomputed integrals display = ... - optional - can be true  or false  (default), to display or not the solutions on the screen as they are computed HINT = ... - optional - a list with the functional form of the infinitesimals of a symmetry generator invariants = ... - optional - a list of lists, where each inner list contains invariants  to be used to compute group invariant solutions numberofsolutions = ... - optional - can be a non-negative integer or the keyword all , to specify the maximum number of invariant solutions desired; default is 10 onlythetransformation = ... - optional - can be true  or false , to return or not just the different transformations used to reduce the number of independent variables of PDESYS reducedsystemsolver = ... - optional - can be pdsolve  or InvariantSolutions , specifies the command to be used to solve the system resulting from reducing the number of independent variables of PDESYS reduceinonego = ... - optional - a positive integer up to infinity , specifies how many independent variables are reduced in one go ; default is as much as possible removeredundant = ... - optional - can be true  or false  (default), to remove or not, redundant solutions which happen to be particular cases of more general solutions being returned together simplifier = ... - optional - indicates the simplifier to be used instead of the default simplify/size specialize_Fn = ... - optional - can be true , false  (default), or a set of procedures, to specialize any arbitrary functions that may appear in the infinitesimals specialize_Cn = ... - optional - can be true  (default) or false , to split or not the list of infinitesimals into cases by specializing the integration constants $\mathrm{\_Cn}$ typeofsymmetry = ... - optional - can be any of pointlike  (default), contact , evolutionary , or general , or any of polynomial  or functionfield , or And(<kind>, <functionality>)  where kind  is any of the first three types mentioned and functionality is any of polynomial  or functionfield , indicates the type of symmetry to be computed usestandardinvariants = ... - optional - can be true  or false  (default), to use the standard Invariants  instead of the CharacteristicQInvariants  for the purpose of computing the invariant solutions Description • Given a PDE problem ( PDESYS ), as an equation or a set or list of them, the InvariantSolutions  command computes the so-called group invariant  solutions of PDESYS ; these are solutions derived from n -dimensional symmetry groups admitted by PDESYS  by reducing the number of independent variables by n  in one go. From a practical point of view, this is an important difference with regards to the so-called similarity  solutions (also implemented in PDEtools  - see SimilaritySolutions ) where the number of independent variables are reduced one at a time and to arrive at a solution requires computing further symmetries for the reduced equation at each step. A "similarity solutions" approach can also be performed using InvariantSolutions  using the optional argument reduceinonego = 1 . InvariantSolutions  also works with anticommutative variables  set using the Physics  package using the approach explained in PerformOnAnticommutativeSystem . • The group-invariant solutions for PDESYS  returned in one step by InvariantSolutions , can also be computed step by step, by first computing the determining PDE  system, next solving this system to obtain the infinitesimals . With the infinitesimals, you compute an invariant transformation  to reduce the number of independent variables of PDESYS . Finally you solve the reduced system to change variables back, resulting in the solutions for PDESYS . You can compute any of these steps directly departing from PDESYS , respectively using the commands: DeterminingPDE , Infinitesimals , InvariantTransformation , dchange  and InvariantSolutions . • If DepVars  is not given, InvariantSolutions  will consider all the differentiated unknown functions in PDESYS  as unknowns of the problems. 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 . • It is possible to indicate the functional form of the symmetries, as shown in the examples, otherwise InvariantSolutions  will compute only the solutions that can be derived from the point  symmetries admitted by PDESYS . • 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. Optional arguments A number of optional arguments can change the default computational behavior of InvariantSolutions  in ways that, depending on the example or your needs, may be preferred. • 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 and everything else derived from it would be meaningless. • closesystem = true | false By default InvariantSolutions  processes PDESYS  as given. If, in PDESYS , all the integrability conditions that can be derived from its equations are explicitly present, then the determining system for the symmetries is assured  to have for solution all  the symmetries of PDESYS . Otherwise, the determining system for the symmetries 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. Correspondingly, some invariant solutions may not be present in the output of InvariantSolutions . 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 also the Library  routine CloseSystem . • degreeofinfinitesimals = nonnegint | set(function = nonnegint) This option is related to the option typeofsymmetry  being polynomial  or functionfield  explained below. By default, an upper bound for the degree of the polynomial dependency of the infinitesimals on the jet variables  is determined by the Library  routine UpperBounds . To override this behavior pass degreeofinfinitesimals = n , where $n$  is a non-negative integer, so that Infinitesimals , called by InvariantSolutions , will search for infinitesimals all of them of degree $n$ . You can also specify the degree of each infinitesimal, say, ${\mathrm{\_\&xi;}}_1$ , ${\mathrm{\_\&xi;}}_2$ , ${\mathrm{\_\&eta;}}_1$ , passing, say, degreeofinfinitesimals = {_xi[1] = n, _xi[2] = m, _eta[1] = p} , where $n$ , $m$ , $p$  are non-negative integers. • dependency = name | nonnegint | range(nonnegint) | set(name) | list(name), or a set or list of any of the previous types The right-hand-side is a positive integer, or a set or list of independent variables of the PDE problem represented by PDESYS , respectively indicating either the minimum number of independent variables on which the solutions are expected to depend, or the actual independent variables expected in this dependency. This option is relevant in connection with the fact that some of the group invariant solutions do not depend on all the independent variables of the problem.  If this option is not indicated, the solutions returned are computed regardless of their dependency, using as many symmetries as possible in one go. • dependencyofinfinitesimals = name | nonnegint | range(nonnegint) | set(name) | list(name), or a set or list of any of the previous types This option is used to specify the dependency of the infinitesimals, typically to restrict it so that their computation is easier. – When dependencyofinfinitesimals = name , all the infinitesimals are assumed to depend only on that name variable. – dependencyofinfinitesimals = n  where $n$  is a non-negative integer generates a search for infinitesimals depending on any of the permutations of $n$ jet variables ; this is useful in contexts where we need only infinitesimals depending on just some ( $n$ ) variables. – Specifying an integer range dependencyofinfinitesimals = n..m  generates a search for infinitesimals depending on any permutation of $k$  variables where n < k < m , so depending on not less than $n$  and no more than $m$  variables. – Specifying a set of lists, say the jet variables are $x$ , $t$ , and $u$ , so for instance as in dependencyofinfinitesimals = {[x], [x, t], [t, u]}  will generate a search for infinitesimals depending on any of the tree lists of variables specified. – Specifying a set or list of names generates a search for infinitesimals depending on any but only one of the indicated variables. • discardinfinitesimalswithintegrals = true | false The right-hand-side is either true  or false , to discard or not infinitesimals that contain uncomputed integrals before attempting the computation invariants with them and the related group invariant solutions. The default value of this option is false. Passing discardinfinitesimalswithintegrals = true  however can remove from the computational process invariants that due to the presence of integrals slow down the solving process significantly even when not help to find any more solutions. • display = true | false The right-hand-side is either true  or false , to display or not the solutions on the screen as they are computed  - that is before the final form of the output, typically involving many  group invariant solutions, is ready to be returned. The default value of this option is false. • HINT = <functional form of the infinitesimals> The functional form of the infinitesimals can be specified as a list of algebraic expressions, there must be as many as the number $n$  of independent variables plus the number $m$  of dependent variables to be correspondingly associated with ${\mathrm{\_\&xi;}}_1$  to ${\mathrm{\_\&xi;}}_n$  and with ${\mathrm{\_\&eta;}}_1$  to ${\mathrm{\_\&eta;}}_m$ . Alternatively this list can contain $n+m$  equations, where the left-hand-sides contain any functions of name variables and the right-hand-sides contains the actual hint being suggested. In both cases, the algebraic expressions (or those in the right-hand- sides) may or not contain functions to be determined so that the problem has solution. If there are no functions, then either the given HINT  exactly solves the determining system for the infinitesimals, and hence is used as given to construct solutions, or it doesn't in which case they are discarded. • invariants = [[algebraic, ...], [algebraic, ...], ...] The right-hand-side is a list of lists, where each inner list contains invariants  represented by algebraic expressions, to be used to compute the group invariant solutions. When this option is used, InvariantSolutions  will not attempt to compute any infinitesimals for PDESYS  and will instead use the given lists of invariants directly. Note there is a connection between the number of invariants (e.g., in one of these lists) and the number of independent variables that can be removed from the problem in one go. For a PDESYS  problem involving n  independent and m  dependent variables, a point symmetry of PDESYS  has n + m -  1 differential invariants of order zero . Recalling: invariants of order j  involve at most derivatives of order j  of the dependent variables. For each additional point symmetry enlarging the symmetry group, this number o differential invariants of order zero, now simultaneously invariant under the whole group, gets reduced  by one. So, when n + m - k  invariants of order zero are passed, InvariantSolutions  will attempt first a reduction of k  variables in one go, then k -  1 and so on until either achieving the numberofsolutions  or exhausting the use of the invariants received for the purpose of computing solutions. • numberofsolutions = nonnegint The right-hand-side is a positive integer indicating how many solutions InvariantSolutions  will attempt to compute before interrupting the process and returning the solutions found. This option is relevant in connection that, depending on the Infinitesimals  and related Invariants  found, the computation of all the related solutions may take more time than desired in a situation where only some solutions are required.  The default behavior is to compute the 10 group invariant solutions obtained by reducing the number of independent variables as much as possible in one go  (that is using many symmetries at the same time). • onlythetransformation = ... The right-hand-side is either true or false, to return or not just the different transformations used to reduce the number of independent variables of PDESYS . This option is used by the InvariantTransformation  command to compute these transformations. By combining this option with the other options you can search in an optimized way for only  the transformations you are interested in. This option is also of use when the transformation to reduce - say by k  - the number of independent variables in one go exist, but the system fails in solving the resulting reduced problem. In these cases it may be that by working with the transformation and the reduced PDEs system you may arrive to the solution to the original PDESYS  after some interactive steps. • reducedsystemsolver = ... The right-hand-side is either pdsolve  or InvariantSolutions , specifies the command to be used to attempt solving the reduced PDESYS. The default command used is a combination of pdsolve , dsolve  and solve  and whenever the solution obtained involves ODESolStruc  or DESol , this solution is discarded as unsuitable for representing PDE solutions. ( Note:  this may change in the future.) • reduceinonego = ... The right-hand-side is a positive integer, say k ; infinity  is allowed meaning as much as possible. This number k  specifies how many independent variables are to be reduced in one go . The default value is infinity . When this option is passed, only  solutions obtained by reducing k  variables in one go are returned. Note also that this option interacts with the numberofsolutions = ...  option in that InvariantSolutions  will not compute more solutions than the value of numberofsolutions . • removeredundant = ... The right-hand-side is either true or false, to remove or not redundant solutions that frequently appear in the set of group invariant solutions; they may be redundant in that they may be particular cases of more general solutions being returned together. • simplifier = ... The right-hand-side indicates the simplifier to be used instead of the default which is simplify/size . This option is of use not just to simplify  the solutions. You can specify any  procedure to be used as simplifier; in this way you can have the solutions computed InvariantSolutions  post-processed entirely in the way you need before being returned (e.g. discard some, do something else with other ones, etc.). • specialize_Cn = ... The right-hand-side is either true or false, to split or not the list of infinitesimals  into cases by specializing the integration constants $\mathrm{\_Cn}$  (see related discussion in Infinitesimals ) before computing invariants in the framework of the process to derive group-invariant solutions from these infinitesimals. The default of this option is true . This option is irrelevant when invariants are passed using the invariants = ...  option described above. • specialize_Fn = true | false | set(procedure) By default the arbitrary functions $\mathrm{Fn}\left(\mathrm{...}\right)$  entering the returned infinitesimals are not specialized. When specialize_Fn = true , these functions are specializing, taking all of them but one equal to 0, the remaining one equal to 1 and to each of its arguments. So if an arbitrary function has - say - three arguments, then there will be four specializations. Alternatively you can specify a set of procedures to be applied to the arguments of the arbitrary functions, each one resulting in a different desired specialization. • typeofsymmetry = pointlike | evolutionary | contact | general | polynomial | functionfield | And(symmetry_type, functionality_type) By default InvariantSolutions  will search for point symmetries to construct the corresponding invariant solutions. Alternatively you can request it to search for symmetries with polynomial  or functionfield  type of functionality, or of evolutionary , contact  or general  symmetry types instead of pointlike , 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. The right-hand-side can also be of the form And(symmetry_type, functionality_type)  where symmetry_type  is any of pointlike , evolutionary , contact  or general , and functionality_type  is one of polynomial  or functionfield , the latter related to the FunctionFieldSolutions  command. In these cases, You can set the degree of the polynomial or functionfield form using the option degreeofinfinitesimals  explained above. • usestandardinvariants = true | false By default, InvariantSolutions  will search for solutions by computing invariants from the characteristicQ function of the symmetries . The same computation can be carried out using the standard Invariants , normally resulting in a more complicated computation, though sometimes preferred for investigation or pedagogical reasons. For that purpose pass usestandardinvariants = true Examples InvariantSolutions  is programmed to automatically return invariant solutions (up to 10) reducing as many independent variables as possible in one go. So invoking it as in InvariantSolutions(PDE)  would suffice for that purpose. Most of the time, however, we are interested in either the most general invariant solutions, or those that can be computed faster, or that depend only on certain variables and so on. The following examples aim at illustrating the different ways of computing all these solutions. Consider the wave equation in four dimensions; to avoid redundant typing on input and on the display use diff_table  and declare > $\mathrm{with}\left(\mathrm{PDEtools}\right)\&colon;$ > $U≔\mathrm{diff\_table}\left(u\left(x\&comma;y\&comma;z\&comma;t\right)\right)\&colon;$ > $\mathrm{declare}\left(U\left[\right]\right)$ $u\left(x\&comma;y\&comma;z\&comma;t\right)\mathrm{will\; now\; be\; displayed\; as}u$ (1) > $\mathrm{pde}\left[1\right]≔U\left[x,x\right]+U\left[y,y\right]+U\left[z,z\right]-U\left[t,t\right]=0$ ${\mathrm{pde}}_1≔u_{x,x}+u_{y,y}+u_{z,z}-u_{t,t}=0$ (2) Compute just two solutions, both depending on the four variables { x, y, z, t }, reducing the number of independent variables only by three in one go  (so that the problem becomes an ODE) > $\mathrm{InvariantSolutions}\left(\mathrm{pde}\left[1\right]\&comma;\mathrm{numberofsolutions}=2\&comma;\mathrm{dependency}=4\&comma;\mathrm{reduceinonego}=3\right)$ $u=\mathrm{\_C1}+\frac{\mathrm{\_C2}}{\sqrt{\frac{t^2-y^2-z^2}{x^2}}},u=\frac{\mathrm{\_C1}\left(t^2-x^2-y^2-z^2\right)+\mathrm{\_C2}}{t^2-x^2-y^2-z^2}$ (3) Compute four solution depending at least on { x, t } (note these may depend on arbitrary functions $\mathrm{\_Fn}$  where n  is an integer); display each of them as soon as it is computed  (before computing the other ones) > $\mathrm{InvariantSolutions}\left(\mathrm{pde}\left[1\right]\&comma;\mathrm{numberofsolutions}=4\&comma;\mathrm{dependency}=\left\{t\&comma;x\right\}\&comma;\mathrm{display}\right)$ $\mathrm{invariants}=\left[t^2-x^2\&comma;u\right]$ $\mathrm{solutions}=\left[\left\{u=\mathrm{\_C2}\ln \left(\left(t-x\right)\left(t+x\right)\right)+\mathrm{\_C1}\right\}\right]$ $\mathrm{\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_}$ $\mathrm{invariants}=\left[\frac{t}{x}\&comma;u\right]$ $\mathrm{solutions}=\left[\left\{u=\mathrm{\_C3}+\frac{\mathrm{\_C4}\ln \left(\frac{t-x}{x}\right)}{2}-\frac{\mathrm{\_C4}\ln \left(\frac{t+x}{x}\right)}{2}\right\}\right]$ $\mathrm{\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_}$ $\mathrm{invariants}=\left[x\&comma;t\&comma;u\right]$ $\mathrm{solutions}=\left[\left\{u=\mathrm{\_F1}\left(t+x\right)+\mathrm{\_F2}\left(t-x\right)\right\}\right]$ $\mathrm{\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_}$ $\mathrm{invariants}=\left[\frac{t}{x}\&comma;ux\right]$ $\mathrm{solutions}=\left[\left\{u=\frac{\mathrm{\_C5}t+\mathrm{\_C6}x}{t^2-x^2}\right\}\right]$ $\mathrm{\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_}$ $u=\mathrm{\_F1}\left(t+x\right)+\mathrm{\_F2}\left(t-x\right),u=\mathrm{\_C2}\ln \left(\left(t-x\right)\left(t+x\right)\right)+\mathrm{\_C1},u=\frac{\mathrm{\_C1}t+\mathrm{\_C2}x}{t^2-x^2},u=\mathrm{\_C1}+\frac{\mathrm{\_C2}\ln \left(\frac{t-x}{x}\right)}{2}-\frac{\mathrm{\_C2}\ln \left(\frac{t+x}{x}\right)}{2}$ (4) This display  option is also useful in that you can select the invariants = [....] with the mouse, copy, and paste it as extra argument to InvariantSolutions , in order to compute only  the solutions it can derive from those invariants. In this sense the invariants you see on the screen associated to each solution are a sort of fingerprint of the solution. For example > $\mathrm{InvariantSolutions}\left(\mathrm{pde}\left[1\right]\&comma;\mathrm{invariants}=\left[-x^2+t^2\&comma;u\right]\right)$ $u=\mathrm{\_C2}\ln \left(\left(t-x\right)\left(t+x\right)\right)+\mathrm{\_C1}$ (5) If the purpose is to compute some solutions rapidly, regardless of their generality, you can use use the options typeofsymmetry = polynomial  or in addition specify that the degreeofinfinitesimals = 0  so that only constant infinitesimals are searched to construct invariant solutions; the same can be obtained specifying that the dependencyofinfinitesimals = 0  (they depend on 0 variables) or that they depend on no variables, for example > $\left[\mathrm{InvariantSolutions}\left(\mathrm{pde}\left[1\right]\&comma;\mathrm{dependencyofinfinitesimals}=\mathrm{none}\right)\right]$ $\left[u=\mathrm{\_C1}y+\mathrm{\_C2}\&comma;u=\mathrm{\_C1}z+\mathrm{\_C2}\&comma;u=\mathrm{\_C1}t+\mathrm{\_C2}\&comma;u=\mathrm{\_C1}x+\mathrm{\_C2}\right]$ (6) > $\left[\mathrm{InvariantSolutions}\left(\mathrm{pde}\left[1\right]\&comma;\mathrm{typeofsymmetry}=\mathrm{polynomial}\&comma;\mathrm{degreeofinfinitesimals}=0\right)\right]$ $\left[u=\mathrm{\_C1}y+\mathrm{\_C2}\&comma;u=\mathrm{\_C1}z+\mathrm{\_C2}\&comma;u=\mathrm{\_C1}t+\mathrm{\_C2}\&comma;u=\mathrm{\_C1}x+\mathrm{\_C2}\right]$ (7) This results in solutions computed fast and at the same time a bit more general than in the above; we can omit typeofsymmetry = polynomial  because it is already implicit when entering degreeofinfinitesimals = ... > $\left\{\mathrm{InvariantSolutions}\left(\mathrm{pde}\left[1\right]\&comma;\mathrm{degreeofinfinitesimals}=1\right)\right\}$ $\left\{u=\frac{\mathrm{\_C1}\left(t^2-x^2-y^2-z^2\right)+\mathrm{\_C2}}{t^2-x^2-y^2-z^2}\&comma;u=\frac{\mathrm{\_F1}\left(\sqrt{t^2-x^2-y^2}+z\right)+\mathrm{\_F2}\left(-\sqrt{t^2-x^2-y^2}+z\right)}{\sqrt{t^2-x^2-y^2}}\&comma;u=\frac{\mathrm{\_F1}\left(\sqrt{t^2-x^2-z^2}+y\right)+\mathrm{\_F2}\left(-\sqrt{t^2-x^2-z^2}+y\right)}{\sqrt{t^2-x^2-z^2}}\&comma;u=\frac{\mathrm{\_F1}\left(\sqrt{t^2-y^2-z^2}+x\right)+\mathrm{\_F2}\left(-\sqrt{t^2-y^2-z^2}+x\right)}{\sqrt{t^2-y^2-z^2}}\&comma;u=\mathrm{\_C1}+\frac{\mathrm{\_C2}}{\sqrt{\frac{t^2-x^2-z^2}{y^2}}}\&comma;u=\mathrm{\_C1}+\frac{\mathrm{\_C2}}{\sqrt{\frac{t^2-y^2-z^2}{x^2}}}\&comma;u=\mathrm{\_C2}+\frac{\mathrm{\_C1}}{\sqrt{\frac{t^2-x^2-y^2}{z^2}}}\&comma;u=\mathrm{\_C2}\ln \left(\frac{\left(t-x\right)\left(t+x\right)}{y^2+z^2}\right)+\mathrm{\_C1}\&comma;u=\mathrm{\_C2}\ln \left(\frac{\left(t-y\right)\left(t+y\right)}{x^2+z^2}\right)+\mathrm{\_C1}\&comma;u=\mathrm{\_C2}\ln \left(\frac{\left(t-z\right)\left(t+z\right)}{x^2+y^2}\right)+\mathrm{\_C1}\right\}$ (8) From where are the group invariant solutions being derived? They are derived from pde [1]'s symmetries, represented by its infinitesimals , that is, a list with the coefficients of the infinitesimal generator  of symmetry transformations  leaving invariant (in form) pde [1]. These infinitesimals can be computed via > $S\left[1\right]≔\mathrm{Infinitesimals}\left(\mathrm{pde}\left[1\right]\&comma;\mathrm{displayfunctionality}=\mathrm{false}\right)$ $S_1≔\left[{\mathrm{\_\&xi;}}_x=0\&comma;{\mathrm{\_\&xi;}}_y=1\&comma;{\mathrm{\_\&xi;}}_z=0\&comma;{\mathrm{\_\&xi;}}_t=0\&comma;{\mathrm{\_\&eta;}}_u=0\right],\left[{\mathrm{\_\&xi;}}_x=0\&comma;{\mathrm{\_\&xi;}}_y=0\&comma;{\mathrm{\_\&xi;}}_z=1\&comma;{\mathrm{\_\&xi;}}_t=0\&comma;{\mathrm{\_\&eta;}}_u=0\right],\left[{\mathrm{\_\&xi;}}_x=0\&comma;{\mathrm{\_\&xi;}}_y=0\&comma;{\mathrm{\_\&xi;}}_z=0\&comma;{\mathrm{\_\&xi;}}_t=1\&comma;{\mathrm{\_\&eta;}}_u=0\right],\left[{\mathrm{\_\&xi;}}_x=1\&comma;{\mathrm{\_\&xi;}}_y=0\&comma;{\mathrm{\_\&xi;}}_z=0\&comma;{\mathrm{\_\&xi;}}_t=0\&comma;{\mathrm{\_\&eta;}}_u=0\right],\left[{\mathrm{\_\&xi;}}_x=0\&comma;{\mathrm{\_\&xi;}}_y=t\&comma;{\mathrm{\_\&xi;}}_z=0\&comma;{\mathrm{\_\&xi;}}_t=y\&comma;{\mathrm{\_\&eta;}}_u=0\right],\left[{\mathrm{\_\&xi;}}_x=0\&comma;{\mathrm{\_\&xi;}}_y=0\&comma;{\mathrm{\_\&xi;}}_z=t\&comma;{\mathrm{\_\&xi;}}_t=z\&comma;{\mathrm{\_\&eta;}}_u=0\right],\left[{\mathrm{\_\&xi;}}_x=t\&comma;{\mathrm{\_\&xi;}}_y=0\&comma;{\mathrm{\_\&xi;}}_z=0\&comma;{\mathrm{\_\&xi;}}_t=x\&comma;{\mathrm{\_\&eta;}}_u=0\right],\left[{\mathrm{\_\&xi;}}_x=x\&comma;{\mathrm{\_\&xi;}}_y=y\&comma;{\mathrm{\_\&xi;}}_z=z\&comma;{\mathrm{\_\&xi;}}_t=t\&comma;{\mathrm{\_\&eta;}}_u=0\right],\left[{\mathrm{\_\&xi;}}_x=0\&comma;{\mathrm{\_\&xi;}}_y=0\&comma;{\mathrm{\_\&xi;}}_z=0\&comma;{\mathrm{\_\&xi;}}_t=0\&comma;{\mathrm{\_\&eta;}}_u=u\right],\left[{\mathrm{\_\&xi;}}_x=0\&comma;{\mathrm{\_\&xi;}}_y=z\&comma;{\mathrm{\_\&xi;}}_z=-y\&comma;{\mathrm{\_\&xi;}}_t=0\&comma;{\mathrm{\_\&eta;}}_u=0\right],\left[{\mathrm{\_\&xi;}}_x=z\&comma;{\mathrm{\_\&xi;}}_y=0\&comma;{\mathrm{\_\&xi;}}_z=-x\&comma;{\mathrm{\_\&xi;}}_t=0\&comma;{\mathrm{\_\&eta;}}_u=0\right],\left[{\mathrm{\_\&xi;}}_x=y\&comma;{\mathrm{\_\&xi;}}_y=-x\&comma;{\mathrm{\_\&xi;}}_z=0\&comma;{\mathrm{\_\&xi;}}_t=0\&comma;{\mathrm{\_\&eta;}}_u=0\right],\left[{\mathrm{\_\&xi;}}_x=xz\&comma;{\mathrm{\_\&xi;}}_y=zy\&comma;{\mathrm{\_\&xi;}}_z=\frac{z^2}{2}+\frac{t^2}{2}-\frac{x^2}{2}-\frac{y^2}{2}\&comma;{\mathrm{\_\&xi;}}_t=zt\&comma;{\mathrm{\_\&eta;}}_u=-uz\right],\left[{\mathrm{\_\&xi;}}_x=xy\&comma;{\mathrm{\_\&xi;}}_y=\frac{y^2}{2}+\frac{t^2}{2}-\frac{x^2}{2}-\frac{z^2}{2}\&comma;{\mathrm{\_\&xi;}}_z=zy\&comma;{\mathrm{\_\&xi;}}_t=yt\&comma;{\mathrm{\_\&eta;}}_u=-uy\right],\left[{\mathrm{\_\&xi;}}_x=xt\&comma;{\mathrm{\_\&xi;}}_y=yt\&comma;{\mathrm{\_\&xi;}}_z=zt\&comma;{\mathrm{\_\&xi;}}_t=\frac{t^2}{2}+\frac{x^2}{2}+\frac{y^2}{2}+\frac{z^2}{2}\&comma;{\mathrm{\_\&eta;}}_u=-ut\right],\left[{\mathrm{\_\&xi;}}_x=-\frac{x^2}{2}-\frac{t^2}{2}+\frac{y^2}{2}+\frac{z^2}{2}\&comma;{\mathrm{\_\&xi;}}_y=-xy\&comma;{\mathrm{\_\&xi;}}_z=-xz\&comma;{\mathrm{\_\&xi;}}_t=-xt\&comma;{\mathrm{\_\&eta;}}_u=ux\right]$ (9) You can test these results for correctness using the SymmetryTest  command > $\mathrm{map}\left(\mathrm{SymmetryTest}\&comma;\left[S\left[1\right]\right]\&comma;\mathrm{pde}\left[1\right]\right)$ $\left[\left\{0\right\}\&comma;\left\{0\right\}\&comma;\left\{0\right\}\&comma;\left\{0\right\}\&comma;\left\{0\right\}\&comma;\left\{0\right\}\&comma;\left\{0\right\}\&comma;\left\{0\right\}\&comma;\left\{0\right\}\&comma;\left\{0\right\}\&comma;\left\{0\right\}\&comma;\left\{0\right\}\&comma;\left\{0\right\}\&comma;\left\{0\right\}\&comma;\left\{0\right\}\&comma;\left\{0\right\}\right]$ (10) Sometimes, instead of these infinitesimals it is preferable to see the corresponding infinitesimal generator ; here is how you get it for example for the infinitesimals of the fourth list of infinitesimals above > $\mathrm{InfinitesimalGenerator}\left(S\left[1\right]\left[9\right]\&comma;U\left[\right]\right)$ $f\&rarr;u\left(\frac{\partial }{\partial u}f\right)$ (11) How are these infinitesimals  used? With the infinitesimals you can compute invariants  under the whole symmetry group (one or more symmetries); with them, transformations reducing the number of independent variables can systematically be computed (provided the computer has enough memory and time for that purpose - this is frequently the case). For example, these are the invariants of order zero  associated to the fifth list of infinitesimals in S[1], and this is how you can pass them to InvariantSolutions  to be used at the time of constructing the solution > $\mathrm{invs}≔\left[\mathrm{Invariants}\left(S\left[1\right]\left[5\right]\&comma;U\left[\right]\&comma;0\right)\right]$ $\mathrm{invs}≔\left[x\&comma;z\&comma;u\&comma;t^2-y^2\right]$ (12) > $\mathrm{InvariantSolutions}\left(\mathrm{pde}\left[1\right]\&comma;U\left[\right]\&comma;\mathrm{invariants}=\mathrm{invs}\right)$ $u=\left(\mathrm{\_C1}\mathrm{\_C3}{\&ExponentialE;}^{2\sqrt{{\mathrm{\_c}}_1}x+2\sqrt{{\mathrm{\_c}}_2}z}+{\&ExponentialE;}^{2\sqrt{{\mathrm{\_c}}_1}x}\mathrm{\_C1}\mathrm{\_C4}+\mathrm{\_C2}\left({\&ExponentialE;}^{2\sqrt{{\mathrm{\_c}}_2}z}\mathrm{\_C3}+\mathrm{\_C4}\right)\right)\left(\mathrm{\_C6}\mathrm{BesselY}\left(0\&comma;\sqrt{-{\mathrm{\_c}}_1-{\mathrm{\_c}}_2}\sqrt{t^2-y^2}\right)+\mathrm{\_C5}\mathrm{BesselJ}\left(0\&comma;\sqrt{-{\mathrm{\_c}}_1-{\mathrm{\_c}}_2}\sqrt{t^2-y^2}\right)\right){\&ExponentialE;}^{-\sqrt{{\mathrm{\_c}}_1}x-\sqrt{{\mathrm{\_c}}_2}z}$ (13) The symmetry used when computing the invariant solutions can also be specified > $\mathrm{InvariantSolutions}\left(\mathrm{pde}\left[1\right]\&comma;U\left[\right]\&comma;S\left[1\right]\left[6\right]\right)$ $u=\left(\mathrm{\_C1}\mathrm{\_C3}{\&ExponentialE;}^{2\sqrt{{\mathrm{\_c}}_1}x+2\sqrt{{\mathrm{\_c}}_2}y}+{\&ExponentialE;}^{2\sqrt{{\mathrm{\_c}}_1}x}\mathrm{\_C1}\mathrm{\_C4}+\mathrm{\_C2}\left({\&ExponentialE;}^{2\sqrt{{\mathrm{\_c}}_2}y}\mathrm{\_C3}+\mathrm{\_C4}\right)\right)\left(\mathrm{\_C6}\mathrm{BesselY}\left(0\&comma;\sqrt{-{\mathrm{\_c}}_1-{\mathrm{\_c}}_2}\sqrt{t^2-z^2}\right)+\mathrm{\_C5}\mathrm{BesselJ}\left(0\&comma;\sqrt{-{\mathrm{\_c}}_1-{\mathrm{\_c}}_2}\sqrt{t^2-z^2}\right)\right){\&ExponentialE;}^{-\sqrt{{\mathrm{\_c}}_1}x-\sqrt{{\mathrm{\_c}}_2}y}$ (14) By default, InvariantSolutions  will specialize arbitrary functions $\mathrm{\_Fn}\left(\mathrm{...}\right)$  found in the infinitesimals  before attempting a reduction of the number of independent variables. This is so because infinitesimals with arbitrary functions are difficult to use, the invariants  necessary may result, in practice, impossible to compute. On the other hand, when those invariants can be computed, the corresponding invariant solution is more general without specializing these arbitrary functions. To handle the latter case use the specialize_Fn = false , or with any procedure of your choice on the right-hand-side. Here are two examples illustrating the two situations. Consider the wave equation in two dimensions: > $U≔\mathrm{diff\_table}\left(u\left(x\&comma;t\right)\right)\&colon;$ > $\mathrm{PDE}≔U\left[x,x\right]-U\left[t,t\right]=0$ $\mathrm{PDE}≔u_{x,x}-u_{t,t}=0$ (15) Infinitesimals  for this PDE  can be computed and involve arbitrary functions > $S≔\mathrm{Infinitesimals}\left(\mathrm{PDE}\right)$ $S≔\left[{\mathrm{\_\&xi;}}_x\left(x\&comma;t\&comma;u\right)=\mathrm{\_F1}\left(t+x\right)+\mathrm{\_F2}\left(t-x\right)\&comma;{\mathrm{\_\&xi;}}_t\left(x\&comma;t\&comma;u\right)=-\mathrm{\_F2}\left(t-x\right)+\mathrm{\_F1}\left(t+x\right)+1\&comma;{\mathrm{\_\&eta;}}_u\left(x\&comma;t\&comma;u\right)=\mathrm{\_F3}\left(t+x\right)+\mathrm{\_F4}\left(t-x\right)\right],\left[{\mathrm{\_\&xi;}}_x\left(x\&comma;t\&comma;u\right)=\mathrm{\_F5}\left(t+x\right)+\mathrm{\_F6}\left(t-x\right)\&comma;{\mathrm{\_\&xi;}}_t\left(x\&comma;t\&comma;u\right)=-\mathrm{\_F6}\left(t-x\right)+\mathrm{\_F5}\left(t+x\right)\&comma;{\mathrm{\_\&eta;}}_u\left(x\&comma;t\&comma;u\right)=u+\mathrm{\_F7}\left(t+x\right)+\mathrm{\_F8}\left(t-x\right)\right]$ (16) It is not possible to compute a SimilarityTransformation  or for the case any symmetry related transformation with them, because of the arbitrary functions, all these commands return NULL: > $\mathrm{SimilarityTransformation}\left(S\left[1\right]\&comma;U\left[\right]\right)$ > $\mathrm{InvariantTransformation}\left(S\left[1\right]\&comma;U\left[\right]\right)$ > $\mathrm{SymmetryTransformation}\left(S\left[1\right]\&comma;U\left[\right]\right)$ In order to be able to use the infinitesimals S , these arbitrary functions are automatically specialized by all the symmetry commands (but for Infinitesimals ), as you could see by interactively using specialize_Fn = true  in the call to Infinitesimals > $S≔\left[\mathrm{Infinitesimals}\left(\mathrm{PDE}\&comma;\mathrm{specialize\_Fn}\right)\right]$ $S≔\left[\left[{\mathrm{\_\&xi;}}_x\left(x\&comma;t\&comma;u\right)=1\&comma;{\mathrm{\_\&xi;}}_t\left(x\&comma;t\&comma;u\right)=0\&comma;{\mathrm{\_\&eta;}}_u\left(x\&comma;t\&comma;u\right)=0\right]\&comma;\left[{\mathrm{\_\&xi;}}_x\left(x\&comma;t\&comma;u\right)=1\&comma;{\mathrm{\_\&xi;}}_t\left(x\&comma;t\&comma;u\right)=2\&comma;{\mathrm{\_\&eta;}}_u\left(x\&comma;t\&comma;u\right)=0\right]\&comma;\left[{\mathrm{\_\&xi;}}_x\left(x\&comma;t\&comma;u\right)=0\&comma;{\mathrm{\_\&xi;}}_t\left(x\&comma;t\&comma;u\right)=1\&comma;{\mathrm{\_\&eta;}}_u\left(x\&comma;t\&comma;u\right)=1\right]\&comma;\left[{\mathrm{\_\&xi;}}_x\left(x\&comma;t\&comma;u\right)=1\&comma;{\mathrm{\_\&xi;}}_t\left(x\&comma;t\&comma;u\right)=1\&comma;{\mathrm{\_\&eta;}}_u\left(x\&comma;t\&comma;u\right)=u\right]\&comma;\left[{\mathrm{\_\&xi;}}_x\left(x\&comma;t\&comma;u\right)=1\&comma;{\mathrm{\_\&xi;}}_t\left(x\&comma;t\&comma;u\right)=\mathrm{-1}\&comma;{\mathrm{\_\&eta;}}_u\left(x\&comma;t\&comma;u\right)=u\right]\&comma;\left[{\mathrm{\_\&xi;}}_x\left(x\&comma;t\&comma;u\right)=0\&comma;{\mathrm{\_\&xi;}}_t\left(x\&comma;t\&comma;u\right)=1\&comma;{\mathrm{\_\&eta;}}_u\left(x\&comma;t\&comma;u\right)=t+x\right]\&comma;\left[{\mathrm{\_\&xi;}}_x\left(x\&comma;t\&comma;u\right)=0\&comma;{\mathrm{\_\&xi;}}_t\left(x\&comma;t\&comma;u\right)=1\&comma;{\mathrm{\_\&eta;}}_u\left(x\&comma;t\&comma;u\right)=t-x\right]\&comma;\left[{\mathrm{\_\&xi;}}_x\left(x\&comma;t\&comma;u\right)=0\&comma;{\mathrm{\_\&xi;}}_t\left(x\&comma;t\&comma;u\right)=0\&comma;{\mathrm{\_\&eta;}}_u\left(x\&comma;t\&comma;u\right)=u+1\right]\&comma;\left[{\mathrm{\_\&xi;}}_x\left(x\&comma;t\&comma;u\right)=0\&comma;{\mathrm{\_\&xi;}}_t\left(x\&comma;t\&comma;u\right)=0\&comma;{\mathrm{\_\&eta;}}_u\left(x\&comma;t\&comma;u\right)=u+t+x\right]\&comma;\left[{\mathrm{\_\&xi;}}_x\left(x\&comma;t\&comma;u\right)=0\&comma;{\mathrm{\_\&xi;}}_t\left(x\&comma;t\&comma;u\right)=0\&comma;{\mathrm{\_\&eta;}}_u\left(x\&comma;t\&comma;u\right)=u+t-x\right]\&comma;\left[{\mathrm{\_\&xi;}}_x\left(x\&comma;t\&comma;u\right)=t+x\&comma;{\mathrm{\_\&xi;}}_t\left(x\&comma;t\&comma;u\right)=1+t+x\&comma;{\mathrm{\_\&eta;}}_u\left(x\&comma;t\&comma;u\right)=0\right]\&comma;\left[{\mathrm{\_\&xi;}}_x\left(x\&comma;t\&comma;u\right)=t-x\&comma;{\mathrm{\_\&xi;}}_t\left(x\&comma;t\&comma;u\right)=-t+x+1\&comma;{\mathrm{\_\&eta;}}_u\left(x\&comma;t\&comma;u\right)=0\right]\&comma;\left[{\mathrm{\_\&xi;}}_x\left(x\&comma;t\&comma;u\right)=t+x\&comma;{\mathrm{\_\&xi;}}_t\left(x\&comma;t\&comma;u\right)=t+x\&comma;{\mathrm{\_\&eta;}}_u\left(x\&comma;t\&comma;u\right)=u\right]\&comma;\left[{\mathrm{\_\&xi;}}_x\left(x\&comma;t\&comma;u\right)=t-x\&comma;{\mathrm{\_\&xi;}}_t\left(x\&comma;t\&comma;u\right)=-t+x\&comma;{\mathrm{\_\&eta;}}_u\left(x\&comma;t\&comma;u\right)=u\right]\right]$ (17) In this way, all these commands can deal with PDE  properly > $\mathrm{SymmetryTransformation}\left(S\left[1\right]\&comma;U\left[\right]\right)$ $\left\{\mathrm{\_t1}=x+\mathrm{\_\&epsilon;}\&comma;\mathrm{\_t2}=t\&comma;\mathrm{\_u1}\left(\mathrm{\_t1}\&comma;\mathrm{\_t2}\right)=u\left(x\&comma;t\right)\right\}$ (18) > $\mathrm{InvariantTransformation}\left(S\left[11\right]\&comma;U\left[\right]\right)$ $\left\{\mathrm{\_t1}={\&ExponentialE;}^{-2t+2x}\left(1+2t+2x\right)\&comma;\mathrm{\_t2}=t\&comma;\mathrm{\_u1}\left(\mathrm{\_t1}\right)=u\left(x\&comma;t\right)\right\},\left\{t=\mathrm{\_t2}\&comma;x=\frac{\mathrm{LambertW}\left(\mathrm{\_t1}{\&ExponentialE;}^{4\mathrm{\_t2}+1}\right)}{2}-\mathrm{\_t2}-\frac{1}{2}\&comma;u\left(x\&comma;t\right)=\mathrm{\_u1}\left(\mathrm{\_t1}\right)\right\}$ (19) > $\left[\mathrm{InvariantSolutions}\left(\mathrm{PDE}\right)\right]$ $\left[u\left(x\&comma;t\right)=\mathrm{\_C1}\sqrt{t+x}\&comma;u\left(x\&comma;t\right)=\frac{\mathrm{\_C1}}{\sqrt{t-x}}\&comma;u\left(x\&comma;t\right)=\mathrm{\_C1}x+\mathrm{\_C2}+t\&comma;u\left(x\&comma;t\right)=\mathrm{\_C1}\left(-2x+t\right)+\mathrm{\_C2}\&comma;u\left(x\&comma;t\right)=\frac{\mathrm{\_C1}{\&ExponentialE;}^{\frac{t}{2}-\frac{x}{2}}}{{\&ExponentialE;}^{-x}}\&comma;u\left(x\&comma;t\right)=\frac{\mathrm{\_C1}{\&ExponentialE;}^{-\frac{t}{2}-\frac{x}{2}}}{{\&ExponentialE;}^{-x}}\&comma;u\left(x\&comma;t\right)=\frac{x^2}{2}+\frac{\left(2t+2\mathrm{\_C1}\right)x}{2}+\frac{t^2}{2}+\mathrm{\_C2}\&comma;u\left(x\&comma;t\right)=\frac{x^2}{2}+\frac{\left(-2t+2\mathrm{\_C2}\right)x}{2}+\frac{t^2}{2}+\mathrm{\_C1}\&comma;u\left(x\&comma;t\right)=\mathrm{\_C2}\ln \left({\&ExponentialE;}^{-2t+2x}\left(1+2t+2x\right)\right)+\mathrm{\_C1}\&comma;u\left(x\&comma;t\right)=\mathrm{\_C2}\ln \left({\&ExponentialE;}^{2t+2x}\left(2t-2x-1\right)\right)+\mathrm{\_C1}\right]$ (20) As an example where this approach leads to less general solutions consider this other PDE2 , also admitting infinitesimals that depend on arbitrary functions > $\mathrm{PDE2}≔\mathrm{diff}\left(\mathrm{diff}\left(\mathrm{diff}\left(u\left(x\&comma;t\right)\&comma;t\right)\&comma;x\right)\&comma;x\right)=-\mathrm{diff}\left(\mathrm{diff}\left(\mathrm{diff}\left(u\left(x\&comma;t\right)\&comma;t\right)\&comma;t\right)\&comma;x\right)$ $\mathrm{PDE2}≔u_{t,x,x}=-u_{t,t,x}$ (21) > $\mathrm{S2}≔\mathrm{Infinitesimals}\left(\mathrm{PDE2}\right)$ $\mathrm{S2}≔\left[{\mathrm{\_\&xi;}}_x\left(x\&comma;t\&comma;u\right)=0\&comma;{\mathrm{\_\&xi;}}_t\left(x\&comma;t\&comma;u\right)=1\&comma;{\mathrm{\_\&eta;}}_u\left(x\&comma;t\&comma;u\right)=\mathrm{\_F1}\left(t\right)+\mathrm{\_F2}\left(x\right)+\mathrm{\_F3}\left(t-x\right)\right],\left[{\mathrm{\_\&xi;}}_x\left(x\&comma;t\&comma;u\right)=1\&comma;{\mathrm{\_\&xi;}}_t\left(x\&comma;t\&comma;u\right)=0\&comma;{\mathrm{\_\&eta;}}_u\left(x\&comma;t\&comma;u\right)=\mathrm{\_F4}\left(t\right)+\mathrm{\_F5}\left(x\right)+\mathrm{\_F6}\left(t-x\right)\right],\left[{\mathrm{\_\&xi;}}_x\left(x\&comma;t\&comma;u\right)=x\&comma;{\mathrm{\_\&xi;}}_t\left(x\&comma;t\&comma;u\right)=t\&comma;{\mathrm{\_\&eta;}}_u\left(x\&comma;t\&comma;u\right)=\mathrm{\_F7}\left(t\right)+\mathrm{\_F8}\left(x\right)+\mathrm{\_F9}\left(t-x\right)\right],\left[{\mathrm{\_\&xi;}}_x\left(x\&comma;t\&comma;u\right)=0\&comma;{\mathrm{\_\&xi;}}_t\left(x\&comma;t\&comma;u\right)=0\&comma;{\mathrm{\_\&eta;}}_u\left(x\&comma;t\&comma;u\right)=u+\mathrm{\_F10}\left(t\right)+\mathrm{\_F11}\left(x\right)+\mathrm{\_F12}\left(t-x\right)\right]$ (22) For these other infinitesimals, despite the arbitrary functions it is possible to compute all kinds of symmetry transformations, although the results involve uncomputed integrals > $\mathrm{SymmetryTransformation}\left(\mathrm{S2}\left[1\right]\&comma;U\left[\right]\right)$ $\left\{\mathrm{\_t1}=x\&comma;\mathrm{\_t2}=t+\mathrm{\_\&epsilon;}\&comma;\mathrm{\_u1}\left(\mathrm{\_t1}\&comma;\mathrm{\_t2}\right)=\mathrm{\_F2}\left(x\right)\mathrm{\_\&epsilon;}+\int \left(-\mathrm{\_F1}\left(t\right)-\mathrm{\_F3}\left(t-x\right)\right)\&DifferentialD;t-\left({\int }_{}^{t+\mathrm{\_\&epsilon;}}\left(-\mathrm{\_F1}\left(\mathrm{\_a}\right)-\mathrm{\_F3}\left(\mathrm{\_a}-x\right)\right)\&DifferentialD;\mathrm{\_a}\right)+u\left(x\&comma;t\right)\right\}$ (23) > $\mathrm{SimilarityTransformation}\left(\mathrm{S2}\left[2\right]\&comma;U\left[\right]\right)$ $\left\{\mathrm{\_\&psi;}=x\&comma;{\mathrm{\_\&phi;}}_1=t\&comma;{\mathrm{\_\&phi;}}_2\left({\mathrm{\_\&phi;}}_1\right)=-\mathrm{\_F4}\left(t\right)x-\left(\int \left(\mathrm{\_F5}\left(x\right)+\mathrm{\_F6}\left(t-x\right)\right)\&DifferentialD;x\right)+u\left(x\&comma;t\right)\right\},\left\{t={\mathrm{\_\&phi;}}_1\&comma;x=\mathrm{\_\&psi;}\&comma;u\left(x\&comma;t\right)={\mathrm{\_\&phi;}}_2\left({\mathrm{\_\&phi;}}_1\right)+\mathrm{\_F4}\left({\mathrm{\_\&phi;}}_1\right)\mathrm{\_\&psi;}+{\int }_{}^{\mathrm{\_\&psi;}}\left(\mathrm{\_F5}\left(\mathrm{\_a}\right)+\mathrm{\_F6}\left({\mathrm{\_\&phi;}}_1-\mathrm{\_a}\right)\right)\&DifferentialD;\mathrm{\_a}\right\}$ (24) Consequently, the default solution returned by InvariantSolutions  is less general than the one you obtain enforcing the use of the infinitesimals with specialize_Fn = false > $\mathrm{InvariantSolutions}\left(\mathrm{PDE2}\right)$ $u\left(x\&comma;t\right)=\left(\mathrm{\_C2}+\mathrm{\_C3}\right)\ln \left(\frac{t-x}{x}\right)-\mathrm{\_C3}\ln \left(\frac{t}{x}\right)+t+\mathrm{\_C1},u\left(x\&comma;t\right)=\left(\mathrm{\_C2}+\mathrm{\_C3}\right)\ln \left(\frac{t-x}{x}\right)-\mathrm{\_C3}\ln \left(\frac{t}{x}\right)+x+\mathrm{\_C1},u\left(x\&comma;t\right)=\left(\mathrm{\_C2}+\mathrm{\_C3}\right)\ln \left(\frac{t-x}{x}\right)-\mathrm{\_C3}\ln \left(\frac{t}{x}\right)+\mathrm{\_C1}+\ln \left(x\right),u\left(x\&comma;t\right)=\left(\mathrm{\_C2}+\mathrm{\_C3}\right)\ln \left(\frac{t-x}{x}\right)-\mathrm{\_C3}\ln \left(\frac{t}{x}\right)+t-x+\mathrm{\_C1}$ (25) > $\mathrm{InvariantSolutions}\left(\mathrm{PDE2}\&comma;\mathrm{specialize\_Fn}=\mathrm{false}\right)$ $u\left(x\&comma;t\right)={\int }_{}^x\frac{\mathrm{\_F1}\left(\frac{t\mathrm{\_a}}{x}\right)+\mathrm{\_F2}\left(\mathrm{\_a}\right)+\mathrm{\_F3}\left(\frac{\mathrm{\_a}\left(t-x\right)}{x}\right)}{\mathrm{\_a}}\&DifferentialD;\mathrm{\_a}+\left(\mathrm{\_C2}+\mathrm{\_C3}\right)\ln \left(\frac{t-x}{x}\right)-\mathrm{\_C3}\ln \left(\frac{t}{x}\right)+\mathrm{\_C1}$ (26) Compatibility • The typeofsymmetry  option was updated in Maple 15. See Also CharacteristicQ CharacteristicQInvariants declare DeterminingPDE diff_table InvariantTransformation PDEtools PerformOnAnticommutativeSystem Physics PolynomialSolutions Setup SimilaritySolutions SimilarityTransformation SymmetrySolutions SymmetryTest SymmetryTransformation

Download Help Document

Thank you for submitting feedback on this help document. Your feedback will be used to improve Maple's help in the future.

Find a Non-negative Scaling Invariant Solution Having the Form

Source: https://www.maplesoft.com/support/help/maple/view.aspx?path=PDEtools%2FInvariantSolutions

Share :