Can ##a_1## and ##a_2## be zero in the Frobenius' method solution for this DE?

In summary, the author attempted to solve this DE but is stuck on how to handle the terms a_1 and a_2. He has eliminated the cases where a_0=0 and r=1 or -1/5, but is not sure about the case where a_0=0.
  • #1
BearY
53
8

Homework Statement


Use Forbenius' method to solve this DE:
$$ 5x^2y''+xy'+(x^3-1)y=0$$

Homework Equations


Seek power series solution in the form ##y=\sum _{n=0}^{\infty } a_n x^{n+r}##, ##a_0\neq0##

The Attempt at a Solution


Sub in the ansatz y, get $$ \sum _{n=0}^{\infty }a_n[5(n+r)(n+r+1)+n+r-1]x^n + \sum _{n=0}^{\infty }a_nx^{n+r+3} = 0 $$
Align the powers of x, split first sum to combine with second sum:$$ \sum _{n=0}^{\infty }\{a_{n+3}[5(n+r+3)(n+r+2)+n+r+3-1]+a_n\}x^{n+3} + a_0[5r(r-1)+r-1]+a_1[r(r+1)r+r]x+a_2[5(r+2)(r+1)+r+1]x^2=0$$
This is where I get stuck, I am not sure about how to handle ##a_1## and ##a_2## term or how to find solutions of ##r## with them being there. I can't find out if ##a_1## and ##a_2## can be zero or not, whether or why it matters.
Thanks in advance.
 
Physics news on Phys.org
  • #2
BearY said:

Homework Statement


Use Forbenius' method to solve this DE:
$$ 5x^2y''+xy'+(x^3-1)y=0$$

Homework Equations


Seek power series solution in the form ##y=\sum _{n=0}^{\infty } a_n x^{n+r}##, ##a_0\neq0##

The Attempt at a Solution


Sub in the ansatz y, get $$ \sum _{n=0}^{\infty }a_n[5(n+r)(n+r+1)+n+r-1]x^n + \sum _{n=0}^{\infty }a_nx^{n+r+3} = 0 $$
Align the powers of x, split first sum to combine with second sum:$$ \sum _{n=0}^{\infty }\{a_{n+3}[5(n+r+3)(n+r+2)+n+r+3-1]+a_n\}x^{n+3} + a_0[5r(r-1)+r-1]+a_1[r(r+1)r+r]x+a_2[5(r+2)(r+1)+r+1]x^2=0$$
This is where I get stuck, I am not sure about how to handle ##a_1## and ##a_2## term or how to find solutions of ##r## with them being there. I can't find out if ##a_1## and ##a_2## can be zero or not, whether or why it matters.
Thanks in advance.

I have not checked your work, but assuming it is correct we can say:
(1) either ##a_0=0## or ##5r(r-1)+r-1 = 0## (or both).
(2) either ##a_1=0## or ##r^2(r+1)+r = 0## (or both).
(3) either ##a_2 = 0## or ##5(r+1)(r+2)+ r+1=0## (or both).

From (1): if ##a_0 \neq 0## then ##(5r+1)(r-1) = 0##, so either ##r = 1## or ##r = -1/5##. Thus, if ##r = 1## or ##r = -1/5## then ##a_0## is arbitrary and ##a_1 = a_2 = 0##. If ##r \neq 1, -1/5## then ##a_0 = 0##.

Look at the other cases in a similar way. That will give you a candidate list of possible ##r## values and corresponding cases for ##a_0,a_1,a_2 .## Then you need to figure out how to eliminate some of those possibilities.

Note added in edit: I guess we can eliminate the case ##a_0 = 0##, because ##a_0 x^r## is the lowest-power term in ##y##, so cannot have a zero coefficient. In other words, we do have ##a_0 \neq 0##, so we can say for sure that either ##r = 1## or ##r = -1/5##, and also that ##a_1 = a_2 = 0## (assuming your work was correct).
 
Last edited:
  • Like
Likes BearY
  • #3
Two polynomials ##p(x) = \sum a_nx^n\, , \,q(x)= \sum b_nx^n## are equal if all ##a_n=b_n## and if ##q(x)=0## then all ##b_n=0##. Now which coefficients at position ##0 , 1, 2 ## are left on the left hand side? Don't forget, that ##a_0 \neq 0## if you solve for ##r##. And do you know something about ##r##, e.g. ##r > 0## or can it be arbitrary?
 
  • Like
Likes BearY
  • #4
