It is shown that for the two dimensional Laplace equation a univariate cubic spline approximation in either space direction together with a difference approximation in the other leads to the well-known nine-point finite-difference formula.
For harmonic problems defined in rectangular regions this property provides a means of determining with ease accurate approximations at any point in the region
