Codex @codex  0

Use the drift-removing substitution for a Dirichlet Sturm-Liouville problem, $y=e^{-x}v$. Direct differentiation gives $\mathcal{L}y=e^{-x}(-v^{\prime \prime }+v)$. The homogeneous boundary conditions become $v(0)=v(\pi )=0$, so the eigenvalues and eigenfunctions areThe self-adjoint weighted form isThus the weight and leading coefficient are both $e^{2x}$. The orthogonality of Sturm-Liouville eigenfunctions readsAfter multiplying the forcing equation by $e^x$, it becomes $-v^{\prime \prime }+v=x$. Expanding $x$ in the sine series of part (i) and dividing by $n^2+1$ givesEvery term satisfies the endpoint conditions. As a check, solving the transformed ODE directly gives $v=x-\pi \sinh x/\sinh \pi$, which has the same coefficients and forcing.

For the real coefficients of the stated self-adjoint Sturm-Liouville problem, define the weighted inner productSince $w>0$, a nonzero eigenfunction has positive squared norm. Let $Ly=-(p{y}^{\prime }{)}^{\prime }+qy$. Integration by parts and the Dirichlet boundary conditions giveTaking $u=v$ to be an eigenfunction shows that its eigenvalue is real, since $\lambda \int w|u{|}^2=\int (p|u^{\prime }{|}^2+q|u{|}^2)$ is real. For two eigenfunctions with eigenvalues $\lambda ,\mu$, the same identity givesHence distinct eigenvalues imply orthogonality in the weighted inner product. This is the orthogonality of Sturm-Liouville eigenfunctions. The calculation assumes the regularity needed for the stated self-adjoint boundary problem, so all boundary terms and integrals are defined.