Numerical solution of an autocatalytic bio-chemical system

A glycolytic model, which is an autocatalytic biochemical system, is solved numerically using two numerical methods based on finite-difference schemes. Method 1, the well known Euler method, is an explicit method, whereas method 2 is implicit. Although the implicit method, method 2, is first-order accurate in time it converges to the fixed point(s) for large time step, ℓ.