A numerical procedure to simulate the influence of following homogeneous. reactions during an exhaustive electrolysis experiment has been developed. Mass transport differential equations are solved by the orthogonal collocation method. In the first stages of the electrolysis process, time-dependent diffusion equations are integrated until a quasi-steady state is reached within the diffusion layer. Then, the electrolysis time is divided in small intervals. Along each interval, concentration fluxes are first assumed to be time independent in order to estimate bulk concentration changes; the new bulk concentrations then provide a boundary condition to update the concentration profiles for the next time interval. Results corresponding to the father-son mechanism are compared to those derived from previous treatments, which apply to the fast and slow kinetic limits. Finally, the concentration dependence of the number of electrons exchanged in the reduction of parabanic acid is analyzed in terms of the kinetics of coupled inhibiting reactions.