code: def fluxfrac(*mags):
      
      Ftot = 0
      for mag in mags:
          Ftot += 10**(-0.4*mag)
      F1 = 10**(-0.4*mags[0])
      return F1/Ftot
docs: Returns fraction of total flux in first argument, assuming all are magnitudes.