$\begin{cases}0 & \text{if $(l\ge \mathrm{l\_0})\land (\mathrm{time\_}< 0.15)$}\\ 1 & \text{otherwise}\end{cases}$ $\begin{cases}0 & \text{if $\mathrm{isotonic}=0$}\\ 1 & \text{if $(\mathrm{isotonic}=1)\land (\mathrm{F\_muscle}\ge \mathrm{F\_afterload})$}\\ 0 & \text{if $(\mathrm{isotonic}=1)\land (l\ge \mathrm{l\_0})\land (\mathrm{flag}=1)$}\end{cases}$ $\begin{cases}\mathrm{q\_1}-\frac{\mathrm{q\_2}v}{\mathrm{v\_max}} & \text{if $v\le 0$}\\ \frac{(\mathrm{q\_4}-\mathrm{q\_3})v}{\mathrm{v\_star}}+\mathrm{q\_3} & \text{if $(v\le \mathrm{v\_star})\land (0< v)$}\\ \frac{\mathrm{q\_4}}{(1+\frac{\mathrm{beta\_Q}(v-\mathrm{v\_star})}{\mathrm{v\_max}})^{\mathrm{alpha\_Q}}} & \text{otherwise}\end{cases}$ $\frac{\mathrm{v\_max}}{10}$ $\mathrm{lambda}\mathrm{p\_v}N$ $\mathrm{F\_XSE}$ $\mathrm{beta\_3}(e^{\mathrm{alpha\_3}\mathrm{l\_3}}-1)$ $\mathrm{beta\_1}(e^{\mathrm{alpha\_1}(\mathrm{l\_2}-\mathrm{l\_1})}-1)$ $\mathrm{beta\_2}(e^{\mathrm{alpha\_2}\mathrm{l\_2}}-1)$ $\mathrm{K\_chi}$ $\mathrm{k\_p\_v}\mathrm{M\_A}\mathrm{n\_1}\mathrm{L\_oz}(1-N)-\mathrm{k\_m\_v}N$ $\frac{A^{\mathrm{mu}}}{A^{\mathrm{mu}}+\mathrm{A\_half}^{\mathrm{mu}}}$ $0.6\mathrm{l\_1}+0.5$ $\frac{\mathrm{l\_1}+\mathrm{S\_0}}{0.46+\mathrm{S\_0}}$ $\mathrm{chi}\mathrm{chi\_0}\mathrm{q\_v}\mathrm{m\_0}\mathrm{G\_star}$ $\begin{cases}\mathrm{chi\_0}\mathrm{q\_v}(1-\mathrm{chi}\mathrm{m\_0}\mathrm{G\_star}) & \text{if $v\le \mathrm{v\_star}$}\\ \mathrm{chi\_0}(\mathrm{q\_4}(1-\mathrm{chi}\mathrm{m\_0}\mathrm{G\_star})+\frac{\mathrm{q\_star}(v-\mathrm{v\_star})}{\mathrm{v\_0}-\mathrm{v\_star}}) & \text{otherwise}\end{cases}$ $v$ $\mathrm{dl\_1\_dt}$ $\begin{cases}\mathrm{phi\_chi\_2} & \text{if $\mathrm{k\_S\_vis}=0$}\\ w & \text{otherwise}\end{cases}$ $\mathrm{dl\_2\_dt}$ $\begin{cases}0 & \text{if $\mathrm{isotonic\_mode}=1$}\\ -\mathrm{phi\_chi\_2} & \text{if $(\mathrm{isotonic\_mode}=0)\land (\mathrm{k\_S\_vis}=0)$}\\ -w & \text{otherwise}\end{cases}$ $\mathrm{dl\_3\_dt}$ $\mathrm{l\_2}+\mathrm{l\_3}$ $\begin{cases}\frac{\mathrm{alpha\_1}\mathrm{beta\_1}e^{\mathrm{alpha\_1}(\mathrm{l\_2}-\mathrm{l\_1})}(\mathrm{phi\_chi\_2}-v)-\mathrm{lambda}\mathrm{K\_chi}\mathrm{p\_v}+\mathrm{alp\_p}\mathrm{k\_P\_vis}v^2}{\mathrm{lambda}N\mathrm{p\_prime\_v}+\mathrm{k\_P\_vis}} & \text{if $\mathrm{k\_S\_vis}=0$}\\ \mathrm{phi\_chi} & \text{otherwise}\end{cases}$ $\begin{cases}\frac{-(\mathrm{lambda}\mathrm{K\_chi}\mathrm{p\_v}+\mathrm{alp\_p}\mathrm{k\_P\_vis}v^2+\mathrm{alpha\_2}\mathrm{beta\_2}e^{\mathrm{alpha\_2}\mathrm{l\_2}}w)}{\mathrm{lambda}N\mathrm{p\_prime\_v}+\mathrm{k\_P\_vis}} & \text{if $\mathrm{isotonic\_mode}=1$}\\ \frac{-(\mathrm{lambda}\mathrm{K\_chi}\mathrm{p\_v}+\mathrm{alp\_p}\mathrm{k\_P\_vis}v^2+(\mathrm{alpha\_2}\mathrm{beta\_2}e^{\mathrm{alpha\_2}\mathrm{l\_2}}+\mathrm{alpha\_3}\mathrm{beta\_3}e^{\mathrm{alpha\_3}\mathrm{l\_3}})w)}{\mathrm{lambda}N\mathrm{p\_prime\_v}+\mathrm{k\_P\_vis}} & \text{otherwise}\end{cases}$ $\begin{cases}\frac{\mathrm{alpha\_1}\mathrm{beta\_1}e^{\mathrm{alpha\_1}(\mathrm{l\_2}-\mathrm{l\_1})}v}{\mathrm{alpha\_1}\mathrm{beta\_1}e^{\mathrm{alpha\_1}(\mathrm{l\_2}-\mathrm{l\_1})}+\mathrm{alpha\_2}\mathrm{beta\_2}e^{\mathrm{alpha\_2}\mathrm{l\_2}}} & \text{if $\mathrm{isotonic\_mode}=1$}\\ \frac{\mathrm{alpha\_1}\mathrm{beta\_1}e^{\mathrm{alpha\_1}(\mathrm{l\_2}-\mathrm{l\_1})}v}{\mathrm{alpha\_1}\mathrm{beta\_1}e^{\mathrm{alpha\_1}(\mathrm{l\_2}-\mathrm{l\_1})}+\mathrm{alpha\_2}\mathrm{beta\_2}e^{\mathrm{alpha\_2}\mathrm{l\_2}}+\mathrm{alpha\_3}\mathrm{beta\_3}e^{\mathrm{alpha\_3}\mathrm{l\_3}}} & \text{otherwise}\end{cases}$ $\begin{cases}\mathrm{beta\_P\_lengthening}e^{\mathrm{alpha\_P\_lengthening}\mathrm{l\_1}} & \text{if $v\le 0$}\\ \mathrm{beta\_P\_shortening}e^{\mathrm{alpha\_P\_shortening}\mathrm{l\_1}} & \text{otherwise}\end{cases}$ $\begin{cases}\mathrm{alpha\_P\_lengthening} & \text{if $v\le 0$}\\ \mathrm{alpha\_P\_shortening} & \text{otherwise}\end{cases}$ $\begin{cases}\frac{\mathrm{k\_S\_vis}(\mathrm{phi\_chi}-\mathrm{alp\_s}(w-v)^2)-\mathrm{alpha\_1}\mathrm{beta\_1}e^{\mathrm{alpha\_1}(\mathrm{l\_2}-\mathrm{l\_1})}(w-v)-\mathrm{alpha\_2}\mathrm{beta\_2}e^{\mathrm{alpha\_2}\mathrm{l\_2}}w}{\mathrm{k\_S\_vis}} & \text{if $\mathrm{isotonic\_mode}=1$}\\ \mathrm{phi\_chi}-\mathrm{alp\_s}(w-v)^2-\frac{\mathrm{alpha\_1}\mathrm{beta\_1}e^{\mathrm{alpha\_1}(\mathrm{l\_2}-\mathrm{l\_1})}(w-v)+(\mathrm{alpha\_2}\mathrm{beta\_2}e^{\mathrm{alpha\_2}\mathrm{l\_2}}+\mathrm{alpha\_3}\mathrm{beta\_3}e^{\mathrm{alpha\_3}\mathrm{l\_3}})w}{\mathrm{k\_S\_vis}} & \text{otherwise}\end{cases}$ $\begin{cases}\mathrm{beta\_S\_lengthening}e^{\mathrm{alpha\_S\_lengthening}(\mathrm{l\_2}-\mathrm{l\_1})} & \text{if $w\le v$}\\ \mathrm{beta\_S\_shortening}e^{\mathrm{alpha\_S\_shortening}(\mathrm{l\_2}-\mathrm{l\_1})} & \text{otherwise}\end{cases}$ $\begin{cases}\mathrm{alpha\_S\_lengthening} & \text{if $w\le v$}\\ \mathrm{alpha\_S\_shortening} & \text{otherwise}\end{cases}$ $\frac{\mathrm{P\_star}}{\mathrm{G\_star}}$ $\begin{cases}\mathrm{case\_1} & \text{if $v\le -\mathrm{v\_max}$}\\ \mathrm{case\_2} & \text{if $(-\mathrm{v\_max}< v)\land (v\le 0)$}\\ \mathrm{case\_3} & \text{if $(0< v)\land (v\le \mathrm{v\_1})$}\\ \mathrm{case\_4} & \text{otherwise}\end{cases}$ $\begin{cases}\frac{a(1+\frac{v}{\mathrm{v\_max}})}{a-\frac{v}{\mathrm{v\_max}}} & \text{if $v\le 0$}\\ 1+\mathrm{d\_h}-\frac{\mathrm{d\_h}^{2}a}{\frac{a\mathrm{d\_h}}{\mathrm{gamma}}\left(\frac{v}{\mathrm{v\_max}}\right)^2+\frac{(a+1)v}{\mathrm{v\_max}}+a\mathrm{d\_h}} & \text{otherwise}\end{cases}$ $\begin{cases}1+\frac{0.6v}{\mathrm{v\_max}} & \text{if $(-\mathrm{v\_max}\le v)\land (v\le 0)$}\\ \frac{\mathrm{P\_star}}{\frac{(0.4a+1)v}{a\mathrm{v\_max}}+1} & \text{if $(0< v)\land (v\le \mathrm{v\_1})$}\\ \frac{\mathrm{P\_star}e^{-\mathrm{alpha\_G}\left(\frac{v-\mathrm{v\_1}}{\mathrm{v\_max}}\right)^{\mathrm{alpha\_P}}}}{\frac{(0.4a+1)v}{a\mathrm{v\_max}}+1} & \text{otherwise}\end{cases}$ $\frac{a\mathrm{d\_h}\left(\frac{\mathrm{v\_1}}{\mathrm{v\_max}}\right)^2}{3a\mathrm{d\_h}-\frac{(a+1)\mathrm{v\_1}}{\mathrm{v\_max}}}$ $\frac{a(0.4+0.4a)}{\mathrm{v\_max}((a+1)\times 0.4)^2}$ $\frac{a\times 1(1+0.4a+\frac{1.2v}{\mathrm{v\_max}}+0.6\left(\frac{v}{\mathrm{v\_max}}\right)^2)}{\mathrm{v\_max}((a-\frac{v}{\mathrm{v\_max}})(1+\frac{0.6v}{\mathrm{v\_max}}))^2}$ $\frac{0.4a+1}{a\mathrm{v\_max}}$ $\frac{1}{\mathrm{v\_max}}e^{-\mathrm{alpha\_G}(\frac{v}{\mathrm{v\_max}}-\frac{\mathrm{v\_1}}{\mathrm{v\_max}})^{\mathrm{alpha\_P}}}(\frac{0.4a+1}{a}+\mathrm{alpha\_G}\mathrm{alpha\_P}(1+\frac{(0.4a+1)v}{a\mathrm{v\_max}})(\frac{v}{\mathrm{v\_max}}-\frac{\mathrm{v\_1}}{\mathrm{v\_max}})^{(\mathrm{alpha\_P}-1)})$ $\mathrm{dA\_dt}$ $\mathrm{a\_on}(\mathrm{A\_tot}-A)\mathrm{Ca\_C}-\mathrm{a\_off}e^{-\mathrm{k\_A}A}\mathrm{pi\_N\_A}A$ $\frac{\mathrm{A\_tot}N}{\mathrm{L\_oz}A}$ $\begin{cases}1 & \text{if $\mathrm{N\_A}\ge 1$}\\ 0.02^{\mathrm{N\_A}} & \text{otherwise}\end{cases}$ $\mathrm{dB\_dt}$ $\mathrm{b\_on}(\mathrm{B\_tot}-B)\mathrm{Ca\_C}-\mathrm{b\_off}B$ $\begin{cases}4\mathrm{a\_c}\mathrm{Ca\_m}\mathrm{time\_}(1-e^{-\mathrm{a\_c}\mathrm{time\_}^2})e^{-\mathrm{a\_c}\mathrm{time\_}^2} & \text{if $\mathrm{time\_}< \mathrm{t\_d}$}\\ -\mathrm{dA\_dt}-\mathrm{dB\_dt}-\mathrm{r\_Ca}e^{-\mathrm{q\_Ca}\mathrm{Ca\_C}}\mathrm{Ca\_C} & \text{otherwise}\end{cases}$