Recall that we can only find a series solution about \({x_0} = 0\) if this point is an ordinary point, or in other words, if the coefficient of the highest derivative term is not zero at \({x_0} = 0\). That is clearly the case here so let’s start with the form of the solutions as well as the derivatives that we’ll need for this solution.

\[y\left( x \right) = \sum\limits_{n = 0}^\infty {{a_n}{x^n}} \hspace{0.25in}\,\,\,\,\,\,\,\,\,\,\,y'\left( x \right) = \sum\limits_{n = 1}^\infty {n{a_n}{x^{n - 1}}} \hspace{0.25in}y'''\left( x \right) = \sum\limits_{n = 3}^\infty {n\left( {n - 1} \right)\left( {n - 2} \right){a_n}{x^{n - 3}}} \]

Plugging into the differential equation gives,

\[\sum\limits_{n = 3}^\infty {n\left( {n - 1} \right)\left( {n - 2} \right){a_n}{x^{n - 3}}} + {x^2}\sum\limits_{n = 1}^\infty {n{a_n}{x^{n - 1}}} + x\sum\limits_{n = 0}^\infty {{a_n}{x^n}} = 0\]

Now, move all the coefficients into the series and do appropriate shifts so that all the series are in terms of \({x^n}\).

\[\begin{align*}\sum\limits_{n = 3}^\infty {n\left( {n - 1} \right)\left( {n - 2} \right){a_n}{x^{n - 3}}} + \sum\limits_{n = 1}^\infty {n{a_n}{x^{n + 1}}} + \sum\limits_{n = 0}^\infty {{a_n}{x^{n + 1}}} & = 0\\ \sum\limits_{n = 0}^\infty {\left( {n + 3} \right)\left( {n + 2} \right)\left( {n + 1} \right){a_{n + 3}}{x^n}} + \sum\limits_{n = 2}^\infty {\left( {n - 1} \right){a_{n - 1}}{x^n}} + \sum\limits_{n = 1}^\infty {{a_{n - 1}}{x^n}} & = 0\end{align*}\]

Next, let’s notice that we can start the second series at \(n = 1\) since that term will be zero. So let’s do that and then we can combine the second and third terms to get,

\[\begin{align*}\sum\limits_{n = 0}^\infty {\left( {n + 3} \right)\left( {n + 2} \right)\left( {n + 1} \right){a_{n + 3}}{x^n}} + \sum\limits_{n = 1}^\infty {\left[ {\left( {n - 1} \right) + 1} \right]{a_{n - 1}}{x^n}} & = 0\\ \sum\limits_{n = 0}^\infty {\left( {n + 3} \right)\left( {n + 2} \right)\left( {n + 1} \right){a_{n + 3}}{x^n}} + \sum\limits_{n = 1}^\infty {n{a_{n - 1}}{x^n}} & = 0\end{align*}\]

So, we got a nice simplification in the new series that will help with some further simplification. The new second series can now be started at \(n = 0\) and then combined with the first series to get,

\[\sum\limits_{n = 0}^\infty {\left[ {\left( {n + 3} \right)\left( {n + 2} \right)\left( {n + 1} \right){a_{n + 3}} + n{a_{n - 1}}} \right]{x^n}} = 0\]

We can now set the coefficients equal to get a fairly simply recurrence relation.

\[\left( {n + 3} \right)\left( {n + 2} \right)\left( {n + 1} \right){a_{n + 3}} + n{a_{n - 1}} = 0\hspace{0.25in}n = 0,1,2, \ldots \]

Solving the recurrence relation gives,

\[{a_{n + 3}} = \frac{{ - n{a_{n - 1}}}}{{\left( {n + 1} \right)\left( {n + 2} \right)\left( {n + 3} \right)}}\hspace{0.25in}n = 0,1,2, \ldots \]

Now we need to start plugging in values of \(n\) and this will be one of the main areas where we can see a somewhat significant increase in the amount of work required when moving into a higher order differential equation.

\( \displaystyle n = 0:\hspace{0.25in}{a_3} = 0\)

\( \displaystyle n = 1:\hspace{0.25in}{a_4} = \frac{{ - {a_0}}}{{\left( 2 \right)\left( 3 \right)\left( 4 \right)}}\)

\( \displaystyle n = 2:\hspace{0.25in}{a_5} = \frac{{ - 2{a_1}}}{{\left( 3 \right)\left( 4 \right)\left( 5 \right)}}\)

