docs: def absolute_magnitude(distance_modulus,g,r,prob=None):
      """ Calculate the absolute magnitude from a set of bands """
      V = g - 0.487*(g - r) - 0.0249
          
      flux = np.sum(10**(-(V-distance_modulus)/2.5))
      Mv = -2.5*np.log10(flux)
      return Mv