Suppose
an
+ ⋯ + a1 dny dxn
+ a0 y ...= ...f ( x ) dy dx
is an nth order linear differential equation (l.d.e.). The equation is homogenous ...if ...f ( x ) = 0 , otherwise it is inhomogenous.
a
+ b d2y dx2
+ c y ...= ...0 ..................a, b, c ...constants. dy dx
Try ...y = emx ...then ...a m2 emx + b m emx + c emx ...= ...0 ......⇒ ......a m2 + b m + c ...= ...0
m ...= ...–b ± √ b² – 4ac 2a
In general there will be two roots for m , ...p and q say. The general solution of the equation is ...y ...= ...A e px + B e qx
The equation, ...a m2 + b m + c ...= ...0 , ...is called the auxiliary equation or indicial equation.
Example 1
– d2y dx2
– 2 y ...= ...0 dy dx
y = emx ...is a solution if ...m2 – m – 2 ...= ...0 ...⇔ ...( m + 1 ) ( m – 2 ) ...= ...0
i.e. the general soultion is ...y ...= ...A e –x + B e 2x
Example 2
– k2 y ...= ...0 d2y dx2
y = emx ...is a solution if ...m2 – k2 ...= ...0 ...⇔ ...m ...= ...± k .
The general soultion is ...y ...= ...A e k x + B e – k x
Example 3
+ k2 y ...= ...0 d2y dx2
y = emx ...is a solution if ...m2 + k2 ...= ...0 ...⇔ ...m ...= ...± i k . ...The general soultion is
y ...= ...A e i k x + B e – i k x ...= ...( A + B ) cos k x ...+ ...i ( A – B ) sin k x
Note that this is a more general version of the solution , we found earlier, (where k = 1). If the initial conditions require that the solution be real, then the general (real) solution is :
y ...= ...C cos k x ...+ ...D sin k x , ......C, D ∊ IR
where ...A = C/2 – i D/2 ...and ...B = C/2 + i D/2 . ...Of course this is an additional constraint on the constants.
If ...b² = 4ac ...then the two roots will be equal: ...p = q.
So the 'general' solution becomes ...y ...= ...A e px + B e px ...= ...C e px , ...i.e. only one constant.
But in this case ...y = x e px ...is also a solution, since ...p = –b ⁄ 2a , ...and
y ...= ...x e px ......⇒ ......y' ...= ...x pe px + e px = ( p x + 1 ) e px
⇒ ............y'' ...= ...p e px + ( p x + 1 ) p e px ...= ...( p2 x + 2 p ) e px
So substituting these values in the differential equation:
a y'' + b y' + c y ...= ...( ( a p2 + b p + c ) x + 2 a p + b ) e px
but ...a p2 + b p + c ...= ...0 , ...and ...2 a p ...= ...2 a × –b ⁄ 2 a ...= ...-b , ...proving that ...y = x e px ...is a solution.
The general solution is therefore: .........y ...= ...A e px + B x e px
Example 4
+ 2 d2y dx2
+ y ...= ...0 dy dx
The auxiliary equation is ...m2 + 2m + 1 ...= ...0 , ...i.e. ...( m – 1 )2 ...= ...0 ......⇒ ......p = q = 1
The general solution is therefore: .........y ...= ...A e x + B x e x
By dividing through by a, a second order equation can be written in the form
+ 2d d2y dx2
+ k y ...= ...0 dy dx
If ...k > 0, ...and ...d 2 < k ...then the roots of the auxiliary equation are imaginary, and the general solution is
y ...= ...A exp { ( –d + i √ k – d ² ...) x } ...+ ...B exp { ( –d – i √ k – d ² ...) x }
......= ...e–dx [ A exp { i x √ k – d ² } ...+ ...B exp { – i x √ k – d ² } ]
......= ...e–dx [ ( A + B ) cos ( x √ k – d ² ) ...+ ...( A – B ) i sin ( x √ k – d ² ) ]
In this case the general real valued solution is
y ...= ...e–dx [ C cos ( x √ k – d ² ) ...+ ...D sin ( x √ k – d ² ) ]
......= ...e–dx ...E ...sin ( x √ k – d ² ...+ ...φ ) , ......where ...E = √ C ² + D ² , ...φ = sin–1 C / E
If ...k > 0, ...and ...d 2 > k ...then the roots of the auxiliary equation are real, and the general solution is
y ...= ...A exp { ( –d + √ d ² – k ...) x } ...+ ...B exp { ( –d – √ d ² – k ...) x }
......= ...e–dx [ A exp { x √ d ² – k } ...+ ...B exp { – x √ d ² – k } ]
If ...k > 0, ...and ...d 2 = k ...then the roots are concurrent, and the general solution is
y ...= ...e–dx [ A + Bx ]
We can generalize the above by stating that the equation
an
+ ⋯ + a1 dny dxn
+ a0 y ...= ...0 dy dx
will have a general solution ......A1 exp ( m1 x ) + ⋯ + An exp ( mn x ) , ......where ...m1, ⋯ mn ...are the separate roots of the auxiliary equation
a1 m12 + ⋯ + an mn2 ...= ...0
If this equation does not have n separate roots, say for example ...mi = mj , ...then the general solution becomes:
A1 exp ( m1 x ) + ⋯ + Ai exp ( mi x ) + ⋯ + Aj x exp ( mi x ) + ⋯ + An exp ( mn x )
and if ...mi = mj = mk :
A1 exp ( m1 x ) + ⋯ + Ai exp ( mi x ) + ⋯ + Aj x exp ( mi x ) + ⋯ + Ak x2 exp ( mi x ) + ⋯ + An exp ( mn x )
etcetera.