the vortex states in mesoscopic superconducting samples and for vortex arrays in bulk superconductors. These profiles are shown to be strongly affected by the Andreev interference effects and the normal scattering at the sample boundaries. In particular, we find that in contrast to the case of isolated vortices far from the boundaries the positions of peaks in the LDOS profiles in vortex configurations in small size samples are shifted from the phase singularity points. Taking a generic example of the two vortex system, we have considered the evolution of the LDOS profile which accompanies the merging of two vortices and appearance of the doubly quantized vortex. In this case the distance between LDOS peaks reduces slower than the intervortex distance $a$. Finally, when $a$ is smaller than the vortex core size the two-peaked structure of LDOS transforms into an anisotropic ring. The anisotropy vanishes when the vortex positions coincide and a doubly quantized vortex is formed.

Earlier (see Ref. ) it was shown that the spectrum of a single vortex placed near the flat surface is transformed analogously to the vortex-antivortex system and when the distance is of the order of the vortex core size, the interlevel spacing in the 