In this paper, we study numerical approximations of a semilinear parabolic problem in one-dimension, of which the nonlinearity appears both in source term and in Neumann boundary condition. By a semidiscretization using finite difference method, we obtain a system of ordinary differential equations which is an approximation of the original problem. We obtain some conditions under which the positive solution of our system quenches in a finite time and estimate its semidiscrete quenching time. Convergence of the numerical quenching time to the theoretical one is established. Next, we show that the quenching rate of the numerical scheme is different from the continuous one. Finally, we give some numerical results to illustrate our analysis.