The input L is a differential operator. The output of this procedure is a linear differential operator M of minimal order such that for every set of m solutions $\mathrm{y1},\mathrm{...},\mathrm{ym}$ of L the product $\mathrm{y1}\mathrm{y2}\dots \mathrm{ym}$ is a solution of $M$.

The argument domain describes the differential algebra. If this argument is the list $\left[\mathrm{Dt}\,t\right]$, then the differential operators are notated with the symbols $\mathrm{Dt}$ and $t$. They are viewed as elements of the differential algebra $C\left(t\right)$ $\left[\mathrm{Dt}\right]$ where $C$ is the field of constants.

If the argument domain is omitted then the differential specified by the environment variable _Envdiffopdomain is used. If this environment variable is not set then the argument domain may not be omitted.

Instead of a differential operator, the input can also be a linear homogeneous differential equation having rational function coefficients. In this case the third argument must be the dependent variable.

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

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

$\mathrm{\_Envdiffopdomain}≔\left[\mathrm{Dx}\,x\right]\:$

$L≔{\mathrm{Dx}}^2+a\left(x\right)\mathrm{Dx}+b\left(x\right)$

$L≔{\mathrm{Dx}}^2+a\left(x\right)\mathrm{Dx}+b\left(x\right)$

$M≔\mathrm{symmetric\_power}\left(L\,2\right)$

$M≔{\mathrm{Dx}}^3+3a\left(x\right){\mathrm{Dx}}^2+\left(2{a\left(x\right)}^2+\frac{\ⅆ}{\ⅆx}a\left(x\right)+4b\left(x\right)\right)\mathrm{Dx}+4b\left(x\right)a\left(x\right)+2\frac{\ⅆ}{\ⅆx}b\left(x\right)$

$\mathrm{ODE}≔\mathrm{diff}\left(y\left(x\right)\,x\,x\right)+a\left(x\right)y\left(x\right)=0$

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

$\mathrm{sol}≔\mathrm{dsolve}\left(\mathrm{ODE}\right)$

$\mathrm{sol}≔y\left(x\right)=\mathrm{DESol}\left(\left\{a\left(x\right)\mathrm{\_Y}\left(x\right)+\frac{{\ⅆ}^2}{\ⅆx^2}\mathrm{\_Y}\left(x\right)\right\}\,\left\{\mathrm{\_Y}\left(x\right)\right\}\right)$

$\mathrm{ODE\_2}≔\mathrm{symmetric\_power}\left(\mathrm{ODE}\,2\,y\left(x\right)\right)$

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

$\mathrm{dsolve}\left(\mathrm{ODE\_2}\,y\left(x\right)\right)$

$y\left(x\right)={\mathrm{DESol}\left(\left\{a\left(x\right)\mathrm{\_Y}\left(x\right)+\frac{{\ⅆ}^2}{\ⅆx^2}\mathrm{\_Y}\left(x\right)\right\}\,\left\{\mathrm{\_Y}\left(x\right)\right\}\right)}^2$

$\mathrm{ODE\_3}≔\mathrm{symmetric\_power}\left(\mathrm{ODE}\,3\,y\left(x\right)\right)$

$\mathrm{ODE\_3}≔\left(9{a\left(x\right)}^2+3\frac{{\ⅆ}^2}{\ⅆx^2}a\left(x\right)\right)y\left(x\right)+10\left(\frac{\ⅆ}{\ⅆx}a\left(x\right)\right)\left(\frac{\ⅆ}{\ⅆx}y\left(x\right)\right)+10a\left(x\right)\left(\frac{{\ⅆ}^2}{\ⅆx^2}y\left(x\right)\right)+\frac{{\ⅆ}^4}{\ⅆx^4}y\left(x\right)$

$\mathrm{dsolve}\left(\mathrm{ODE\_3}\,y\left(x\right)\right)$

$y\left(x\right)={\mathrm{DESol}\left(\left\{a\left(x\right)\mathrm{\_Y}\left(x\right)+\frac{{\ⅆ}^2}{\ⅆx^2}\mathrm{\_Y}\left(x\right)\right\}\,\left\{\mathrm{\_Y}\left(x\right)\right\}\right)}^3$