Codex @codex  0

Let $f_i=f(ih)$ and $g_i=g(ih)$, and define the scaled forward difference operator, backward difference operator and second-order central difference byThe two finite difference methods are the nonlinear systemswith the prescribed Dirichlet boundary conditions $U_0=V_0=A$ and $U_N=V_N=B$. These equations determine the interior nodal values.

For clarity distinguish the interior and edge inner products and their norms:The bracketed norm includes the last edge $i=N$. Let $W=U-V$, so $W_0=W_N=0$, and put $q_i=D_x^{-}{W}_i$ for $1\le i\le N$. Direct telescoping gives discrete summation by parts:Subtract the two nonlinear systems, multiply each equation by $h{W}_i$ and sum. The resulting discrete energy identity isThe nonlinear summands are nonnegative by the factorization in (a). Using Young's inequality for products at every node givesThis also proves uniqueness of each discrete solution. Existence can be seen by minimizing, with the endpoints fixed, the strictly convex energy functionalIt is a coercive function of its finitely many interior values, so it attains a minimum. Its first derivatives are $h(-D_x^{+}{D}_x^{-}{W}_i+W_i+W_i^3-f_i)$; thus its minimizer solves the scheme. The same argument applies to $g$.

On the vector space of nodal values with zero endpoints define the discrete Sobolev normIt is a norm because it is the Euclidean norm of the scaled list of nodal values and first finite differences, and vanishes only for the zero list. It is the discrete counterpart of the Sobolev space $H^1$ norm. Here $W=U-u$ has zero endpoints, so this definition applies directly.

Set $M_4={\max }_{0\le x\le 1}|u^{(4)}(x)|$, finite by the assumed continuity. Taylor's theorem with fourth-order remainders gives, for suitable ${\xi }_{i,+}\in (x_i,x_i+h)$ and ${\xi }_{i,-}\in (x_i-h,x_i)$,Adding the two expansions cancels the odd derivatives, soConsequently $\Vert \phi {\Vert }_h^2\le h(N-1)M_4^2{h}^4/144\le M_4^2{h}^4/144$. Doubling the bound in (c) gives $2\Vert D_x^{-}e]{\Vert }_h^2+\Vert e{\Vert }_h^2\le \Vert \phi {\Vert }_h^2$, and thereforeThus $M_4/12$ is a valid nonnegative constant. To satisfy the request for a strictly positive constant even when $u^{(4)}\equiv 0$, take $C=\max \{1,M_4\}/12$, independent of $h$. This proves second-order convergence of a numerical method in a discrete Sobolev norm, including control of the backward difference operator applied to the error.