A finite element method software FEMALB has been implemented to the kinetic analysis of complex reaction systems. The established protocol provides fast solutions to the 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 in the simulation of the 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 matrixes 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.