%
% Reaction equation
%
function [rv] = reaction_v(u,v)

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

  rv = gg*(alpha*(bb-v)-rho*u.*v./(1+u+K*u.^2));


