be confirmed by the detailed numerical calculations.

This paper is organized as follows. In Sec. \[sec:theory\] we give an overview of the theoretical framework which is employed in this work, namely the Bogoliubov–de Gennes theory and the quasiclassical Eilenberger approach. We discuss the spectrum and LDOS patterns for multi vortex configurations in Sec. \[subsec:molecules\] and address the case of a vortex in a small mesoscopic cylinder in Sec. \[subsec:cylinder\]. We give our conclusions in Sec. \[sec:summary\].

Model and basic equations {#sec:theory}
=========================

Our further consideration is based on two approaches: (i) Bogoliubov–de Gennes (BdG) equations; (ii) quasiclassical Eilenberger theory. The first one appears to be more transparent for the qualitative analysis of the interference effects and convenient for numerical calculation of LDOS profiles around the vortex placed in a superconducting disk. The Bogoliubov–de Gennes (BdG) equations for particle– ($u$) and hole–like ($v$) parts of the wave function have the following form: $$\begin{aligned}
\label{BdG}
  \nonumber
  -\frac{1}{2m}\left(\hat {\bf p}-\frac{e}{c}{\bf A}\right)^2u+\Delta v
  &=&(\varepsilon+\varepsilon_F) u\,,\\
  \frac{1}{2m}\left(\hat {\bf p}+\frac{e}{c}{\bf A}\right)^2v
  +\Delta^* u&
  =&(\varepsilon-\varepsilon_F) v\,.\end{aligned}$$ Here $\Delta$ is a gap function, ${\bf A}$ is a vector potential, $\hat
{\bf p}=-i\hbar(\partial/\partial x,\partial/\partial y)$, and ${\bf
r}=(x,y)$ is a radius vector in the plane perpendicular to the vortex 