function sum = Calculate_Sum_Of_Avals()

% Variables
global mRNAa;
global A;
global A2;
global Pz;
global PzA2;
global PzA4;
global mRNAz;

global l1;

% Kinetic Coefficients
global k_xlate;
global k_xcribe;
global k_dec_mrna_a;
global k_dec_mrna_z;
global k_dec_a;
global k_dec_a2;
global k_dimer_a;
global k_sngl_a;
global k_rprs_a2;
global k_rprs_a4;
global k_dissoc_a2;
global k_dissoc_a4;
global k_l1;

% DEFINITION
global REGRESSION;
global END_TIME;
global RECORD_STEP;
global EQ_NUMBERS;
global AVALS;
global PLASMID_COPIES;
global rRNA;
global RNAp;

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%% Function Definition
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

sum = 0.0;

% compute the actual reaction rate
AVALS(1) = k_xlate       * mRNAa * rRNA;
AVALS(2) = k_dec_mrna_a  * mRNAa;
AVALS(3) = k_dimer_a     * A     * A;
AVALS(4) = k_dec_a       * A;
AVALS(5) = k_sngl_a      * A2;
AVALS(6) = k_dec_a2      * A2;
AVALS(7) = k_rprs_a2     * Pz    * A2;
AVALS(8) = k_dissoc_a2   * PzA2;
AVALS(9) = k_dec_a2      * PzA2;
AVALS(10) = k_dec_a2     * PzA4;
AVALS(11) = k_rprs_a4    * PzA2  * A2;
AVALS(12) = k_dissoc_a4  * PzA4;
AVALS(13) = k_xcribe     * Pz    * RNAp;
AVALS(14) = k_dec_mrna_z * mRNAz;

AVALS(15) = k_l1 * l1;

% compute sum of all reaction rates
for i = 1 : EQ_NUMBERS,
    sum = sum + AVALS(i);
end