The general form of
2
nd
{\displaystyle 2^{\text{nd}}}
order equation is
y
″
=
f
(
t
,
y
,
y
′
)
.
{\displaystyle {\begin{aligned}y''=f(t,y,y').\end{aligned}}}
We call them
linear non-homogeneous if the equation can be written in the form
y
″
+
p
(
t
)
y
′
+
q
(
t
)
y
=
g
(
t
)
{\displaystyle y''+p(t)y'+q(t)y=g(t)}
and
linear homogeneous if, in addition to being linear non-homogeneous,
g
(
t
)
=
0
{\displaystyle g(t)=0}
y
″
+
p
(
t
)
y
′
+
q
(
t
)
y
=
0.
{\displaystyle y''+p(t)y'+q(t)y=0.}
The method of characteristic equations is for homogeneous equations and the methods of undetermined coefficients and of variation of parameters for homogeneous equations.
If the equation is linear homogeneous and further
p
(
t
)
,
q
(
t
)
{\displaystyle p(t),q(t)}
are constant, then the equation is referred to as a constant-coefficients equation:
a
y
″
+
b
y
′
+
c
y
=
0
{\displaystyle ay''+by'+cy=0}
and we can apply the method of characteristic equations to solve such an equation. Note that
a
{\displaystyle a}
is assumed to be non-zero since we are working with a second order equation.
We assume that the solution is of the form
y
(
t
)
=
e
r
t
{\displaystyle y(t)=e^{rt}}
(this is called making an ansatz). This gives
(
a
r
2
+
b
r
+
c
)
e
r
t
=
0
⟹
a
r
2
+
b
r
+
c
=
0
,
{\displaystyle (ar^{2}+br+c)e^{rt}=0\Longrightarrow \,ar^{2}+br+c=0,}
which equation is called the characteristic equation .
So to solve the above ODE, it suffices to find the two roots
r
1
,
r
2
{\displaystyle r_{1},r_{2}}
.
Then the general solution is of the form:
y
(
t
)
=
c
1
e
r
1
t
+
c
2
e
r
2
t
.
{\displaystyle y(t)=c_{1}e^{r_{1}t}+c_{2}e^{r_{2}t}.}
Consider a mass
m
{\displaystyle m}
hanging at rest on the end of a vertical spring of length
l
{\displaystyle l}
, spring constant
k
{\displaystyle k}
and damping constant
γ
{\displaystyle \gamma }
.
Let
u
(
t
)
{\displaystyle u\left(t\right)}
denote the displacement, in units of feet, from the equilibrium position. Note that since
u
(
t
)
{\displaystyle u(t)}
represents the amount of displacement from the spring's equilibrium position (the position obtained when the downward force of gravity is matched by the will of the spring to not allow the mass to stretch the spring further) then
u
(
t
)
{\displaystyle u(t)}
should increase downward. Then by Newton's Third Law one can obtain the equation
m
u
″
(
t
)
+
γ
u
′
(
t
)
+
k
u
(
t
)
=
F
(
t
)
,
{\displaystyle mu''(t)+\gamma u'(t)+ku(t)=F(t),}
where
F
(
t
)
{\displaystyle F\left(t\right)}
is any external force, which for simplicity we will assume to be zero.
First we obtain the characteristic equation:
m
r
2
+
γ
r
+
k
=
0.
{\displaystyle mr^{2}+\gamma r+k=0.}
Suppose that
m
=
1
lb
,
γ
=
5
lb
/
ft
/
s
{\displaystyle m=1{\text{lb}},\gamma =5{\text{lb}}/{\text{ft}}/{\text{s}}}
and
k
=
6
lb
/
ft
{\displaystyle k=6{\text{lb}}/{\text{ft}}}
then we obtain the roots
r
1
=
−
2
{\displaystyle r_{1}=-2}
,
r
2
=
−
3
{\displaystyle r_{2}=-3}
.
Therefore, the general solution will be
u
(
t
)
=
c
1
e
−
2
t
+
c
2
e
−
3
t
.
{\displaystyle u(t)=c_{1}e^{-2t}+c_{2}e^{-3t}.}
Further if
u
(
0
)
=
0
,
u
′
(
0
)
=
1
{\displaystyle u(0)=0,\,u'(0)=1}
we obtain
c
1
=
1
,
c
2
=
−
1
{\displaystyle c_{1}=1,c_{2}=-1}
:
u
(
t
)
=
e
−
2
t
−
e
−
3
t
.
{\displaystyle u(t)=e^{-2t}-e^{-3t}.}
4
y
″
−
y
=
0
,
y
(
−
2
)
=
1
,
y
′
(
−
2
)
=
−
1.
{\displaystyle 4y''-y=0,\,\,y(-2)=1,\,\,y'(-2)=-1.}
We obtain the characteristic equation
4
r
2
−
1
=
0
⇒
r
=
±
1
2
{\displaystyle 4r^{2}-1=0\Rightarrow r=\pm {\frac {1}{2}}}
and so the general solution will be
y
(
t
)
=
c
1
e
t
2
+
c
2
e
−
t
2
.
{\displaystyle y(t)=c_{1}e^{\frac {t}{2}}+c_{2}e^{-{\frac {t}{2}}}.}
Using the initial conditions we obtain:
1
=
c
1
e
−
1
+
c
2
e
and
−
1
=
1
2
(
c
1
e
−
1
−
c
2
e
)
.
{\displaystyle 1=c_{1}e^{-1}+c_{2}e\,{\text{ and }}\,-1={\frac {1}{2}}(c_{1}e^{-1}-c_{2}e).}
Solving these two equations gives:
c
1
=
−
1
2
e
,
c
2
=
3
2
e
−
1
{\displaystyle c_{1}={\frac {-1}{2}}e,c_{2}={\frac {3}{2}}e^{-1}}
and so the solution for our IVP is:
y
(
t
)
=
−
1
2
e
1
+
t
2
+
3
2
e
−
t
2
−
1
.
{\displaystyle y\left(t\right)=-{\frac {1}{2}}e^{1+{\frac {t}{2}}}+{\frac {3}{2}}e^{-{\frac {t}{2}}-1}.}
Therefore, as
t
→
+
∞
{\displaystyle t\to +\infty }
we obtain
y
→
−
∞
{\displaystyle y\to -\infty }
.
y
″
+
5
y
′
+
6
y
=
0
,
y
(
0
)
=
2
,
y
′
(
0
)
=
β
{\displaystyle y''+5y'+6y=0,\,\,y\left(0\right)=2,\,\,y'\left(0\right)=\beta }
The characteristic equation is
r
2
+
5
r
+
6
=
0
⇒
r
=
−
2
,
−
3
{\displaystyle r^{2}+5r+6=0\Rightarrow r=-2,-3}
and so the general solution will be:
y
(
t
)
=
c
1
e
−
2
t
+
c
2
e
−
3
t
{\displaystyle y(t)=c_{1}e^{-2t}+c_{2}e^{-3t}}
Using the initial conditions we obtain:
2
=
c
1
+
c
2
and
β
=
−
2
c
1
−
3
c
2
.
{\displaystyle 2=c_{1}+c_{2}\,{\text{ and }}\,\beta =-2c_{1}-3c_{2}.}
Solving these two equations gives:
c
1
=
(
6
+
β
)
,
c
2
=
−
(
4
+
β
)
{\displaystyle c_{1}=(6+\beta ),\,c_{2}=-\left(4+\beta \right)}
and so the solution for our IVP is:
y
(
t
)
=
(
6
+
β
)
e
−
2
t
−
(
4
+
β
)
e
−
3
t
.
{\displaystyle y(t)=(6+\beta )e^{-2t}-(4+\beta )e^{-3t}.}
Therefore, as
t
→
+
∞
{\displaystyle t\to +\infty }
we obtain
y
→
0
{\displaystyle y\to 0}
.