Azimuthal-correlation functions of $\Dzero$ mesons with associated particles, after the subtraction of the baseline, in pp collisions at $\sqrt{s}$ = 13~TeV, for 3 $ \ptD $ 5~\GeVc, 5 $ \ptD $ 8~\GeVc, 8 $ \ptD $ 16~\GeVc, and 16 $\ptD$ 24~\GeVc (from left to right) and $\ptass>$ 0.3~\GeVc, 0.3 $\ptass$ 1~\GeVc, 1 $\ptass$ 2~\GeVc, 2 $\ptass$ 3~\GeVc (from top to bottom) for different multiplicity classes estimated with V0M. The four multiplicity classes are shown with different marker styles. Statistical and $\Delta\varphi$-dependent systematic uncertainties are shown as vertical error bars and boxes, respectively, and $\Delta\varphi$-independent uncertainties are written as text. The uncertainties from the subtraction of the baseline are displayed as boxes at $\Delta\varphi>\pi$.