Solving Second Order ODEs by Computing Integrating Factors for Them - Maple Programming Help

Online Help

All Products    Maple    MapleSim


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

Solving Second Order ODEs by Computing Integrating Factors for Them

 

Description

Examples

Description

  

Integrating factors - depending on two variables - for second order ODEs

• 

For second order ODEs integrating factors depending on three variables, (x,y,y') always exist - their determination however may be as difficult as solving the ODEs themselves. However, integrating factors of the form mu(x,y), mu(x,y'), mu(y,y'), that is, depending on just two variables, when they exist, can be determined systematically (see E.S. Cheb-Terrab and A.D. Roche, "Integrating factors for second order ODEs", Journal of Symbolic Computation V. 27, No. 5 (1999) 501). This method is implemented in dsolve.

• 

The main idea of this method consists of algorithmically determining whether or not a given second order ODE is member of one of the three reducible ODE families respectively admitting integrating factors of the forms mu(x,y), mu(x,y'), mu(y,y') (see redode):

with(DEtools,redode):

PDEtools[declare](y(x), prime=x);

yxwill now be displayed asy

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

(1)

subs(_F1 = G, redode(mu(x,y(x)),y(x),2) );            # ODE family admitting mu(x,y)

y''=D2μx,yy'2+D1μx,yy'+D2Gx,yy'+D1Gx,yμx,y

(2)

convert( subs(diff(y(x),x,x) = `y''`, diff(y(x),x) = `y'`,y(x)=y,(2)), diff);

y''=y'2μy+y'Gy+y'μx+Gxμx,y

(3)

subs(_F1 = G, redode(mu(x,diff(y(x),x)),y(x),2) );    # ODE family admitting mu(x,y')

y''=∫y'μxⅆ_a+D1Gx,y+D2Gx,yy'μx,y'

(4)

convert( subs(diff(y(x),x,x) = `y''`, diff(y(x),x) = `y'`,y(x)=y,(4)), diff);

y''=∫y'μxⅆ_a+Gx+Gyy'μx,y'

(5)

subs(_F1 = G, redode(mu(y(x),diff(y(x),x)),y(x),2) ); # ODE family admitting mu(y,y')

y''=∫y'D1μy,_ay'ⅆ_a+D1Gx,y+D2Gx,yy'μy,y'

(6)

convert( subs(diff(y(x),x,x) = `y''`, diff(y(x),x) = `y'`,y(x)=y,(6)), diff);

y''=∫y'μyy'ⅆ_a+Gx+Gyy'μy,y'

(7)
  

and in doing that, determining mu. Concretely, this means that members of these three ODE families can now be systematically reduced to first order ODEs - for arbitrary mu and G - as if they were simple exact equations (total derivatives). This has enlarged in a noticeable way the decision procedures available for tackling second order ODEs in the nonlinear case.

Examples

Although the relevant thing here is that the new algorithms don't require solving any auxiliary differential equation, it is also worth mentioning that the ODE families above include members not having point symmetries and nonetheless reducible to quadratures; this fact was the main result of an interesting article in 1988 (see A. Gonzalez-Lopez, Phys. Lett. A, (1988) 190). As an ODE example, consider

withDEtools,firint,intfactor,mutest,odeadvisor:

odeⅆ2ⅆx2yxⅆⅆxyx2yx+sinxyxⅆⅆxyx+cosxyx2=0

ode:=y''y'2y+sinxyy'+cosxy2=0

(8)

This ODE has the following integrating factor (see intfactor)

Μintfactorode,yx

Μ:=1y

(9)

This integrating factor can be tested using mutest

mutestΜ,ode

0

(10)

The product Mu*ode is an exact ODE (a total derivative) and hence one can take advantage of this fact to reduce the second order ODE to a first order one, in this case of Bernoulli type (in turn fully solvable in the general case):

first_integralfirintΜode,yx

first_integral:=sinxy+y'y+_C1=0

(11)

odeadvisorfirst_integral,yx

_Bernoulli

(12)

From where the solution to the original ode follows straightforwardly:

dsolvefirst_integral

y=_C12+1ⅇ_C1x_C12_C2+ⅇ_C1x_C2sinx_C1cosx

(13)

odetest,ode

0

(14)

(for more details on how to build solutions from the knowledge of integrating factors see dsolve,education, subsection Constructing solutions using integrating factors.) All these steps are performed automatically by dsolve by calling

dsolveode,yx

It is also possible to invoke the use of "just" this method for reducible ODEs or for each of the three forms of integrating factors via

dsolveode,yx,reducible

dsolveode,yx,mu_xy

dsolveode,yx,mu_y_y1

dsolveode,yx,mu_x_y1

For the specific example under discussion, the first three callings above will solve the problem.

As a concluding remark, an ODE not having point symmetries and having such an integrating factor of the form 1/y is still a simple problem if compared with the three ODE families shown on top, related to arbitrary integrating factors depending on two variables.

That the ode under discussion has no point symmetries can be checked with ease as follows. First we generate the determining system for the infinitesimals of the symmetry generator (see gensys):

sysDEtools[gensys]ode,ξx,y,ηx,y

sys:=y2ξy,yyξy,2ξysinxy3+ηy,yy22ξx,yy2ηyy+ηx,y,3ξycosxy4+ξx,ycosxy3+ξxsinxy3+ηx,ysinxy2+2ηx,yy2ξx,xy22ηxy,ξx,ysinxy4+2ξxcosxy4ηycosxy4+ηxsinxy3+2ηx,ycosxy3+ηx,xy2

(15)

nopssys

4

(16)

Then we simplify this system with respect to its integrability conditions using casesplit (in turn using DifferentialAlgebra or DEtools,Rif):

PDEtools[casesplit]sys

ηx,y=0,ξx,y=0 &where

(17)

verifying that the only solution to this system is the trivial (useless) one. There are in fact infinitely many cases like the one just mentioned, so that the integrating factor and the symmetry approaches end up complementing each other in a quite useful manner.

See Also

DEtools

dsolve/education

firint

firtest

gensys

intfactor

muchange

mutest

PDEtools

PDEtools[casesplit]

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