$\frac{d \mathrm{q\_Ca\_o}}{d \mathrm{time}}=-\mathrm{v\_Re\_depol}\frac{d \mathrm{q\_Ca\_i}}{d \mathrm{time}}=\mathrm{v\_Re\_depol}\frac{d \mathrm{q\_R\_Mg}}{d \mathrm{time}}=-\mathrm{v\_Re\_bind}\frac{d \mathrm{q\_R\_Mg\_open}}{d \mathrm{time}}=\mathrm{v\_Re\_bind}-\mathrm{v\_Re\_depol}\frac{d \mathrm{q\_R\_open}}{d \mathrm{time}}=\mathrm{v\_Re\_depol}\frac{d \mathrm{q\_Gly}}{d \mathrm{time}}=-\mathrm{v\_Re\_bind}\mathrm{nGly}\frac{d \mathrm{q\_Glu}}{d \mathrm{time}}=-\mathrm{v\_Re\_bind}\mathrm{nGlu}\frac{d \mathrm{q\_Mg\_o}}{d \mathrm{time}}=\mathrm{v\_Re\_depol}$ $\mathrm{mu\_Ca\_o}=RT\ln (\mathrm{K\_Ca\_o}\mathrm{q\_Ca\_o})\mathrm{mu\_Ca\_i}=RT\ln (\mathrm{K\_Ca\_i}\mathrm{q\_Ca\_i})\mathrm{mu\_R\_Mg}=RT\ln (\mathrm{K\_R\_Mg}\mathrm{q\_R\_Mg})\mathrm{mu\_R\_Mg\_open}=RT\ln (\mathrm{K\_R\_Mg\_open}\mathrm{q\_R\_Mg\_open})\mathrm{mu\_R\_open}=RT\ln (\mathrm{K\_R\_open}\mathrm{q\_R\_open})\mathrm{mu\_Gly}=RT\ln (\mathrm{K\_Gly}\mathrm{q\_Gly})\mathrm{mu\_Glu}=RT\ln (\mathrm{K\_Glu}\mathrm{q\_Glu})\mathrm{mu\_Mg\_o}=RT\ln (\mathrm{K\_Mg\_o}\mathrm{q\_Mg\_o})\mathrm{v\_Re\_bind}=\mathrm{kappa\_Re\_bind}(e^{\frac{\mathrm{mu\_R\_Mg}+\mathrm{nGly}\mathrm{mu\_Gly}+\mathrm{nGlu}\mathrm{mu\_Glu}}{RT}}-e^{\frac{\mathrm{mu\_R\_Mg\_open}}{RT}})\mathrm{v\_Re\_depol}=\mathrm{kappa\_Re\_depol}(e^{\frac{\mathrm{mu\_Ca\_o}+\mathrm{mu\_R\_Mg\_open}+\mathrm{mu\_Glu}}{RT}}-e^{\frac{\mathrm{mu\_R\_open}+\mathrm{mu\_Mg\_o}+\mathrm{mu\_Ca\_i}}{RT}})$