Empirical modeling of the number densities of star-forming and quiescent galaxies as well as their bulge-dominated subpopulations. The symbols correspond to the measurements, and the curves show the best-fit models as described by the set of differential equations. The bottom panels show the best-fit parametric forms of the quenching rates, integrated SFE, and bulge formation function.