Where does the general solution for second order linear ODEs come from?

Join the discussion
Ask a follow-up here, or get your own question answered by working scientists, mathematicians and engineers — people, not an autocomplete.
Real named experts · corrections over time · the nuance an AI answer skips
2 replies · 2K views
sandy.bridge
Messages
797
Reaction score
1
If
[tex]ay+b\int^y_0ydy+cy'=0[/tex]
then
[tex]ay'+by+cy''=0[/tex]
now, let
[tex]y=e^{sx}[/tex]
thus,
[tex]s^2+a/cs+b/c=0[/tex]
and then one solves for s. It is then plugged into what sources are deeming a "general solution"
[tex]y=C_1e^{s_1x}+C_2e^{s_2x}[/tex]
however, none of these texbooks explain or derive where this comes from, and I have not taken any linear algebra or differential equations. Can someone explain where [tex]y=C_1e^{s_1x}+C_2e^{s_2x}[/tex] comes from?
Thanks!
 
Physics news on Phys.org
sandy.bridge said:
If
[tex]ay+b\int^y_0ydy+cy'=0[/tex]
then
[tex]ay'+by+cy''=0[/tex]
now, let
[tex]y=e^{sx}[/tex]
thus,
[tex]s^2+a/cs+b/c=0[/tex]
and then one solves for s. It is then plugged into what sources are deeming a "general solution"
[tex]y=C_1e^{s_1x}+C_2e^{s_2x}[/tex]
however, none of these texbooks explain or derive where this comes from, and I have not taken any linear algebra or differential equations. Can someone explain where [tex]y=C_1e^{s_1x}+C_2e^{s_2x}[/tex] comes from?
Thanks!

Hey sandy.bridge.

Lets take it from the step if y = e^mx, y' = me^mx, y'' = m^2 e^mx then you end up getting an equation of the form:

(am^2 + bm + c) x e^(mx) = 0. Now divide by e^(mx) which gives you the quadratic term:

am^2 + bm + c = 0.

Now if you have two solutions for m, what does that imply for the final solution?

The real 'guts' of proving this general property for linear ODE's or systems of ODE's is based on the superposition principle. The idea uses linear algebra, in particular spectral theory to show the solution in terms of superpositions. If you want to prove it yourself, you need to setup the ODE in terms of a matrix system, and then do some eigen-analysis to decompose the matrix.

I didn't have to do any proving when I took my prior DE course, but if you want to understand why you get this kind of formula in a general sense, look at the superposition principle with reference to ODE's and a detailed proof of it.

Also I didn't link to Wikipedia since they are having a blackout for SOPA awareness.
 
Note that while "trying" [itex]e^{rx}[/itex] will reduce a linear differential equation with constant coefficients to its "characteristic equation", a polynomial, that does not necessarily mean the solution will involve exponentials.

If the solution to the characteristic equation is a complex number, a+ bi, then (assuming all coefficients are real) a- bi is also a root and [itex]e^{(a+ bi)x}= e^{ax}\left(cos(bx)+ sin(bx)\right)[/itex]. Also, if r is a multiple root (If the differential equation is y'''+ 3y''+ 3y'+ y= 0, the characteristic equation is [itex]r^3+ 3r^2+ 3r+ 1= (r+1)^3= 0[/itex] which has r= -1 as a triple root. The general solution to the differential equation is [tex]y= Ae^{-x}+ Bxe^{-x}+ Dx^2e^{-x}[/tex].