© 2021, SBMAC - Sociedade Brasileira de Matemática Aplicada e Computacional.This article concerns with the numerical investigations of the reaction–diffusion systems (RDSs) arising in the study of pattern formation in biological and chemical systems with the employment of the quartic-trigonometric B-spline functions. The computationally numerical scheme uses collocation method which is established by a relatively new B-splines for the spatial discretizations and, for time integration Crank–Nicolson technique is adapted. Therefore, solutions of the RDSs are assembled by the wholly discretized space-time scheme. A matrix stability analysis is performed for the numerical scheme after linearization process. Experimental cases include Brusselator model, Gray–Scott model, Schnakenberg model as well as a linear problem in one-dimensional domain. Numerical solutions are compared to the existing studies. Spatial pattern formation is demonstrated by present computational algorithm.