eta_k receives a pair of infinitesimals; k, the order of the required prolongation; and the dependent variable y(x), and returns the k-prolongation of eta (see infgen).

This command also works with dynamical symmetries, in which case the ODE that is assumed to be invariant under the given infinitesimals is also required as an argument. The right hand side of the given nth order ODE is then used to replace the nth order derivatives of the dependent variable appearing in the prolongation of eta.

If the infinitesimal generator of a Lie group is written as

X = (F -> xi(x,y)*diff(F,x) + eta(x,y)*diff(F,y));

$X=\left(F\→\mathrm{\ξ}\left(x\,y\right)\left(\frac{\partial }{\partial x}F\right)\+\mathrm{\η}\left(x\,y\right)\left(\frac{\partial }{\partial y}F\right)\right)$

where y(x) plays the role of the dependent variable, the k prolongation of eta is obtained as follows.

1) The linear operator A corresponding to the given ODE is built

A = (F -> diff(F,x) + _y1*diff(F,y) + _y2*diff(F,y,y) + `...` + rhs_ODE * diff(y(x),x$(n-1)));

$A=\left(F\→\frac{\partial }{\partial x}F\+\mathrm{\_y1}\left(\frac{\partial }{\partial y}F\right)\+\mathrm{\_y2}\left(\frac{{\partial }^2}{\partial y^2}F\right)\+\mathrm{...}\+\mathrm{rhs\_ODE}\left(\frac{{\ⅆ}^{n-1}}{\ⅆx^{n-1}}y\left(x\right)\right)\right)$

(rhs_ODE means the right hand side of the ODE).

2) The expression

eta[n] = A(eta[n-1])- diff(y(x),x$n)*A(xi);

${\mathrm{\eta }}_n=A\left({\mathrm{\eta }}_{n-1}\right)-\left(\frac{{\ⅆ}^n}{\ⅆx^n}y\left(x\right)\right)A\left(\mathrm{\xi }\right)$

is evaluated recursively, and in the case of dynamical symmetries, the highest order derivative is replaced by the right hand side of the ODE.

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

$\mathrm{with}\left(\mathrm{DEtools}\,\mathrm{eta\_k}\,\mathrm{infgen}\right)$

$\left[\mathrm{eta\_k}\,\mathrm{infgen}\right]$

$X≔\left[-y\,x\right]$

$X≔\left[-y\,x\right]$

$\mathrm{\eta }\left[1\right]=\mathrm{eta\_k}\left(X\,1\,y\left(x\right)\right)$

${\mathrm{\eta }}_1={\mathrm{\_y1}}^2+1$

$\mathrm{\eta }\left[2\right]=\mathrm{eta\_k}\left(X\,2\,y\left(x\right)\right)$

${\mathrm{\eta }}_2=3\mathrm{\_y2}\mathrm{\_y1}$

$\mathrm{\eta }\left[3\right]=\mathrm{eta\_k}\left(X\,3\,y\left(x\right)\right)$

${\mathrm{\eta }}_3=4\mathrm{\_y3}\mathrm{\_y1}+3{\mathrm{\_y2}}^2$

$\mathrm{infgen}\left(X\,3\,y\left(x\right)\right)$

$\mathrm{\_F1}\→-y\left(\frac{\partial }{\partial x}\mathrm{\_F1}\right)\+x\left(\frac{\partial }{\partial y}\mathrm{\_F1}\right)\+\left({\mathrm{\_y1}}^2\+1\right)\left(\frac{\partial }{\partial \mathrm{\_y1}}\mathrm{\_F1}\right)\+3\mathrm{\_y2}\mathrm{\_y1}\left(\frac{\partial }{\partial \mathrm{\_y2}}\mathrm{\_F1}\right)\+\left(4\mathrm{\_y1}\mathrm{\_y3}\+3{\mathrm{\_y2}}^2\right)\left(\frac{\partial }{\partial \mathrm{\_y3}}\mathrm{\_F1}\right)$

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

$\mathrm{ODE}≔\frac{{\ⅆ}^2}{\ⅆx^2}y\left(x\right)=F\left(x\,y\left(x\right)\,\frac{\ⅆ}{\ⅆx}y\left(x\right)\right)$

