The previously suggested method-oriented approach to the construction of algorithms of the numerical solution of initial boundary value problems occurring in electrochemical kinetics is extended to the class of kinetic problems characterized by the simultaneous presence of bulk species (distributed in the electrolyte) and interfacial species (localized at the electrodes). Such problems are described by coupled sets of partial differential and differential-algebraic equations, which in particular cases may also take the form of the sets of partial differential and algebraic equations, or partial differential and ordinary differential equations. Seven algorithms based on finite-difference discretizations are formulated, and their accuracy, efficiency and numerical stability are examined in calculations using typical examples of kinetic problems. An improvement of the coding of text input formulae and run-time interpretation of the codes, based on the idea of the three-address code in quadruple form, is also described.