Nowadays, numerical models have great importance in every field of science, especially for solving nonlinear differential equations, partial differential equations, biochemical reactions, etc. The total time evolution of the reactant species which interacts with other species is simulated by the Runge-Kutta of order four (RK4) and by the Non-Standard finite difference (NSFD) method. A NSFD model has been constructed for the biochemical reaction problems and numerical experiments were performed for different values of discretization parameter 'h'. The results were compared with a well-known numerical scheme, i.e., RK4. The developed scheme NSFD gave better results than RK4.