The dchange command realizes changes of variables in mathematical expressions (PDEs, multiple integrals, integro-differential equations, limits, and so on) and in procedures. Only the global variables of a procedure can be changed using the dchange command, not the local ones.

This command was meant as a unified command for changing variables in arbitrary mathematical expressions, according to a given transformation of variables. In addition, it can use any specific indication for the inverse transformation, in which case a different algorithm for changing variables in derivatives is used. This is particularly useful when the transformation involves special functions not known to the system, or the "target" involves derivatives and the transformation is not one-to-one (see the examples).

These are the conventions for the arguments:

a) The first and second arguments have a fixed position. The first argument is the transformation, in the form old_var = f(new_vars) or as a set of equations with this form, and the second is the target, in the form of an algebraic expression, an equation, or a set or list of them.  The other (all optional) arguments may appear in any order after the second one.

b) The old and new variables involved in the transformation equations must be of type name or of type function. This second version of the dchange command also accepts transformation equations for the dependent variables (functions), allowing for a direct study of the invariance groups of a given PDE. The extension of the transformation to any required differential order is automatic.

c) The sets of transformation equations must contain only old (new) variables on the left hand side and only new (old) variables on the right hand side, as in

{x = xi+eta, y = xi-eta}; # the transformation

d) An indication of the new variables is required when it is not obvious, and it must then be given as a list. For example, in the transformation above, there are only two variables on the right hand sides of a set of two transformation equations; dchange assumes that the new variables are those appearing on the right hand side. If, however, these numbers are not equal, one must indicate which are the new variables:

[xi, eta] (or) [eta, xi]

Since in Maple $f\left(\mathrm{\xi }\,\mathrm{\eta }\right)\ne f\left(\mathrm{\eta }\,\mathrm{\xi }\right)$, the order given to the new variables in the list above may be relevant.

e) If only one set of transformation equations is given, then it must contain the new variables on the right hand side.

The dchange command has two different criteria for changing variables. The first one, applied by default to "known" functions (functions with a differentiation rule known to Maple), changes the functional dependence using the explicit relation between the old and new variables. For example, if the transformation is given by

{x = xi+eta, y = xi-eta},

then dchange will transform known functions as follows:

cos(x*y) -> cos( (xi+eta)*(xi-eta) )

The second criterion, applied by default to "unknown" functions (functions whose differentiation rule is unknown to Maple), changes the functional dependence taking into account only the implicit relation between the old and new variables. For example, with the transformation above, the dchange command will transform any "unknown" function (say f) as in:

f(x, y) -> f(xi, eta)

It is possible to change these defaults and force dchange to apply a different criterion to a given function, by specifying it using the 'known' or 'unknown' options.

When a function depends on more variables than those which are being changed, the variables not entering the process are all put to the right. For example, the transformation above would take f(x, y, z, t) -> f(xi, eta, z, t).

In addition, you can indicate the use of any standard or special procedure for the automatic simplification of the resulting expression. This is useful when nonobvious mathematical equivalencies must be taken into account, or just to make the resulting expression more compact. The simplification is applied not only to the final result but also to all inner intermediate results.

Special rules are applied to multiple integrals, since for definite ones, the mapping of the old integration domain into the new one may present additional complications. The command could fail to determine the (perhaps functional) range for the new variables defining the new integration domain. In these cases, a special warning message is displayed and only the functional dependence of the new ranges is determined.

Optional arguments can be given alone or in conjunction, and in any order.

This function is part of the PDEtools package, and so it can be used in the form dchange(..) only after executing the command with(PDEtools). However, it can always be accessed through the long form of the command by using PDEtools[dchange](..).

$\mathrm{with}\left(\mathrm{PDEtools}\right)\:$

$\mathrm{PDE}≔\mathrm{diff}\left(f\left(x\,y\right)\,x\right)+g\left(x\,y\right)+\mathrm{diff}\left(f\left(x\,y\right)\,y\right)=0$

$\mathrm{PDE}≔\frac{\partial }{\partial x}f\left(x\,y\right)+g\left(x\,y\right)+\frac{\partial }{\partial y}f\left(x\,y\right)=0$

$\mathrm{tr}≔\left\{x=r+s\,y=r-s\right\}\:$

$\mathrm{dchange}\left(\mathrm{tr}\,\mathrm{PDE}\right)$

