function dF_dV=odeP6_16(V,F)
global k1 k2 k3 k4 k5 Cto
%
%  Differential PFR design equation: dFj_dV=rj
%  First, set up the rate expressions in terms of molar flows
%
Ft=F(1)+F(2)+F(3)+F(4)+F(5)+F(6)+F(7)+F(8);
r1=k1*(Cto^1.5)*(Ft^(-1.5))*(F(1))*(F(7)^0.5);
r2=k2*(Cto^1.5)*(Ft^(-1.5))*(F(2))*(F(7)^0.5);
r3=k3*(Cto^1.5)*(Ft^(-1.5))*(F(3))*(F(7)^0.5);
r4=k4*(Cto^1.5)*(Ft^(-1.5))*(F(4))*(F(7)^0.5);
r5=k5*(Cto^1.5)*(Ft^(-1.5))*(F(5))*(F(7)^0.5);
%
%  Then write differential mole balances for all species
%  using the appropriate stoichiometry 
%
dF_dV(1) = -r1;  % this is mole balance for C11
dF_dV(2) = r1-r2;  % this is mole balance for C10
dF_dV(3) = r2-r3;  % this is mole balance for C9
dF_dV(4) = r3-r4;  % this is mole balance for C8
dF_dV(5) = r4-r5;  % this is mole balance for C7
dF_dV(6) = r5;  % this is mole balance for C6
dF_dV(7) = -r1-r2-r3-r4-r5;  % this is mole balance for H2
dF_dV(8) = r1+r2+r3+r4+r5;     % this is mole balance for M
dF_dV = dF_dV';
