code: def compute_blackbody_response(self, Teffs=None):
          
  
          if Teffs is None:
              log10Teffs = np.linspace(2.5, 5.7, 97) # this corresponds to the 316K-501187K range.
              Teffs = 10**log10Teffs
  
          # Energy-weighted intensities:
          log10ints_energy = np.array([np.log10(self._bb_intensity(Teff, photon_weighted=False)) for Teff in Teffs])
          self._bb_func_energy = interpolate.splrep(Teffs, log10ints_energy, s=0)
          self._log10_Inorm_bb_energy = lambda Teff: interpolate.splev(Teff, self._bb_func_energy)
  
          # Photon-weighted intensities:
          log10ints_photon = np.array([np.log10(self._bb_intensity(Teff, photon_weighted=True)) for Teff in Teffs])
          self._bb_func_photon = interpolate.splrep(Teffs, log10ints_photon, s=0)
          self._log10_Inorm_bb_photon = lambda Teff: interpolate.splev(Teff, self._bb_func_photon)
  
          self.content.append('blackbody')
          self.atmlist.append('blackbody')
docs: Computes blackbody intensities across the entire range of
          effective temperatures. It does this for two regimes, energy-weighted
          and photon-weighted. It then fits a cubic spline to the log(I)-Teff
          values and exports the interpolation functions _log10_Inorm_bb_energy
          and _log10_Inorm_bb_photon.
  
          @Teffs: an array of effective temperatures. If None, a default
          array from ~300K to ~500000K with 97 steps is used. The default
          array is uniform in log10 scale.
  
          Returns: n/a