DEtools - Maple Programming Help

Online Help

All Products    Maple    MapleSim


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

DEtools

  

invariants

  

calculate the differential invariants of a given one-parameter Lie group

 

Calling Sequence

Parameters

Description

Examples

Calling Sequence

invariants([xi, eta], y(x))

invariants([xi, eta], k, y(x))

invariants([xi, eta], k, y(x), ODE, scheme)

Parameters

[xi, eta]

-

list of the coefficients of the symmetry generator (infinitesimals)

k

-

positive integer indicating the order of the differential invariants required

y(x)

-

'dependent variable'; it can be any indeterminate function of one variable

ODE

-

(optional) ODE invariant under the given infinitesimals; required only if they represent dynamical symmetries

scheme

-

(optional); can be can, dif, or generic

Description

• 

The invariants command receives a pair of infinitesimals; an indication of the dependent variable, say y(x); and optionally k, the order of the differential invariants required. It returns a sequence of differential invariants, starting with the one of order zero and finishing with the one of order k. If k is not given, the command returns just the differential invariants of order zero and one.

  

The differential invariants are the solutions of XInvariant=0, where X represents the infinitesimal operator,

X = (f -> xi*diff(f,x) + eta*diff(f,y)+ eta[1]*diff(f,_y1)+`...`);

X=f→ξxf+ηyf+η1_y1f+...

(1)
  

and xi, eta, eta[1], `...` are the infinitesimals and the extensions of eta (see eta_k).

• 

The invariants command uses three different approaches. To indicate the use of one of these, use the optional argument dif, can, or generic.

  

The generic approach works with point and dynamical symmetries, requires giving the ODE (as the ODE parameter) that admits the given symmetries, and consists of three steps.

1. 

Set up a first order PDE for the invariants using DEtools[infgen] replacing the highest derivative by the right-hand side of the given ODE.

2. 

Solve that PDE using the characteristic strip method by calling pdsolve.

3. 

The PDE solution depends on an arbitrary function _Fn...; the arguments of that function are the differential invariants returned.

  

The dif approach is equivalent to the generic method and includes some optimizations possible in this case, such that, for example, it is not necessary to give the ODE supposed to admit the given symmetry. If the ODE is given this method only works with ODEs of differential order two.

  

The can approach only works with point symmetries, does not require giving the ODE admitting the symmetry, and consists of computing the canonical coordinates associated to the given symmetry, and then extracting the differential invariants using standard techniques (see for instance Stephani's book in dsolve,references).

• 

By default, invariants computes the differential invariants according to the following.

  

* If the symmetry is of dynamical type, the ODE parameter is required. If the ODE is of order 2, use the dif approach. Otherwise, use generic.

  

* Otherwise, the symmetry is a point-like symmetry. Use dif. If the algorithm fails, use can.

• 

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

Examples

The infinitesimals xi and eta of the one-parameter rotation group and the first extension of the related infinitesimal generator

withDEtools:

Xy,x

X:=y,x

(2)

invariantsX,yx

_I0=x2+y2,_I1=_y1xy_y1y+x

(3)

The differential invariants of order 0, 1, and 2

invariantsX,2,yx

_I0=x2+y2,_I1=_y1xy_y1y+x,_I2=12_y2x2+_y13+_y1x_y2y2+_y121y_y1y+x3

(4)

Here is a case of a dynamical symmetry in the context of second order ODEs. When working with dynamical symmetries, the ODE itself is required as an argument.

X0,_y12

X:=0,_y12

(5)

The general second order ODE invariant under X

general_ODEequinvX,yx,2

general_ODE:=ⅆ2ⅆx2yx=ⅆⅆxyxⅆⅆxyx_F1ⅆⅆxyx,12xⅆⅆxyx+yx+x

(6)

A particular case of the above

ODEevalsubs_F1=a,b→ab,general_ODE

ODE:=ⅆ2ⅆx2yx=ⅆⅆxyxⅆⅆxyx212xⅆⅆxyx+yx+x

(7)

The related differential invariants of order 0, 1

invariantsX,1,ODE

_I0=x,_I1=12_y1x2yⅇ16_y13

(8)

See Also

DEtools

DEtools[canoni]

DEtools[equinv]

DEtools[eta_k]

DEtools[infgen]

DEtools[odepde]

DEtools[symgen]

dsolve[Lie]

PDEtools

 


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