Power series solution for 1st-order ODE

Click For Summary
SUMMARY

The power series solution for the first-order ordinary differential equation (ODE) given by y' = x^2y is derived using the series expansion y = ∑(n=0 to ∞) a_n x^n. The differentiation yields y' = ∑(n=1 to ∞) n a_n x^(n-1), which leads to the equation ∑(n=0 to ∞) (n + 1)a_(n + 1)x^n = ∑(n=2 to ∞) a_(n - 2)x^n. The coefficients a_1 and a_2 are determined to be zero, resulting in the solution y = a_0 exp(x^3/3). This confirms that the initial conditions y'(0) = 0 and y''(0) = 0 are satisfied.

PREREQUISITES
  • Understanding of power series and their convergence
  • Familiarity with ordinary differential equations (ODEs)
  • Knowledge of Taylor series expansion
  • Basic calculus, including differentiation and series manipulation
NEXT STEPS
  • Study the method of power series solutions for higher-order ODEs
  • Learn about the Frobenius method for solving differential equations
  • Explore the implications of initial conditions on series solutions
  • Investigate the relationship between power series and special functions, such as exponential functions
USEFUL FOR

Mathematicians, physics students, and engineers looking to deepen their understanding of differential equations and power series methods in solving ODEs.

Luminous Blob
Messages
50
Reaction score
0
I am trying to find the power series solution of

[tex]y' = x^2y[/tex]

but don't know how to arrive at the answer of [tex]y = a_0exp(x^3/3)[/tex]. [I know that it's an easily solved separable equation, I'm just trying to figure out how to find the power series solution]

My solution so far:

Assume
[tex]y = \sum_{n=0}^\infty a_nx^n[/tex]

then
[tex]y' = \sum_{n=1}^\infty na_nx^{n-1}[/tex]

giving:

[tex]\sum_{n=1}^\infty na_nx^{n-1} = x^2 \sum_{n=0}^\infty a_nx^n[/tex]
[tex]\sum_{n=1}^\infty na_nx^{n-1} = \sum_{n=0}^\infty a_nx^{n+2}[/tex]

changing the index for the LHS to give [tex]x^n[/tex]

[tex]\sum_{n=0}^\infty (n + 1)a_{n + 1}x^n[/tex]

changing the index for the RHS to give [tex]x^n[/tex]:

[tex]\sum_{n=2}^\infty a_{n - 2}x^{n}[/tex]

Then taking the first two terms out of the LHS sum, so that both sums start from the same point:

[tex]a_1 + 2a_2x + \sum_{n=2}^\infty (n + 1)a_{n + 1}x^n = \sum_{n=2}^\infty a_{n - 2}x^{n}[/tex]

I don't know what to do after this (I'm not entirely sure if what I've done so far is right, either).

If the y term in the initial equation didn't have the [tex]x^2[/tex] in front of it, it would be easy to equate the coefficients of [tex]x^n[/tex] to get the recursion formula. But having the terms [tex]a_1[/tex] and [tex]2a_2x[/tex] in front of the sum on the LHS throws me - can anyone explain clearly to me the correct steps required to solve the problem?
 
Last edited:
Physics news on Phys.org
We know that [tex]y'(0) = a_1[/tex] and [tex]y''(0) = 2 a_2[/tex]

We also know that [tex]y'(x) = x^2 \ y(x)[/tex] so that [tex]y'(0) = 0 = a_1[/tex]

Differentiating both sides we get [tex]y''(x) = 2x \ y(x) + x^2 \ y'(x)[/tex] and it is seen that [tex]y''(0) = 0 = 2 a_2[/tex] so now we know that [tex]a_1=a_2=0[/tex].

Hope this helps!
 
MathNerd said:
We know that [tex]y'(0) = a_1[/tex] and [tex]y''(0) = 2 a_2[/tex]

We also know that [tex]y'(x) = x^2 \ y(x)[/tex] so that [tex]y'(0) = 0 = a_1[/tex]

Differentiating both sides we get [tex]y''(x) = 2x \ y(x) + x^2 \ y'(x)[/tex] and it is seen that [tex]y''(0) = 0 = 2 a_2[/tex] so now we know that [tex]a_1=a_2=0[/tex].

Hope this helps!

Thanks, that was a great help.
 

Similar threads

  • · Replies 1 ·
Replies
1
Views
4K
  • · Replies 1 ·
Replies
1
Views
1K
  • · Replies 4 ·
Replies
4
Views
2K
  • · Replies 7 ·
Replies
7
Views
4K
  • · Replies 0 ·
Replies
0
Views
5K
  • · Replies 1 ·
Replies
1
Views
3K
  • · Replies 2 ·
Replies
2
Views
2K
  • · Replies 4 ·
Replies
4
Views
3K
  • · Replies 3 ·
Replies
3
Views
2K
  • · Replies 4 ·
Replies
4
Views
2K