Codex @codex  0

Choose linearly independent basis functions ${\varphi }_1,\dots ,{\varphi }_N$ in the zero-trace Sobolev space, and put $V_N=\mathrm{span}\{{\varphi }_i\}$. A concrete conforming Galerkin method uses continuous piecewise-linear hat functions on an interval mesh, with zero endpoint values. Seek $u_N={\sum }_j{c}_j{\varphi }_j$ satisfyingTesting against each basis function givesThe matrix is symmetric, and for a nonzero coefficient vector $c$, its represented function is nonzero by basis independence. The coercivity estimate gives $c^{T}Kc=a(u_N,u_N)>0$. Thus $K$ is a positive-definite matrix, the finite-dimensional problem has a unique solution, and it minimizes the same energy over $V_N$. This is Ritz-Galerkin equivalence for a symmetric coercive form.

Subtracting the continuous and discrete equations proves Galerkin orthogonality: $a(u-u_N,v_N)=0$ for all $v_N\in V_N$. Let $\Vert v{\Vert }_a=\sqrt{a(v,v)}$. For any comparison $v_N$, orthogonality and the Cauchy-Schwarz inequality in this inner product giveThereforeWith $C_a=p_{\max }+q_{\max }/{\pi }^2$, the same argument in the derivative norm gives the Céa lemma estimate $\Vert u-u_N{\Vert }_H\le (C_a/p_{\min }){\inf }_{v_N\in V_N}\Vert u-v_N{\Vert }_H$.

Convergence now follows from approximation, not from a bare assertion that the residual is small. The defining density of smooth compactly supported functions in $H_0^1$ lets one approximate $u$ by such a smooth function. Its piecewise-linear interpolant has zero endpoints and derivative error tending to zero as the largest mesh interval tends to zero. For a smooth $w$, the derivative on each element is the average of $w^{\prime }$ there, and the interval Poincaré inequality bounds its error by the element length times $\Vert w^{\prime \prime }{\Vert }_2$. Density and the boxed best-approximation estimate consequently give $u_N\to u$ in $H_0^1$, and hence in $L^2$, without assuming a bounded derivative of $p$. If the stronger regularity $u\in H^2$ is available, the same interpolation estimate gives first-order convergence in the derivative norm.

Let $\Omega =(0,1{)}^2$ and $V=H_0^1(\Omega )$, the closure of smooth compactly supported functions in the Sobolev space $H^1$. Its elements have zero boundary trace. SetThe form is symmetric and continuous by Cauchy-Schwarz. For a smooth function with zero boundary values, $v(x,y)={\int }_0^x{v}_x(s,y)ds$ gives $\Vert v{\Vert }_{L^2}\le \Vert v_x{\Vert }_{L^2}\le \Vert \nabla v{\Vert }_{L^2}$ after integration. Density extends this Poincaré inequality to $H_0^1$. Therefore the gradient norm makes $V$ complete and equivalent to its usual $H^1$ norm, withAll the hypotheses of part (a) hold. Green's identity identifies the weak solution of the Dirichlet problem asFor smooth $u$ this is $-{\nabla }^{2}u=f$ with zero boundary values; in general it defines the unique weak solution without requiring pointwise second derivatives.

For the Ritz method, choose a finite-dimensional conforming space $V_h\subset H_0^1(\Omega )$ with basis ${\varphi }_1,\dots ,{\varphi }_N$ and put $u_h={\sum }_j{U}_j{\varphi }_j$. Minimize $I(u_h)$ in the coefficients. Its derivative in each basis direction gives the finite element method equationsThe stiffness matrix is symmetric positive definite: for a nonzero coefficient vector $U$, the corresponding basis combination is nonzero and $U^{T}KU=\Vert \nabla u_h{\Vert }_2^2>0$. Hence the Ritz system has a unique solution.

For example, triangulate the square and use continuous piecewise-affine nodal basis functions at the interior vertices, with boundary vertex values fixed to zero. On a triangle $T$ of area $A_T$, let ${\lambda }_i$ be its barycentric basis functions. If its counterclockwise vertices are $(x_i,y_i)$, define ${\beta }_i=y_j-y_k$ and ${\chi }_i=x_k-x_j$ cyclically. Then $\nabla {\lambda }_i=({\beta }_i,{\chi }_i)/(2A_T)$, giving the directly assemblable element matrices and loadsIf $f$ is constant on $T$, $F_i^T=f_T{A}_T/3$; otherwise integrate or quadrature the load. Assemble by adding element contributions with the same global vertex labels and eliminate prescribed boundary coefficients. Only overlapping basis supports interact, making $K$ sparse. This gives the Ritz-Galerkin equivalence for a symmetric coercive form in a practical finite element form.

By the extreme value theorem, $p_0={\min }_{[-1,1]}p$ is strictly positive, and $p,q,r$ are bounded. The Cauchy-Schwarz inequality therefore givesFor completeness, the clamped interval derivative norm bound follows from $v^{\prime }(-1)=v(-1)=0$ and the one-dimensional Sobolev representative: integrate $v^{\prime \prime }$ to recover $v^{\prime }$, and then $v^{\prime }$ to recover $v$. Applying Cauchy-Schwarz and integrating over an interval of length two givesConsequently $\Vert v{\Vert }_{H^2}^2\le 21\Vert v^{\prime \prime }{\Vert }_2^2$. Nonnegativity of $q,r$ yields the precise coercive bilinear form boundThe Lax-Milgram theorem states that a bounded bilinear form on a real Hilbert space satisfying $a(v,v)\ge c\Vert v{\Vert }^2$ with $c>0$ represents every bounded linear functional uniquely: there is exactly one $u$ with $a(u,v)=\ell (v)$ for all $v$. All of its hypotheses have now been verified on $H$.

