On the monotonicity and discrete maximum principle of the finite difference implementation of C⁰-Q² finite element method

We show that the fourth order accurate finite difference implementation of continuous finite element method with tensor product of quadratic polynomial basis is monotone thus satisfies the discrete maximum principle for solving a scalar variable coefficient equation $-\nabla\cdot(a\nabla u)+cu=f$ under a suitable mesh constraint.