The canoni command tries to find a set of transformations from the original coordinates to the canonical coordinates by knowing the coefficients of the symmetry generator (infinitesimals) of a one-parameter Lie group.

If there is more than one derivative in the ODE, canoni requires an extra argument (anywhere in the calling sequence), say y(x), specifying the dependent variable.

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

$\mathrm{with}\left(\mathrm{DEtools}\,\mathrm{symgen}\,\mathrm{canoni}\right)$

$\left[\mathrm{symgen}\,\mathrm{canoni}\right]$

$\mathrm{with}\left(\mathrm{PDEtools}\,\mathrm{dchange}\right)$

$\left[\mathrm{dchange}\right]$

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

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

$\mathrm{infinitesimals}≔\mathrm{symgen}\left(\mathrm{ODE}\right)$

$\mathrm{infinitesimals}≔\left[\mathrm{\_\ξ}=y{x}^2\,\mathrm{\_\η}=-2a\right]$

$\mathrm{tr}≔\mathrm{canoni}\left(\mathrm{infinitesimals}\,y\left(x\right)\,s\left(r\right)\right)$

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

$\mathrm{itr}≔\mathrm{op}\left(1\,\left[\mathrm{solve}\left(\mathrm{tr}\,\left\{x\,y\left(x\right)\right\}\right)\right]\right)$

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

$\mathrm{ODE2}≔\mathrm{dchange}\left(\mathrm{itr}\,\mathrm{ODE}\,\left[r\,s\left(r\right)\right]\,\mathrm{simplify}\right)\:$

$\mathrm{ODE1}≔\mathrm{op}\left(\mathrm{solve}\left(\mathrm{ODE2}\,\left\{\mathrm{diff}\left(s\left(r\right)\,r\right)\right\}\right)\right)$

$\mathrm{ODE1}≔\frac{\ⅆ}{\ⅆr}s\left(r\right)=-\frac{1}{8a^{2}F\left(r\right)}$