We introduce two numerical schemes for Solving a system of ordinary differential equations which characterizes several kinds of linear reactions and diffusion from biochemistry, physiology, etc. The methods consist of sequential applications of the simple exact solver for a reversible reaction. We prove absolute stability and convergence of the proposed explicit methods. One is of first order and the other is of second order. Numerical results are included.