Boundary Value Problems for Second Order Differential Equations

[fluidistic]

Homework Statement

1)Calculate the general solution of y''+\frac{y'x}{1+x}-\frac{y}{1+x}=1+x.
2)What behavior do the solutions have in x=-1?
3)Solve the boundary problem to the DE in the interval [0,1] with y(0)=0 and y(1)+y'(1)=0.
4)Write down the DE under Sturm-Liouville's form and find a Green function for the DE with the boundary conditions in 3) over the interval [0,1]. Then solve the DE using Green's function.
Hint:Use solutions y_1 and y_2 such that y_1(0)=0 \neq y_2 (0) and y_2 (1)+y _2 '(1)=0 \neq y_1 (1)+y_1 ' (1).

Homework Equations

For part 1, not sure.

The Attempt at a Solution

I tried to apply Frobenius method because of the singularity of the DE in x=-1.
I reached \sum _{n=0}^ \infty (n+c)(n+c-1)a_n(x-1)^{n+c+2}+ \frac{x}{1+x}\sum _{n=0}^ \infty (n+c)a_n (x-1)^{n+c-1}- \frac{1}{1+x} \sum _{n=0}^ \infty a_n (x-1)^{n+c}=1+x.
Noting really simplifies well, if I continue further I think I'll have problems in finding a recurrence relation that could be useful. I wonder if I'm not doing things wrong so far and if I should continue this way.

Edit: this is not going well at all, my right hand side is different from 0.

 

[DryRun]

For part (1), i think the forcing function or free term should be zero.

The first step would be to write up the auxiliary equation and solve.

I get: m=-1 \; or\; m=\frac{1}{x+1}

 

[fluidistic]

For part (1), the forcing function or free term should be zero.

Write up the auxiliary equation and solve.

Thanks for helping me. I never met such terms before. If google helped me well, you mean that there's a function of only the variable x that should be 0?
What I reached is (x-1)^{c-2}a_0c(c-1)+(x-1)^{c-1}a_1 (c-1)c +\frac{x}{1+x} (x-1)^{c-1} c a_0 + \sum _{n=0}^{\infty } (x-1)^{n+c} [a_{n+2}(n+c+2)(n+1+c)+\frac{x}{1+x} (n+c+1)a_{n+1}- \frac{a_n}{x+1}] = 1+x. I don't know where is the free term nor why a term should be 0.

 

[fluidistic]

