A numerical method is used to characterize the steady-state voltammetry at microdisk electrodes of a new electrocatalytic reaction. This reaction, which occurs for the oxidation of N,N-dimethylphenylenediamine (DMPD, A) in the presence of H2S (X) is believed to proceed via the following route: A - 2e(-) reversible arrow B, B + X --> B-X (k(2)), B-X - 2e(-) reversible arrow [B-X](2+). Due to the presence of a reagent restricted homogeneous kinetic step, the reaction is labeled EC2XE. The numerical method for simulating this reaction scheme is based on the finite-difference formulation of coupled mass transport and kinetic equations in oblate spherical coordinates. The method is illustrated for not only the EC2XE but also the EC' reaction and is applicable to the simulation of steady-state limiting currents at microdisk electrodes. Iterative solutions are calculated using a Gauss-Newton scheme to overcome nonlinear homogeneous kinetic terms. The spatial convergence of the simulation for both reactions is investigated by considering the form of the concentration function describing the species. Via the comparison of working surfaces generated from simulated results, measuring the steady-state limiting current is shown to be insensitive to the resolution of EC', ECE, and EC2XE reactions. Experimental steady-state limiting current data is reported for the DMPD/H2S system at microelectrodes of 7.3, 19.5, and 25.0 mum diameter to verify the theory behind the EC2XE reaction. These results are shown to closely fit experimental data using a working surface interpolation method. Specifically, this method correctly predicts the variation of the steady-state limiting current with the concentration of H2S for a 19.5 mum diameter microelectrode to a relative standard deviation of 1.9%. Similar analysis for the 7.3 and 25.0 mum electrodes results in a mean value of 1.4 x 10(7) mol(-1) cm(3) s(-1) for the rate constant k(2) in the DMPD/H2S system.