this LDOS transformation for small intervortex distances was studied in Ref.  within a perturbation approach. A nonperturbative approach which allows to describe the spectrum transformation accompanying the decay of the multiquantum vortex has been suggested in Refs. [@MelnikovSilaev; @MelnikovPrb2008].

The second factor which is crucial in mesoscopic superconducting samples is a quasiparticle scattering at the sample boundaries which comes into play when vortices approach a superconductor surface or a size of superconducting sample is small enough. For vortices positioned rather close to the sample surface the effect of quasiparticle reflection at the boundary on the spectrum and total DOS was investigated recently in Refs. . The early stage of transformation of the LDOS profiles for vortices approaching a flat boundary has been studied numerically in Ref. .

In the present paper we analyze the effect of both the hybridization of vortex core states and boundary scattering on the LDOS structures and for this purpose consider three generic examples: (i) vortex-vortex pair; (ii) vortex positioned near a flat boundary; (iii) vortex positioned in the center of a superconducting disk. We will show that the LDOS profiles can be obtained qualitatively from the spectra found in Refs. . Our qualitative considerations will 