%
% Reaction equation
%
function [r] = reaction(u)


   r = 2*u;

