DEtools - Maple Help

Online Help

All Products    Maple    MapleSim


Home : Support : Online Help : Mathematics : Differential Equations : Lie Symmetry Method : Commands for ODEs : DEtools/symgen

DEtools

  

symgen

  

look for a symmetry generator for a given ODE

 

Calling Sequence

Parameters

Description

Examples

References

Calling Sequence

symgen(ODE, y(x), way=xxx, HINT=...)

Parameters

ODE

-

ordinary differential equation

y(x)

-

(optional) dependent variable

way=xxx

-

(optional) algorithm where xxx is one of 2, 3, 4, 5, 6, abaco1, abaco2, all, exp_sym, formal, mu, patterns, and pdsolve; forces the use of only one xxx of the nine internal algorithms for determining the coefficients of the symmetry generator (infinitesimals)

HINT = [e1, e2]

-

(optional) e1 and e2 indicate a possible functional form for the infinitesimals

HINT = [[e1, e2], [e3, e4], ...]

-

(optional) list of possible functional forms for the infinitesimals

Description

• 

The symgen command looks for the infinitesimals ξx,y and ηx,y of a one-parameter Lie group that make the given ODE invariant (see odepde). For instance, in the case of a first order ODE of the form y`' `=Φx,y, symgen looks for a pair of functions η,ξ that satisfies the following:

with(DEtools):

odepde(diff(y(x),x)=Phi(x,y(x)));

x_ηx,y+y_ηx,yx_ξx,yΦx,yy_ξx,yΦx,y2_ξx,yxΦx,y_ηx,yyΦx,y

(1)
• 

The symgen command is the engine that implements symmetry methods in dsolve (see dsolve,Lie). All nine of symgen's schemes assume that the highest derivative present in the given ODE can be isolated. By default, symgen looks for as many pairs of infinitesimals as the differential order of the ODE. However, you can change the number of symmetries by assigning the environment variable _Env_ODE_syms to a positive integer (the number of desired symmetries); symgen then returns an answer as soon as these symmetries are found, or when it has tried all nine schemes.

• 

Although in the framework of symgen, the definition of a symmetry includes dynamical symmetries -- that is: derivatives of the dependent variable can be present in the infinitesimals -- when using the optional argument way = formal, symgen will compute only point symmetries -- that is: infinitesimals depending only on the independent and dependent variables.

• 

The options available for this command are the same way=xxx and HINT=xxx options of dsolve,Lie. The HINT option is particularly useful since it enables you to try your own heuristics for finding the infinitesimals if none of symgen's schemes are successful.

• 

A brief summary of symgen's schemes is as follows:

• 

The first algorithms, called formal, formulates a linear PDE system for the infinitesimals ξ,η, then formally triangularize this system - using differential algebra techniques - finally tacking the resulting uncoupled system. When successful, this algorithm returns the complete set of point symmetries admitted by a given ODE of order 2 or higher. NOTE: this algorithm is advantageous mainly for 2nd and higher order ODEs. The algorithm works as well in the case of first order ODEs, but in this case the subproblems it will need to solve to find the symmetries are as difficult to solve as the first order ODE itself.

• 

The second algorithm, called abaco1, consists of two sets of four schemes each. The first four schemes look for xi, eta taking one of the infinitesimals as 0, and the other as a function of only one variable, namely ξ=0,η=fx, ξ=0,η=fy, ξ=fx,η=0, ξ=fy,η=0.

  

The possible success of each case is determined by algebraic factorization of the PDE above, and when this possibility is confirmed, an auxiliary ODE is built and sent to dsolve to determine the explicit resulting form for f.

  

The second sequence of schemes looks for xi, eta taking one of the infinitesimals as 0, and the other one as an expression containing both x and y, namely ξ=0,η=fxgy,ξ=fxgy,η=0, where g is an algebraic expression built by extracting factors from the numerator and denominator of the ODE to be solved, and f is an unknown function determined by solving auxiliary ODEs as in the first four schemes.

• 

