return the nth order ODE having a given integrating factor - Maple Help

Online Help

All Products    Maple    MapleSim


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

DEtools[redode] - return the nth order ODE having a given integrating factor

Calling Sequence

redode(mu, n, y(x), R)

Parameters

mu

-

integrating factor - one or a list of them - depending on x,y,...,_ym, m=n1

n

-

order of the requested ODE

y(x)

-

dependent variable

R

-

(optional) expected reduced ODE depending on x,y,...,_yk 

Description

• 

Given an integrating factor (x,y,..._yn), where _yn represents the nth order derivative, or a list of them, redode returns the ODE of order n having all the values of mu as integrating factors.

• 

This command is useful to identify the general ODE problem related to a given mu, and to understand the possible links between the integrating factor scheme for reducing the order and other reduction schemes (for example, symmetries).

• 

When the expected reduced ODE, here called R, is also given as argument, its differential order, say m, must satisfy m<n. The routine then proceeds as follows. First, it redefines R by differentiating it such that its differential order becomes n1. Then it performs a test to see if the problem is solvable. The test for solvability is

mu(x,y,`...`,_ym) = nu(x,y,`...`,_y||(m-1))*diff(R(x,y,`...`,_ym),_ym);

&mu;x&comma;y&comma;...&comma;_ym&equals;&nu;x&comma;y&comma;...&comma;_y||m1_ymRx&comma;y&comma;...&comma;_ym

(1)
  

for some function nux,y,...,_ym. If the problem is solvable, redode then returns an nth order ODE satisfying

mu(x,y,`...`,diff(y(x),x$m)) * ODE=Diff(R(x,y(x),`...`,diff(y(x),x$m)),x);

&mu;x&comma;y&comma;...&comma;&DifferentialD;m&DifferentialD;xmyxODE&equals;xRx&comma;yx&comma;...&comma;&DifferentialD;m&DifferentialD;xmyx

(2)
  

Note that the reduced ODE being differentiated on the right-hand side is defined up to a constant. For example, R&plus;_C1 also satisfies the above equation.

• 

When the given mu does not depend on _ym and R is nonlinear in _ym, the requested nth order ODE exists if R can be solved for _ym; then nu can be determined as the ratio μd_ymdR.

• 

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

Examples

The redode command is interesting mainly as a tool for generating solving schemes for given ODE families; we illustrate this with two examples.

withDEtools&colon;

1. Consider the second order ODE family having an integrating factor μ=Fx -- an arbitrary function -- such that the reduced ODE has the same integrating factor. We want to set up an algorithm such that, given a second order linear ODE,

ode_psi:=&DifferentialD;2&DifferentialD;x2yx&equals;&psi;1x&DifferentialD;&DifferentialD;xyx&plus;&psi;2xyx&plus;&psi;3x

ode_psi:=&DifferentialD;2&DifferentialD;x2yx&equals;&psi;1x&DifferentialD;&DifferentialD;xyx&plus;&psi;2xyx&plus;&psi;3x

(3)

where there are no restrictions on psi1x&comma;psi2x, or psi3x, the scheme determines if the ODE belongs to the family just described, and if so it also determines Fx. The knowledge of Fx is enough to completely solve the given ODE.

To start with, we obtain the general first order ODE having Fx as integrating factor using:

ode_1:=redodeFx&comma;yx

ode_1:=&DifferentialD;&DifferentialD;xyx&equals;&DifferentialD;&DifferentialD;xFxyx&plus;_F1xFx

(4)

where _F1x is an arbitrary function. To obtain the second order ODE mentioned above, we pass ode_1 as an argument (playing the role of the reduced ODE) with the integrating factor Fx to obtain:

ode_2:=redodeFx&comma;yx&comma;ode_1

ode_2:=&DifferentialD;2&DifferentialD;x2yx&equals;2&DifferentialD;&DifferentialD;xFx&DifferentialD;&DifferentialD;xyx&plus;&DifferentialD;2&DifferentialD;x2Fxyx&plus;&DifferentialD;&DifferentialD;x_F1xFx

(5)

Taking this general ODE pattern as our starting point, we set up the required solving scheme by comparing coefficients in ode_2 and ode_psi obtaining the following expressions.

e1:=2&DifferentialD;&DifferentialD;xFxFx&equals;&psi;1x

e1:=2&DifferentialD;&DifferentialD;xFxFx&equals;&psi;1x

(6)

e2:=1&DifferentialD;2&DifferentialD;x2FxFx&equals;&psi;2x

e2:=&DifferentialD;2&DifferentialD;x2FxFx&equals;&psi;2x

(7)

By solving e1, we get Fx as:

ans_F:=dsolvee1

ans_F:=Fx&equals;_C1&ExponentialE;&int;12&psi;1x&DifferentialD;x

(8)

and by substituting this result into the second one we get the pattern identifying the ODE family.

ode_pattern:=expandevale2&comma;ans_F

ode_pattern:=12&DifferentialD;&DifferentialD;x&psi;1x14&psi;1x2&equals;&psi;2x

(9)

2. Consider the second order ODE family having an integrating factor μ=Fx -- an arbitrary function -- and the symmetry ξ=0&comma;η=Fx (see infgen), and such that the reduced ODE is the most general first order linear ODE:

ode_1:=&DifferentialD;&DifferentialD;xyx&equals;Axyx&plus;Bx

ode_1:=&DifferentialD;&DifferentialD;xyx&equals;Axyx&plus;Bx

