“…One of the program options enables solution of electrochemical reaction models. The numerical algorithm is based on an implementation of the method of orthogonal collocation on finite elements (OCFE) (Berninger et al, 1991;Gardini et al, 1985;Kaczmarski, 1996;Kaczmarski et al, 1997;Ma & Guiochon, 1991;Villadsen & Michelsen, 1978;Yu & Wang, 1989). After discretization of the space derivatives due to the method of OCFE, the obtained set of ordinary differential equations is solved with the backward differentiation formulae, implemented in the VODE procedure (Brown et al, n.d.).…”