# Help Understanding Sturm-Liouville Boundary Value Problem

Consider the following Sturm-Liouville boundary value problem:

$$\text{Given parameters }c > 0 \text{ and } \beta > 0, \text{ let } y=y(x) \text{ for } 0 \leq x \leq c. \text{ We have }$$
$$y''+ \lambda y=0$$ $$\text{ with boundary conditions }$$ $$y'(0)=\beta y(0),$$ $$y'(c)=\beta y(c),$$

To start, I want to show that this has exactly one negative eigenvalue $$\lambda_0$$. Furthermore this eigenvalue will be independent on the parameter $$c$$ > 0.

Considering the negative eigenvalues $$\lambda < 0$$, let $$\lambda = -\alpha^2$$
This is where I run into trouble. I am not sure what shape we should assume the general solution $$y(x)$$ to have. I know of Sturm-Liouville cases where the general solution we use in the case of the negative eigenvalues is assumed to be $$y(x)=A\cosh(\alpha x)+B\sinh(\alpha x)$$. I also know of Sturm-Liouville problems where the general we use in the case of the negative eigenvalues is assumed to be $$y(x)=Ae^{\alpha x}+Be^{-\alpha x}$$ (I can provide specific details upon request). So it seems that, depending on the boundary conditions, we have to assume $$y(x)$$ to have a certain shape. (I know of more "shapes" that I listed here, I could provide more details if interested.)

So how can I figure out what I should assume $$y(x)$$ to be in this case? After I have done so and found $$y'(0), y'(c), y(0), y(c)$$, where do I go after applying this to the boundary conditions? What exactly am I looking to do to show that there is exactly one negative eigenvalue $$\lambda_0$$ independent of $$c>0$$? After that, how can I go about finding an associated eigenfunction $$y_0(x)$$?

In a similar manner, I am interested in other eigenvalues. How can I determine whether or not $$\lambda = 0$$ is an eigenvalue? How do I determine its associated eigenfunction? The same goes for positive eigenvalues and associated eigenfunctions.

The boundary conditions in this problem are what make it particularly difficult, nothing really "cancels" out nicely with what I tried for $$y(x)$$. Do I have to assume a different form of $$y(x)$$ for every case $$\lambda < 0$$, $$\lambda = 0$$, and $$\lambda > 0$$? Why is this? Does this change the manner in which the associated eigenfunctions are found?

In general, for simple linear ordinary differential equations, a standard ansatz(together with the uniquess theorem for ODEs), tell us that $$y(x)=e^{kx}$$ for some complex $$k$$ (associated with the eigenvalues) span the space of solutions. So all you have to do is to plug it in, and see what that tells you about $$k$$. You get $$(k^2+\lambda)y=0$$ or $$k^2=-\lambda$$.

Now, apply the boundary conditions to $$y(x)=e^{\pm i\sqrt{\lambda}\,x}$$, to obtain $$\pm \sqrt{\lambda}=-i\beta \hspace{3mm}$$ (both boundary conditions yield this equality). You conclude that $$\lambda=-\beta^2$$ is the unique (and negative) eigenvalue for this problem.

As for other eigenvalues, you're free to investigate them but I don't see the point, since you already demonstrated that there can be only one. However, for the fun of it, investigate the zero eigenvalue: $$y''=0\implies y(x)=ax+b$$ for real $$a$$ and $$b$$. The first boundary condition gives us $$a=\beta b$$, and the second one $$a=\beta(ac+b)=\beta b (\beta c+1)$$. So unless $$\beta c+1=1$$ you can have no solution; but $$c$$ and $$\beta$$ are both strictly positive, so zero is not an eigenvalue

In general trying to guess the eigenfunctions and the eigenvalues is really hard. There is a robust branch of functional analysis that dedicates itself, among other things, to that task.

• How do I determine the eigenfunction associated with that negative eigenvalue? Also are you saying that there are no positive eigenvalues (and thus no associated eigenfunctions)? To show this, do we assume y(x)=exp(kx)? Or is there some other form? My collaboration with peers suggests there are positive eigenvalues and associated eigenfunctions. Were they all incorrect? Commented Nov 16, 2023 at 7:11
• @RungSoup indeed the uniquess theorem guarantess that if you find a solution it is unique. The arrived solution, which is of the form $$y(x)=a e^{i\beta x}+b e^{-i\beta x}$$is then the unique solution. This is the desired eigenfunction. In principle, there are no more eigenvalues. What made your colleagues think there were? Commented Nov 16, 2023 at 17:32
• I have a few more questions. First how did you come up with that for the eigenfunction? We assumed initially $y(x)=e^{\pm i\sqrt{\lambda}x}$ but if $\pm\sqrt{\lambda}=-i\beta$, then $\pm i\sqrt{\lambda}=\beta$, so y(x) would be $y(x)=e^{\beta x}$ Why is there a coefficient of i in front? And in a more general sense, why did you initially choose $y(x)=e^{kx}$? As I said in my question, I've seen many different "choosing" for y(x), depending on the problem. Commented Nov 16, 2023 at 23:07
• (2) Why would we not choose $y(x)=A\cosh(kx)+B\sinh(kx)$, for example? Also could you send me some references on this uniqueness theorem. How would I go about finding the normalized eigenfunction? Furthermore I believe the thing that was throwing off my colleagues was the wording of what we were originally going over. That in combination with misunderstanding a strong accent. Commented Nov 16, 2023 at 23:17

The advantage of using the hyperbolic basis functions $$y(x)=A\cosh(αx)+B\sinh(αx)$$ is that the relations to value and derivative at $$x=0$$ are simple, $$y(0)=A,~~~y'(0)=αB.$$ So the first boundary condition removes one degree of freedom by the condition $$αB=βA$$.

At the other boundary one gets of course the full mix of parameters, $$0=y'(c)-βy(c)=α(A\sinh(αc)+B\cosh(αc))-β(A\cosh(αc)+B\sinh(αc)) \\=(αA-βB)\sinh(αc)$$ This indeed is only possible if $$α^2=β^2$$, as neither $$A$$ nor $$B$$ are zero for a non-zero solution. The scale of the solution remains a free parameter, as is usually the case for eigenvectors or eigenfunctions.