Can we construct Sturm Liouville problems from an orthogonal basis of functions? Given a sequence of functions orthogonal over some interval, which satisfy Dirichlet boundary conditions at that Interval, can we construct a Sturm Liouville problem that gives these as its eigenfunctions?

For example, if we take the sequence of scaled Bessel's functions $J_n (\zeta_i x)$ for all positive integral values of $i$ where $\zeta_i$s are roots of the Bessel's function. Note that here, $n$ is fixed and non - negative.    
We already know that these form an orthogonal basis over the weight function $x$ and range $[0,1]$ such that $$\int_0^1 J_n (\zeta x)J_n (\zeta_j x) x dx =\frac12 \delta_{ij}(J_n' (\zeta_i ))^2$$
As the $\zeta_i$s are roots (and the Bessel's functions of the first kind are zero at the origin for positive $n$), we already have $J_n(\zeta_i 0)=J_n(\zeta_i 1)=0$, which is a Dirichlet boundary condition on $[0,1]$.
To me, this seems to be the setup for a Sturm Liouville problem. Can we find $p(x)$ and $q(x)$ such that the functions $y_i(x)=J_n(\zeta_i x)$ satisfy $$\frac{\mathrm d}{\mathrm d x}\left[p(x)\frac{\mathrm dy}{\mathrm d x}\right] + q(x)y(x)=\lambda_i xy(x)$$
Update: Turns out we can, for the example given above. The equation is the radial part of the circular membrane problem, i.e.:
$$x^2y''(x) +xy'(x) +(\lambda^2x^2-n^2)y(x)=0$$
with eigenvalue $\lambda$, where the eigenvalues turn out to be successive roots of $J_n$ and the eigenfunctions are $J_n(\zeta_i x)$ for root $\zeta_i$.
However, can this be done in the general case? If I have a series of functions which:


*

*Form an orthogonal basis over an interval with a certain weight function

*Have Dirichlet BCs on the same interval
can I always construct a Sturm-Liouville problem/ODE for the same wight function and boundary that gives rise to these as eigenfunctions? If not, when is this possible?
 A: The theory of 'integral equations' as described in chapter 3 of Hildebrand's 'Methods of Applied Mathematics' or Hilbert/Courant's book, shows you that the eigenfunctions derive from a Fredholm integral equation 
$$y_n(x) = \lambda_n \int_a^b K(x,z)y_n(z)dz$$
and since the $y_n$ are a basis we have
$$f(x) = \sum_n A_n y_n(x) = \sum_n A_n \lambda_n \int_a^b K(x,z)y_n(z)dz= \int_a^b K(x,z)F(z)dz$$
One can naturally convert a Fredholm integral equation with a symmetric kernel back into a boundary value Sturm-Liouville problem, it's in Hildebrand for boundary conditions more general than just Dirichlet.
To understand this though, we should reverse the process. Given a linear IVP second order ode we can convert it into an equivalent Volterra (Variable endpoint) integral equation (3.2 of Hildebrand). That's a great way to motivate integral equations. If we then try to do the same thing for a boundary value problem, the Fixed endpoint problem gives a Fredholm integral equation, but for a general ode $y'' + ay' + by = f$ we end up with a discontinuous integral kernel. The method motivates the necessity for turning $y'' + ay' + by = f$ into Sturm-Liouville form as the natural way of forming a symmetric integral Kernel that is continuous! This is in Hildebrand 3.3 - 3.4. I think it's more natural to think of Sturm theory as being part of the theory of integral equations rather than differential equations for this reason, as well as the fact all the orthogonality proofs invoke integrals randomly, yet here they arise naturally :)
