Time series expansion to find solutions of nonlinear systems: an application to enzymatic reactions