DEtools
gensys
return the determining PDE system for either the coefficients of the symmetry generator or the integrating factor
Calling Sequence
Parameters
Description
Examples
gensys(ODE, [_xi = f(x, y, ..), _eta = g(x, y, ..)], y(x))
gensys(ODE, _mu = f(x, y, ..), y(x))
ODE
-
ODE (ordinary differential equation)
[_xi=.., _eta=..]
(optional) list containing the infinitesimal; right hand sides may be functions or algebraic expressions
_mu =..]
(optional) equation where the right hand side may be a function or an algebraic expression
y(x)
(optional) required if the ODE has derivatives of more than one function
Given an ODE, gensys returns the determining PDE system related to either its symmetries or its integrating factors, depending on the second argument. This command is typically used together with other commands to determine symmetries and integrating factors for ODEs. gensys first calls odepde to get the determining PDE and then splits it -- when possible -- by taking coefficients with respect to the derivatives of the dependent variable.
Concerning the optional second argument, gensys works like odepde; this second argument may be either an explicit form for the symmetry or the integrating factor to be taken into account at the time of calculating the determining PDE.
This function is part of the DEtools package, and so it can be used in the form gensys(..) only after executing the command with(DEtools). However, it can always be accessed through the long form of the command by using DEtools[gensys](..).
with⁡DEtools:
with⁡PDEtools:
declare⁡y⁡x,prime=x
y⁡x⁢will now be displayed as⁢y
derivatives with respect to⁢x⁢of functions of one variable will now be displayed with '
1. Determine if an integrating factor of the form μ⁡x,y exists for the second order nonlinear ODE
ODE≔diff⁡y⁡x,x,x=−exp⁡x⁢diff⁡y⁡x,x⁢x⁢y⁡x2+y⁡x2−diff⁡y⁡x,x⁢xx⁢y⁡x2⁢exp⁡x
ODE≔y''=−ⅇx⁢y'⁢x⁢y2+y2−y'⁢xx⁢y2⁢ⅇx
and if so, determine the integrating factor itself.
Set up the determining PDE system for these integrating factors (gensys enters in this step), then try to simplify it (perhaps using casesplit), then try to solve it. The determining PDE system is set up as follows:
gensys⁡ODE,_μ=μ⁡x,y⁡x
μy⁢ⅇx⁢x,μy,y⁢ⅇx⁢x,μx,y⁢ⅇx⁢x,−μ⁡x,y⁢xy2−ⅇx⁢y2−1⁢μx⁢xy2+μx,x⁢ⅇx⁢x+μy
This system (each equation above is assumed to be equal to zero) has four equations,
nops⁡
4
and by taking into account their integrability conditions it is simplified to
casesplit⁡
μx=μ⁡x,y,μy=0where
Thus, mu does not depend on y and the remaining equation leads to
μ=exp⁡x
μ=ⅇx
Note that this result can be tested using mutest
mutest⁡exp⁡x,ODE
0
and can be used to obtain a first integral (that is, to reduce the order of ODE).
firint⁡exp⁡x⁢ODE
1y+ⅇx⁢y'+ln⁡x+_C1=0
First integrals can be tested using firtest.
2. Determine if an integrating factor of the form mu⁡y,y'⁢ exists for the same second order ODE.
gensys⁡ODE,_μ=μ⁡y,_y1
−μ_y1,_y1⁢ⅇx⁢_y1⁢x2⁢y3+μ_y1,y⁢ⅇx⁢_y1⁢x2⁢y3−2⁢μ_y1⁢ⅇx⁢x2⁢y3+2⁢μy⁢ⅇx⁢x2⁢y3+μ_y1,_y1⁢_y1⁢x2⁢y−μ_y1,_y1⁢x⁢y3+2⁢μ_y1⁢x2⁢y,−μ_y1,y⁢ⅇx⁢_y12⁢x2⁢y3+_y12⁢μy,y⁢ⅇx⁢x2⁢y3+μ_y1,y⁢_y12⁢x2⁢y−μ_y1,y⁢_y1⁢x⁢y3−2⁢μ_y1⁢_y12⁢x2−μ_y1⁢_y1⁢x2⁢y+μ_y1⁢x⁢y3+μy⁢x⁢y3+μ_y1⁢y3−μ⁡y,_y1⁢y⁢x2
2
μ⁡y,_y1=0where
So this ODE has no integrating factor of the form mu⁡y,y'⁢.
3. Determine whether or not this ODE has point symmetries.
Set the determining system for the coefficients of the infinitesimal symmetry generator (so-called "infinitesimals") as follows:
gensys⁡ODE,ξ⁡x,y,η⁡x,y
ξy,y⁢ⅇx⁢x2⁢y3,2⁢ⅇx⁢ξy⁢x2⁢y3+ηy,y⁢ⅇx⁢x2⁢y3−2⁢ξx,y⁢ⅇx⁢x2⁢y3−2⁢ξy⁢x2⁢y,ⅇx⁢ξx⁢x2⁢y3+2⁢ηx,y⁢ⅇx⁢x2⁢y3−ξx,x⁢ⅇx⁢x2⁢y3+3⁢ξy⁢x⁢y3−ξx⁢x2⁢y+ξ⁡x,y⁢x2⁢y+2⁢η⁡x,y⁢x2,ⅇx⁢ηx⁢x2⁢y3+ηx,x⁢ⅇx⁢x2⁢y3+2⁢ξx⁢x⁢y3−ηy⁢x⁢y3−ξ⁡x,y⁢x⁢y3−ξ⁡x,y⁢y3−ηx⁢x2⁢y
(in above, [_xi = xi(x,y), _eta = eta(x,y)] also works). This system consists of 4 equations,
and in simplifying these equations by taking into account their integrability conditions one obtains
η⁡x,y=0,ξ⁡x,y=0where
showing that this ODE has no point symmetries.
See Also
DEtools[buildsym]
DEtools[equinv]
DEtools[firint]
DEtools[firtest]
DEtools[intfactor]
DEtools[Lie]
DEtools[mutest]
DEtools[odepde]
DEtools[symgen]
dsolve[Lie]
PDEtools
PDEtools[casesplit]
pdsolve
Download Help Document