Take the 2-minute tour ×
Mathematics Stack Exchange is a question and answer site for people studying math at any level and professionals in related fields. It's 100% free, no registration required.

My friend and I have conflicting answers and since his phone is off, I can't get his full solution and I don't understand his argument.

Consider this ODE

$$(x+1)y''+\frac{1}{x}y' + (x+3)y= 0$$

Basically what I did was divide out that $(x+1)$ on $y''$ and got

$$y''+\frac{1}{x(x+1)}y' + \frac{(x+3)}{(x+1)}y= 0$$

The singularities are x = -1, and 0 (both are regular)

My friend said we had to get our ODE in the form of

$$a(x-x_0)^2y'' + b(x-x_0)y' + (x-x_0)y =0$$ otherwise we cant' do anything. and he got x = 3 as an irregular singular point which i have no idea how even got this. Sorry if this is too vague, but my first source was also as vague.

what was my friend doing and who's right?

share|improve this question
    
This might help:math.washington.edu/~toro/Courses/99-00/135/regular.pdf –  Lyapunov Jun 26 '12 at 20:02
    
It's basically what I did with different notations. Do you know how he came up with x = 3? Sorry if I am asking you to "read his mind" –  jip Jun 26 '12 at 20:15
    
I see nothing singular in this equation at $x=3$. –  user31373 Jun 26 '12 at 20:17
1  
In fact en.wikipedia.org/wiki/… , mathworld.wolfram.com/RegularSingularPoint.html , mathworld.wolfram.com/IrregularSingularity.html and mathworld.wolfram.com/… had already tell us all the necessaries for finding singular points of any 2nd order linear ODE. –  doraemonpaul Jun 27 '12 at 0:13

1 Answer 1

up vote 1 down vote accepted

$(x+1)y''+\dfrac{1}{x}y'+(x+3)y=0$

$y''+\dfrac{1}{x(x+1)}y'+\dfrac{x+3}{x+1}y=0$

It should be noted that the positions of finite singular points are always appear at the positions that makes one of the coefficients diverge.

$\therefore$ the positions of finite singular points in this question are $x=0$ and $x=-1$ only, absolute not like your friend said $x=3$ also is the finite singular point.

The next step is to determine whether the finite singular points are regular or irregular.

$\lim_{x\to 0}\left(x\times\dfrac{1}{x(x+1)}\right)=\lim_{x\to 0}\dfrac{1}{x+1}=0$

$\lim_{x\to 0}\left(x^2\times\dfrac{x+3}{x+1}\right)=\lim_{x\to 0}\dfrac{x^2(x+3)}{x+1}=0$

$\therefore$ the finite singular point $x=0$ is regular.

$\lim_{x\to-1}\left((x+1)\times\dfrac{1}{x(x+1)}\right)=\lim_{x\to-1}\dfrac{1}{x}=-1$

$\lim_{x\to-1}\left((x+1)^2\times\dfrac{x+3}{x+1}\right)=\lim_{x\to-1}(x+1)(x+3)=0$

$\therefore$ the finite singular point $x=-1$ is regular.

That's not the end. We should also check the singularities at infinity. Because these also act as one of the key points of distinguishing the ODE type.

Let $u=\dfrac{1}{x}$ ,

Then $\dfrac{dy}{dx}=\dfrac{dy}{du}\dfrac{du}{dx}=-\dfrac{1}{x^2}\dfrac{dy}{du}=-u^2\dfrac{dy}{du}$

$\dfrac{d^2y}{dx^2}=\dfrac{d}{dx}\left(-u^2\dfrac{dy}{du}\right)=\dfrac{d}{du}\left(-u^2\dfrac{dy}{du}\right)\dfrac{du}{dx}=\left(-u^2\dfrac{d^2y}{du^2}-2u\dfrac{dy}{du}\right)\left(-\dfrac{1}{x^2}\right)=\left(-u^2\dfrac{d^2y}{du^2}-2u\dfrac{dy}{du}\right)(-u^2)=u^4\dfrac{d^2y}{du^2}+2u^3\dfrac{dy}{du}$

$\therefore u^4\dfrac{d^2y}{du^2}+2u^3\dfrac{dy}{du}+\dfrac{1}{\dfrac{1}{u}\left(\dfrac{1}{u}+1\right)}\left(-u^2\dfrac{dy}{du}\right)+\dfrac{\dfrac{1}{u}+3}{\dfrac{1}{u}+1}y=0$

$u^4\dfrac{d^2y}{du^2}+2u^3\dfrac{dy}{du}-\dfrac{u^4}{u+1}\dfrac{dy}{du}+\dfrac{3u+1}{u+1}y=0$

$u^4\dfrac{d^2y}{du^2}+\dfrac{u^4+2u^3}{u+1}\dfrac{dy}{du}+\dfrac{3u+1}{u+1}y=0$

$\dfrac{d^2y}{du^2}+\dfrac{u+2}{u(u+1)}\dfrac{dy}{du}+\dfrac{3u+1}{u^4(u+1)}y=0$

$\lim_{u\to 0}\left(u\times\dfrac{u+2}{u(u+1)}\right)=\lim_{u\to 0}\dfrac{u+2}{u+1}=2$

$\lim_{u\to 0}\left(u^2\times\dfrac{3u+1}{u^4(u+1)}\right)=\lim_{u\to 0}\dfrac{3u+1}{u^2(u+1)}=\infty$

$\therefore$ the singularities at $x=\pm\infty$ are irregular.

So $(x+1)y''+\dfrac{1}{x}y'+(x+3)y=0$ belongs to Heun's Confluent type ODE. If you solve it in MATLAB, MATLAB will express its general solution in terms of Heun's Confluent function.

share|improve this answer

Your Answer

 
discard

By posting your answer, you agree to the privacy policy and terms of service.

Not the answer you're looking for? Browse other questions tagged or ask your own question.