The reactions are:
1 | A->B | -r1A = k1*cA; r1A = -k1*cA |
2 | B->A | r2A = k2*cB |
3 | B->C | -r3B = k3*cB |
4 | 2A->D | -r4A = k4*cA^2 |
FA - FA0 = (S4i=1 riA)*V = (r1A + r2A + r3A + r4A)*V
FA = v0*cA; FB = v0*cB
cA - cA0 = (r1A + r2A + r3A + r4A)*V/v0 = (r1A + r2A + r3A + r4A)*t
where t = V/v0
cA - cA0 = (-k1cA+k2cB-k4cA^2)*t cB - cB0 = t*S4i=1 riB = (k1cA-k2cB-k3cB)*t
Return