The second, third and fourth algorithms consist of bivariate ansatze which are polynomial and rational in x,y. You can set the maximum degree of the polynomials by assigning the global 'dgun' variable; otherwise, it is set automatically by an internal routine.

• 

The fifth algorithm looks for a solution for xi and eta as follows:

  

1) a basis of functions and algebraic objects is built by taking, from the given ODE, all the known functions (their derivative rules are known by the Maple system) and composite algebraic objects, together with their derivatives, as well as all the unknown functions;

  

2) a polynomial of degree two in such objects is built; its coefficients, in turn, are polynomials of degree d (the dgun variable mentioned above) in x,y, with undetermined coefficients.

  

This ansatz is introduced into a PDE equivalent to the PDE for xi, eta above, resulting in a system of algebraic equations for the undetermined coefficients mentioned above.

• 

The sixth algorithm looks for a solution for xi and eta when the ODE is of first order and of Riccati type. This algorithm consists of three different schemes for finding the infinitesimals, exploring changes of variables, and tackling the PDE for the infinitesimals using separation of variables.

• 

The seventh algorithm, called pdsolve, tries to solve a PDE equivalent to the PDE for xi, eta using the pdsolve command from the PDEtools package.

• 

The eighth algorithm, called abaco2, consists of four schemes for determining the infinitesimals xi, eta, making use of the HINT option of symgen, with two indeterminate functions of only one variable each, namely

  

ξ=gx,η=fx,ξ=gx,η=fy,ξ=gy,η=fx,ξ=gy,η=fy.

• 

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

Examples

First order ODEs

Kamke's ODE 120; for better displaying of formulas we use here the PDEtools[declare] command:

withDEtools:

PDEtools[declare]yx,prime=x

yxwill now be displayed asy

derivatives with respect toxof functions of one variable will now be displayed with '

(2)

If infolevel is set to a greater integer (possible settings are 1 through 5), more detailed information about the computation method is displayed.

infolevelsymgen2

infolevelsymgen:=2

(3)

ode1xⅆⅆxyxyxxlnx2yx+2=0

ode1:=xy'yxlnx2y+2=0

(4)

symgenode1

   -> Computing symmetries using: way = 2-> Computing symmetries using: way = 3
   -> Computing symmetries using: way = 4

   <- successful computation of symmetries.

1&comma;2yx

_&xi;&equals;1&comma;_&eta;&equals;2yx

(5)

Kamke's ODE 236

ode2xyx&plus;4&DifferentialD;&DifferentialD;xyxyx22yx2x&equals;0

ode2:=xy&plus;4y'y22y2x&equals;0

(6)

odeadvisorode2

_rational&comma;_Abel&comma;2nd type&comma;class B

(7)

sym2symgenode2

   -> Computing symmetries using: way = 2

   <- successful computation of symmetries.

0&comma;xy2&plus;y34x2&plus;4xyxy&plus;4&comma;0&comma;x33x2y&plus;2xy24x2&plus;4xyxy&plus;4

sym2:=_&xi;&equals;0&comma;_&eta;&equals;xyx2y4y&plus;4&comma;_&xi;&equals;0&comma;_&eta;&equals;y2&plus;4xxyxy&plus;4

(8)

These symmetries can be tested using symtest.

mapsymtest&comma;sym2&comma;ode2

0&comma;0

(9)

Kamke's ODE 357 (no classification)

ode3x&DifferentialD;&DifferentialD;xyxlnxsinyx&plus;cosyx1xcosyx&equals;0

ode3:=xy'lnxsiny&plus;cosy1xcosy&equals;0

(10)

sym3symgenode3&comma;way&equals;5

   -> Computing symmetries using: way = 5

   <- successful computation of symmetries.

0&comma;cos2y&plus;1lnxsiny&comma;0&comma;xcos2y2cosylnx&plus;xlnxsiny

sym3:=_&xi;&equals;0&comma;_&eta;&equals;cos2y&plus;1lnxsiny&comma;_&xi;&equals;0&comma;_&eta;&equals;xcos2y2cosylnx&plus;xlnxsiny

