In this thesis, an efficient numerical method is proposed for the numerical solution of the chemical reaction-diffusion model governed by a non-linear system of partial differential equations known as a Brusselator model. The method proposed is based on a combination of higher-order Compact Finite Difference schemes and stable time integrator known as an adaptive step-size Runge-Kutta method. The performance of adaptive step-size Runge-Kutta formula of fifth-order accurate in time and Compact Finite Difference scheme of sixth-order in space are investigated. The method is implemented to solve three test problems and reveals that the method is capable of achieving high efficiency, accuracy and reliability. The results obtained are sufficiently accurate compared to some available results in the literature.