\int_{-\infty}^{+\infty}\frac{N({\bf r},\varepsilon)}{N_0}
  \frac{\partial f(\varepsilon-eV)}{\partial V}\,d\varepsilon\,,$$ where $(dI/dV)_{N}$ is a conductance of the normal metal junction and $f(E)=[1+\exp(E/T)]^{-1}$ is a Fermi distribution function. The typical plots of the differential conductance are shown in Fig. \[dIdV\_M=1\] for different disk radii.

![\[dIdV\_M=1\] (Color online) Three-dimensional plot of the normalized local differential conductance $dI/dV$ as a function of voltage ($eV$) and distance from the vortex center ($r$) in disks with a large radius $R=7\xi$ (a) and small radius $R=2\xi$ (b) for $T/\Delta_0=0.005$. We choose here $k_F\xi=200$.](dIdV.eps){width="1.1\linewidth"}

The zero bias peak at the vortex center is clearly seen for a large disk radius when the normal scattering at the surface can be neglected. The decrease in the disk radius results in the shift of this conductance peak to higher voltages and formation of ring structure of local zero-bias conductance with finite radius $r\sim r_0$.

Note, that normal scattering of quasiparticles at the boundaries can be also important for multiquanta vortex configurations, e.g., vortex pairs, trapped in rather small samples. In our consideration of the LDOS profiles for a vortex-vortex pair in the previous subsection we have neglected the normal scattering at the sample surface assuming the vortex centers to be situated rather far from the 