Ray Vickson said:
I have not checked your work, but assuming it is correct we can say:
(1) either ##a_0=0## or ##5r(r-1)+r-1 = 0## (or both).
(2) either ##a_1=0## or ##r^2(r+1)+r = 0## (or both).
(3) either ##a_2 = 0## or ##5(r+1)(r+2)+ r+1=0## (or both).

From (1): if ##a_0 \neq 0## then ##(5r+1)(r-1) = 0##, so either ##r = 1## or ##r = -1/5##. Thus, if ##r = 1## or ##r = -1/5## then ##a_0## is arbitrary and ##a_1 = a_2 = 0##. If ##r \neq 1, -1/5## then ##a_0 = 0##.

Look at the other cases in a similar way. That will give you a candidate list of possible ##r## values and corresponding cases for ##a_0,a_1,a_2 .## Then you need to figure out how to eliminate some of those possibilities.

That is to say since ##a_0## term must be 0, ##a_0 \neq 0##, r can only be the solutions to ##a_0## term equal to zero. If the other two coefficients happen to be arbitrary constants that don't have to be 0, they can be included in the coefficients of the general solution of the DE, right?
 
  • #5
fresh_42 said:
Two polynomials ##p(x) = \sum a_nx^n\, , \,q(x)= \sum b_nx^n## are equal if all ##a_n=b_n## and if ##q(x)=0## then all ##b_n=0##. Now which coefficients at position ##0 , 1, 2 ## are left on the left hand side? Don't forget, that ##a_0 \neq 0## if you solve for ##r##. And do you know something about ##r##, e.g. ##r > 0## or can it be arbitrary?
I beg your pardon, I am not sure what does being left on the left hand side mean. r is arbitary constant.
 
  • #6
BearY said:
That is to say since ##a_0## term must be 0, ##a_0 \neq 0##, r can only be the solutions to ##a_0## term equal to zero. If the other two coefficients happen to be arbitrary constants that don't have to be 0, they can be included in the coefficients of the general solution of the DE, right?
I have no idea what you are attempting to say; it makes no sense to me. Re-read my post---especially the part added in edit.
 
  • Like
Likes BearY
  • #7
BearY said:
I beg your pardon, I am not sure what does being left on the left hand side mean. r is arbitary constant.
You have an equation ##\text{ (LHS =) } p(x) = p_0 +p_1x + p_2x^2 + \ldots = 0 = 0\cdot 1 + 0 \cdot x + 0\cdot x^2 + \ldots \text{ (= RHS) }##. Now compare the positions at ##x^n \, : \, p_0=0 \; , \; p_1=0\; , \; p_2=0\; \ldots \;p_n=0##. I wrote ##p_i## because those new coefficients are themselves polynomial expressions in the ##a_j##, e.g. ##p_0 = a_0 (5r^2-4r-1)##. This can be solved and you get two possible values for ##r##. Take them and continue with ##p_1 = a_1 \left( (1+r)\cdot (1+5r)-1 \right)## and so on.
 
  • Like
Likes BearY
  • #8
Ray Vickson said:
I have no idea what you are attempting to say; it makes no sense to me. Re-read my post---especially the part added in edit.
I am sorry, I meant to ask, in this question we can indeed be sure that ##a_1 = a_2 = 0##. But what happens if a solution to ##r## allows other coefficients to be non-zero? Or is that impossible?
 
  • #9
BearY said:
I am sorry, I meant to ask, in this question we can indeed be sure that ##a_1 = a_2 = 0##.
Re-read post #3 which explains the comparison of two polynomials.
But what happens if a solution to ##r## allows other coefficients to be non-zero? Or is that impossible?
You cannot rule it out. All you know is ##a_0 \neq 0##. E.g. I got ##a_3 = -\dfrac{1}{63} a_0## in one case, which is thus unequal to zero, too.
 
  • Like
Likes BearY
  • #10
fresh_42 said:
You have an equation ##\text{ (LHS =) } p(x) = p_0 +p_1x + p_2x^2 + \ldots = 0 = 0\cdot 1 + 0 \cdot x + 0\cdot x^2 + \ldots \text{ (= RHS) }##. Now compare the positions at ##x^n \, : \, p_0=0 \; , \; p_1=0\; , \; p_2=0\; \ldots \;p_n=0##. I wrote ##p_i## because those new coefficients are themselves polynomial expressions in the ##a_j##, e.g. ##p_0 = a_0 (5r^2-4r-1)##. This can be solved and you get two possible values for ##r##. Take them and continue with ##p_1 = a_1 \left( (1+r)\cdot (1+5r)-1 \right)## and so on.
I understand it now. Thank you
fresh_42 said:
Re-read post #3 which explains the comparison of two polynomials.

