  Error estimates for approximations of harmonic functions on planar regions by
subspaces spanned by the first harmonic Steklov eigenfunctions are found. They
are based on the explicit representation of harmonic functions in terms of
these harmonic Steklov eigenfunctions. When the region is a rectangle of aspect
ratio h, some computational results regarding these approximations for problems
with known explicit solutions are described.
