This note is inspired by the paper of Bialecki, Fairweather, and Bennett [1] which studies a method for the numerical solution of u00 = f on an interval, and u = f on a rectangle, with zero boundary data; their scheme calculates that C piecewise cubic U for which U 00 (or U) agrees with f at the two (or four) Gauss quadrature points of each subdivision of the original domain. They observe that ...