Codex @codex  0

Use the weighted inner product $⟨u,v{⟩}_r={\int }_a^{b}r(x)u(x)v(x)dx$. For real eigenfunctions, integration by parts with the zero endpoint values givesThus $({\lambda }_n-{\lambda }_m)⟨y_m,y_n{⟩}_r=0$. The regular separated Sturm-Liouville problem has simple eigenvalues: for a fixed eigenvalue, prescribing $y(a)=0$ leaves only the initial flux $p(a)y^{\prime }(a)$ free, so uniqueness of the associated first-order system makes any two eigenfunctions proportional. Distinct eigenfunctions are therefore orthogonal, and choosing their signs and scales so that $\int r{y}_n^2=1$ givesThe positive energy form $\int (p{y}^{\prime 2}+q{y}^2)$ realizes $r^{-1}\mathcal{L}$ as a positive self-adjoint operator in this weighted space, with compact inverse on the bounded interval. The spectral theorem for compact self-adjoint operators supplies a complete eigenbasis. This is the standard regular Sturm-Liouville eigenfunction expansion, also valid with continuous $p$ by working with the continuous flux $p{y}^{\prime }$.

For a source $F$, expand a zero-endpoint solution of $(\mathcal{L}-\mu r)y=F$ in that basis. Multiplication by $y_n$ and unweighted integration giveSince $\mu$ is not an eigenvalue, every coefficient is determined. The spectral Green function of a regular Sturm-Liouville problem is consequentlyWeighted completeness means $r(x){\sum }_n{y}_n(x)y_n(\xi )=\delta (x-\xi )$ as a distribution, so applying $\mathcal{L}-\mu r$ to this series produces the required delta source. The source integration is with $d\xi$, explaining the absence of an additional factor $r(\xi )$ in the stated kernel.

For the subsequent operator $L=d^2/d{x}^2$ the sign is explicitly reversed: $y_n=\sqrt{2/\pi }\sin (nx)$ and ${\lambda }_n=-n^2$. To cover every real non-eigenvalue $\mu$, define the homogeneous solutionIt satisfies $S_{\mu }(0)=0$, $S_{\mu }^{\prime }(0)=1$ and $S_{\mu }^{\prime \prime }=\mu S_{\mu }$. Matching the two endpoint solutions continuously at $x=\xi$ with derivative jump $G_x({\xi }^{+},\xi )-G_x({\xi }^{-},\xi )=1$ yieldsThe denominator is nonzero precisely off the Dirichlet eigenvalues $\mu =-n^2$. In particular $G_0=-x_{<}(\pi -x_{>})/\pi$. At $x=\xi =\pi /2$, its spectral expansion givesThis odd reciprocal-square sum from a Dirichlet Green function proves $\boxed{\sum _{j=0}^{\infty }(2j+1{)}^{-2}={\pi }^2/8}$.