Another series solution ODE problem

Click For Summary
SUMMARY

The discussion focuses on solving the second-order ordinary differential equation (ODE) given by (1-x^2)y''+(1-x)y'-3y = 0 using a power series approach. The user encountered a situation where three coefficients were derived for the same power of x, which is typical in second-order differential equations. The solution can be expressed in terms of two arbitrary constants, a_0 and a_1, leading to a general solution format similar to y(t) = a_0cos(ωt) + a_1sin(ωt). The user successfully constructed a series solution, although it did not fully align with results obtained from Mathematica, indicating potential algebraic simplifications needed.

PREREQUISITES
  • Understanding of second-order ordinary differential equations (ODEs)
  • Familiarity with power series solutions
  • Knowledge of recursive relations in series expansions
  • Basic algebraic manipulation skills
NEXT STEPS
  • Study the method of power series solutions for second-order ODEs
  • Learn about recursive relations in series expansions
  • Explore the use of Mathematica for solving differential equations
  • Investigate simplification techniques for power series
USEFUL FOR

Students and educators in mathematics, particularly those focusing on differential equations, as well as anyone interested in applying power series methods to solve ODEs.

TheFerruccio
Messages
216
Reaction score
0

Homework Statement



Find a basis of solutions.

Homework Equations



[tex](1-x^2)y''+(1-x)y'-3y = 0[/tex]

The Attempt at a Solution



Using the series approach, having:

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

I ended up with an equation representing the coefficients for [tex]x^0[/tex]

[tex]2a_2+a_1-3a_0 = 0[/tex]

I'm confused because this is the first time I have encountered three coefficients for the same power of x. I just want to verify that I didn't make an algebra mistake, because I am used to finding that certain coefficients are linked to even powers, while others are linked to odd powers.
 
Physics news on Phys.org
Hi TheFerruccio,

Your work looks fine, I just went through the problem and obtain the same result you got. I have seen a lot of power series solutions that end up with three coefficients in their recursive relations, so it is no problem. This happens when you have second order differential equations. Given that a second order differential equation will admit two arbitrary constants for the general solution, your single equations in three unknowns is a way of saying there are two 'free parameters' (if I may borrow the term from linear algebra). Your final solution will be able to be expressed in terms of two solutions that depend each on one of the arbitrary constants: [tex]a_0[/tex] and [tex]a_1[/tex]. If you had initial/boundary conditions, you would be able to find these constants.

For instance, the differential equation:

[tex]y'' + \omega^2y = 0[/tex]

If you used a power series solution to find the result, you would get a recursive relation in three coefficients, and you could compute several terms. In the end you would be able to reduce all the coefficients to be put in terms the first two coefficients: [tex]a_0[/tex] and [tex]a_1[/tex], such that your final solution would look like:

[tex]y(t) = a_0\left(1 - \frac{(\omega t)^2}{2!} + \frac{(\omega t)^4}{4!} - \ldots\right) + a_1\left(\omega t - \frac{(\omega t)^3}{3!} + \frac{(\omega t)^5}{5!} - \ldots \right)[/tex]

Which you would recognize as:

[tex]y(t) = a_0\cos\omega t + a_1\sin\omega t[/tex]

That only happens though if your recursive relation would have three parameters in it to start with. It is a common situation though, so no worries.
 
Last edited:
Thanks for the help! I managed to put together some kind of solution, though it does not really match what was found using mathematica. I'm probably missing something with the algebra or the simplification. Here is the solution I have:

[tex]y = a_0(1 + \frac{3}{2}x^2 + \frac{1}{2}x^3 + \frac{3}{4}x^4 + \frac{3}{4}x^5 + \ldots) + a_1(x - \frac{1}{2}x^2 + \frac{1}{2}x^3 + \frac{1}{4}x^4 + \frac{1}{4}x^5 + \ldots)[/tex]

EDIT: It looks like this works just fine. There is likely a simplification for some of the series, but I don't know it.
 
Last edited:

Similar threads

Replies
8
Views
3K
  • · Replies 16 ·
Replies
16
Views
2K
  • · Replies 1 ·
Replies
1
Views
2K
  • · Replies 1 ·
Replies
1
Views
2K
Replies
3
Views
2K
  • · Replies 2 ·
Replies
2
Views
2K
  • · Replies 1 ·
Replies
1
Views
2K
  • · Replies 3 ·
Replies
3
Views
2K
  • · Replies 5 ·
Replies
5
Views
2K
  • · Replies 3 ·
Replies
3
Views
2K