We construct necessary and sufficient conditions for the existence of solution of seminonlinear boundary value problem for a parametric excitation system of difference equations. The convergent iteration algorithms for the construction of the solutions of the semi-nonlinear boundary value problem for a parametric excitation system difference equations in the critical case have been found. The investigation of periodic and Noetherian boundary-value problems in the critical cases is traditionally performed under the assumption that the differential equation and boundary conditions are known and fixed. As a rule, the study of periodic problems in the case of parametric resonance is reduced to the investigation of the problems of stability. At the same time, due to numerous applications in electronics, geodesy, plasma theory, nonlinear optics, mechanics, and machine-building, the analysis of periodic boundary-value problems in the case of parametric resonance requires not only to find the solutions but also to determine the eigenfunctions of the corresponding differenсе equation. The investigation of autonomous Noetherian boundary-value problems is also reduced to the study of Noetherian boundary-value problems in the case of parametric resonance because the change of the independent variable in the critical case gives a nonautonomous boundary-value problem with an additional unknown quantity. The aim of the present paper is to construct the solutions of Noetherian boundary-value problems in the case of parametric resonance whose solvability is guaranteed by the corresponding choice of the eigenfunction of the analyzed boundary-value problem. The applied classification of Noetherian boundary-value prob\-lems in the case of parametric resonance depending on the simplicity or multiplicity of roots of the equation for generating constants noticeably differs from a similar classification of periodic problems in the case of parametric resonance and corresponds to the general classification of periodic and Noetherian boundary-value problems. The equation for generating constants obtained for the Noetherian boundary-value problems in the case of parametric resonance strongly differs from the conventional equation for generating constants in the absence of parametric resonance by the dependence both of the equation and of its roots on a small parameter, which leads to noticeable corrections of the approximate solutions as compared with the approximations obtained by the Poincare method. Using the convergent iteration algorithms we expand solution of seminonlinear two-point boundary value problem for a parametric excitation Mathieu type difference equation in the neighborhood of the generating solution. Estimates for the value of residual of the solutions of the seminonlinear two-point boundary value problem for a parametric excitation Mathieu type difference equation are found.