(10)

where A(x) and B(x) are arbitrary functions. To start with, we obtain the second order ODE mentioned above which has the integrating factor Fx as in example 1.

ode_2:=redodeFx&comma;yx&comma;ode_1

ode_2:=&DifferentialD;2&DifferentialD;x2yx&equals;Ax&DifferentialD;&DifferentialD;xyxFx&plus;Axyx&DifferentialD;&DifferentialD;xFx&plus;&DifferentialD;&DifferentialD;xAxyxFx&plus;Bx&DifferentialD;&DifferentialD;xFx&plus;&DifferentialD;&DifferentialD;xBxFx&DifferentialD;&DifferentialD;xFx&DifferentialD;&DifferentialD;xyxFx

(11)

In the previous step, ode_2 is in fact the most general second order linear ODE we can imagine. If we now impose the symmetry condition Xode_2&equals;0 (see odepde), where X&equals;0&comma;Fx, we arrive at the following restriction on Ax.

constraint:=expandodepdeode_2&comma;yx&comma;0&comma;Fx

constraint:=2Ax&DifferentialD;&DifferentialD;xFx&DifferentialD;&DifferentialD;xAxFx&plus;&DifferentialD;&DifferentialD;xFx2Fx&plus;&DifferentialD;2&DifferentialD;x2Fx

(12)

We can solve this ODE for Ax.

ans_A:=dsolveconstraint&comma;Ax

ans_A:=Ax&equals;Fx&DifferentialD;&DifferentialD;xFx&plus;_C1Fx2

(13)

Introducing the result into ode_2 and disregarding the nonhomogeneous term (irrelevant in the solving scheme), we obtain the homogeneous ODE family pattern:

ode_2_H:=&DifferentialD;2&DifferentialD;x2yx&equals;collectselecthas&comma;expandrhsode_2&comma;y&comma;yx&comma;&DifferentialD;&DifferentialD;xyx

ode_2_H:=&DifferentialD;2&DifferentialD;x2yx&equals;Ax&DifferentialD;&DifferentialD;xFxFx&DifferentialD;&DifferentialD;xyx&plus;Ax&DifferentialD;&DifferentialD;xFxFx&plus;&DifferentialD;&DifferentialD;xAxyx

(14)

ode_pattern:=subs_C1&equals;1&comma;collectexpandsubsans_A&comma;ode_2_H&comma;y

ode_pattern:=&DifferentialD;2&DifferentialD;x2yx&equals;&DifferentialD;&DifferentialD;xFxFx3&plus;&DifferentialD;2&DifferentialD;x2FxFxyx&plus;&DifferentialD;&DifferentialD;xyxFx2

(15)

where we have set the constant _C1 equal to 1 without loss of generality (μ=Fx is defined up to a constant factor). To set up a solving scheme here, check if the coefficient of dydx in a given ODE is related to the coefficient of y as in the equation above, in which case the integrating factor is just Fx.

3. Consider the following three integrating factors of some (unknown) third order ODE:

&mu;:=yxx&comma;x14&plus;129124yx&comma;x14129124yx

&mu;:=yxx&comma;x14&plus;1429yx&comma;x141429yx

(16)

The most general third order ODE simultaneously admitting all these integrating factors depends on an arbitrary function of _F1x and can be computed as follows.

redode&mu;&comma;3&comma;yx

&DifferentialD;3&DifferentialD;x3yx&equals;138yxx3&plus;52&DifferentialD;2&DifferentialD;x2yxx134&DifferentialD;&DifferentialD;xyxx2&plus;3&DifferentialD;&DifferentialD;xyx&DifferentialD;2&DifferentialD;x2yx&plus;52&DifferentialD;&DifferentialD;xyx2x&plus;_F1xyx

(17)

mapDEtools&lsqb;mutest&rsqb;&comma;&mu;&comma;

0&comma;0&comma;0

(18)

The test performed above verifies - returning zero - that the given integrating factors (the list mu) indeed turn the returned ODE exact - see mutest.

The number of given integrating factors can be different than the differential order. For example, taking just the first two integrating factors of the previous example, the most general third order ODE turned exact by them actually depends - in addition - on an arbitrary function of two variables, _F3x&comma;y.

&mu;1..2

yxx&comma;x14&plus;1429yx

(19)

redode&comma;3&comma;yx

&DifferentialD;3&DifferentialD;x3yx&equals;D2_F3x&comma;yx&DifferentialD;&DifferentialD;xyx2&plus;182x_F3x&comma;yx29&plus;16D1_F3x&comma;yxx22x_F3x&comma;yx529&plus;19&DifferentialD;&DifferentialD;xyxx2&plus;_F3x&comma;yx&DifferentialD;2&DifferentialD;x2yx&plus;_F3x&comma;yx&DifferentialD;&DifferentialD;xyx23&DifferentialD;&DifferentialD;xyx&DifferentialD;2&DifferentialD;x2yx&plus;&int;yx18229x_F3x&comma;_ax2&plus;8x32x2_F3x&comma;_a2x_F3x&comma;_ax2&plus;52919_ax3&DifferentialD;_a&plus;_F4xyx

(20)

mapDEtools&lsqb;mutest&rsqb;&comma;&mu;1..2&comma;

0&comma;0

(21)

See Also

DEtools, DEtools[equinv], DEtools[intfactor], DEtools[odeadvisor], DEtools[symgen], DEtools[symtest], dsolve, odetest, PDEtools, PDEtools[dchange]


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