6.3: Series Solutions and Convergence
- Page ID
- 392
\( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)
\( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)
\( \newcommand{\id}{\mathrm{id}}\) \( \newcommand{\Span}{\mathrm{span}}\)
( \newcommand{\kernel}{\mathrm{null}\,}\) \( \newcommand{\range}{\mathrm{range}\,}\)
\( \newcommand{\RealPart}{\mathrm{Re}}\) \( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)
\( \newcommand{\Argument}{\mathrm{Arg}}\) \( \newcommand{\norm}[1]{\| #1 \|}\)
\( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)
\( \newcommand{\Span}{\mathrm{span}}\)
\( \newcommand{\id}{\mathrm{id}}\)
\( \newcommand{\Span}{\mathrm{span}}\)
\( \newcommand{\kernel}{\mathrm{null}\,}\)
\( \newcommand{\range}{\mathrm{range}\,}\)
\( \newcommand{\RealPart}{\mathrm{Re}}\)
\( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)
\( \newcommand{\Argument}{\mathrm{Arg}}\)
\( \newcommand{\norm}[1]{\| #1 \|}\)
\( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)
\( \newcommand{\Span}{\mathrm{span}}\) \( \newcommand{\AA}{\unicode[.8,0]{x212B}}\)
\( \newcommand{\vectorA}[1]{\vec{#1}} % arrow\)
\( \newcommand{\vectorAt}[1]{\vec{\text{#1}}} % arrow\)
\( \newcommand{\vectorB}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)
\( \newcommand{\vectorC}[1]{\textbf{#1}} \)
\( \newcommand{\vectorD}[1]{\overrightarrow{#1}} \)
\( \newcommand{\vectorDt}[1]{\overrightarrow{\text{#1}}} \)
\( \newcommand{\vectE}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash{\mathbf {#1}}}} \)
\( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)
\( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)
\(\newcommand{\avec}{\mathbf a}\) \(\newcommand{\bvec}{\mathbf b}\) \(\newcommand{\cvec}{\mathbf c}\) \(\newcommand{\dvec}{\mathbf d}\) \(\newcommand{\dtil}{\widetilde{\mathbf d}}\) \(\newcommand{\evec}{\mathbf e}\) \(\newcommand{\fvec}{\mathbf f}\) \(\newcommand{\nvec}{\mathbf n}\) \(\newcommand{\pvec}{\mathbf p}\) \(\newcommand{\qvec}{\mathbf q}\) \(\newcommand{\svec}{\mathbf s}\) \(\newcommand{\tvec}{\mathbf t}\) \(\newcommand{\uvec}{\mathbf u}\) \(\newcommand{\vvec}{\mathbf v}\) \(\newcommand{\wvec}{\mathbf w}\) \(\newcommand{\xvec}{\mathbf x}\) \(\newcommand{\yvec}{\mathbf y}\) \(\newcommand{\zvec}{\mathbf z}\) \(\newcommand{\rvec}{\mathbf r}\) \(\newcommand{\mvec}{\mathbf m}\) \(\newcommand{\zerovec}{\mathbf 0}\) \(\newcommand{\onevec}{\mathbf 1}\) \(\newcommand{\real}{\mathbb R}\) \(\newcommand{\twovec}[2]{\left[\begin{array}{r}#1 \\ #2 \end{array}\right]}\) \(\newcommand{\ctwovec}[2]{\left[\begin{array}{c}#1 \\ #2 \end{array}\right]}\) \(\newcommand{\threevec}[3]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \end{array}\right]}\) \(\newcommand{\cthreevec}[3]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \end{array}\right]}\) \(\newcommand{\fourvec}[4]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \\ #4 \end{array}\right]}\) \(\newcommand{\cfourvec}[4]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \\ #4 \end{array}\right]}\) \(\newcommand{\fivevec}[5]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \\ #4 \\ #5 \\ \end{array}\right]}\) \(\newcommand{\cfivevec}[5]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \\ #4 \\ #5 \\ \end{array}\right]}\) \(\newcommand{\mattwo}[4]{\left[\begin{array}{rr}#1 \amp #2 \\ #3 \amp #4 \\ \end{array}\right]}\) \(\newcommand{\laspan}[1]{\text{Span}\{#1\}}\) \(\newcommand{\bcal}{\cal B}\) \(\newcommand{\ccal}{\cal C}\) \(\newcommand{\scal}{\cal S}\) \(\newcommand{\wcal}{\cal W}\) \(\newcommand{\ecal}{\cal E}\) \(\newcommand{\coords}[2]{\left\{#1\right\}_{#2}}\) \(\newcommand{\gray}[1]{\color{gray}{#1}}\) \(\newcommand{\lgray}[1]{\color{lightgray}{#1}}\) \(\newcommand{\rank}{\operatorname{rank}}\) \(\newcommand{\row}{\text{Row}}\) \(\newcommand{\col}{\text{Col}}\) \(\renewcommand{\row}{\text{Row}}\) \(\newcommand{\nul}{\text{Nul}}\) \(\newcommand{\var}{\text{Var}}\) \(\newcommand{\corr}{\text{corr}}\) \(\newcommand{\len}[1]{\left|#1\right|}\) \(\newcommand{\bbar}{\overline{\bvec}}\) \(\newcommand{\bhat}{\widehat{\bvec}}\) \(\newcommand{\bperp}{\bvec^\perp}\) \(\newcommand{\xhat}{\widehat{\xvec}}\) \(\newcommand{\vhat}{\widehat{\vvec}}\) \(\newcommand{\uhat}{\widehat{\uvec}}\) \(\newcommand{\what}{\widehat{\wvec}}\) \(\newcommand{\Sighat}{\widehat{\Sigma}}\) \(\newcommand{\lt}{<}\) \(\newcommand{\gt}{>}\) \(\newcommand{\amp}{&}\) \(\definecolor{fillinmathshade}{gray}{0.9}\)In the last section, we saw how to find series solutions to second order linear differential equations. We did not investigate the convergence of these series. In this discussion, we will derive an alternate method to find series solutions. We will also learn how to determine the radius of convergence of the solutions just by taking a quick glance of the differential equation.
Consider the differential equation
\[ y'' + y' + ty = 0. \nonumber\]
As before we seek a series solution
\[ y = a_0 + a_1t + a_2t^2 + a_3t^3 + a_4t^4 + ... \;. \nonumber\]
The theory for Taylor Series states that
\[ n!\; a_n = y^{(n)}(0). \nonumber \]
We have
\[ y'' = -y' -ty. \nonumber \]
Plugging in 0 gives
\[ 2!\, a_2 = y''(0) = -y'(0) + 0 = -a_1 \nonumber \]
\[ a_2 = -\dfrac{a_1}{2}. \nonumber\]
Taking the derivative of the differential equation gives
\[ (y'' + y' + ty)' = y''' + y'' + ty' + y = 0 \nonumber \]
or
\[ y''' = -y'' - ty' - y. \nonumber\]
Plugging in zero gives
\[ 3!\, a_3 = a_1 - a_0 \nonumber\]
\[ a_3 = \dfrac{a_1}{6} - \dfrac{a_0}{6}. \nonumber\]
Taking another derivative gives
\[ (y''' + y'' + ty' + y)' = y^{(iv)} + y''' + ty'' + 2y' = 0 \nonumber \]
or
\[ y^{(iv)} = -y''' - ty'' - 2y' . \nonumber\]
Plugging in zero gives
\[ 4! \,a_4 = -a_1 + a_0 - 2a_1 \nonumber \]
\[ a_4 = -\dfrac{49}{24} a_1 + \dfrac{a_0}{24}. \nonumber\]
The important thing to note here is that all of the coefficients can be written in terms of the first two. To come up with a theorem regarding this, we first need a definition.
A function \(f(x)\) is called analytic at \(x_0\) if \(f(x)\) is equal to its power series.
It turns out that if \(p(x)\) and \(q(x)\) are analytic then there always exists a power series solution to the corresponding differential equation. We state this fact below without proof. If \(x_0\) is a point such that \(p(x)\) and \(p(x)\) are analytic, then \(x_0\) is called an ordinary point of the differential equation.
Let \(x_0\) be an ordinary point of the differential equation
\[ L(y) = y'' + p(t)y' + q(t)y = 0. \]
Then the general solution can be represented by the power series
\[ y= \sum_{n=0}^\infty a_n(x-x_0)^n = a_0\,y_1(x) + a_1\,y_2(x). \]
where \(a_0\) and \(a_1\) are arbitrary constants and \(y_1\) and \(y_2\) are analytic at \(x_0\). The radii of convergence for \(y_1\) and \(y_2\) are at least as large as the minimum radii of convergences for \(p\) and \(q\).
Remark: The easiest way of find the radii of convergence of most functions us by using the following fact
If \(f(x)\) is an analytic function for all \(x\), then the radius of convergence for \(1/f(x)\) is the distance from the center of convergence to the closest root (possibly complex) of \(f(x)\).
Find a lower bound for the radius of convergence of series solutions about \(x = 1\) for the differential equation
\[ (x^2 + 4)\; y'' + \text{sin} \; (x)y' + e^xy = 0. \nonumber \]
Solution
We have
\[ p(x) = \dfrac{\sin x}{x^2 + 4} \nonumber\]
\[ q(x) = \dfrac{e^x}{x^2 + 4} . \nonumber\]
Both of these are quotient of analytic functions. The roots of \(x^2 + 4\) are
\[2i \;\;\; \text{and} \;\;\; -2i. \nonumber \]
The distance from \(1\) to \(2i\) is the same as the distance from \((1,0)\) to \((0,2)\) which is \( \sqrt{5} \). We get the same distance from \(1\) to \(-2i\). Hence the radii of convergence of the solutions are both at least \( \sqrt{5} \).
Contributors and Attributions
- Larry Green (Lake Tahoe Community College)
Integrated by Justin Marshall.