The main purpose of this paper is to investigate the solution of general reaction–diffusion glycolysis system numerically. Glycolysis model demonstrates the positive solution as the unknown variables show concentration of chemical substances. Three numerical methods are used to solve glycolysis model. Two methods are well-known finite difference (FD) schemes and one is proposed FD scheme. The proposed scheme is explicit in nature. The main feature of the proposed FD scheme is to preserve the property of positivity retained by the glycolysis model. Results are compared with forward Euler explicit scheme and Crank Nicolson implicit scheme. All the attributes are verified by simulations.
