Take a smooth variation $u+\epsilon v$ with $v=0$ on the boundary. Differentiating under the integral gives the first variationIntegration by parts removes the derivatives of $v$, and the boundary contribution vanishes. The fundamental lemma of the calculus of variations therefore givesFor a regular constraint $J[u]=c$, the Lagrange multiplier rule says that an extremal satisfies $\delta I=\lambda \delta J$ for a constant $\lambda$. Equivalently, apply the preceding Euler–Lagrange equation to $F-\lambda G$. Regularity means that $\delta J$ is not identically zero as a linear functional on the admissible variations. For the quadratic constraint below this holds at every nonzero admissible $u$; a zero constraint leaves only $u=0$.

In polar coordinates the constrained integrand isIts Euler–Lagrange equation is $-2\lambda ru-{\partial }_r(2r{u}_r)-{\partial }_{\theta }(2u_{\theta }/r)=0$, orThe field is single-valued and periodic in $\theta$. Substitution $u=f(r)\cos (n\theta )$ with $n\in \{0,1,2,\dots \}$ givesThe centre is part of the disk, so one must also select the regular radial solution: for $n\ge 1$, smoothness gives $f(r)=O(r^n)$ and the regular derivative behaviour; for $n=0$, $f$ is bounded and $f^{\prime }(0)=0$. A singular solution at zero does not define an admissible smooth field on the disk.

For fixed $n$, let $c_n={\int }_0^{2\pi }{\cos }^2(n\theta )d\theta$, so $c_0=2\pi$ and $c_n=\pi$ for $n\ge 1$. The angular derivative term has integral $n^2{c}_n$, including zero when $n=0$. ConsequentlyThe one-dimensional Euler–Lagrange equation for $\hat{I}-\lambda \hat{J}$ isexactly the same radial equation. The common factor $c_n$ cancels.

For two regular solutions with the same $n$ and different multipliers, multiply the equation for $f_1$ by $f_2$, multiply that for $f_2$ by $f_1$, subtract and integrate. The result isAt $a$, both functions vanish; at zero, regularity makes the boundary term tend to zero. Since the multipliers are distinct,This is weighted Sturm–Liouville orthogonality with weight $r$.

Finally, multiply one radial equation by $f$ and integrate by parts. The same endpoint conditions yieldMultiplication by $c_n$ provesFor any nonzero admissible solution the energy on the left before subtraction is positive, so its multiplier is the positive Rayleigh quotient $\lambda =\hat{I}/\hat{J}$.

For a continuously differentiable diffeomorphism $\varphi :J\to I$ between intervals, the integration by substitution statesIt holds as an identity element of extended integrals for nonnegative measurable $h$. For signed or complex-valued measurable $h$, integrability on one side is equivalent to integrability after substitution, and the identity element then holds. The absolute derivative accounts for a decreasing substitution; on oriented finite intervals the equivalent formula uses ${\varphi }^{\prime }$ and oriented endpoints.

Put $C={\int }_0^{\infty }(1+x^4{)}^{-1}dx=\pi /(2\sqrt{2})$. The map $x=y^{-1/4}$ is a decreasing diffeomorphism of the positive half-line onto itself andApplying the nonnegative version before assuming finiteness proves both existence and value:Tonelli's theorem permits either iterated integral of a nonnegative measurable function on a product of sigma-finite measure spaces and identifies both with its extended product integral. Fubini's theorem states that for an absolutely integrable measurable function on such a product, almost every section is integrable, the section-integral functions are integrable, and the two iterated integrals equal the product integral.

The given two-variable function is continuous and positive on the positive quadrant. For fixed $y>0$, substituting $t=x{y}^{1/4}$ givesTonelli now givesThis proves absolute integrability, so Fubini is applicable.

For $x\ne 1$, integrate the displayed partial fraction identity element in $y$ up to $R$:At $x=1$ the section integral is ${\int }_0^{\infty }(1+y{)}^{-2}dy=1$, agreeing with the continuous limit of this expression. Therefore Fubini identifies the remaining integral and givesThe apparent singularity at one is removable, with limiting integrand value $1/4$. Positivity and the finite product integral already establish convergence at both ends.

The monotone convergence theorem states that if measurable functions $0\le f_n\uparrow f$ pointwise on a measure space, then $\int f_n\uparrow \int f$, allowing the limit to be infinite.

Extend each restriction to $I_n$ by zero on $I\setminus I_n$. The functions $f{1}_{I_n}$ are nonnegative and increase pointwise to $f$ because the intervals are increasing and exhaust $I$. HenceFor nonnegative $f$, integrability means finiteness of this integral, soNo boundedness or closedness of the intervals is required.