Solving Linear ODEs Which Are Fully Exact - Maple Help

Online Help

All Products    Maple    MapleSim


Home : Support : Online Help : Mathematics : Differential Equations : Classifying ODEs : First Order : odeadvisor/fully_exact_linear

Solving Linear ODEs Which Are Fully Exact

 

Description

Examples

Description

• 

A fully exact nth-order linear ODE is one which can be obtained by directly differentiating a first order linear ODE n-1 times; for example, a fourth order ODE of this type is of the form:

fully_exact_linear_4th_order_ODE := Diff( A(x)*diff(y(x),x) + B(x)*y(x) + C(x) = 0, x$3);

fully_exact_linear_4th_order_ODE:=ⅆ3ⅆx3Axⅆⅆxyx+Bxyx+Cx=0

(1)
• 

Linear ODEs (LODEs) of this form can be recognized by the fact that the adjoint of their homogeneous part is always a first order ODE in disguised form (see "How it works" below). This type of LODE is systematically solved by dsolve.

Examples

withDEtools,adjoint,firint,odeadvisor

adjoint,firint,odeadvisor

(2)

PDEtools[declare]yx,prime=x

yxwill now be displayed asy

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

(3)

Consider the following 4th order ODE:

ODE4x2+ⅇxxⅆ4ⅆx4yx+4ⅇx+6x+3ⅇxxⅆ3ⅆx3yx+6+3ⅇxx+9ⅇxⅆ2ⅆx2yx+ⅇxx+6ⅇxⅆⅆxyx+ⅇxyx=0

ODE4:=x2+ⅇxxy''''+4ⅇx+6x+3ⅇxxy'''+6+3ⅇxx+9ⅇxy''+ⅇxx+6ⅇxy'+ⅇxy=0

(4)

odeadvisorODE4,yx

_high_order,_fully,_exact,_linear

(5)

Once we know ODE4 is fully exact we can reduce it to a first order linear ODE straightforwardly; this reduction can be performed interactively, for instance, using the firint command as follows:

ODE3firintODE4,yx

ODE3:=ⅇxy+ⅇxx+4ⅇx+2y'+2ⅇxx+3ⅇx+4xy''+x2+ⅇxxy'''+_C1=0

(6)

ODE2firintODE3,yx

ODE2:=ⅇxy+ⅇxx+2ⅇx+2xy'+x2+ⅇxxy''+_C1x+_C2=0

(7)

ODE1firintODE2,yx

ODE1:=ⅇxy+x2+ⅇxxy'+12_C1x2+_C2x+_C3=0

(8)

and the first order LODE above - ODE1 - can be solved using dsolve leading to the solution to ODE4

sol1dsolveODE1,yx

sol1:=y=∫12ⅇ∫ⅇxxⅇx+xⅆx_C1x2+2_C2x+2_C3xⅇx+xⅆx+_C4ⅇ∫ⅇxxⅇx+xⅆx

(9)

odetestsol1,ODE4

0

(10)

The solution above is essentially the same one obtains by calling dsolve directly:

dsolveODE4,yx

y=_C4+∫3_C1x2+2_C2x+_C3ⅇ∫ⅇxx2+ⅇxxⅆxx2+ⅇxxⅆxⅇ∫ⅇxx2+ⅇxxⅆx

(11)

How it works

To detect that ODE4 is fully exact, we first look at the equation's adjoint:

ADJadjointODE4,yx

ADJ:=ⅇxx+2xy'''+x2+ⅇxxy''''

(12)

This adjoint contains only third and fourth order derivatives, and hence can be seen as a first order ODE in y'''. This means that the original fourth order ODE4 can be obtained by differentiating 3 times a certain first order linear ODE1 (not the one obtained with firint lines above). In turn, this ODE1 can be obtained directly from the adjoint equation ADJ above by replacing y''' -> y and calculating the adjoint again:

subsⅆ3ⅆx3yx=yx,ADJ

ⅇxx+2xy+x2+ⅇxxy'

(13)

adjoint,yx

ⅇxy+x2+ⅇxxy'

(14)

Actually, differentiating 3 times the equation above one reobtains the original fourth order ODE:

3x3=0

ⅇxy+3y'ⅇx+3y''ⅇx+y'''ⅇx+ⅇxx+3ⅇxy'+32+ⅇxx+2ⅇxy''+32x+ⅇxx+ⅇxy'''+x2+ⅇxxy''''=0

(15)

normalODE4

0=0

(16)

The method just explained works with ODEs of arbitrary order, and - as in the homogeneous case - if the given fully exact LODE is non-homogeneous, then it can be obtained by differentiating a non-homogeneous linear first order ODE, and this fact can be detected in the same way by analyzing the adjoint of the homogeneous part of the given LODE.

One of the interesting properties of this method is that a closed form "solution-formula" can be written directly, for all differential orders and for the "non-homogeneous case" (the homogeneous one is just a particular case) - thus shortcutting all the steps and the need for first solving the homogeneous part to afterwards finding a particular solution, for instance, through variation of parameters.

What is at the base of this method is the fact that integrating factors for a LODE are solutions of its adjoint (see for instance Murphy's book and others in dsolve, references).

See Also

adjoint

dsolve

odeadvisor

firint

firtest

intfactor

mutest

odetest

quadrature

missing

reducible

linear_ODEs

exact_linear

exact_nonlinear

odeadvisor,types

 


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