code: def get_total_flux(self, energies):
          
  
          fluxes = []
  
          for src in self._point_sources:
  
              fluxes.append(self._point_sources[src](energies))
  
          return np.sum(fluxes, axis=0)
docs: Returns the total differential flux at the provided energies from all *point* sources
  
          :return: