function [out] = ChemReac23(t, y)%%  Evaluate the right hand side function for the%  chemical reaction system%%    y_1'(x) = -k1*y_1(x) %    y_2'(x) =  k1*y_1(x) - k2*y_2(x)%    y_3'(x) =              k2*y_2(x)%%%global k1 k2out  = [ -k1   0   0;          k1  -k2  0;          0    k2  0 ]* y;