Infrared luminosity function found in this work at different redshift bins, as indicated in the plot. The faint-end data (Table 1), determined through stacking (see Section 3.1.1 for details), are depicted by black points with horizontal error bars. The remaining bright-end data (Table 3) found using the ALMA follow-up data of the S2CLS UKIDSS UDS sources (AS2UDS; Stach et al. 2019) are also shown (Section 3.1.2). The best-fit Schechter functions of Equation (6) are represented by dashed black lines, with the gray area representing 1σ uncertainties. The faint-end slope was determined at two low redshift bins, and the more accurate result, α = −0.26 ± 0.11, was adopted at all the remaining redshifts, as explained in Section 3.1.3.