$X≔\left[\mathrm{\xi }\left(x\,y\right)\,\mathrm{\eta }\left(x\,y\right)\right]$

$X≔\left[\mathrm{\xi }\left(x\,y\right)\,\mathrm{\eta }\left(x\,y\right)\right]$

$\mathrm{\eta }\left[0\right]=\mathrm{eta\_k}\left(X\,0\,\mathrm{ODE}\right)$

${\mathrm{\eta }}_0=\mathrm{\eta }\left(x\,y\right)$

$\mathrm{\eta }\left[1\right]=\mathrm{eta\_k}\left(X\,1\,\mathrm{ODE}\right)$

${\mathrm{\eta }}_1=-\left(\frac{\partial }{\partial y}\mathrm{\xi }\left(x\,y\right)\right){\mathrm{\_y1}}^2+\left(\frac{\partial }{\partial y}\mathrm{\eta }\left(x\,y\right)-\frac{\partial }{\partial x}\mathrm{\xi }\left(x\,y\right)\right)\mathrm{\_y1}+\frac{\partial }{\partial x}\mathrm{\eta }\left(x\,y\right)$

$\mathrm{\eta }\left[2\right]=\mathrm{eta\_k}\left(X\,2\,\mathrm{ODE}\right)$

${\mathrm{\eta }}_2=-\left(\frac{{\partial }^2}{\partial y^2}\mathrm{\xi }\left(x\,y\right)\right){\mathrm{\_y1}}^3+\left(-2\frac{{\partial }^2}{\partial x\partial y}\mathrm{\xi }\left(x\,y\right)+\frac{{\partial }^2}{\partial y^2}\mathrm{\eta }\left(x\,y\right)\right){\mathrm{\_y1}}^2+\left(-3F\left(x\,y\,\mathrm{\_y1}\right)\left(\frac{\partial }{\partial y}\mathrm{\xi }\left(x\,y\right)\right)+2\frac{{\partial }^2}{\partial x\partial y}\mathrm{\eta }\left(x\,y\right)-\frac{{\partial }^2}{\partial x^2}\mathrm{\xi }\left(x\,y\right)\right)\mathrm{\_y1}+F\left(x\,y\,\mathrm{\_y1}\right)\left(\frac{\partial }{\partial y}\mathrm{\eta }\left(x\,y\right)\right)-2F\left(x\,y\,\mathrm{\_y1}\right)\left(\frac{\partial }{\partial x}\mathrm{\xi }\left(x\,y\right)\right)+\frac{{\partial }^2}{\partial x^2}\mathrm{\eta }\left(x\,y\right)$

$X≔\left[\mathrm{\xi }\left(x\,y\,\mathrm{\_y1}\right)\,\mathrm{\eta }\left(x\,y\,\mathrm{\_y1}\right)\right]$

$X≔\left[\mathrm{\xi }\left(x\,y\,\mathrm{\_y1}\right)\,\mathrm{\eta }\left(x\,y\,\mathrm{\_y1}\right)\right]$

$\mathrm{\eta }\left[1\right]=\mathrm{eta\_k}\left(X\,1\,\mathrm{ODE}\right)$

${\mathrm{\eta }}_1=-\left(\frac{\partial }{\partial y}\mathrm{\xi }\left(x\,y\,\mathrm{\_y1}\right)\right){\mathrm{\_y1}}^2+\left(-F\left(x\,y\,\mathrm{\_y1}\right)\left(\frac{\partial }{\partial \mathrm{\_y1}}\mathrm{\xi }\left(x\,y\,\mathrm{\_y1}\right)\right)+\frac{\partial }{\partial y}\mathrm{\eta }\left(x\,y\,\mathrm{\_y1}\right)-\frac{\partial }{\partial x}\mathrm{\xi }\left(x\,y\,\mathrm{\_y1}\right)\right)\mathrm{\_y1}+F\left(x\,y\,\mathrm{\_y1}\right)\left(\frac{\partial }{\partial \mathrm{\_y1}}\mathrm{\eta }\left(x\,y\,\mathrm{\_y1}\right)\right)+\frac{\partial }{\partial x}\mathrm{\eta }\left(x\,y\,\mathrm{\_y1}\right)$