(11)

mapsymtest&comma;sym3&comma;ode3

0&comma;0

(12)

The answer is returned by dsolve using the first of these symmetries:

dsolveode3

   -> Computing symmetries using: way = 3-> Computing symmetries using: way = 4
   -> Computing symmetries using: way = 5

   <- successful computation of symmetries.

0&comma;cos2y&plus;1lnxsiny&comma;0&comma;xcos2y2cosylnx&plus;xlnxsiny

y&equals;arccoslnx_C1&plus;x

(13)

Here is an ODE with an arbitrary function Fx&comma;y (no classification):

ode4&DifferentialD;&DifferentialD;xyx&equals;Fyx323&ExponentialE;x2&ExponentialE;xyx12

ode4:=y'&equals;Fy3&sol;232&ExponentialE;x&ExponentialE;xy

(14)

symgenode4&comma;way&equals;abaco2

   -> Computing symmetries using: way = abaco2

   <- successful computation of symmetries.

&ExponentialE;x&comma;1y

_&xi;&equals;&ExponentialE;x&comma;_&eta;&equals;1y

(15)

Second order ODEs

Kamke's nonlinear second order example number 33; in this example, only formal methods succeed in computing a 2-D symmetry group depending on an arbitrary function f(x):

k33&DifferentialD;2&DifferentialD;x2yx&plus;yx&plus;3fx&DifferentialD;&DifferentialD;xyxyx3&plus;yx2fx&plus;yx&DifferentialD;&DifferentialD;xfx&plus;2fx2&equals;0

k33:=y''&plus;y&plus;3fxy'y3&plus;y2fx&plus;yf'&plus;2fx2&equals;0

(16)

sym33symgenk33&comma;yx&comma;way&equals;formal

   -> Computing symmetries using: way = formal

   <- successful computation of symmetries.

1&ExponentialE;&int;fx&DifferentialD;x&comma;fxy&ExponentialE;&int;fx&DifferentialD;x&comma;&int;&ExponentialE;&int;fx&DifferentialD;x&DifferentialD;x&ExponentialE;&int;fx&DifferentialD;x&comma;&int;&ExponentialE;&int;fx&DifferentialD;x&DifferentialD;xfx&plus;&ExponentialE;&int;fx&DifferentialD;xy&ExponentialE;&int;fx&DifferentialD;x

sym33:=_&xi;&equals;&ExponentialE;&int;fx&DifferentialD;x&comma;_&eta;&equals;&ExponentialE;&int;fx&DifferentialD;xfxy&comma;_&xi;&equals;&int;&ExponentialE;&int;fx&DifferentialD;x&DifferentialD;x&ExponentialE;&int;fx&DifferentialD;x&comma;_&eta;&equals;yfx&ExponentialE;&int;fx&DifferentialD;x&int;&ExponentialE;&int;fx&DifferentialD;x&DifferentialD;x&plus;1

(17)

mapsymtest&comma;sym33&comma;k33&comma;yx

0&comma;0

(18)

This example shows an ODE containing an arbitrary function G:

ode5&DifferentialD;2&DifferentialD;x2yx&equals;4&DifferentialD;&DifferentialD;xyx24yxG&DifferentialD;&DifferentialD;xyxx1yx2

ode5:=y''&equals;4y'24yGxy'12y

(19)

odeadvisorode5

_2nd_order&comma;_with_linear_symmetries

(20)

This one has two pairs of infinitesimals for ode5 (a point and a dynamical symmetry):

sym5symgenode5

   -> Computing symmetries using: way = 3

   <- successful computation of symmetries.-> Computing symmetries using: way = 2
   -> Computing symmetries using: way = abaco1

   <- successful computation of symmetries.

x&comma;0

x&comma;0&comma;0&comma;1_y1

sym5:=_&xi;&equals;x&comma;_&eta;&equals;0&comma;_&xi;&equals;0&comma;_&eta;&equals;1_y1

(21)

mapsymtest&comma;sym5&comma;ode5

0&comma;0

