Kansas State University

STEP 1: Plug in $\displaystyle y(x) = \sum_{n=0}^{\infty}a_n x^n $ and compute all the different terms in the equation $$ \begin{align} (5x - 4)y &= \sum_{ n = 0 }^{ \infty }\,5 a_{ n } x^{ n+1 } + \sum_{ n = 0 }^{ \infty }\,-4 a_{ n } x^{ n }\\ 7y' &= \sum_{ n = 1 }^{ \infty }\,7n a_{ n } x^{ n-1 }\\ (6x - 1)y'' &= \sum_{ n = 2 }^{ \infty }\,6n(n-1) a_{ n } x^{ n-1 } + \sum_{ n = 2 }^{ \infty }\,-n(n-1) a_{ n } x^{ n-2 } \end{align} $$ STEP 2: Make the substitutions $ k = n-2 $, $ j = n-1 $ and $ p = n+1 $ to make all terms of the form $ x^{\text{index}} $ rather than the $ x^{\text{index}-1} $ or $ x^{\text{index}+2} $ or whatever. $$ \begin{align} (5x - 4)y&= \sum_{ p = 1 }^{ \infty }\,5 a_{ p-1 } x^{ p } + \sum_{ n = 0 }^{ \infty }\,-4 a_{ n } x^{ n }\\ 7y'&= \sum_{ j = 0 }^{ \infty }\,7(j+1) a_{ j+1 } x^{ j }\\ (6x - 1)y''&= \sum_{ j = 1 }^{ \infty }\,6(j+1)j a_{ j+1 } x^{ j } + \sum_{ k = 0 }^{ \infty }\,-(k+2)(k+1) a_{ k+2 } x^{ k }\end{align} $$ STEP 3: Change all the indices to the same letter (I use $ m $) and plug into the equation. $$ \begin{align} (6x - 1)y'' + 7y' + (5x - 4)y &= \sum_{ m = 1 }^{ \infty }\,6(m+1)m a_{ m+1 } x^{ m } + \sum_{ m = 0 }^{ \infty }\,-(m+2)(m+1) a_{ m+2 } x^{ m } \\ &+ \sum_{ m = 0 }^{ \infty }\,7(m+1) a_{ m+1 } x^{ m } \\ &+ \sum_{ m = 1 }^{ \infty }\,5 a_{ m-1 } x^{ m } + \sum_{ m = 0 }^{ \infty }\,-4 a_{ m } x^{ m } \end{align} $$ STEP 4: Collect like terms. (While the general term starts at m=1, it doesn't hurt to separate out the x term as we do below). $$ \begin{align} &(-2a_2 + 7a_1 - 4a_0) + (-6a_{3} + 12a_{2} + 14a_{2} - 4a_{1} + 5a_{0})x \\ &+ \sum_{ m = 2 }^{ \infty }\,(-(m+2)(m+1)a_{m+2} + 6(m+1)ma_{m+1} + 7(m+1)a_{m+1} - 4a_m + 5a_{m-1})x^m = 0 \end{align} $$ STEP 5: Equate coefficients to 0.

Equating the constant term to 0 we get $$ a_2 = \frac{-7a_1 + 4a_0}{-2} $$ Equating the linear term to 0 we get $$ a_3 = \frac{-26a_2 + 4a_1 - 5a_0}{-6} $$ Finally, equating the general term to 0, we find that for $ m \ge 2,$ $$ a_{m+2} = \frac{-(6(m+1)m + 7(m+1))a_{m+1} + 4a_m - 5a_{m-1}}{-(m+2)(m+1)} $$ STEP 6: We know that $ a_0 = y(0) = 5 $ and $ a_1 = y'(0) = 1.$ We then plug these values into the formulas found in step 5 to compute the coefficients of the solution.

From the equation for the constant term we get $$ a_2 = \frac{-7(1) + 4(5)}{-2} = -13/2 $$ From the equation for the linear term we get $$ a_3 = \frac{-26(-13/2) + 4(1) - 5(5)}{-6} = -74/3 $$ Finally, using the recurrence equation with $ m = 2 $ we get $$ a_4 = \frac{-(6(2+1)2 + 7(2+1))(-74/3) + 4(-13/2) - 5(1)}{-(4)(3)} = -1375/12 $$

So our solution is $$ y(x) = 5 + x - (13/2)x^2 - (74/3)x^3 - (1375/12)x^4 + \cdots $$

You may reload this page to generate additional examples.