look for the most general ODE invariant under a given symmetry
equinv([xi, eta], y(x))
equinv([[xi1, eta1], [xi2, eta2],...], y(x))
equinv([xi, eta], y(x), n)
equinv([[xi1, eta1], [xi2, eta2],...], y(x), n)
list of the coefficients of the symmetry generator (infinitesimals)
indeterminate function to be used in the result
positive integer, indicating that an nth order ODE is to be returned
equinv takes as first argument either a list of two algebraic expressions, to be seen as the infinitesimals of a one-parameter Lie group, or a list of such pairs representing many different symmetry groups. It returns, within the possibilities of the system, the most general ODE simultaneously invariant under those symmetry groups. The second argument is the indeterminate function of that ODE.
If the number n is given as extra argument, equinv looks for an invariant ODE of order n; otherwise, it will look for a first order ODE.
As a rule, if we represent the infinitesimal generator as
where x is the independent variable and y is the dependent one; equinv assumes that the order in which the infinitesimals xi and eta appear in the received list is ξ,η.
In principle, there are no restrictions on the number of pairs of infinitesimals you can give at one time to equinv. However, the complexity of the problem increases with the number of pairs of infinitesimals given, and the problem may not have a solution. Also, when an nth order ODE is required (extra argument n) the given infinitesimals may contain _y1, _y2, ... which, as a convention, represent the derivative of the "indeterminate function" of the problem. In short, equinv also accepts dynamical symmetries as arguments.
equinv and buildsym may be useful in connection with the odeadvisor command and the HINT=e1,e2 option of dsolve, in order to study the relationship between symmetry patterns and ODE patterns (see dsolve/Lie).
This function is part of the DEtools package, and so it can be used in the form equinv(..) only after executing the command with(DEtools). However, it can always be accessed through the long form of the command by using DEtools[equinv](..).
Three pairs of infinitesimals not containing a 2-D subalgebra
sym := [[1,1],[x,y],[x^2,y^2]];
The most general second order ODE simultaneously invariant under these three point symmetries
ODE2 := equinv(sym,y(x),2);
The fact that ODE2 is invariant under the three pairs of infinitesimals above can be tested using the symtest command (see symtest):
Two pairs of dynamical symmetries
dyn_sym := [[_y1,0],[0,1/_y1]];
The most general second order ODE simultaneously invariant under these two dynamical symmetries
ODE2 := equinv(dyn_sym,y(x),2);
Download Help Document
What kind of issue would you like to report? (Optional)