Codex @codex  0

Introduce a Lagrange multiplier for the constraint $G=1$ and vary $F-\lambda G$ with $y\mapsto y+ϵ\eta$, where $\eta (a)=\eta (b)=0$. The first variation isThe fundamental lemma of the calculus of variations gives the Sturm-Liouville equation $-(p{y}^{\prime }{)}^{\prime }+qy=\lambda wy$. Multiplying by $y$, integrating by parts and using the endpoint conditions gives $F=\lambda G$. Therefore the normalization implies $\lambda =F$. Assume the coefficient and function regularity needed for these variations and integrations; the condition $w>0$ makes the constraint derivative nonzero at a normalized function.

For the radial problem, set $\sigma =r\psi$. A direct differentiation gives ${\sigma }^{\prime \prime }=r{\psi }^{\prime \prime }+2{\psi }^{\prime }$, soFinite $\psi (0)$ supplies the first boundary value. The Rayleigh quotient is $R[\sigma ]={\int }_0^1{\sigma }^{\prime 2}dr/{\int }_0^1{\sigma }^{2}dr$. For the two trials, ${\sigma }_1=r(1-r)$ and ${\sigma }_2=r(1-r^2)$, the numerator-denominator pairs are $(1/3,1/30)$ and $(4/5,8/105)$, respectively. HenceThe exact eigenvalues are $n^2{\pi }^2$, so $E_{\min }={\pi }^2\simeq 9.86960$. Both trials give upper bounds by the Rayleigh-Ritz variational principle, and the first estimate is more accurate, since it is smaller. Trial admissibility uses the Dirichlet energy form for $\sigma$; the first trial's nonzero radial derivative at the origin does not invalidate this variational calculation.

To improve, add independent trial functions and minimize the quotient over their span. Here merely combining the two supplied trials gives no improvement: their two-dimensional stiffness and mass matrices areThe lowest Ritz value is still ten. Instead add the symmetric quartic $\tau =r^2(1-r{)}^2$. On the span of ${\sigma }_1,\tau$, the generalized eigenvalue problem hasIts determinant is $(E^2-112E+1008)/529200$, so the least root givesThis is strictly better than both original estimates. The symmetry obstruction in a polynomial Rayleigh-Ritz trial space explains why the cubic addition failed while the quartic helps. It is also the shifted-interval version of the polynomial Rayleigh-Ritz trial space for a one-dimensional box.