Title: An O(h6) cubic spline interpolating procedure for harmonic functions