%
% Reaction equation
%
function [ru] = reaction_u(u,v)

 global D1 D2 allt allx ally allu allv ran;
 global aa bb gg alpha rho K

 ru = gg*(aa-u-rho*u.*v./(1+u+K*u.^2));

