Summary: This article studies the pattern formation of reaction-diffusion Brusselator model along with Neumann boundary conditions arising in chemical processes. To accomplish this work, a new modified trigonometric cubic B-spline functions based differential quadrature algorithm is developed which is more general than [R. C. Mittal and R. Jiwari, Appl. Math. Comput. 217, No. 12, 5404–5415 (2011; Zbl 1209.65102); R. Jiwari and J. Yuan, J. Math. Chem. 52, No. 6, 1535–1551 (2014; Zbl 1296.92252)]. The reaction-diffusion model arises in enzymatic reactions, in the formation of ozone by atomic oxygen via a triple collision, and in laser and plasma physics in multiple couplings between modes. The algorithm converts the model into a system of ordinary differential equations and the obtained system is solved by Runge-Kutta method. To check the precision and performance of the proposed algorithm four numerical problems are contemplated and computed results are compared with the existing methods. The computed results pamper the theory of Brusselator model that for small values of diffusion coefficient, the steady state solution converges to equilibrium point \(( {\mu, \lambda /\mu})\) if \(1-\lambda +\mu^{2}>0\).