  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.