I haven't made much progress on this DE. I tried several ways, even with a friend. He told me apparently that the solution can only be kept as an integral though the integral is so horrible that even the basic special functions wouldn't help.
On my own I've attempted to solve the homogeneous DE y''+\frac{y'x}{1+x}-\frac{y}{1+x}=0 via infinite series expansions.
I reached that y_1 (z)= \sum _{n=0} ^{\infty } \frac{(-1)^nz^n}{n!} and y_2 (z)= \sum _{n=0} ^{\infty } \frac{(-1)^n z^{n+2}}{(n+1)!}. I do not recognize any function I know in these 2 linearly independent solutions. By the way here z=x+1.
Now in order to find the particular solution to the general DE y''+\frac{y'x}{1+x}-\frac{y}{1+x}=1+x I must calculate the wronskian of y_1 and y_2, which seems way too aweful to even be attempted.
However according to wolfram alpha, the general solution is of the form y(x)=Ae^{-x}+Bx+Cx^2+D.
Hmm now that I think about it, my solution y_1(z)=e^{-z}. But I'm still unsure about my y_2(z). Any help is immensely appreciated on how to proceed further.

 

[DryRun]

For part (1): The general solution of the reduced equation is: y_c = Ae^{x/(1+x)}+Be^{-x}
The Wronskian is: \frac{x^2+2x+2}{(1+x)^2}

 

[fluidistic]

Ok thanks a lot, it means I've got one solution to the homogeneous problem right; the other wrong. I've spotted one error but even when corrected, the other solution to the DE doesn't match Ae^{x/(1+x)}.
Thus I fear I have to show every single steps I've made.
Let z=x+1 so that the orignal DE becomes y''+ \frac{(z-1)y'}{z}-\frac{y}{z}=0. I propose a solution of the form y(z)= \sum _{n=0}^{\infty } a_n z^{n+c} \Rightarrow y'(z)=y''(z)=\sum _{n=0} ^{\infty } (n+c) a_n z^{n+c-1} \Rightarrow \sum _{n=0} ^{\infty } (n+c) (n+c-1) a_n z^{n+c-2}.
Now I replace those series into the DE:
\sum _{n=0} ^{\infty } (n+c) (n+c-1) a_n z^{n+c-2}+\frac{(z-1)}{z}\sum _{n=0} ^{\infty } (n+c) a_n z^{n+c-1}- \frac{1}{z} \sum _{n=0}^{\infty } a_n z^{n+c}=0.
Developping a bit:
\sum _{n=0} ^{\infty } (n+c) (n+c-1) a_n z^{n+c-2}+\sum _{n=0}^{\infty } (n+c) a_n z^{n+c-1} - \sum _{n=0}^{\infty } (n+c)a_nz^{n+c-2}- \sum _{n=0} ^{\infty } a_n z^{n+c-1}=0.
Now I set every series in powers of z^{n+c-2}.
\sum _{n=0} ^{\infty } z^{n+c-2} a_n [(n+c)(n+c-1)-(n+c)] + \sum _{n=1}^{\infty } z^{n+c-2} (n+c-1)a_{n-1}-\sum _{n=1}^{\infty } z^{n+c-2}a_{n-1}=0
So \sum _{n=0}^{\infty } z^{n+c-2}a_n (n+c)(n+c-2)+ \sum {n=1} ^{\infty } z^{n+c-2} a_{n-1}(n+c-2)=0.
The secular equation is thus z^{c-2} a_0 c(c-2)=0 with roots c=0 or c=2.
Continuing the work on the infinite series, I get that \sum _{n=1}^{\infty } z^{n+c-2} [a_n (n+c)(n+c-2)+a_{n-1}(n+c-2)]=0.
Hence the reccurrence relation I get is a_n (n+c)(n+c-2)+a_{n-1}(n+c-2)=0.
First case: c=0:
a_n=\frac{-a_{n-1}}{n}.
Choosing a_0=1, I get a_1=-1, a_2=\frac{1}{2}, a_3=-\frac{1}{2 \cdot 3}. In general I notice the pattern a_n = \frac{(-1)^n}{n!} so that y_1 (z)= \frac{(-1)^n z^n}{n!}=e^{-z}. This is a good sign because you obtained the same solution.
Second case: c=2:
a_n=- \frac{a_{n-1}}{n+2}.
Choosing a_0=1, I get a_1=- \frac{1}{3}, a_2 = \frac{1}{3 \cdot 4 }, a_3=- \frac{1}{3 \cdot 4 \cdot 5 }.
In general I notice the pattern a_n = \frac{(-1)^n \cdot 2 }{(n+2)!} \Rightarrow y_2 (z) = \sum _{n=0} ^{\infty }\frac { 2 \cdot (-1)^n z^{n+2}}{(n+1)!}.
Hmm now that I think about it, I don't really see how it's possible to reach a solution of the form e^{-\frac{1}{z}} (this is what I should reach according to your solution) because this would mean a different form than powers of z^{n+c}, in fact this would require infinitely many negative powers of z. I'm at a loss at understanding what's going on.

 

[fluidistic]

I just read some theory and I've "learned" something.
If the indicial equation leads to two roots which differ by an integer then I cannot assume that (assuming that c_1<c_2) y_2(z)= \sum _{n=0}^{\infty } a_n z^{n+c_2}. This explains why my solution for y_1(z) is right but not the one for y_2(z).
In order to find y_2(z), I used the method which starts by assuming that y_2(z)=v(z)y_1(z).
After a lot of algebra I reached that y_2(z)=c_1(z-1)+c_2e^{-z}.
Which implies that the complentary solution to the homogeneous DE is y_c (z)=c_1e^{-z}+c_2(z-1).
From this fact, I had some fun and checked out the validity of Abel's theorem.
It only remains to find the particular solution to the DE y''+ \left ( \frac{z-1}{z} \right ) y' - \frac{y}{z}=z.
I used variation of parameters but still haven't found out the solution (I think I should reach y_p(z)=z^2-2z+2 which is equivalent to x^2+1). I reach an extremely complicated function given by wolfram; I'll recheck my algebra.
I think your solution y_c = Ae^{x/(1+x)}+Be^{-x} is equivalen to Ae^{-\frac{1}{z}}+Be^{-z} which isn't equivalent to my solution. Is something I've said wrong?