Nowadays, numerical models have great importance in every field of science, especially for solving the nonlinear differential equations, partial differential equations, biochemical reactions, etc. The total time evolution of the reactant concentrations in the basic enzyme-substrate reaction is simulated by the Runge-Kutta of order four (RK4) and by nonstandard finite difference (NSFD) method. A NSFD model has been constructed for the biochemical reaction problem and numerical experiments are performed for different values of discretization parameter ‘h’. The results are compared with the well-known numerical scheme, i.e. RK4. Unlike RK4 which fails for large time steps, the developed scheme gives results that converge to true steady states for any time step used.