Solving Linear Ordinary Differential Equations (LODEs) by computing integrating factors for them - Maple Help

Online Help

All Products    Maple    MapleSim


Home : Support : Online Help : Mathematics : Differential Equations : dsolve : dsolve/integrating_factors_for_LODEs

Solving Linear Ordinary Differential Equations (LODEs) by computing integrating factors for them

Description

  

Integrating factors for second and higher order linear ODEs

• 

For linear ODEs (LODEs) of order 2 or greater, it is possible to calculate integrating factors by solving the adjoint of the LODE. This could be as difficult as the original problem, or much easier, depending on the example. This method is implemented in dsolve.

Examples

withDEtools,adjoint,firint,intfactor:

ode:=ⅆ2ⅆx2yx=yxx+yxlnxx2+ⅆⅆxyxlnxx2ⅇxlnxx2yxlnxx2

ode:=ⅆ2ⅆx2yx=yxx+yxlnxx2+ⅆⅆxyxlnxx2ⅇxlnxx2yxlnxx2

(1)

This ODE has the following adjoint

adjointode,yx

ⅇxx1yxlnxx2+ⅇxⅆⅆxyx+ⅆ2ⅆx2yx

(2)

This adjoint equation is in turn solvable by dsolve

dsolve

yx=∫ⅇⅇxlnx2ⅆx_C1+_C2lnx

(3)

Now the solutions to the adjoint equation are integrating factors of the original LODE, so the two independent solutions implied in the general solution above

rhs_C1=0,_C2=1|rhs_C1=0,_C2=1

lnx

(4)

rhs_C1=1,_C2=2|rhs_C1=1,_C2=2

∫ⅇⅇxlnx2ⅆx+2lnx

(5)

are integrating factors of ode. These integrating factors could also be found using the intfactor directly

Μ:=intfactorode

Μ:=lnx,∫ⅇⅇxlnx2ⅆxlnx

(6)

Constructing solutions using integrating factors

How are these integrating factors transformed into a solution to the original problem? By using them to construct two first integrals; that is: two ODEs of lower order (in this case two first order ODEs). For that purpose it is provided the firint command which receives an exact ODE and returns a first integral. The idea is simple: an exact ODE is a total derivative - say dR/dx; firint returns the R + _C1:

firintΜ1ode

lnxx2ⅇx+xyxx2+lnxⅆⅆxyx+_C1=0

(7)

firintΜ2ode

∫ⅇⅇxlnx2ⅆxlnx2ⅇxx2+∫ⅇⅇxlnx2ⅆxlnxx+ⅇⅇxx2yxlnxx2+∫ⅇⅇxlnx2ⅆxlnxⅆⅆxyx+_C1=0

(8)

Eliminating y' from these two first integrals (and replacing _C1 by _C2 in one of them) leads to the solution f(x,y(x),_C1,_C2) = 0 to this ode. So this process could be run interactively, as shown, or in one step

dsolveode

yx=_C1ⅇⅇxlnx+_C2ⅇⅇx∫ⅇⅇxlnx2ⅆxlnx

(9)

See Also

DEtools, dsolve,education, firint, firtest, intfactor, muchange, mutest, PDEtools, redode


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