This paper deals with the design of a numerical method for the Schnakenberg model. The patterns are nicely captured by way of the parameter values of the Schnakenberg model. The spatial integration of the equation is achieved by using a finite element method setting up the trigonometric quadratic B-spline collocation method over the subelements of the problem domain. The Crank–Nicolson technique is employed to get fully integrated Schnakenberg model. Numerical examples are given to show the good agreement with the Schnakenberg patterns.