Mathematical methods of Physics, ODE

Join the discussion
Ask a follow-up here, or get your own question answered by working scientists, mathematicians and engineers — people, not an autocomplete.
Real named experts · corrections over time · the nuance an AI answer skips
10 replies · 3K views
fluidistic
Gold Member
Messages
3,934
Reaction score
286

Homework Statement


For what values of K does the DE [itex]xy''-2xy'+(K-3x)y=0[/itex] (1) has a bounded solution in [itex](0, \infty)[/itex]?

Homework Equations


Not so sure, Frobenius method maybe.

The Attempt at a Solution


First, I check what happens when x tends to infinity. I see that the DE behaves like [itex]\phi ''-2 \phi '-3 \phi =0[/itex]. Where [itex]\phi (x)=c_1e^{3x}+c_2e^{-x}[/itex]. I checked out and this indeed satisfies the DE (2).
I propose a solution to (1) of the form [itex]y(x)=f(x) \phi (x)[/itex]. So I calculated y' and y'' and replaced all this into the DE (1) and after a lot of algebra I reached that the DE (1) is equivalent to the DE [itex]f''(x)+4f'(x)+f(x) \left ( \frac{Kc_1e^{3x}}{x} \right ) \cdot \left ( \frac{1}{c_1e^{3x}+c_2e^{-x}} \right )=0[/itex] (3) where these [itex]c_1[/itex] and [itex]c_2[/itex] differs from the ones given in the expression for [itex]\phi (x)[/itex] but this doesn't matter, they are constants. Now I must say I'm not really eager to use Frobenius's method on the DE (3). Mainly due to the denominator. I don't really know how to proceed... any help will be appreciated as usual.
 
Physics news on Phys.org
Not sure if it will help, but surely you wouldn't want to include the e3x in your calculations because it isn't bounded. I would try just the e-x part, but no guarantees; I haven't tried it myself.
 
LCKurtz said:
Not sure if it will help, but surely you wouldn't want to include the e3x in your calculations because it isn't bounded. I would try just the e-x part, but no guarantees; I haven't tried it myself.
Whoops, you're absolutely right! I forgot about this. :redface:
 
I progressed on the problem.
I skip the algebra, and I must solve [itex]f''(x)-4f'(x)+\frac{K}{x}f(x)=0[/itex]. I use Frobenius's method on it and this gave me the indicial equation s=0 (that I discard later) and s=1.
I reach the recurrence relation [itex]a_{n+1}=\frac{-a_n [4(n+s)+K]}{(n+2)(n+1)}[/itex].
The values of K they ask for is when [itex]a_{n+1}=0[/itex] and [itex]a_n \neq 0[/itex]. So [itex]K=-4n-4[/itex] with n going from 0 to infinity. (This is, I hope, the answer to the problem).
However say I want to find the general terms [itex]a_n[/itex] in [itex]\sum _{n=0}^{\infty} a_nx^{n+1}[/itex].
I'm currently trying to find a pattern, by setting [itex]a_0=1[/itex] between the [itex]a_n[/itex]'s but I'm not getting successful in it.
[itex]a_0=1[/itex], [itex]a_1=-\frac{4+K}{2}[/itex], [itex]a_2=\left ( \frac{4+K}{2} \right )\left ( \frac{8+K}{3\cdot 2} \right )[/itex], [itex]a_3=-\frac{(4+K)}{2} \left ( \frac{8+K}{3\cdot 2} \right ) \left ( \frac{12+K}{4\cdot 3} \right )[/itex].

So far I realize that [itex]a_n=\frac{(-1)^n}{n!(n+1)!}[/itex] multiplied by something that I'm not sure (and don't find yet). Maybe something similar to [itex]\frac{(4n+K)!}{(4n+K-1)(4n+K-2)(4n+K-3)}[/itex]. Can someone help me on this task?
 
vela said:
Did you perhaps make a sign error? I find the recursion relation to be
[tex]a_{n+1} = -\frac{K-4(n+1)}{(n+1)(n+2)} a_n.[/tex]

Yes I think I did, though now I get a different answer from ours.
From [itex]f''(x)-4f'(x)+\frac{K}{x}f(x)=0[/itex] and [itex]f(x)=\sum_{n=0}^{\infty}a_n x^{n+s}[/itex], I reach that it's equivalent to [itex]\sum_{n=-1}^{\infty}a_{n+1}(n+s)(n+s+1) x^{n+s-2}-4\sum_{n=0}^{\infty}a_n(n+s) x^{n+s-1}-K \sum_{n=0}^{\infty}a_n x^{n+s-1}=0[/itex]
[itex]\Rightarrow a_0s(s-1)x^{s-2}+\sum_{n=0}^{\infty}a_n x^{n+s-1} [-4(n+s)-K]+a_{n+1}(n+s+1)(n+s)x^{n+s-1}=0[/itex].
So [itex]-a_n [4(n+s)+K]+a_{n+1}(n+s+1)(n+s)=0\Rightarrow a_{n+1}=\frac{a_n[4(n+s)+K]}{(n+s+1)(n+s)}[/itex].

Mod note: fixed LaTeX
 
Last edited by a moderator:
fluidistic said:
Yes I think I did, though now I get a different answer from ours.
From [itex]f''(x)-4f'(x)+\frac{K}{x}f(x)=0[/itex] and [itex]f(x)=\sum_{n=0}^{\infty}a_n x^{n+s}[/itex], I reach that it's equivalent to [tex]\sum_{n=-1}^{\infty}a_{n+1}(n+s)(n+s+1) x^{n+s-2}-4\sum_{n=0}^{\infty}a_n(n+s) x^{n+s-1}-K \sum_{n=0}^{\infty}a_n x^{n+s-1}=0[/tex]
You changed the sign on the K term, so you should get
[tex]a_0s(s-1)x^{s-2}+\sum_{n=0}^{\infty}a_n x^{n+s-1} [K-4(n+s)]+a_{n+1}(n+s+1)(n+s)x^{n+s-1}=0.[/tex]Then you made one more sign error which got rid of the alternating sign.
 
You are right, I made a typo when I copied the DE when I turned the page of my draft for the sign of K. I now reach your result.
I'm trying to find a pattern by calculating the first terms, assuming [itex]a_0=1[/itex].
[itex]a_1=-\frac{(K-4)}{2}, a_2=\frac{(K-4)(K-8)}{2\cdot 3\cdot 2}[/itex], etc.
My attempt is [itex]a_n=\frac{(-1)^n(K-4)(K-8)...(K-4n)}{n!(n+1)!}[/itex] but this doesn't work for [itex]a_0[/itex] since I get K instead of 1.
 
vela said:
If you require that the series terminates, you know K has the form K=4m where m=1, 2, 3, ...

Yes indeed. [itex]K-4(n+1)=0 \Rightarrow k=4(n+1)[/itex] with n=0, 1, etc. which is equivalent to K=4m for m=1, 2, etc.
I was just trying to do more than what they asked me. I wanted to express [itex]a_n[/itex] in terms of n only rather than [itex]a_{n+1}[/itex] in terms of [itex]a_n[/itex]. I don't think it really matters anyway, but since I've a weakness in doing it I wanted to eradicate it. :smile: