$${\mathrm{d}I_f =-\frac{hfB_{21}}{c} (N_1-N_2 ) I_f \mathrm{d}z}$$