 Maximum-norm convergence of the corrected nine-point Poisson scheme

On the unit square with exact Dirichlet boundary data, let $D_{9}U=B_{h}f$ and suppose the normalized residual on the exact smooth solution has supremum norm at most $\epsilon$. The error satisfies $D_{9}e=-\tau$ and vanishes on the boundary. The barrier $q=[x(1-x)+y(1-y)]/4$ has $D_{9}q=-1$, $q\ge 0$ on the boundary, and maximum $1/8$. Apply the discrete maximum principle to $e-\epsilon q$ and $-e-\epsilon q$, obtaining $|e|\le \epsilon q$. Thus $\Vert e{\Vert }_{\infty }\le \epsilon /8$, and fourth-order residual gives fourth-order supremum norm nodal convergence. The proof needs a uniform smoothness bound and a compatible boundary discretization.