\( \displaystyle n = 3:\hspace{0.25in}{a_6} = \frac{{ - 3{a_2}}}{{\left( 4 \right)\left( 5 \right)\left( 6 \right)}}\)

\( \displaystyle n = 4:\hspace{0.25in}{a_7} = \frac{{ - 4{a_3}}}{{\left( 5 \right)\left( 6 \right)\left( 7 \right)}} = 0\)

\( \displaystyle n = 5:\hspace{0.25in}{a_8} = \frac{{ - 5{a_4}}}{{\left( 6 \right)\left( 7 \right)\left( 8 \right)}} = \frac{{5{a_0}}}{{\left( 2 \right)\left( 3 \right)\left( 4 \right)\left( 6 \right)\left( 7 \right)\left( 8 \right)}}\)

\( \displaystyle n = 6:\hspace{0.25in}{a_9} = \frac{{ - 6{a_5}}}{{\left( 7 \right)\left( 8 \right)\left( 9 \right)}} = \frac{{\left( 2 \right)\left( 6 \right){a_1}}}{{\left( 3 \right)\left( 4 \right)\left( 5 \right)\left( 7 \right)\left( 8 \right)\left( 9 \right)}}\)

\( \displaystyle n = 7:\hspace{0.25in}{a_{10}} = \frac{{ - 7{a_6}}}{{\left( 8 \right)\left( 9 \right)\left( {10} \right)}} = \frac{{\left( 3 \right)\left( 7 \right){a_2}}}{{\left( 4 \right)\left( 5 \right)\left( 6 \right)\left( 8 \right)\left( 9 \right)\left( {10} \right)}}\)

\( \displaystyle n = 8:\hspace{0.25in}{a_{11}} = \frac{{ - 8{a_7}}}{{\left( 9 \right)\left( {10} \right)\left( {11} \right)}} = 0\)

\( \displaystyle n = 9:\hspace{0.25in}{a_{12}} = \frac{{ - 9{a_8}}}{{\left( {10} \right)\left( {11} \right)\left( {12} \right)}} = \frac{{ - \left( 5 \right)\left( 9 \right){a_0}}}{{\left( 2 \right)\left( 3 \right)\left( 4 \right)\left( 6 \right)\left( 7 \right)\left( 8 \right)\left( {10} \right)\left( {11} \right)\left( {12} \right)}}\)

\( \displaystyle n = 10:\hspace{0.25in}{a_{13}} = \frac{{ - 10{a_9}}}{{\left( {11} \right)\left( {12} \right)\left( {13} \right)}} = \frac{{ - \left( 2 \right)\left( 6 \right)\left( {10} \right){a_1}}}{{\left( 3 \right)\left( 4 \right)\left( 5 \right)\left( 7 \right)\left( 8 \right)\left( 9 \right)\left( {11} \right)\left( {12} \right)\left( {13} \right)}}\)

\( \displaystyle n = 11:\hspace{0.25in}{a_{14}} = \frac{{ - 11{a_{10}}}}{{\left( {12} \right)\left( {13} \right)\left( {14} \right)}} = \frac{{ - \left( 3 \right)\left( 7 \right)\left( {11} \right){a_2}}}{{\left( 4 \right)\left( 5 \right)\left( 6 \right)\left( 8 \right)\left( 9 \right)\left( {10} \right)\left( {12} \right)\left( {13} \right)\left( {14} \right)}}\)

Okay, we can now break the coefficients down into 4 sub cases given by \({a_{4k}}\), \({a_{4k + 1}}\), \({a_{4k + 2}}\) and \({a_{4k + 3}}\) for \(k = 0,1,2,3, \ldots \) We’ll give a semi-detailed derivation for \({a_{4k}}\) and then leave the rest to you with only couple of comments as they are nearly identical derivations.

First notice that all the \({a_{4k}}\) terms have \({a_0}\) in them and they will alternate in sign. Next notice that we can turn the denominator into a factorial, \(\left( {4k} \right)!\) to be exact, if we multiply top and bottom by the numbers that are already in the numerator and so this will turn these numbers into squares. Next notice that the product in the top will start at 1 and increase by 4 until we reach \(4k - 3\). So, taking all of this into account we get,

