The Brusselator equations, originally studied as a model for chemical dynamics, exhibit a wide range of behavior. In this paper, a powerful analytical method, called Homotopy analysis method (HAM) is used to solve the system of nonlinear differential equations. Furthermore, in this work the numerical simulation of the problem is also reported using Scilab/Matlab program. Our analytical results are compared with simulation results. A good agreement between analytical and numerical results is noted. Keywords: Chemical dynamics; Mathematical modeling; Brusselator equation; Non-linear reaction- diffusion equation; Homotopy analysis method.