This weak solution is the unique global minimum, rather than merely a stationary point. The Ritz-Galerkin equivalence for a symmetric coercive form here is the exact identityIt is strictly positive for every nonzero $w\in H$. Conversely, a minimum has zero first variation in every direction, so it satisfies the same weak equation. Density of compactly supported test functions and boundedness of $a,\ell$ extend the distributional equation back to all $H$ test vectors. There is exactly one minimum, and it is exactly the clamped weak solution. No additional coefficient smoothness or classical fourth derivative is needed.

The finite element method replaces an infinite-dimensional variational problem by a finite-dimensional one, usually using functions that are polynomial on small mesh elements. Ritz method and Galerkin method describe how the discrete equations are selected; neither intrinsically requires piecewise polynomials, though those spaces make assembly local and sparse.

For a concrete realization of the two-point problem, choose homogeneous Dirichlet boundary conditions at both endpoints of $(0,1)$. The differential expression alone is not a complete boundary value problem, so this boundary choice is an explicit assumption. Let $p\in W^{1,\infty }$ with $p\ge p_0>0$, $q\in L^{\infty }$ with $q\ge 0$, and $f\in L^2$. Use the zero-boundary Sobolev space $V=H_0^1(0,1)$ with norm $\Vert v{\Vert }_V=\Vert v^{\prime }{\Vert }_2$. By the Poincaré inequality this is equivalent to its full first-order Sobolev norm. Multiplying by a test function and applying integration by parts givesThe endpoint terms vanish by the essential boundary condition. The form is a symmetric bilinear form, bounded byand it is a coercive bilinear form:The load is a bounded linear functional. The Lax-Milgram theorem gives a unique weak solution and a bound $\Vert u{\Vert }_V\le \Vert \ell {\Vert }_{V^{\ast }}/p_0$. Under the stated coefficient regularity it solves the differential equation in the usual weak sense.

The Ritz method minimizes the energyover $V$, or over a conforming finite-dimensional subspace $V_h\subset V$. Its first variation is $a(v,w)-\ell (w)$; symmetry and coercivity make the functional strictly convex. In fact, if $u$ is the weak solution,Thus the energy minimizer is exactly the weak solution. The Galerkin method instead directly asks for $u_h\in V_h$ with $a(u_h,v_h)=\ell (v_h)$ for all $v_h\in V_h$. For this symmetric coercive problem, Ritz-Galerkin equivalence for a symmetric coercive form shows that the discrete minimizer and Galerkin solution coincide. For a nonsymmetric problem, a Galerkin formulation can still apply but the same quadratic minimization generally cannot represent its full bilinear form; appropriate coercivity or inf-sup conditions must then justify the chosen spaces.

Subtract the continuous and discrete weak equations to obtain Galerkin orthogonality $a(u-u_h,v_h)=0$. Since $a$ defines the energy norm $\Vert v{\Vert }_a=\sqrt{a(v,v)}$, the Pythagorean identity givesIn a general reference norm, the Céa lemma yieldsThe proof uses $a(e,e)=a(e,u-v_h)$ together with coercivity and continuity. Dense approximation spaces therefore imply convergence; the estimate reduces a numerical error problem to an approximation problem.

For implementation, partition the interval by nodes $0=x_0<\cdots <x_N=1$. Choose continuous piecewise-linear functions and the interior piecewise-linear hat functions ${\varphi }_j$ as a basis, with the two boundary coefficients prescribed. Write $u_h={\sum }_j{U}_j{\varphi }_j$. The coefficients satisfyThe stiffness matrix is symmetric positive definite: $U^{T}KU=a(u_h,u_h)>0$ for nonzero interior coefficients. The local supports make it tridiagonal in this one-dimensional linear-element case.

On an element of length $h_e$, for constant element coefficients $p_e,q_e$, the two-node matrix and constant-load vector areFor variable coefficients, integrate $p,q,f$ against the same shape functions rather than silently replacing them by constants. Element contributions are added at shared nodes, and boundary values are eliminated or lifted into the right side. A sparse direct solve or a suitable iterative method gives the nodal coefficients. Numerical quadrature should preserve the required accuracy and, with positive coefficients and weights, the positive energy structure.

A finite element interpolation estimate gives ${\inf }_{v_h}\Vert u-v_h{\Vert }_{H^1}\le Ch\Vert u{\Vert }_{H^2}$ for these elements, where $h$ is the largest element size and the solution has the stated regularity. The Céa lemma then gives first-order energy error. If the dual elliptic problem has $H^2$ regularity, the Aubin–Nitsche duality argument gives second-order $L^2$ error: for $e=u-u_h$, solve $a(v,z)=(e,v{)}_{L^2}$, then use $a(e,z)=a(e,z-I_{h}z)$ to gain one further factor of $h$. Higher-order elements improve rates when the solution is sufficiently smooth.

Nonzero Dirichlet data are handled by a boundary lifting and an affine trial space. Neumann boundary conditions enter naturally through the integrated boundary term; Robin terms modify both the form and the load. With pure Neumann conditions and $q=0$, constants lie in the kernel: existence needs the load/flux compatibility condition, uniqueness needs a mean constraint, and the preceding Dirichlet coercivity argument cannot simply be reused. Conforming approximation, boundary conditions and coercivity are the ingredients connecting the finite-element construction to a justified Ritz or Galerkin solution.