Rest-frame 35 μm luminosity functions of IR-selected galaxies in redshift bins, in comoving volume units. The black datapoints are from a Spitzer-MIPS sample selected at the effective wavelength of 70 μm by Magnelli et al. (2009), and computed at 35 μm to minimize the K-corrections. Red datapoints are from Gruppioni et al. (2013). The units in ordinate and abscissa are as in Fig. 4. The lines show our multi-population fit to the data (see caption to Fig. 4). The upper black continuous line is the total predicted function best-fitting the luminosity function data. We limit our analysis to z = 1.3 because higher-redshift sources do not contribute significantly to the EBL.