\[{a_{4k}} = \frac{{{{\left( { - 1} \right)}^k}{{\left( 1 \right)}^2}{{\left( 5 \right)}^2} \cdots {{\left( {4k - 3} \right)}^2}{a_0}}}{{\left( {4k} \right)!}}\hspace{0.25in}k = 1,2,3, \ldots \]

and notice that this will only work starting with \(k = 1\) as we won’t get \({a_0}\) out of this formula as we should by plugging in \(k = 0\).

Now, for \({a_{4k + 1}}\) the derivation is almost identical and so the formula is,

\[{a_{4k + 1}} = \frac{{{{\left( { - 1} \right)}^k}{{\left( 2 \right)}^2}{{\left( 6 \right)}^2} \cdots {{\left( {4k - 2} \right)}^2}{a_1}}}{{\left( {4k + 1} \right)!}}\hspace{0.25in}k = 1,2,3, \ldots \]

and again notice that this won’t work for \(k = 0\)

The formula for \({a_{4k + 2}}\) is again nearly identical except for this one note that we also need to multiply top and bottom by a 2 in order to get the factorial to appear in the denominator and so the formula here is,

\[{a_{4k + 2}} = \frac{{2{{\left( { - 1} \right)}^k}{{\left( 3 \right)}^2}{{\left( 7 \right)}^2} \cdots {{\left( {4k - 1} \right)}^2}{a_2}}}{{\left( {4k + 2} \right)!}}\hspace{0.25in}k = 1,2,3, \ldots \]

noticing yet one more time that this won’t work for \(k = 0\).

Finally, we have \({a_{4k + 3}} = 0\) for \(k = 0,1,2,3, \ldots \)

Now that we have all the coefficients let’s get the solution,

\[\begin{align*}y\left( x \right) & = {a_0} + {a_1}x + {a_2}{x^2} + {a_3}{x^3} + \cdots + {a_{4k}}{x^{4k}} + {a_{4k + 1}}{x^{4k + 1}} + {a_{4k + 3}}{x^{4k + 3}} + {a_{4k + 3}}{x^{4k + 3}} + \cdots \\ & = {a_0} + {a_1}x + {a_2}{x^2} \cdots + \frac{{{{\left( { - 1} \right)}^k}{{\left( 1 \right)}^2}{{\left( 5 \right)}^2} \cdots {{\left( {4k - 3} \right)}^2}{a_0}\,}}{{\left( {4k} \right)!}}{x^{4k}} + \\ & \hspace{2.05in}\frac{{{{\left( { - 1} \right)}^k}{{\left( 2 \right)}^2}{{\left( 6 \right)}^2} \cdots {{\left( {4k - 2} \right)}^2}{a_1}}}{{\left( {4k + 1} \right)!}}{x^{4k + 1}} + \\ & \hspace{2.05in}\frac{{2{{\left( { - 1} \right)}^k}{{\left( 3 \right)}^2}{{\left( 7 \right)}^2} \cdots {{\left( {4k - 1} \right)}^2}{a_2}\,}}{{\left( {4k + 2} \right)!}}{x^{4k + 2}} + \cdots \end{align*}\]

Collecting up the terms that contain the same coefficient (except for the first one in each case since the formula won’t work for those) and writing everything as a set of series gives us our solution,

\[\begin{align*}y\left( x \right) & = {a_0}\left\{ {1 + \sum\limits_{k = 1}^\infty {\frac{{{{\left( { - 1} \right)}^k}{{\left( 1 \right)}^2}{{\left( 5 \right)}^2} \cdots {{\left( {4k - 3} \right)}^2}{x^{4k}}\,}}{{\left( {4k} \right)!}}} } \right\} + \\ & \hspace{0.5in}{a_1}\left\{ {x + \sum\limits_{k = 1}^\infty {\frac{{{{\left( { - 1} \right)}^k}{{\left( 2 \right)}^2}{{\left( 6 \right)}^2} \cdots {{\left( {4k - 2} \right)}^2}{x^{4k + 1}}}}{{\left( {4k + 1} \right)!}}} } \right\} + \\ & \hspace{0.5in}{a_2}\left\{ {{x^2} + \sum\limits_{k = 1}^\infty {\frac{{2{{\left( { - 1} \right)}^k}{{\left( 3 \right)}^2}{{\left( 7 \right)}^2} \cdots {{\left( {4k - 1} \right)}^2}{x^{4k + 2}}}}{{\left( {4k + 2} \right)!}}} } \right\}\end{align*}\]