expand_by_r.m

％The Metabolic Networks Toolbox contains functions to create, ％modify, display, and simulate bioche

function delta_y=expand_by_R(RY_1,RY_2,delta_par);

%delta_y = expand_by_R(RY_1, RY_2, delta_par);
%
% expand steady-state quantity y using the response coefficients 
%
% delta_par: column vector of parameter differences

delta_y = RY_1*delta_par + 0.5*tensor_product(RY_2,delta_par)*delta_par;