Given S-L equation $\dfrac{1}{x}[\dfrac{d}{dx}(xy')+(\dfrac{-m^2}{x})y]=-\lambda y$
Say $\mathcal{L}$ is the Sturm-Liouville operator, $y_k$ is eigenfunction $J_m(j_{mk}x)$ where $J_m$ is Bessel function and $j_{mk}$ is kth zero of $J_m(x)$. Associated eigenvalues $\lambda_k=(j_{mk})^2$. Let $y(x,\lambda)=J_m(\sqrt{\lambda}x)$.
How can I evaluate $(\mathcal{L}y_k)y-y_k(\mathcal{L}y)$ in two different ways, then form a suitable integral to show $\int^1_0 xJ_m(\sqrt{\lambda}x)J_m(\sqrt{\lambda_k}x)dx=\dfrac{\sqrt{\lambda_k}J'_m(\sqrt{\lambda_k})J_m(\sqrt{\lambda})}{\lambda-\lambda_k}$
One way is just a straight computation. i.e $$( \mathcal{L} y_k )y - y_k ( \mathcal{L}(y)) = \frac{y}{x} \left [ (xy_k')' - \frac{m^2 y_k}{x} \right ] -\frac{y_k}{x} \left [ (xy')' - \frac{m^2 y}{x} \right ] =\frac{W_{1/x}[y,y_k]'}{x}$$ where $W_p[y_1,y_2]$ is the $1/p$ weighted Wronskian. The other is to notice that $$ ( \mathcal{L} y_k )y - y_k ( \mathcal{L}(y)) = (\lambda - \lambda_k) y y_{k} $$ Thus, the identities together give us that $$ x y y_k = \frac{ W_{1/x} [y,y_k]' }{\lambda- \lambda_k} $$