A power series solution to a differential equation is a function with infinitely many terms, each term containing a different power of the dependent variable. The general solution has the form $y=a_{0}+a_{1}x+a_{2}x^{2}+a_{3}x^{3}+\cdots$ .

To illustrate how to use the power series method to arrive at a solution, let's take the equation $y''+x^{2}y'+xy=0$ .

Start with the general form of the power series: $y=a_{0}+a_{1}x+a_{2}x^{2}+a_{3}x^{3}+\cdots$

Take derivatives of the general power series until an equation exists for every order of derivative in the equation. You can shift the nth term so that all the terms end with $x^{n}$ . For example,

$y'=a_{1}+2a_{2}x+\cdots +na_{n}x^{n-1}$ or in alternative form, $y'=a_{1}+2a_{2}x+\cdots +(n+1)a_{n+1}x^{n}$

To find a "recursive" solution, compare the general nth terms of each equation and substitute them into the original differential equation. If there are dependent variables in the differential equation, you might have to alter the nth term to get all the powers of the dependent variable the same. For example: