finding another particular solution of linear ODE


Consider the homogeneousPlanetmathPlanetmathPlanetmathPlanetmath (http://planetmath.org/HomogeneousLinearDifferentialEquation) second-order linear ordinary differential equationMathworldPlanetmath

y′′+P⁢(x)⁢y′+Q⁢(x)⁢y= 0. (1)

If one knows one particular solution (http://planetmath.org/SolutionsOfOrdinaryDifferentialEquation)  y=y1⁢(x)≢0  of (1), it’s possible to derive from it via two quadratures another solution  y2⁢(x),  linearly independentMathworldPlanetmath on  y1⁢(x);  thus one can write the general solution

y=C1⁢y1⁢(x)+C2⁢y2⁢(x)

of that homogeneous differential equation.

We will now show the derivation procedure.

We put

y=u⁢v (2)

which renders (1) to

(v′′+P⁢v′+Q⁢v)⁢u+(2⁢v′+P⁢v)⁢u′+u′′⁢v= 0. (3)

Here one can choose  v:=y1⁢(x), whence the first addend vanishes, and (3) gets the form

(2⁢y1′+P⁢y1)⁢u′+y1⁢u′′= 0. (4)

This equation may be written as  u′′u′=-2⁢y1′y1-P, which is integrated to

ln⁡|d⁢ud⁢x|=ln⁡1y12-∫P⁢𝑑x+constant,

i.e.

d⁢ud⁢x=Cy12⁢e-∫P⁢𝑑x.

A new integration results from this the general solution of (4):

u=C⁢∫e-∫P⁢𝑑xy12⁢𝑑x+C′.

Thus by (2), we have obtained the wanted other solution

y2⁢(x)=y1⁢(x)⁢∫e-∫P⁢𝑑xy12⁢𝑑x

which is clearly linearly independent on y_1(x).

Consequently, we can express the general solution of the differential equation (1) as

y=y1⁢(x)⁢u=C1⁢y1⁢(x)+C2⁢y1⁢(x)⁢∫e-∫P⁢𝑑xy12⁢𝑑x,

where C1 and C2 are arbitrary constants.

Remark.  The substitution

y:=e-12⁢∫P⁢(x)⁢𝑑x⁢u

converts the equation (1) into the form

d2⁢ud⁢x2+(Q-P24-P′2)⁢u= 0

not containing the derivativePlanetmathPlanetmath d⁢ud⁢x.

References

  • 1 Ernst Lindelöf: Differentiali- ja integralilasku ja sen sovellutukset III.1.  Mercatorin Kirjapaino Osakeyhtiö, Helsinki (1935).
Title finding another particular solution of linear ODE
Canonical name FindingAnotherParticularSolutionOfLinearODE
Date of creation 2014-02-28 14:31:42
Last modified on 2014-02-28 14:31:42
Owner pahio (2872)
Last modified by pahio (2872)
Numerical id 10
Author pahio (2872)
Entry type AlgorithmMathworldPlanetmath
Classification msc 34A05