Dust mass function derived using the V_MAX method for dust masses of the galaxies computed assuming T_D = 25K (cyan squares and black circles with errors, computed following Gehrels 1986, upper panel), T_D = 35K (magenta squares and black circles with errors, central panel) and using the dust mass from CIGALE (red squares and black circles with errors, lower panel). The best fit, computed using the A³COSMOS data points, is displayed as a solid red line, with shaded error bands of the same color. For this fit, M_D^*, Φ_D^*, and α were free to vary. For comparison, different estimates from the literature are reported. The light blue circles and squares are the values obtained by Pozzi et al. (2020), and the dashed cyan lines correspond to the best fit. The dashed gray curves are from Dunne et al. (2003). The blue circle show the estimate by Magnelli et al. (2019), and the dotted line shows the ALPINE dust mass function by Pozzi et al. (2021), assuming T_D = 25K.