Location: A computational human atrial myofibroblast electrophysiological model @ 203520b0e61a / atrialmyofibroblast.cellml

Author:
Jorge Sanchez <jorge.arciniegas@kit.edu>
Date:
2021-10-21 09:59:06+02:00
Desc:
First version of the model added
Permanent Source URI:
https://models.physiomeproject.org/workspace/785/rawfile/203520b0e61a344db31677d3f9f950d85471e8f1/atrialmyofibroblast.cellml

<?xml version='1.0' encoding='UTF-8'?>
<model xmlns="http://www.cellml.org/cellml/2.0#" xmlns:cellml="http://www.cellml.org/cellml/2.0#" name="atrialmyofibroblast">
  <component name="cell">
    <variable name="B" units="dimensionless"/>
    <variable name="Ca_c" units="dimensionless"/>
    <variable name="Ca_d" units="dimensionless" initial_value="7.1e-05"/>
    <variable name="Ca_i" units="dimensionless" initial_value="6.5e-05"/>
    <variable name="Ca_o" units="dimensionless" initial_value="1.8"/>
    <variable name="Ca_rel" units="dimensionless" initial_value="0.632613"/>
    <variable name="Ca_up" units="dimensionless" initial_value="0.649195"/>
    <variable name="Cm" units="dimensionless" initial_value="6.3"/>
    <variable name="E_Ca" units="dimensionless"/>
    <variable name="E_Ca_app" units="dimensionless" initial_value="60.0"/>
    <variable name="E_K" units="dimensionless"/>
    <variable name="E_Na" units="dimensionless"/>
    <variable name="Eca_app" units="dimensionless" initial_value="60.0"/>
    <variable name="F" units="dimensionless" initial_value="96487.0"/>
    <variable name="F1" units="dimensionless" initial_value="0.470055"/>
    <variable name="F2" units="dimensionless" initial_value="0.002814"/>
    <variable name="IK1" units="dimensionless"/>
    <variable name="IKv" units="dimensionless"/>
    <variable name="INKA" units="dimensionless"/>
    <variable name="INa" units="dimensionless"/>
    <variable name="INab" units="dimensionless"/>
    <variable name="I_stim" units="dimensionless" interface="public"/>
    <variable name="I_up_max" units="dimensionless" initial_value="2800.0"/>
    <variable name="Iion" units="dimensionless"/>
    <variable name="J_O" units="dimensionless"/>
    <variable name="J_O_C" units="dimensionless"/>
    <variable name="J_O_Calse" units="dimensionless"/>
    <variable name="J_O_TC" units="dimensionless"/>
    <variable name="J_O_TMgC" units="dimensionless"/>
    <variable name="J_O_TMgMg" units="dimensionless"/>
    <variable name="K_NaCa" units="dimensionless" initial_value="0.0374842"/>
    <variable name="K_i" units="dimensionless" initial_value="124.4059"/>
    <variable name="K_mK" units="dimensionless" initial_value="1.0"/>
    <variable name="K_mNa" units="dimensionless" initial_value="11.0"/>
    <variable name="K_o" units="dimensionless" initial_value="5.3581"/>
    <variable name="Mg_i" units="dimensionless" initial_value="2.5"/>
    <variable name="Na_h_Vhalf" units="dimensionless" initial_value="64.7"/>
    <variable name="Na_h_hill" units="dimensionless" initial_value="6.9"/>
    <variable name="Na_i" units="dimensionless" initial_value="13.3789"/>
    <variable name="Na_m_Vhalf" units="dimensionless" initial_value="29.9"/>
    <variable name="Na_m_hill" units="dimensionless"/>
    <variable name="Na_o" units="dimensionless" initial_value="130.011"/>
    <variable name="O_C" units="dimensionless" initial_value="0.026766"/>
    <variable name="O_Calse" units="dimensionless" initial_value="0.431547"/>
    <variable name="O_TC" units="dimensionless" initial_value="0.012922"/>
    <variable name="O_TMgC" units="dimensionless" initial_value="0.190369"/>
    <variable name="O_TMgMg" units="dimensionless" initial_value="0.714463"/>
    <variable name="P_Na" units="dimensionless" initial_value="2.25e-05"/>
    <variable name="R" units="dimensionless" initial_value="8314.0"/>
    <variable name="T" units="dimensionless" initial_value="306.15"/>
    <variable name="V_rev" units="dimensionless"/>
    <variable name="Vm" units="dimensionless" initial_value="-26.7096"/>
    <variable name="Vol_d" units="dimensionless" initial_value="1.17679999999999998e-04"/>
    <variable name="Vol_i" units="dimensionless" initial_value="0.00137"/>
    <variable name="Vol_rel" units="dimensionless" initial_value="4.41e-05"/>
    <variable name="Vol_up" units="dimensionless" initial_value="0.0003969"/>
    <variable name="alpha_rel" units="dimensionless" initial_value="200000.0"/>
    <variable name="d_L" units="dimensionless" initial_value="0.0"/>
    <variable name="d_L_factor" units="dimensionless"/>
    <variable name="d_L_infinity" units="dimensionless"/>
    <variable name="d_NaCa" units="dimensionless" initial_value="0.0003"/>
    <variable name="diff_Ca_d" units="dimensionless"/>
    <variable name="diff_Ca_i" units="dimensionless"/>
    <variable name="diff_Ca_rel" units="dimensionless"/>
    <variable name="diff_Ca_up" units="dimensionless"/>
    <variable name="diff_F1" units="dimensionless"/>
    <variable name="diff_F2" units="dimensionless"/>
    <variable name="diff_K_i" units="dimensionless"/>
    <variable name="diff_Na_i" units="dimensionless"/>
    <variable name="diff_O_C" units="dimensionless"/>
    <variable name="diff_O_Calse" units="dimensionless"/>
    <variable name="diff_O_TC" units="dimensionless"/>
    <variable name="diff_O_TMgC" units="dimensionless"/>
    <variable name="diff_O_TMgMg" units="dimensionless"/>
    <variable name="diff_f_L1" units="dimensionless"/>
    <variable name="diff_f_L2" units="dimensionless"/>
    <variable name="diff_h1" units="dimensionless"/>
    <variable name="diff_h2" units="dimensionless"/>
    <variable name="f_Ca" units="dimensionless"/>
    <variable name="f_L1" units="dimensionless" initial_value="0.998597"/>
    <variable name="f_L2" units="dimensionless" initial_value="0.998586"/>
    <variable name="f_L_factor" units="dimensionless"/>
    <variable name="f_L_infinity" units="dimensionless"/>
    <variable name="g_B_Ca" units="dimensionless" initial_value="0.078681"/>
    <variable name="g_Ca_L" units="dimensionless" initial_value="6.75"/>
    <variable name="g_K1" units="dimensionless"/>
    <variable name="g_Kv" units="dimensionless"/>
    <variable name="g_NKA" units="dimensionless"/>
    <variable name="g_Nab" units="dimensionless"/>
    <variable name="gamma_Na" units="dimensionless" initial_value="0.45"/>
    <variable name="h1" units="dimensionless" initial_value="0.0157"/>
    <variable name="h2" units="dimensionless" initial_value="0.0157"/>
    <variable name="h_inf" units="dimensionless"/>
    <variable name="i_B_Ca" units="dimensionless"/>
    <variable name="i_CaP" units="dimensionless"/>
    <variable name="i_CaP_max" units="dimensionless" initial_value="4.0"/>
    <variable name="i_Ca_L" units="dimensionless"/>
    <variable name="i_NaCa" units="dimensionless"/>
    <variable name="i_di" units="dimensionless"/>
    <variable name="i_rel" units="dimensionless"/>
    <variable name="i_rel_f2" units="dimensionless"/>
    <variable name="i_rel_factor" units="dimensionless"/>
    <variable name="i_tr" units="dimensionless"/>
    <variable name="i_up" units="dimensionless"/>
    <variable name="k_Ca" units="dimensionless" initial_value="0.025"/>
    <variable name="k_CaP" units="dimensionless" initial_value="0.0002"/>
    <variable name="k_cyca" units="dimensionless" initial_value="0.0003"/>
    <variable name="k_rel_d" units="dimensionless" initial_value="0.003"/>
    <variable name="k_rel_i" units="dimensionless" initial_value="0.0003"/>
    <variable name="k_srca" units="dimensionless" initial_value="0.5"/>
    <variable name="k_xcs" units="dimensionless" initial_value="0.4"/>
    <variable name="m" units="dimensionless" initial_value="0.3163"/>
    <variable name="m_inf" units="dimensionless"/>
    <variable name="q10IK1" units="dimensionless" initial_value="1.4"/>
    <variable name="q10NKA" units="dimensionless" initial_value="1.63"/>
    <variable name="q10PNa" units="dimensionless" initial_value="1.45"/>
    <variable name="q10exp" units="dimensionless"/>
    <variable name="q10tauNa" units="dimensionless" initial_value="3.0"/>
    <variable name="r" units="dimensionless" initial_value="0.2614"/>
    <variable name="r_Ca_d_factor" units="dimensionless"/>
    <variable name="r_Ca_d_term" units="dimensionless"/>
    <variable name="r_Ca_i_factor" units="dimensionless"/>
    <variable name="r_Ca_i_term" units="dimensionless"/>
    <variable name="r_act" units="dimensionless"/>
    <variable name="r_inact" units="dimensionless"/>
    <variable name="r_inf" units="dimensionless"/>
    <variable name="r_recov" units="dimensionless" initial_value="0.815"/>
    <variable name="rkv_mod" units="dimensionless" initial_value="0.0"/>
    <variable name="rr" units="dimensionless" initial_value="0.0"/>
    <variable name="rr_infinity" units="dimensionless"/>
    <variable name="s" units="dimensionless" initial_value="0.7655"/>
    <variable name="s_factor" units="dimensionless"/>
    <variable name="s_inf" units="dimensionless"/>
    <variable name="skv_mod" units="dimensionless" initial_value="0.0"/>
    <variable name="ss" units="dimensionless" initial_value="0.0"/>
    <variable name="ss_infinity" units="dimensionless"/>
    <variable name="t" units="dimensionless" interface="public"/>
    <variable name="tau_d_L" units="dimensionless"/>
    <variable name="tau_di" units="dimensionless" initial_value="0.01"/>
    <variable name="tau_f_L1" units="dimensionless"/>
    <variable name="tau_f_L2" units="dimensionless"/>
    <variable name="tau_h1" units="dimensionless"/>
    <variable name="tau_h2" units="dimensionless"/>
    <variable name="tau_m" units="dimensionless"/>
    <variable name="tau_r" units="dimensionless"/>
    <variable name="tau_rr" units="dimensionless"/>
    <variable name="tau_s" units="dimensionless"/>
    <variable name="tau_ss" units="dimensionless"/>
    <variable name="tau_tr" units="dimensionless" initial_value="0.01"/>
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply>
        <eq/>
        <ci>B</ci>
        <apply>
          <minus/>
          <cn cellml:units="dimensionless">200.0</cn>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>Ca_c</ci>
        <ci>Ca_o</ci>
      </apply>
      <apply>
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci>t</ci>
          </bvar>
          <ci>Ca_d</ci>
        </apply>
        <ci>diff_Ca_d</ci>
      </apply>
      <apply>
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci>t</ci>
          </bvar>
          <ci>Ca_i</ci>
        </apply>
        <ci>diff_Ca_i</ci>
      </apply>
      <apply>
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci>t</ci>
          </bvar>
          <ci>Ca_rel</ci>
        </apply>
        <ci>diff_Ca_rel</ci>
      </apply>
      <apply>
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci>t</ci>
          </bvar>
          <ci>Ca_up</ci>
        </apply>
        <ci>diff_Ca_up</ci>
      </apply>
      <apply>
        <eq/>
        <ci>E_Ca</ci>
        <apply>
          <times/>
          <apply>
            <divide/>
            <apply>
              <times/>
              <ci>R</ci>
              <ci>T</ci>
            </apply>
            <apply>
              <times/>
              <cn cellml:units="dimensionless">2.0</cn>
              <ci>F</ci>
            </apply>
          </apply>
          <apply>
            <ln/>
            <apply>
              <divide/>
              <ci>Ca_c</ci>
              <ci>Ca_i</ci>
            </apply>
          </apply>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>E_K</ci>
        <apply>
          <times/>
          <apply>
            <divide/>
            <apply>
              <times/>
              <ci>R</ci>
              <ci>T</ci>
            </apply>
            <ci>F</ci>
          </apply>
          <apply>
            <ln/>
            <apply>
              <divide/>
              <ci>K_o</ci>
              <ci>K_i</ci>
            </apply>
          </apply>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>E_Na</ci>
        <apply>
          <times/>
          <apply>
            <divide/>
            <apply>
              <times/>
              <ci>R</ci>
              <ci>T</ci>
            </apply>
            <ci>F</ci>
          </apply>
          <apply>
            <ln/>
            <apply>
              <divide/>
              <ci>Na_o</ci>
              <ci>Na_i</ci>
            </apply>
          </apply>
        </apply>
      </apply>
      <apply>
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci>t</ci>
          </bvar>
          <ci>F1</ci>
        </apply>
        <ci>diff_F1</ci>
      </apply>
      <apply>
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci>t</ci>
          </bvar>
          <ci>F2</ci>
        </apply>
        <ci>diff_F2</ci>
      </apply>
      <apply>
        <eq/>
        <ci>IK1</ci>
        <apply>
          <divide/>
          <apply>
            <times/>
            <apply>
              <times/>
              <apply>
                <times/>
                <apply>
                  <power/>
                  <ci>q10IK1</ci>
                  <ci>q10exp</ci>
                </apply>
                <ci>g_K1</ci>
              </apply>
              <apply>
                <power/>
                <ci>K_o</ci>
                <cn cellml:units="dimensionless">0.4457</cn>
              </apply>
            </apply>
            <apply>
              <minus/>
              <ci>Vm</ci>
              <ci>E_K</ci>
            </apply>
          </apply>
          <apply>
            <plus/>
            <cn cellml:units="dimensionless">1.0</cn>
            <apply>
              <exp/>
              <apply>
                <divide/>
                <apply>
                  <divide/>
                  <apply>
                    <times/>
                    <apply>
                      <times/>
                      <cn cellml:units="dimensionless">1.5</cn>
                      <apply>
                        <plus/>
                        <apply>
                          <minus/>
                          <ci>Vm</ci>
                          <ci>E_K</ci>
                        </apply>
                        <cn cellml:units="dimensionless">3.6</cn>
                      </apply>
                    </apply>
                    <ci>F</ci>
                  </apply>
                  <ci>R</ci>
                </apply>
                <ci>T</ci>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>IKv</ci>
        <apply>
          <times/>
          <apply>
            <times/>
            <apply>
              <times/>
              <ci>g_Kv</ci>
              <ci>r</ci>
            </apply>
            <ci>s</ci>
          </apply>
          <apply>
            <minus/>
            <ci>Vm</ci>
            <ci>E_K</ci>
          </apply>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>INKA</ci>
        <apply>
          <divide/>
          <apply>
            <times/>
            <apply>
              <times/>
              <apply>
                <times/>
                <apply>
                  <times/>
                  <apply>
                    <power/>
                    <ci>q10NKA</ci>
                    <ci>q10exp</ci>
                  </apply>
                  <ci>g_NKA</ci>
                </apply>
                <apply>
                  <divide/>
                  <ci>K_o</ci>
                  <apply>
                    <plus/>
                    <ci>K_o</ci>
                    <ci>K_mK</ci>
                  </apply>
                </apply>
              </apply>
              <apply>
                <divide/>
                <apply>
                  <power/>
                  <ci>Na_i</ci>
                  <cn cellml:units="dimensionless">1.5</cn>
                </apply>
                <apply>
                  <plus/>
                  <apply>
                    <power/>
                    <ci>Na_i</ci>
                    <cn cellml:units="dimensionless">1.5</cn>
                  </apply>
                  <apply>
                    <power/>
                    <ci>K_mNa</ci>
                    <cn cellml:units="dimensionless">1.5</cn>
                  </apply>
                </apply>
              </apply>
            </apply>
            <apply>
              <minus/>
              <ci>Vm</ci>
              <ci>V_rev</ci>
            </apply>
          </apply>
          <apply>
            <minus/>
            <ci>Vm</ci>
            <ci>B</ci>
          </apply>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>INa</ci>
        <apply>
          <divide/>
          <apply>
            <times/>
            <apply>
              <divide/>
              <apply>
                <times/>
                <apply>
                  <times/>
                  <apply>
                    <times/>
                    <apply>
                      <times/>
                      <apply>
                        <times/>
                        <apply>
                          <times/>
                          <apply>
                            <power/>
                            <ci>q10PNa</ci>
                            <ci>q10exp</ci>
                          </apply>
                          <ci>P_Na</ci>
                        </apply>
                        <apply>
                          <power/>
                          <ci>m</ci>
                          <cn cellml:units="dimensionless">3.0</cn>
                        </apply>
                      </apply>
                      <apply>
                        <plus/>
                        <apply>
                          <times/>
                          <cn cellml:units="dimensionless">0.9</cn>
                          <ci>h1</ci>
                        </apply>
                        <apply>
                          <times/>
                          <cn cellml:units="dimensionless">0.1</cn>
                          <ci>h2</ci>
                        </apply>
                      </apply>
                    </apply>
                    <ci>Na_o</ci>
                  </apply>
                  <ci>Vm</ci>
                </apply>
                <apply>
                  <power/>
                  <ci>F</ci>
                  <cn cellml:units="dimensionless">2.0</cn>
                </apply>
              </apply>
              <apply>
                <times/>
                <ci>R</ci>
                <ci>T</ci>
              </apply>
            </apply>
            <apply>
              <minus/>
              <apply>
                <exp/>
                <apply>
                  <divide/>
                  <apply>
                    <divide/>
                    <apply>
                      <times/>
                      <apply>
                        <minus/>
                        <ci>Vm</ci>
                        <ci>E_Na</ci>
                      </apply>
                      <ci>F</ci>
                    </apply>
                    <ci>R</ci>
                  </apply>
                  <ci>T</ci>
                </apply>
              </apply>
              <cn cellml:units="dimensionless">1.0</cn>
            </apply>
          </apply>
          <apply>
            <minus/>
            <apply>
              <exp/>
              <apply>
                <divide/>
                <apply>
                  <divide/>
                  <apply>
                    <times/>
                    <ci>Vm</ci>
                    <ci>F</ci>
                  </apply>
                  <ci>R</ci>
                </apply>
                <ci>T</ci>
              </apply>
            </apply>
            <cn cellml:units="dimensionless">1.0</cn>
          </apply>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>INab</ci>
        <apply>
          <times/>
          <ci>g_Nab</ci>
          <apply>
            <minus/>
            <ci>Vm</ci>
            <ci>E_Na</ci>
          </apply>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>Iion</ci>
        <apply>
          <divide/>
          <apply>
            <plus/>
            <apply>
              <plus/>
              <apply>
                <plus/>
                <apply>
                  <plus/>
                  <apply>
                    <plus/>
                    <apply>
                      <plus/>
                      <apply>
                        <plus/>
                        <apply>
                          <plus/>
                          <ci>INa</ci>
                          <ci>IKv</ci>
                        </apply>
                        <ci>IK1</ci>
                      </apply>
                      <ci>INab</ci>
                    </apply>
                    <ci>INKA</ci>
                  </apply>
                  <ci>i_CaP</ci>
                </apply>
                <ci>i_NaCa</ci>
              </apply>
              <ci>i_Ca_L</ci>
            </apply>
            <ci>i_B_Ca</ci>
          </apply>
          <ci>Cm</ci>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>J_O</ci>
        <apply>
          <plus/>
          <apply>
            <plus/>
            <apply>
              <times/>
              <cn cellml:units="dimensionless">0.08</cn>
              <ci>J_O_TC</ci>
            </apply>
            <apply>
              <times/>
              <cn cellml:units="dimensionless">0.16</cn>
              <ci>J_O_TMgC</ci>
            </apply>
          </apply>
          <apply>
            <times/>
            <cn cellml:units="dimensionless">0.045</cn>
            <ci>J_O_C</ci>
          </apply>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>J_O_C</ci>
        <apply>
          <minus/>
          <apply>
            <times/>
            <apply>
              <times/>
              <cn cellml:units="dimensionless">200000.0</cn>
              <ci>Ca_i</ci>
            </apply>
            <apply>
              <minus/>
              <cn cellml:units="dimensionless">1.0</cn>
              <ci>O_C</ci>
            </apply>
          </apply>
          <apply>
            <times/>
            <cn cellml:units="dimensionless">476.0</cn>
            <ci>O_C</ci>
          </apply>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>J_O_Calse</ci>
        <apply>
          <minus/>
          <apply>
            <times/>
            <apply>
              <times/>
              <cn cellml:units="dimensionless">480.0</cn>
              <ci>Ca_rel</ci>
            </apply>
            <apply>
              <minus/>
              <cn cellml:units="dimensionless">1.0</cn>
              <ci>O_Calse</ci>
            </apply>
          </apply>
          <apply>
            <times/>
            <cn cellml:units="dimensionless">400.0</cn>
            <ci>O_Calse</ci>
          </apply>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>J_O_TC</ci>
        <apply>
          <minus/>
          <apply>
            <times/>
            <apply>
              <times/>
              <cn cellml:units="dimensionless">78400.0</cn>
              <ci>Ca_i</ci>
            </apply>
            <apply>
              <minus/>
              <cn cellml:units="dimensionless">1.0</cn>
              <ci>O_TC</ci>
            </apply>
          </apply>
          <apply>
            <times/>
            <cn cellml:units="dimensionless">392.0</cn>
            <ci>O_TC</ci>
          </apply>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>J_O_TMgC</ci>
        <apply>
          <minus/>
          <apply>
            <times/>
            <apply>
              <times/>
              <cn cellml:units="dimensionless">200000.0</cn>
              <ci>Ca_i</ci>
            </apply>
            <apply>
              <minus/>
              <apply>
                <minus/>
                <cn cellml:units="dimensionless">1.0</cn>
                <ci>O_TMgC</ci>
              </apply>
              <ci>O_TMgMg</ci>
            </apply>
          </apply>
          <apply>
            <times/>
            <cn cellml:units="dimensionless">6.6</cn>
            <ci>O_TMgC</ci>
          </apply>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>J_O_TMgMg</ci>
        <apply>
          <minus/>
          <apply>
            <times/>
            <apply>
              <times/>
              <cn cellml:units="dimensionless">2000.0</cn>
              <ci>Mg_i</ci>
            </apply>
            <apply>
              <minus/>
              <apply>
                <minus/>
                <cn cellml:units="dimensionless">1.0</cn>
                <ci>O_TMgC</ci>
              </apply>
              <ci>O_TMgMg</ci>
            </apply>
          </apply>
          <apply>
            <times/>
            <cn cellml:units="dimensionless">666.0</cn>
            <ci>O_TMgMg</ci>
          </apply>
        </apply>
      </apply>
      <apply>
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci>t</ci>
          </bvar>
          <ci>K_i</ci>
        </apply>
        <ci>diff_K_i</ci>
      </apply>
      <apply>
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci>t</ci>
          </bvar>
          <ci>Na_i</ci>
        </apply>
        <ci>diff_Na_i</ci>
      </apply>
      <apply>
        <eq/>
        <ci>Na_m_hill</ci>
        <apply>
          <minus/>
          <cn cellml:units="dimensionless">3.1</cn>
        </apply>
      </apply>
      <apply>
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci>t</ci>
          </bvar>
          <ci>O_C</ci>
        </apply>
        <ci>diff_O_C</ci>
      </apply>
      <apply>
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci>t</ci>
          </bvar>
          <ci>O_Calse</ci>
        </apply>
        <ci>diff_O_Calse</ci>
      </apply>
      <apply>
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci>t</ci>
          </bvar>
          <ci>O_TC</ci>
        </apply>
        <ci>diff_O_TC</ci>
      </apply>
      <apply>
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci>t</ci>
          </bvar>
          <ci>O_TMgC</ci>
        </apply>
        <ci>diff_O_TMgC</ci>
      </apply>
      <apply>
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci>t</ci>
          </bvar>
          <ci>O_TMgMg</ci>
        </apply>
        <ci>diff_O_TMgMg</ci>
      </apply>
      <apply>
        <eq/>
        <ci>V_rev</ci>
        <apply>
          <minus/>
          <cn cellml:units="dimensionless">150.0</cn>
        </apply>
      </apply>
      <apply>
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci>t</ci>
          </bvar>
          <ci>Vm</ci>
        </apply>
        <apply>
          <minus/>
          <apply>
            <minus/>
            <ci>Iion</ci>
          </apply>
          <ci>I_stim</ci>
        </apply>
      </apply>
      <apply>
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci>t</ci>
          </bvar>
          <ci>d_L</ci>
        </apply>
        <apply>
          <divide/>
          <apply>
            <minus/>
            <ci>d_L_infinity</ci>
            <ci>d_L</ci>
          </apply>
          <ci>tau_d_L</ci>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>d_L_factor</ci>
        <apply>
          <divide/>
          <apply>
            <plus/>
            <ci>Vm</ci>
            <cn cellml:units="dimensionless">35.0</cn>
          </apply>
          <cn cellml:units="dimensionless">30.0</cn>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>d_L_infinity</ci>
        <apply>
          <divide/>
          <cn cellml:units="dimensionless">1.0</cn>
          <apply>
            <plus/>
            <cn cellml:units="dimensionless">1.0</cn>
            <apply>
              <exp/>
              <apply>
                <divide/>
                <apply>
                  <plus/>
                  <ci>Vm</ci>
                  <cn cellml:units="dimensionless">9.0</cn>
                </apply>
                <apply>
                  <minus/>
                  <cn cellml:units="dimensionless">5.8</cn>
                </apply>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>diff_Ca_d</ci>
        <apply>
          <divide/>
          <apply>
            <minus/>
            <apply>
              <plus/>
              <ci>i_Ca_L</ci>
              <ci>i_di</ci>
            </apply>
          </apply>
          <apply>
            <times/>
            <apply>
              <times/>
              <cn cellml:units="dimensionless">2.0</cn>
              <ci>Vol_d</ci>
            </apply>
            <ci>F</ci>
          </apply>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>diff_Ca_i</ci>
        <apply>
          <minus/>
          <apply>
            <divide/>
            <apply>
              <minus/>
              <apply>
                <minus/>
                <apply>
                  <plus/>
                  <apply>
                    <plus/>
                    <ci>i_B_Ca</ci>
                    <ci>i_CaP</ci>
                  </apply>
                  <ci>i_up</ci>
                </apply>
                <apply>
                  <plus/>
                  <apply>
                    <plus/>
                    <ci>i_di</ci>
                    <ci>i_rel</ci>
                  </apply>
                  <apply>
                    <times/>
                    <cn cellml:units="dimensionless">2.0</cn>
                    <ci>i_NaCa</ci>
                  </apply>
                </apply>
              </apply>
            </apply>
            <apply>
              <times/>
              <apply>
                <times/>
                <cn cellml:units="dimensionless">2.0</cn>
                <ci>Vol_i</ci>
              </apply>
              <ci>F</ci>
            </apply>
          </apply>
          <apply>
            <times/>
            <cn cellml:units="dimensionless">1.0</cn>
            <ci>J_O</ci>
          </apply>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>diff_Ca_rel</ci>
        <apply>
          <minus/>
          <apply>
            <divide/>
            <apply>
              <minus/>
              <ci>i_tr</ci>
              <ci>i_rel</ci>
            </apply>
            <apply>
              <times/>
              <apply>
                <times/>
                <cn cellml:units="dimensionless">2.0</cn>
                <ci>Vol_rel</ci>
              </apply>
              <ci>F</ci>
            </apply>
          </apply>
          <apply>
            <times/>
            <cn cellml:units="dimensionless">31.0</cn>
            <ci>J_O_Calse</ci>
          </apply>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>diff_Ca_up</ci>
        <apply>
          <divide/>
          <apply>
            <minus/>
            <ci>i_up</ci>
            <ci>i_tr</ci>
          </apply>
          <apply>
            <times/>
            <apply>
              <times/>
              <cn cellml:units="dimensionless">2.0</cn>
              <ci>Vol_up</ci>
            </apply>
            <ci>F</ci>
          </apply>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>diff_F1</ci>
        <apply>
          <minus/>
          <apply>
            <times/>
            <ci>r_recov</ci>
            <apply>
              <minus/>
              <apply>
                <minus/>
                <cn cellml:units="dimensionless">1.0</cn>
                <ci>F1</ci>
              </apply>
              <ci>F2</ci>
            </apply>
          </apply>
          <apply>
            <times/>
            <ci>r_act</ci>
            <ci>F1</ci>
          </apply>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>diff_F2</ci>
        <apply>
          <minus/>
          <apply>
            <times/>
            <ci>r_act</ci>
            <ci>F1</ci>
          </apply>
          <apply>
            <times/>
            <ci>r_inact</ci>
            <ci>F2</ci>
          </apply>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>diff_K_i</ci>
        <apply>
          <divide/>
          <apply>
            <minus/>
            <apply>
              <minus/>
              <apply>
                <plus/>
                <ci>IK1</ci>
                <ci>IKv</ci>
              </apply>
              <apply>
                <times/>
                <cn cellml:units="dimensionless">2.0</cn>
                <ci>INKA</ci>
              </apply>
            </apply>
          </apply>
          <apply>
            <times/>
            <ci>Vol_i</ci>
            <ci>F</ci>
          </apply>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>diff_Na_i</ci>
        <apply>
          <divide/>
          <apply>
            <minus/>
            <apply>
              <plus/>
              <apply>
                <plus/>
                <ci>INa</ci>
                <ci>INab</ci>
              </apply>
              <apply>
                <times/>
                <cn cellml:units="dimensionless">3.0</cn>
                <ci>INKA</ci>
              </apply>
            </apply>
          </apply>
          <apply>
            <times/>
            <ci>Vol_i</ci>
            <ci>F</ci>
          </apply>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>diff_O_C</ci>
        <ci>J_O_C</ci>
      </apply>
      <apply>
        <eq/>
        <ci>diff_O_Calse</ci>
        <ci>J_O_Calse</ci>
      </apply>
      <apply>
        <eq/>
        <ci>diff_O_TC</ci>
        <ci>J_O_TC</ci>
      </apply>
      <apply>
        <eq/>
        <ci>diff_O_TMgC</ci>
        <ci>J_O_TMgC</ci>
      </apply>
      <apply>
        <eq/>
        <ci>diff_O_TMgMg</ci>
        <ci>J_O_TMgMg</ci>
      </apply>
      <apply>
        <eq/>
        <ci>diff_f_L1</ci>
        <apply>
          <divide/>
          <apply>
            <minus/>
            <ci>f_L_infinity</ci>
            <ci>f_L1</ci>
          </apply>
          <ci>tau_f_L1</ci>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>diff_f_L2</ci>
        <apply>
          <divide/>
          <apply>
            <minus/>
            <ci>f_L_infinity</ci>
            <ci>f_L2</ci>
          </apply>
          <ci>tau_f_L2</ci>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>diff_h1</ci>
        <apply>
          <divide/>
          <apply>
            <minus/>
            <ci>h_inf</ci>
            <ci>h1</ci>
          </apply>
          <ci>tau_h1</ci>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>diff_h2</ci>
        <apply>
          <divide/>
          <apply>
            <minus/>
            <ci>h_inf</ci>
            <ci>h2</ci>
          </apply>
          <ci>tau_h2</ci>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>f_Ca</ci>
        <apply>
          <divide/>
          <ci>Ca_d</ci>
          <apply>
            <plus/>
            <ci>Ca_d</ci>
            <ci>k_Ca</ci>
          </apply>
        </apply>
      </apply>
      <apply>
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci>t</ci>
          </bvar>
          <ci>f_L1</ci>
        </apply>
        <ci>diff_f_L1</ci>
      </apply>
      <apply>
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci>t</ci>
          </bvar>
          <ci>f_L2</ci>
        </apply>
        <ci>diff_f_L2</ci>
      </apply>
      <apply>
        <eq/>
        <ci>f_L_factor</ci>
        <apply>
          <plus/>
          <ci>Vm</ci>
          <cn cellml:units="dimensionless">40.0</cn>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>f_L_infinity</ci>
        <apply>
          <divide/>
          <cn cellml:units="dimensionless">1.0</cn>
          <apply>
            <plus/>
            <cn cellml:units="dimensionless">1.0</cn>
            <apply>
              <exp/>
              <apply>
                <divide/>
                <apply>
                  <plus/>
                  <ci>Vm</ci>
                  <cn cellml:units="dimensionless">27.4</cn>
                </apply>
                <cn cellml:units="dimensionless">7.1</cn>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>g_K1</ci>
        <apply>
          <times/>
          <cn cellml:units="dimensionless">0.25</cn>
          <ci>Cm</ci>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>g_Kv</ci>
        <apply>
          <times/>
          <cn cellml:units="dimensionless">0.37</cn>
          <ci>Cm</ci>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>g_NKA</ci>
        <apply>
          <times/>
          <cn cellml:units="dimensionless">5.39</cn>
          <ci>Cm</ci>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>g_Nab</ci>
        <apply>
          <times/>
          <cn cellml:units="dimensionless">0.0217</cn>
          <ci>Cm</ci>
        </apply>
      </apply>
      <apply>
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci>t</ci>
          </bvar>
          <ci>h1</ci>
        </apply>
        <ci>diff_h1</ci>
      </apply>
      <apply>
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci>t</ci>
          </bvar>
          <ci>h2</ci>
        </apply>
        <ci>diff_h2</ci>
      </apply>
      <apply>
        <eq/>
        <ci>h_inf</ci>
        <apply>
          <divide/>
          <cn cellml:units="dimensionless">1.0</cn>
          <apply>
            <plus/>
            <cn cellml:units="dimensionless">1.0</cn>
            <apply>
              <exp/>
              <apply>
                <divide/>
                <apply>
                  <minus/>
                  <apply>
                    <plus/>
                    <ci>Vm</ci>
                    <ci>Na_h_Vhalf</ci>
                  </apply>
                  <apply>
                    <times/>
                    <apply>
                      <divide/>
                      <apply>
                        <plus/>
                        <cn cellml:units="dimensionless">4.7</cn>
                        <cn cellml:units="dimensionless">11.0</cn>
                      </apply>
                      <cn cellml:units="dimensionless">20.0</cn>
                    </apply>
                    <apply>
                      <minus/>
                      <ci>T</ci>
                      <cn cellml:units="dimensionless">294.15</cn>
                    </apply>
                  </apply>
                </apply>
                <ci>Na_h_hill</ci>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>i_B_Ca</ci>
        <apply>
          <times/>
          <ci>g_B_Ca</ci>
          <apply>
            <minus/>
            <ci>Vm</ci>
            <ci>E_Ca</ci>
          </apply>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>i_CaP</ci>
        <apply>
          <divide/>
          <apply>
            <times/>
            <ci>i_CaP_max</ci>
            <ci>Ca_i</ci>
          </apply>
          <apply>
            <plus/>
            <ci>Ca_i</ci>
            <ci>k_CaP</ci>
          </apply>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>i_Ca_L</ci>
        <apply>
          <times/>
          <apply>
            <times/>
            <apply>
              <times/>
              <ci>g_Ca_L</ci>
              <ci>d_L</ci>
            </apply>
            <apply>
              <plus/>
              <apply>
                <times/>
                <ci>f_Ca</ci>
                <ci>f_L1</ci>
              </apply>
              <apply>
                <times/>
                <apply>
                  <minus/>
                  <cn cellml:units="dimensionless">1.0</cn>
                  <ci>f_Ca</ci>
                </apply>
                <ci>f_L2</ci>
              </apply>
            </apply>
          </apply>
          <apply>
            <minus/>
            <ci>Vm</ci>
            <ci>E_Ca_app</ci>
          </apply>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>i_NaCa</ci>
        <apply>
          <divide/>
          <apply>
            <times/>
            <ci>K_NaCa</ci>
            <apply>
              <minus/>
              <apply>
                <times/>
                <apply>
                  <times/>
                  <apply>
                    <times/>
                    <apply>
                      <times/>
                      <ci>Na_i</ci>
                      <ci>Na_i</ci>
                    </apply>
                    <ci>Na_i</ci>
                  </apply>
                  <ci>Ca_c</ci>
                </apply>
                <apply>
                  <exp/>
                  <apply>
                    <divide/>
                    <apply>
                      <times/>
                      <apply>
                        <times/>
                        <ci>F</ci>
                        <ci>Vm</ci>
                      </apply>
                      <ci>gamma_Na</ci>
                    </apply>
                    <apply>
                      <times/>
                      <ci>R</ci>
                      <ci>T</ci>
                    </apply>
                  </apply>
                </apply>
              </apply>
              <apply>
                <times/>
                <apply>
                  <times/>
                  <apply>
                    <times/>
                    <apply>
                      <times/>
                      <ci>Na_o</ci>
                      <ci>Na_o</ci>
                    </apply>
                    <ci>Na_o</ci>
                  </apply>
                  <ci>Ca_i</ci>
                </apply>
                <apply>
                  <exp/>
                  <apply>
                    <divide/>
                    <apply>
                      <times/>
                      <apply>
                        <times/>
                        <apply>
                          <minus/>
                          <ci>gamma_Na</ci>
                          <cn cellml:units="dimensionless">1.0</cn>
                        </apply>
                        <ci>Vm</ci>
                      </apply>
                      <ci>F</ci>
                    </apply>
                    <apply>
                      <times/>
                      <ci>R</ci>
                      <ci>T</ci>
                    </apply>
                  </apply>
                </apply>
              </apply>
            </apply>
          </apply>
          <apply>
            <plus/>
            <cn cellml:units="dimensionless">1.0</cn>
            <apply>
              <times/>
              <ci>d_NaCa</ci>
              <apply>
                <plus/>
                <apply>
                  <times/>
                  <apply>
                    <times/>
                    <apply>
                      <times/>
                      <ci>Na_o</ci>
                      <ci>Na_o</ci>
                    </apply>
                    <ci>Na_o</ci>
                  </apply>
                  <ci>Ca_i</ci>
                </apply>
                <apply>
                  <times/>
                  <apply>
                    <times/>
                    <apply>
                      <times/>
                      <ci>Na_i</ci>
                      <ci>Na_i</ci>
                    </apply>
                    <ci>Na_i</ci>
                  </apply>
                  <ci>Ca_c</ci>
                </apply>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>i_di</ci>
        <apply>
          <divide/>
          <apply>
            <times/>
            <apply>
              <times/>
              <apply>
                <times/>
                <apply>
                  <minus/>
                  <ci>Ca_d</ci>
                  <ci>Ca_i</ci>
                </apply>
                <cn cellml:units="dimensionless">2.0</cn>
              </apply>
              <ci>Vol_d</ci>
            </apply>
            <ci>F</ci>
          </apply>
          <ci>tau_di</ci>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>i_rel</ci>
        <apply>
          <times/>
          <apply>
            <times/>
            <ci>alpha_rel</ci>
            <ci>i_rel_factor</ci>
          </apply>
          <apply>
            <minus/>
            <ci>Ca_rel</ci>
            <ci>Ca_i</ci>
          </apply>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>i_rel_f2</ci>
        <apply>
          <divide/>
          <ci>F2</ci>
          <apply>
            <plus/>
            <ci>F2</ci>
            <cn cellml:units="dimensionless">0.25</cn>
          </apply>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>i_rel_factor</ci>
        <apply>
          <times/>
          <ci>i_rel_f2</ci>
          <ci>i_rel_f2</ci>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>i_tr</ci>
        <apply>
          <divide/>
          <apply>
            <times/>
            <apply>
              <times/>
              <apply>
                <times/>
                <apply>
                  <minus/>
                  <ci>Ca_up</ci>
                  <ci>Ca_rel</ci>
                </apply>
                <cn cellml:units="dimensionless">2.0</cn>
              </apply>
              <ci>Vol_rel</ci>
            </apply>
            <ci>F</ci>
          </apply>
          <ci>tau_tr</ci>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>i_up</ci>
        <apply>
          <divide/>
          <apply>
            <times/>
            <ci>I_up_max</ci>
            <apply>
              <minus/>
              <apply>
                <divide/>
                <ci>Ca_i</ci>
                <ci>k_cyca</ci>
              </apply>
              <apply>
                <divide/>
                <apply>
                  <times/>
                  <apply>
                    <times/>
                    <ci>k_xcs</ci>
                    <ci>k_xcs</ci>
                  </apply>
                  <ci>Ca_up</ci>
                </apply>
                <ci>k_srca</ci>
              </apply>
            </apply>
          </apply>
          <apply>
            <plus/>
            <apply>
              <divide/>
              <apply>
                <plus/>
                <ci>Ca_i</ci>
                <ci>k_cyca</ci>
              </apply>
              <ci>k_cyca</ci>
            </apply>
            <apply>
              <divide/>
              <apply>
                <times/>
                <ci>k_xcs</ci>
                <apply>
                  <plus/>
                  <ci>Ca_up</ci>
                  <ci>k_srca</ci>
                </apply>
              </apply>
              <ci>k_srca</ci>
            </apply>
          </apply>
        </apply>
      </apply>
      <apply>
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci>t</ci>
          </bvar>
          <ci>m</ci>
        </apply>
        <apply>
          <divide/>
          <apply>
            <minus/>
            <ci>m_inf</ci>
            <ci>m</ci>
          </apply>
          <ci>tau_m</ci>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>m_inf</ci>
        <apply>
          <divide/>
          <cn cellml:units="dimensionless">1.0</cn>
          <apply>
            <plus/>
            <cn cellml:units="dimensionless">1.0</cn>
            <apply>
              <exp/>
              <apply>
                <divide/>
                <apply>
                  <minus/>
                  <apply>
                    <plus/>
                    <ci>Vm</ci>
                    <ci>Na_m_Vhalf</ci>
                  </apply>
                  <apply>
                    <times/>
                    <apply>
                      <divide/>
                      <apply>
                        <plus/>
                        <cn cellml:units="dimensionless">4.3</cn>
                        <cn cellml:units="dimensionless">5.0</cn>
                      </apply>
                      <cn cellml:units="dimensionless">20.0</cn>
                    </apply>
                    <apply>
                      <minus/>
                      <ci>T</ci>
                      <cn cellml:units="dimensionless">294.15</cn>
                    </apply>
                  </apply>
                </apply>
                <ci>Na_m_hill</ci>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>q10exp</ci>
        <apply>
          <divide/>
          <apply>
            <minus/>
            <ci>T</ci>
            <cn cellml:units="dimensionless">306.15</cn>
          </apply>
          <cn cellml:units="dimensionless">10.0</cn>
        </apply>
      </apply>
      <apply>
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci>t</ci>
          </bvar>
          <ci>r</ci>
        </apply>
        <apply>
          <divide/>
          <apply>
            <minus/>
            <ci>r_inf</ci>
            <ci>r</ci>
          </apply>
          <ci>tau_r</ci>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>r_Ca_d_factor</ci>
        <apply>
          <times/>
          <apply>
            <times/>
            <apply>
              <times/>
              <ci>r_Ca_d_term</ci>
              <ci>r_Ca_d_term</ci>
            </apply>
            <ci>r_Ca_d_term</ci>
          </apply>
          <ci>r_Ca_d_term</ci>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>r_Ca_d_term</ci>
        <apply>
          <divide/>
          <ci>Ca_d</ci>
          <apply>
            <plus/>
            <ci>Ca_d</ci>
            <ci>k_rel_d</ci>
          </apply>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>r_Ca_i_factor</ci>
        <apply>
          <times/>
          <apply>
            <times/>
            <apply>
              <times/>
              <ci>r_Ca_i_term</ci>
              <ci>r_Ca_i_term</ci>
            </apply>
            <ci>r_Ca_i_term</ci>
          </apply>
          <ci>r_Ca_i_term</ci>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>r_Ca_i_term</ci>
        <apply>
          <divide/>
          <ci>Ca_i</ci>
          <apply>
            <plus/>
            <ci>Ca_i</ci>
            <ci>k_rel_i</ci>
          </apply>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>r_act</ci>
        <apply>
          <times/>
          <cn cellml:units="dimensionless">203.8</cn>
          <apply>
            <plus/>
            <ci>r_Ca_i_factor</ci>
            <ci>r_Ca_d_factor</ci>
          </apply>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>r_inact</ci>
        <apply>
          <plus/>
          <cn cellml:units="dimensionless">33.96</cn>
          <apply>
            <times/>
            <cn cellml:units="dimensionless">339.6</cn>
            <ci>r_Ca_i_factor</ci>
          </apply>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>r_inf</ci>
        <apply>
          <divide/>
          <cn cellml:units="dimensionless">1.0</cn>
          <apply>
            <plus/>
            <cn cellml:units="dimensionless">1.0</cn>
            <apply>
              <exp/>
              <apply>
                <divide/>
                <apply>
                  <minus/>
                  <apply>
                    <plus/>
                    <ci>Vm</ci>
                    <ci>rkv_mod</ci>
                  </apply>
                </apply>
                <cn cellml:units="dimensionless">11.0</cn>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      <apply>
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci>t</ci>
          </bvar>
          <ci>rr</ci>
        </apply>
        <apply>
          <divide/>
          <apply>
            <minus/>
            <ci>rr_infinity</ci>
            <ci>rr</ci>
          </apply>
          <ci>tau_rr</ci>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>rr_infinity</ci>
        <apply>
          <divide/>
          <cn cellml:units="dimensionless">1.0</cn>
          <apply>
            <plus/>
            <cn cellml:units="dimensionless">1.0</cn>
            <apply>
              <exp/>
              <apply>
                <divide/>
                <apply>
                  <minus/>
                  <ci>Vm</ci>
                  <cn cellml:units="dimensionless">1.0</cn>
                </apply>
                <apply>
                  <minus/>
                  <cn cellml:units="dimensionless">11.0</cn>
                </apply>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      <apply>
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci>t</ci>
          </bvar>
          <ci>s</ci>
        </apply>
        <apply>
          <divide/>
          <apply>
            <minus/>
            <ci>s_inf</ci>
            <ci>s</ci>
          </apply>
          <ci>tau_s</ci>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>s_factor</ci>
        <apply>
          <divide/>
          <apply>
            <plus/>
            <ci>Vm</ci>
            <cn cellml:units="dimensionless">52.45</cn>
          </apply>
          <cn cellml:units="dimensionless">15.8827</cn>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>s_inf</ci>
        <apply>
          <divide/>
          <cn cellml:units="dimensionless">1.0</cn>
          <apply>
            <plus/>
            <cn cellml:units="dimensionless">1.0</cn>
            <apply>
              <exp/>
              <apply>
                <divide/>
                <apply>
                  <plus/>
                  <ci>Vm</ci>
                  <ci>skv_mod</ci>
                </apply>
                <cn cellml:units="dimensionless">7.0</cn>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      <apply>
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci>t</ci>
          </bvar>
          <ci>ss</ci>
        </apply>
        <apply>
          <divide/>
          <apply>
            <minus/>
            <ci>ss_infinity</ci>
            <ci>ss</ci>
          </apply>
          <ci>tau_ss</ci>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>ss_infinity</ci>
        <apply>
          <divide/>
          <cn cellml:units="dimensionless">1.0</cn>
          <apply>
            <plus/>
            <cn cellml:units="dimensionless">1.0</cn>
            <apply>
              <exp/>
              <apply>
                <divide/>
                <apply>
                  <plus/>
                  <ci>Vm</ci>
                  <cn cellml:units="dimensionless">40.5</cn>
                </apply>
                <cn cellml:units="dimensionless">11.5</cn>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>tau_d_L</ci>
        <apply>
          <plus/>
          <apply>
            <times/>
            <cn cellml:units="dimensionless">0.0027</cn>
            <apply>
              <exp/>
              <apply>
                <times/>
                <apply>
                  <minus/>
                  <ci>d_L_factor</ci>
                </apply>
                <ci>d_L_factor</ci>
              </apply>
            </apply>
          </apply>
          <cn cellml:units="dimensionless">0.002</cn>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>tau_f_L1</ci>
        <apply>
          <plus/>
          <apply>
            <times/>
            <cn cellml:units="dimensionless">0.161</cn>
            <apply>
              <exp/>
              <apply>
                <divide/>
                <apply>
                  <divide/>
                  <apply>
                    <times/>
                    <apply>
                      <minus/>
                      <ci>f_L_factor</ci>
                    </apply>
                    <ci>f_L_factor</ci>
                  </apply>
                  <cn cellml:units="dimensionless">14.4</cn>
                </apply>
                <cn cellml:units="dimensionless">14.4</cn>
              </apply>
            </apply>
          </apply>
          <cn cellml:units="dimensionless">0.01</cn>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>tau_f_L2</ci>
        <apply>
          <plus/>
          <apply>
            <times/>
            <cn cellml:units="dimensionless">1.3323</cn>
            <apply>
              <exp/>
              <apply>
                <divide/>
                <apply>
                  <divide/>
                  <apply>
                    <times/>
                    <apply>
                      <minus/>
                      <ci>f_L_factor</ci>
                    </apply>
                    <ci>f_L_factor</ci>
                  </apply>
                  <cn cellml:units="dimensionless">14.2</cn>
                </apply>
                <cn cellml:units="dimensionless">14.2</cn>
              </apply>
            </apply>
          </apply>
          <cn cellml:units="dimensionless">0.0626</cn>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>tau_h1</ci>
        <apply>
          <times/>
          <apply>
            <divide/>
            <cn cellml:units="dimensionless">1.0</cn>
            <apply>
              <power/>
              <ci>q10tauNa</ci>
              <ci>q10exp</ci>
            </apply>
          </apply>
          <apply>
            <plus/>
            <apply>
              <divide/>
              <cn cellml:units="dimensionless">0.03</cn>
              <apply>
                <plus/>
                <cn cellml:units="dimensionless">1.0</cn>
                <apply>
                  <exp/>
                  <apply>
                    <divide/>
                    <apply>
                      <plus/>
                      <ci>Vm</ci>
                      <cn cellml:units="dimensionless">35.1</cn>
                    </apply>
                    <cn cellml:units="dimensionless">3.2</cn>
                  </apply>
                </apply>
              </apply>
            </apply>
            <cn cellml:units="dimensionless">0.0003</cn>
          </apply>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>tau_h2</ci>
        <apply>
          <times/>
          <apply>
            <divide/>
            <cn cellml:units="dimensionless">1.0</cn>
            <apply>
              <power/>
              <ci>q10tauNa</ci>
              <ci>q10exp</ci>
            </apply>
          </apply>
          <apply>
            <plus/>
            <apply>
              <divide/>
              <cn cellml:units="dimensionless">0.12</cn>
              <apply>
                <plus/>
                <cn cellml:units="dimensionless">1.0</cn>
                <apply>
                  <exp/>
                  <apply>
                    <divide/>
                    <apply>
                      <plus/>
                      <ci>Vm</ci>
                      <cn cellml:units="dimensionless">35.1</cn>
                    </apply>
                    <cn cellml:units="dimensionless">3.2</cn>
                  </apply>
                </apply>
              </apply>
            </apply>
            <cn cellml:units="dimensionless">0.003</cn>
          </apply>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>tau_m</ci>
        <apply>
          <times/>
          <apply>
            <divide/>
            <cn cellml:units="dimensionless">1.0</cn>
            <apply>
              <power/>
              <ci>q10tauNa</ci>
              <ci>q10exp</ci>
            </apply>
          </apply>
          <apply>
            <plus/>
            <apply>
              <times/>
              <cn cellml:units="dimensionless">4.2e-05</cn>
              <apply>
                <exp/>
                <apply>
                  <times/>
                  <apply>
                    <minus/>
                    <cn cellml:units="dimensionless">1.0</cn>
                  </apply>
                  <apply>
                    <power/>
                    <apply>
                      <divide/>
                      <apply>
                        <plus/>
                        <ci>Vm</ci>
                        <cn cellml:units="dimensionless">25.57</cn>
                      </apply>
                      <cn cellml:units="dimensionless">28.8</cn>
                    </apply>
                    <cn cellml:units="dimensionless">2.0</cn>
                  </apply>
                </apply>
              </apply>
            </apply>
            <cn cellml:units="dimensionless">2.4e-05</cn>
          </apply>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>tau_r</ci>
        <apply>
          <plus/>
          <cn cellml:units="dimensionless">0.0203</cn>
          <apply>
            <times/>
            <cn cellml:units="dimensionless">0.138</cn>
            <apply>
              <exp/>
              <apply>
                <times/>
                <apply>
                  <minus/>
                  <apply>
                    <divide/>
                    <apply>
                      <plus/>
                      <ci>Vm</ci>
                      <cn cellml:units="dimensionless">20.0</cn>
                    </apply>
                    <cn cellml:units="dimensionless">25.9</cn>
                  </apply>
                </apply>
                <apply>
                  <divide/>
                  <apply>
                    <plus/>
                    <ci>Vm</ci>
                    <cn cellml:units="dimensionless">20.0</cn>
                  </apply>
                  <cn cellml:units="dimensionless">25.9</cn>
                </apply>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>tau_rr</ci>
        <apply>
          <plus/>
          <apply>
            <times/>
            <cn cellml:units="dimensionless">0.0035</cn>
            <apply>
              <exp/>
              <apply>
                <divide/>
                <apply>
                  <divide/>
                  <apply>
                    <times/>
                    <apply>
                      <minus/>
                      <ci>Vm</ci>
                    </apply>
                    <ci>Vm</ci>
                  </apply>
                  <cn cellml:units="dimensionless">30.0</cn>
                </apply>
                <cn cellml:units="dimensionless">30.0</cn>
              </apply>
            </apply>
          </apply>
          <cn cellml:units="dimensionless">0.0015</cn>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>tau_s</ci>
        <apply>
          <plus/>
          <cn cellml:units="dimensionless">1.574</cn>
          <apply>
            <times/>
            <cn cellml:units="dimensionless">5.268</cn>
            <apply>
              <exp/>
              <apply>
                <times/>
                <apply>
                  <minus/>
                  <apply>
                    <divide/>
                    <apply>
                      <plus/>
                      <ci>Vm</ci>
                      <cn cellml:units="dimensionless">23.0</cn>
                    </apply>
                    <cn cellml:units="dimensionless">22.7</cn>
                  </apply>
                </apply>
                <apply>
                  <divide/>
                  <apply>
                    <plus/>
                    <ci>Vm</ci>
                    <cn cellml:units="dimensionless">23.0</cn>
                  </apply>
                  <cn cellml:units="dimensionless">22.7</cn>
                </apply>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>tau_ss</ci>
        <apply>
          <plus/>
          <apply>
            <times/>
            <cn cellml:units="dimensionless">0.025635</cn>
            <apply>
              <exp/>
              <apply>
                <times/>
                <apply>
                  <minus/>
                  <ci>s_factor</ci>
                </apply>
                <ci>s_factor</ci>
              </apply>
            </apply>
          </apply>
          <cn cellml:units="dimensionless">0.01414</cn>
        </apply>
      </apply>
    </math>
  </component>
  <component name="stimulus">
    <variable name="I_stim" units="dimensionless" interface="public"/>
    <variable name="pace" units="dimensionless" initial_value="0.0"/>
    <variable name="stimulus_amplitude" units="dimensionless"/>
    <variable name="t" units="dimensionless" interface="public"/>
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply>
        <eq/>
        <ci>I_stim</ci>
        <apply>
          <times/>
          <ci>pace</ci>
          <ci>stimulus_amplitude</ci>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>stimulus_amplitude</ci>
        <apply>
          <minus/>
          <cn cellml:units="dimensionless">80.0</cn>
        </apply>
      </apply>
    </math>
  </component>
  <connection component_1="cell" component_2="stimulus">
    <map_variables variable_1="I_stim" variable_2="I_stim"/>
    <map_variables variable_1="t" variable_2="t"/>
  </connection>
</model>