The parameter values used in the numerical simulation of equation (2) are as follows: r 1 =0.04 min −1, ... Time-variable control over chemical fluxes by pulse-density modulation.