We construct a nonstandard finite difference scheme for the two coupled
ODE's that model glycolysis. The primary emphasis is having the scheme satisfy a
positivity condition and also retain the limit-cycle behavior for certain values of the
parameters. We show that this is possible and give a full discussion of the scheme
along with some of its numerical properties.