function [f,g] = reactionsTuring(u,v,a,b,alpha)
% compute kinetic terms for Schnackenberg reaction system  
  temp = u.*u.*v;
  f = alpha*(a-u+temp);
  g = alpha*(b-temp);