A finite element method software, FEMALB, has been implemented for the kinetic analysis of complex reaction systems. The established protocol provides fast solutions to coupled differential-algebraic equations. It shows significant advantages over the conventional coding process with the standard implicit Runge-Kutta (IRK) method. The accuracy and high efficiency have been demonstrated by the simulation of reaction processes such as glucose/fructose hydrogenation and catalytic cracking of gasoil. As model validation, the numerical results showed satisfactory agreement with the exact solutions. With the powerful capability of solving large matrices of differential equations (both ODE and PDE) with nonlinear algebraic constrains, such an algorithm has greatly reduced the coding labor in reaction mechanistic studies and provided a unique tool in reactor design and optimization.