A class of conservative methods is developed in the more general framework of cell-centered upwind differences to approximate numerically the solution of one-dimensional non-linear conservation laws with (possibly) stiff reaction source terms. These methods are based on a non-oscillatory...