On the inevitable error of the method of nets

It is proved that no matter what the solution of an arbitrary boundary-value problem for the two-dimensional Laplace equation, unless it is a special fourth-degree harmonic polynomial, the rate of convergence of the method of square nets using the operator for computation of the four-point arithmetic mean can never be better than $h^2$ (where $h$ is the spacing of the net).