DEtools
transinv
look for the set of transformations of variables which leave the ODE invariant
Calling Sequence
Parameters
Description
Examples
transinv([xi(x, y), eta(x, y)], y(x), s(r))
[xi(x, y), eta(x, y)]
-
list of the coefficients of the infinitesimal symmetry generator (infinitesimals)
y(x)
dependent variable
s(r)
new dependent variable
transinv looks for the set of transformations of variables which leave an ODE invariant, by using the coefficients of a symmetry generator (infinitesimals) for it. These transformations are actually the finite form of the one-parameter Lie group of invariance of the ODE.
This function is part of the DEtools package, and so it can be used in the form transinv(..) only after executing the command with(DEtools). However, it can always be accessed through the long form of the command by using DEtools[transinv](..).
An ODE with an arbitrary function F
with⁡DEtools:
with⁡PDEtools:
ODE_y≔diff⁡y⁡x,x=2⁢a−x2⁢y⁡x+x2⋅2⁢F⁡x⁢y⁡x2−4⁢ax⁢a
ODE_y≔ⅆⅆxy⁡x=2⁢a−x2⁢y⁡x+2⁢x2⁢F⁡x⁢y⁡x2−4⁢ax⁢a
odeadvisor⁡ODE_y
_1st_order,_with_symmetry_[F(x),G(y)]
A pair of infinitesimals for it
infinitesimals≔symgen⁡ODE_y
infinitesimals≔_ξ=x2⁢y,_η=−2⁢a
The transformation of variables which leaves ODE_y invariant
tr≔transinv⁡infinitesimals,y⁡x,s⁡r
tr≔r=x_α2⁢a⁢x−_α⁢x⁢y⁡x+1,s⁡r=−2⁢_α⁢a+y⁡x
Note the introduction of _alpha, representing the parameter of the Lie group. Now, to check the invariance of ODE_y under this group, you can change variables as follows:
itr≔solve⁡tr,x,y⁡x
itr≔x=r1+_α2⁢a⁢r+s⁡r⁢_α⁢r,y⁡x=s⁡r+2⁢_α⁢a
The change of variables
ODE_s≔dchange⁡itr,ODE_y,r,s⁡r:
diff⁡s⁡r,r=solve⁡ODE_s,diff⁡s⁡r,r
ⅆⅆrs⁡r=2⁢ar2⁢2⁢F⁡s⁡r2⁢r−4⁢ar⁢a−s⁡r
As can be seen above, we arrived at the original ODE_y just changing x, y by r, s (this is the meaning of "leaving the ODE invariant").
See Also
canoni
dchange
equinv
odeadvisor
PDEtools
symgen
Download Help Document