In this paper, we address the problem of complete description of coexisting stationary localized modes for the system of n nonlinearly coupled Schrödinger-type equations. Our approach is based on the observation that if the nonlinearity is defocusing (repulsive), then generic solutions of the corresponding stationary system have singularities at finite points of the real axis and, consequently, cannot describe profiles of stationary modes. Then the stationary localized modes can be found by a procedure of “filtering out” solutions with singularities. We have formulated and proved three rigorous statements that form the mathematical background of the approach. The first of them, Theorem on parametrization, states that there exists a homeomorphism between S and Rn, where S is the set of solutions of the stationary system that decay to zero componentwise at +∞. This allows to code the solutions from S by n-component vectors. The second statement, Theorem on singularity, states that, under certain conditions, the singular solutions are ubiquitous in the system and gives the conditions for detection of singularities. This statement allows to seek for localized modes by properly arranged numerical scanning in S operating with vectors from Rn. The third theorem, Theorem on collapsing solutions, allows to reduce the scanning procedure to a finite domain in Rn. As a result, if the conditions of the three theorems hold, the scanning procedure yields the complete set of localized modes for the system. We illustrate this method by an examples of two and three nonlinearly coupled equations that describe the steady states of two and three- component Bose-Einstein Condensate. As the main outcome, in the chosen domain of chemical potentials, we have found all coexisting nonlinear modes.