(22)

ode6&DifferentialD;2&DifferentialD;x2yx&equals;1&DifferentialD;&DifferentialD;xyx2x22yx&DifferentialD;&DifferentialD;xyxx&plus;yx2x2

ode6:=y''&equals;x2y'22xyy'&plus;y2x2

(23)

symgenode6

   -> Computing symmetries using: way = 3

   <- successful computation of symmetries.

0&comma;x&comma;x&comma;0

_&xi;&equals;0&comma;_&eta;&equals;x&comma;_&xi;&equals;x&comma;_&eta;&equals;0

(24)

Below, G is again an arbitrary function of its arguments (the ODE does not match any recognizable pattern).

ode7&DifferentialD;2&DifferentialD;x2yx&equals;1&DifferentialD;&DifferentialD;xyx2yxlnyx&DifferentialD;&DifferentialD;xyx&plus;yxlnyx2&ExponentialE;2xG&DifferentialD;&DifferentialD;xyx&ExponentialE;xlnyxlnyxyx

ode7:=y''&equals;y'2ylnyy'&plus;ylny2&ExponentialE;2xGy'&ExponentialE;xlnylnyy

(25)

odeadvisorode7

NONE

(26)

Two point symmetries for ode7. Although symgen succeeds in finding these symmetries without extra arguments, specifying a way might save some time.

symgenode7&comma;way&equals;abaco1

   -> Computing symmetries using: way = abaco1

   <- successful computation of symmetries.

0&comma;lny&comma;&ExponentialE;x&comma;0

_&xi;&equals;0&comma;_&eta;&equals;lny&comma;_&xi;&equals;&ExponentialE;x&comma;_&eta;&equals;0

(27)

High order ODEs

A third order example

ode8&DifferentialD;3&DifferentialD;x3yx&equals;13yx8x3120&DifferentialD;2&DifferentialD;x2yxx2&plus;26&DifferentialD;&DifferentialD;xyxx8x3124&DifferentialD;&DifferentialD;xyxx3&DifferentialD;2&DifferentialD;x2yx20&DifferentialD;&DifferentialD;xyx2x28x3yx

ode8:=y'''&equals;138yx31820x2y''&plus;26xy'x31824x3y'y''20x2y'2x3y

(28)

odeadvisorode8

_3rd_order&comma;_with_linear_symmetries

(29)

sym8symgenode8

   -> Computing symmetries using: way = 2

   <- successful computation of symmetries.-> Computing symmetries using: way = 3

   <- successful computation of symmetries.

0&comma;xy

0&comma;xy&comma;0&comma;y&comma;x&comma;0

sym8:=_&xi;&equals;0&comma;_&eta;&equals;y&comma;_&xi;&equals;x&comma;_&eta;&equals;0&comma;_&xi;&equals;0&comma;_&eta;&equals;xy

(30)

mapsymtest&comma;sym8&comma;ode8

0&comma;0&comma;0

(31)

Here is how to find a solution for ode8 using these symmetries:

ans8dsolveode8

ans8:=y&equals;229145x94&plus;142929_C3145x941429_C229841x94&plus;1429_C3841x941429_C2&plus;1682_C1x&comma;y&equals;229145x94&plus;142929_C3145x941429_C229841x94&plus;1429_C3841x941429_C2&plus;1682_C1x

(32)

References

  

Cheb-Terrab, E.S. "A Computational Approach for the Exact Solving of Systems of Partial Differential Equations." Submitted to Computer Physics Communications, 2001.

  

Cheb-Terrab, E.S.; Duarte, L.G.S.; and da Mota, L.A.C.P. Untitled Article. Computer Physics Communications. Vol. 101. (1997): 254.

See Also

buildsym

canoni

DEtools

dsolve

dsolve,Lie

dsolve,references

equinv

infolevel

odeadvisor

odepde

PDEtools

pdsolve

symtest

 


Download Help Document

Was this information helpful?



Please add your Comment (Optional)
E-mail Address (Optional)
What is ? This question helps us to combat spam