$\mathrm{epsilon}=e^{\frac{F\mathrm{u\_e\_m}}{\mathrm{RT}}}$ $\mathrm{v\_HCO3\_m}=\mathrm{kappa\_HCO3\_m}(\mathrm{Q\_HCO3\_o}\mathrm{epsilon}-\mathrm{Q\_HCO3\_i})$ $\mathrm{v\_H\_m}=\mathrm{kappa\_H\_m}(\mathrm{Q\_H\_o}-\mathrm{Q\_H\_i}\mathrm{epsilon})$ $\mathrm{v\_H2CO3\_o}=\mathrm{kappa\_H2CO3}(\mathrm{Q\_CO2\_o}-\mathrm{Q\_HCO3\_o}\mathrm{Q\_H\_o})\mathrm{v\_H2CO3\_i}=\mathrm{kappa\_H2CO3}(\mathrm{Q\_CO2\_i}-\mathrm{Q\_HCO3\_i}\mathrm{Q\_H\_i})\frac{d \mathrm{q\_CO2\_i}}{d t}=\mathrm{v\_CO2\_m}-\mathrm{v\_H2CO3\_i}\frac{d \mathrm{q\_HCO3\_o}}{d t}=-\mathrm{v\_HCO3\_m}+\mathrm{v\_H2CO3\_o}\frac{d \mathrm{q\_HCO3\_i}}{d t}=\mathrm{v\_HCO3\_m}+\mathrm{v\_H2CO3\_i}\frac{d \mathrm{q\_H\_i}}{d t}=\mathrm{v\_H\_m}+\mathrm{v\_H2CO3\_i}$