 Solution

Write $N[\psi ]={\int }_a^b(p{\psi }^{\prime 2}+q{\psi }^2)dx$ and $D[\psi ]={\int }_a^{b}w{\psi }^{2}dx>0$. For a real variation, with $\Lambda =N/D$,Suitable endpoint conditions are homogeneous Dirichlet boundary conditions, homogeneous Neumann boundary conditions $p{\psi }^{\prime }=0$, or one of each. They remove the variation boundary term and the product $[p\psi {\psi }^{\prime }{]}_a^b$. Stationarity for all interior variations is equivalent toConversely, multiplication of an eigenfunction equation by $\psi$ and integration by parts gives $N=\lambda D$, hence $\boxed{\Lambda [\psi ]=\lambda }$. General Robin data require their boundary-energy contributions unless these endpoint products happen to vanish.

Use an eigenbasis orthonormal with respect to the weight: ${\int }_a^{b}w{\psi }_i^{\ast }{\psi }_{j}dx={\delta }_{ij}$. For complex coefficients use the Hermitian version of the quotient, with absolute squares. Self-adjointness givesWhen $c_0\ne 0$, dividing numerator and denominator by $|c_0{|}^2$ yieldsIf $c_0=0$, these ratios are undefined, but the unscaled formula and upper-bound conclusion remain valid. The assumptions are a problem with a self-adjoint operator bounded below and a complete eigenbasis, with the trial function in its energy-form domain. The Rayleigh-Ritz method minimizes this quotient over an admissible trial family, obtaining an upper bound on ${\lambda }_0$. Higher eigenfunctions are stationary directions and need not be unrestricted extrema.