We prove regularity estimates up to the boundary for solutions of elliptic systems of finite difference equations. The regularity estimates, obtained for boundary-fitted coordinate systems on domains with smooth boundary, involve discrete Sobolev norms and are proved using pseudo-difference operators to treat systems with variable coefficients. The elliptic systems of difference equations and the boundary conditions which are considered are very general in form. We prove that regularity of a regular elliptic system of difference equations is equivalent to the nonexistence of "eigensolutions". The regularity estimates obtained are analogous to those in the theory of elliptic systems of partial differential equations, and to the results of Gustafsson, Kreiss, and SundstrSm [1972] and others for hyperbolic difference equations.