$\frac{\partial }{\partial r}f\left(r\,s\right)+g\left(r\,s\right)=0$

$\mathrm{dchange}\left(\mathrm{tr}\,\mathrm{PDE}\,\mathrm{known}=g\right)$

$\frac{\partial }{\partial r}f\left(r\,s\right)+g\left(r+s\,r-s\right)=0$

$L≔f\→x\left(\frac{\partial }{\partial y}f\right)-y\left(\frac{\partial }{\partial x}f\right)$

$L≔f\→x\left(\frac{\partial }{\partial y}f\right)-y\left(\frac{\partial }{\partial x}f\right)$

$\mathrm{tr}≔\left\{x=s{r}^{\frac{1}{2}}\,y=s{\left(1-r\right)}^{\frac{1}{2}}\right\}$

$\mathrm{tr}≔\left\{x=s\sqrt{r}\,y=s\sqrt{1-r}\right\}$

$\mathrm{`L\text{'}`}=\mathrm{dchange}\left(\mathrm{tr}\,L\,\mathrm{simplify}\right)$

$\mathrm{L\text{'}}=\left(f\→-2\sqrt{r}\sqrt{1-r}\left(\frac{\partial }{\partial r}f\right)\right)$

$d\left[y\right]≔f\→\frac{\partial }{\partial y}f$

$d_y≔f\→\frac{\partial }{\partial y}f$

$\mathrm{dchange}\left(\mathrm{tr}\,d\left[y\right]\,\mathrm{normal}\right)$

$f\→-\frac{\sqrt{1-r}\left(2\left(\frac{\partial }{\partial r}f\right)r-\left(\frac{\partial }{\partial s}f\right)s\right)}{s}$

$\mathrm{macro}\left(\mathrm{ys}=\mathrm{`y*`}\,\mathrm{xs}=\mathrm{`x*`}\right)$

$\mathrm{ys},\mathrm{xs}$

$\mathrm{ODE}≔\frac{-y\left(x\right)+x\mathrm{diff}\left(y\left(x\right)\,x\right)}{x+y\left(x\right)\mathrm{diff}\left(y\left(x\right)\,x\right)}=H\left(\mathrm{sqrt}\left(x^2+{y\left(x\right)}^2\right)\right)$

$\mathrm{ODE}≔\frac{-y\left(x\right)+x\left(\frac{\ⅆ}{\ⅆx}y\left(x\right)\right)}{x+y\left(x\right)\left(\frac{\ⅆ}{\ⅆx}y\left(x\right)\right)}=H\left(\sqrt{x^2+{y\left(x\right)}^2}\right)$

$\mathrm{tr}≔\left\{x=\mathrm{ys}\left(\mathrm{xs}\right)\sin \left(\mathrm{\alpha }\right)+\mathrm{xs}\cos \left(\mathrm{\alpha }\right)\,y\left(x\right)=\cos \left(\mathrm{\alpha }\right)\mathrm{ys}\left(\mathrm{xs}\right)-\sin \left(\mathrm{\alpha }\right)\mathrm{xs}\right\}\:$

$\mathrm{itr}≔\mathrm{simplify}\left(\mathrm{solve}\left(\mathrm{tr}\,\left\{\mathrm{xs}\,\mathrm{ys}\left(\mathrm{xs}\right)\right\}\right)\,\mathrm{trig}\right)$

$\mathrm{itr}≔\left\{\mathrm{x*}=-y\left(x\right)\sin \left(\mathrm{\alpha }\right)+\cos \left(\mathrm{\alpha }\right)x\,\mathrm{y*}\left(\mathrm{x*}\right)=y\left(x\right)\cos \left(\mathrm{\alpha }\right)+\sin \left(\mathrm{\alpha }\right)x\right\}$

$\mathrm{go}≔\mathrm{dchange}\left(\mathrm{tr}\,\mathrm{ODE}\,\left[\mathrm{ys}\left(\mathrm{xs}\right)\,\mathrm{xs}\right]\,u\mapsto \mathrm{simplify}\left(u\,\mathrm{trig}\right)\right)$

$\mathrm{go}≔\frac{\mathrm{x*}\left(\frac{\ⅆ}{\ⅆ\mathrm{x*}}\mathrm{y*}\left(\mathrm{x*}\right)\right)-\mathrm{y*}\left(\mathrm{x*}\right)}{\mathrm{y*}\left(\mathrm{x*}\right)\left(\frac{\ⅆ}{\ⅆ\mathrm{x*}}\mathrm{y*}\left(\mathrm{x*}\right)\right)+\mathrm{x*}}=H\left(\sqrt{{\mathrm{y*}\left(\mathrm{x*}\right)}^2+{\mathrm{x*}}^2}\right)$

$\mathrm{back}≔\mathrm{dchange}\left(\mathrm{itr}\,\mathrm{go}\,\left[x\,y\left(x\right)\right]\,u\mapsto \mathrm{simplify}\left(u\,\mathrm{trig}\right)\right)$

$\mathrm{back}≔\frac{-y\left(x\right)+x\left(\frac{\ⅆ}{\ⅆx}y\left(x\right)\right)}{x+y\left(x\right)\left(\frac{\ⅆ}{\ⅆx}y\left(x\right)\right)}=H\left(\sqrt{x^2+{y\left(x\right)}^2}\right)$

$\mathrm{II}≔\mathrm{Int}\left(\mathrm{Int}\left(\mathrm{diff}\left(f\left(x\,y\right)\,x\right)h\left({\left(x^2+y^2\right)}^{\frac{1}{2}}\right)\,x=0..1\right)\,y=0..1\right)$

$\mathrm{II}≔{\int }_0^1\left({\int }_0^1\left(\frac{\partial }{\partial x}f\left(x\,y\right)\right)h\left(\sqrt{x^2+y^2}\right)\ⅆx\right)\ⅆy$

$\mathrm{tr}≔\left\{x=s{r}^{\frac{1}{2}}\,y=s{\left(1-r\right)}^{\frac{1}{2}}\right\}\:$

$\mathrm{dchange}\left(\mathrm{tr}\,\mathrm{II}\,\mathrm{simplify}\right)$

$-\frac{\left({\int }_{\mathrm{\_\α}}^{\mathrm{\_\β}}\frac{h\left(s\right)\left|s\right|\left({\int }_{{\mathrm{\_l}}_r\left(s\right)}^{{\mathrm{\_u}}_r\left(s\right)}\frac{\left(2\left(\frac{\partial }{\partial r}f\left(r\,s\right)\right)r-\left(\frac{\partial }{\partial s}f\left(r\,s\right)\right)s-2\frac{\partial }{\partial r}f\left(r\,s\right)\right)\sqrt{r}}{\sqrt{\left|r\right|}\sqrt{\left|-1+r\right|}}\ⅆr\right)}{s}\ⅆs\right)}{2}$

$\mathrm{PDE}≔\mathrm{sqrt}\left(1-{\sin \left(s\right)}^2{k}^2\right)\mathrm{diff}\left(f\left(u\,v\,r\,s\right)\,s\right)=G\left(u\,v\,r\,s\right)$

$\mathrm{PDE}≔\sqrt{1-{\sin \left(s\right)}^2{k}^2}\left(\frac{\partial }{\partial s}f\left(u\,v\,r\,s\right)\right)=G\left(u\,v\,r\,s\right)$

$\mathrm{itr}≔\left\{w=\mathrm{Int}\left(\frac{1}{\mathrm{sqrt}\left(1-{\sin \left(r\right)}^2{k}^2\right)}\,r=0..s\right)\right\}$

$\mathrm{itr}≔\left\{w={\int }_0^s\frac{1}{\sqrt{1-{\sin \left(r\right)}^2{k}^2}}\ⅆr\right\}$

$\mathrm{tr}≔\left\{s=\mathrm{am}\left(w\right)\right\}\:$

$\mathrm{dchange}\left(\mathrm{tr}\,\mathrm{PDE}\,\left[w\right]\right)$

$\frac{\sqrt{1-{\sin \left(\mathrm{am}\left(w\right)\right)}^2{k}^2}\left(\frac{\partial }{\partial w}f\left(w\,u\,v\,r\right)\right)}{\frac{\ⅆ}{\ⅆw}\mathrm{am}\left(w\right)}=G\left(w\,u\,v\,r\right)$

$\mathrm{dchange}\left(\mathrm{tr}\,\mathrm{PDE}\,\mathrm{itr}\,\mathrm{simplify}\right)$

$\frac{\partial }{\partial w}f\left(w\,u\,v\,r\right)=G\left(w\,u\,v\,r\right)$