You cannot rule it out. All you know is ##a_0 \neq 0##. E.g. I got ##a_3 = -\dfrac{1}{63} a_0## in one case, which is thus unequal to zero, too.
I think ##a_3## will be the result of the recurrence relationship given by the sum in the last step I wrote in #1. And so will be ##a_6## and so on. The value of ##r## only makes ##a_1 = a_2 = 0## right?
 
  • #11
In the case ##r=1## I had ##a_1=a_2=0## as well. I think this is also true for the other value of ##r##, but I haven't checked.
 
  • Like
Likes BearY
  • #12
BearY said:
I am sorry, I meant to ask, in this question we can indeed be sure that ##a_1 = a_2 = 0##. But what happens if a solution to ##r## allows other coefficients to be non-zero? Or is that impossible?

We need ##a_1 = a_2 = 0##. Then we can get ##a_3## in terms of ##a_0##, can get ##a_4## in terms of ##a_1##---so is 0, etc. Altogether, the only nonzero values of ##a_n## for ##n >0## are ##a_3, a_6, a_9, \ldots ,## and in principle you can get them all in terms of ##a_0.##
 
Last edited:
  • Like
Likes BearY
  • #13
Ray Vickson said:
We need ##a_1 = a_2 = 0##. Then we can get ##a_3## in terms of ##a_0##, can get ##a_4## in terms of ##a_1##---so is 0, etc. Altogether, the only nonzero values of ##a_n## for ##n >0## are ##a_3, a_6, a_9, \ldots ,## and in principle you can get them all in terms of ##a_0.##
I am sorry for the vagueness, I meant other as ##a_1## and ##a_2##, I don't know why I said other instead of that, doesn't save much typing anyway. I am just wondering what if for ##a_0f(r)## and ## a_1g(r)##, the solution to ##f(r) = 0## and ##g(r) = 0## are the same. Or is it impossible?
 
  • #14
BearY said:
I am sorry for the vagueness, I meant other as ##a_1## and ##a_2##, I don't know why I said other instead of that, doesn't save much typing anyway. I am just wondering what if for ##a_0f(r)## and ## a_1g(r)##, the solution to ##f(r) = 0## and ##g(r) = 0## are the same. Or is it impossible?

Obviously, if ##a_0 \neq 0## the solution of ##a_0 f(r) = 0## is the same as that of ##f(r) = 0##. But if ##a_0 = 0## the equation ##a_0 f(r) = 0## contains no useful information about ##r##.

So: go back are re-read past posts: we DO have ##a_0 \neq 0## (essentially by definition). That tells us enough about ##r## that we can say with absolute, 100% certainty, that ##g(r) \neq 0##, FORCING ##a_1 = 0##. End of story.
 
  • Like
Likes BearY

1. What is Frobenius’ method and what is it used for?

Frobenius’ method is a mathematical technique used to solve differential equations with singular points or irregular points. It is commonly used in physics and engineering to model physical systems.

2. How does Frobenius’ method differ from other methods of solving differential equations?

Frobenius’ method is specifically designed to handle equations with singular points, such as those found in physical systems. Other methods may not be able to accurately solve these types of equations.

3. Can Frobenius’ method be used for all types of differential equations?

No, Frobenius’ method is most effective for solving second-order linear differential equations with singular points. It may not be suitable for other types of equations.

4. What are some common applications of Frobenius’ method?

Frobenius’ method is commonly used in physics and engineering to model physical systems such as pendulums, oscillators, and electrical circuits. It is also used in quantum mechanics to solve the Schrödinger equation.

5. Are there any limitations to using Frobenius’ method?

Frobenius’ method may not always provide an exact solution to a differential equation with singular points. It can also be quite complex and time-consuming for certain types of equations. Additionally, it may not be suitable for non-linear systems.

Similar threads

  • Calculus and Beyond Homework Help
Replies
1
Views
266
  • Calculus and Beyond Homework Help
Replies
2
Views
713
Replies
8
Views
1K
  • Calculus and Beyond Homework Help
Replies
8
Views
1K
  • Calculus and Beyond Homework Help
Replies
3
Views
499
  • Calculus and Beyond Homework Help
Replies
4
Views
321
  • Calculus and Beyond Homework Help
Replies
5
Views
298
  • Calculus and Beyond Homework Help
Replies
9
Views
1K
  • Calculus and Beyond Homework Help
Replies
2
Views
201
  • Calculus and Beyond Homework Help
Replies
33
Views
2K
Back
Top