Ordinary differential equations. Series method.

LagrangeEuler
Messages
711
Reaction score
22
Question:
Why equations
x(1-x)\frac{d^2y}{dx^2}+[\gamma-(\alpha+\beta+1)x]\frac{dy}{dx}-\alpha \beta y(x)=0
should be solved by choosing
##y(x)=\sum^{\infty}_{m=0}a_mx^{m+k}##
and not
##y(x)=\sum^{\infty}_{m=0}a_mx^{m}##?
How to know when we need to choose one of the forms.
Also when I sum over ##m##, then ##\sum^{\infty}_{m=0}a_mx^{m+k}=y(x,k)##. Right?
 
Physics news on Phys.org
LagrangeEuler said:
Question:
Why equations
x(1-x)\frac{d^2y}{dx^2}+[\gamma-(\alpha+\beta+1)x]\frac{dy}{dx}-\alpha \beta y(x)=0
should be solved by choosing
##y(x)=\sum^{\infty}_{m=0}a_mx^{m+k}##
and not
##y(x)=\sum^{\infty}_{m=0}a_mx^{m}##?
How to know when we need to choose one of the forms.
Also when I sum over ##m##, then ##\sum^{\infty}_{m=0}a_mx^{m+k}=y(x,k)##. Right?

Ultimately, you are looking for polynomial solutions. A polynomial need not begin at x^0. So you allow it to begin at x^k.

You could begin by substituting into the equation the series
##y(x)=\sum^{\infty}_{m=0}a_mx^{m}##
But then, if the solution polynomial does begin with the x^k term, your first set of coefficients
##a_0, a_1, ..., a^{k-1}##
would all work out to be zero. It is thus less work to allow for this eventuality by starting your series at
##x^k##
You can the solve, among other things, for the value of k.

Technically, the notation you suggest,
##\sum^{\infty}_{m=0}a_mx^{m+k}=y(x,k)##
is correct, but redundant.
 
Ok. Fair enough. I will try to solve tomorrow this with both methods. Could you just tell me how do you know when the first few will be zero or how you know that here when you see equation? Or you always use ## \sum^{\infty}_{m=0}a_mx^{m+k}##?
 
LagrangeEuler said:
Ok. Fair enough. I will try to solve tomorrow this with both methods. Could you just tell me how do you know when the first few will be zero or how you know that here when you see equation? Or you always use ## \sum^{\infty}_{m=0}a_mx^{m+k}##?

Always use ## \sum^{\infty}_{m=0}a_mx^{m+k}##
 
LagrangeEuler said:
Question:
Why equations
x(1-x)\frac{d^2y}{dx^2}+[\gamma-(\alpha+\beta+1)x]\frac{dy}{dx}-\alpha \beta y(x)=0
should be solved by choosing
##y(x)=\sum^{\infty}_{m=0}a_mx^{m+k}##
and not
##y(x)=\sum^{\infty}_{m=0}a_mx^{m}##?
How to know when we need to choose one of the forms.
Also when I sum over ##m##, then ##\sum^{\infty}_{m=0}a_mx^{m+k}=y(x,k)##. Right?
See http://mathworld.wolfram.com/FrobeniusMethod.html

x=0 is a singular point of the differential equation. That's why you have to use the first form and not a plain old Taylor series.
 
  • Like
Likes BvU
The method you are inquiring about is the Method of Frobenius. You have to use the Method of Frobenius whenever there is a singular point in the differential equation.
A singular point can be broken up into 2 categories. They are: ordinary singular point and irregular singular point (not sure if this is the correct name). I suggest reading the section action in your DE book. In other words, no taylor series expansion exist at a singular point. So we change where the Taylor Series is centered, however it gets really messy, so we employ the Method of Frobenius.

There is also a method of getting the values of K, before you even start massaging the problem.
 
Are there any good visualization tutorials, written or video, that show graphically how separation of variables works? I particularly have the time-independent Schrodinger Equation in mind. There are hundreds of demonstrations out there which essentially distill to copies of one another. However I am trying to visualize in my mind how this process looks graphically - for example plotting t on one axis and x on the other for f(x,t). I have seen other good visual representations of...
Back
Top