This paper is concerned with the development of a numerical approach based on the Meshless
Local Petrov-Galerkin (MLPG) method for the approximate solutions of the two dimensional
nonlinear reaction-diffusion Brusselator systems. The method uses finite differences for
discretizing the time variable and the moving least squares (MLS) approximation for field
variables. The application of the weak formulation with the Heaviside type test functions
supported on local subdomains (around the nodes used in MLS approximation) to semidiscretized
partial differential equations yields the finite-volume local weak formulation. A
predictor-corrector scheme is used to handle the nonlinearity of the problem within each time
step. Numerical test problems are given to verify the accuracy of the proposed method. Under
particular conditions, this system exhibits Turing instability which results in a pattern forming
instability. This concept is studied and a test problem is given.