Location: Noble, Varghese, Kohl, Noble, 1998 @ 6e56ebe38a73 / noble_varghese_kohl_noble_1998_c.cellml

Author:
Hanne <Hanne@hanne-nielsens-macbook.local>
Date:
2010-04-12 15:58:30+12:00
Desc:
Added clickable svg for session b
Permanent Source URI:
https://staging.physiomeproject.org/workspace/noble_varghese_kohl_noble_1998/rawfile/6e56ebe38a7367b6271bc30f854555d22459912b/noble_varghese_kohl_noble_1998_c.cellml

<?xml version='1.0' encoding='utf-8'?>
<!--
This CellML file was generated on 02/10/2007 at 16:53:35 using:

COR (0.9.31.793)
Copyright 2002-2007 Dr Alan Garny
http://COR.physiol.ox.ac.uk/ - COR@physiol.ox.ac.uk

CellML 1.0 was used to generate this cellular model
http://www.CellML.org/
--><model xmlns="http://www.cellml.org/cellml/1.0#" xmlns:cmeta="http://www.cellml.org/metadata/1.0#" xmlns:dc="http://purl.org/dc/elements/1.1/" xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#" xmlns:bqs="http://www.cellml.org/bqs/1.0#" xmlns:cellml="http://www.cellml.org/cellml/1.0#" xmlns:dcterms="http://purl.org/dc/terms/" xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" xmlns:ns7="http://www.cellml.org/metadata/simulation/1.0#" cmeta:id="noble_varghese_kohl_noble_1998_stretch" name="noble_varghese_kohl_noble_1998_stretch">

<documentation xmlns="http://cellml.org/tmp-documentation">
<article>
  <articleinfo>
  <title>Improved guinea-pig ventricular cell model incorporating a diadic space, IKr and IKs, and length- and tension-dependent processes8</title>
  <author>
    <firstname>Catherine</firstname>
          <surname>Lloyd</surname>
    <affiliation>
      <shortaffil>Auckland Bioengineering Institute, The University of Auckland</shortaffil>
    </affiliation>
  </author>
</articleinfo>
  <section id="sec_status">
    <title>Model Status</title>
    <para>
This version has been curated and unit-checked by Penny Noble and is known to run in COR and PCEnv. This variant is the STRETCH model, which includes acetylcholine dependent modulation and the electromechanical processes simulated by the extended basic model, and also includes simulation of stretch and fibroblasts. This model is parameterised the same as the BASIC model, with all additional variables set to zero.
    </para>
  </section>
  <sect1 id="sec_structure">
<title>Model Structure</title>

<para>ABSTRACT: The guinea-pig ventricular cell model, originally developed by Noble et al in 1991, has been greatly extended to include accumulation and depletion of calcium in a diadic space between the sarcolemma and the sarcoplasmic reticulum where, according to contempory understanding, the majority of calcium-induced calcium release is triggered. The calcium in this space is also assumed to play the major role in calcium-induced inactivation of the calcium current. Delayed potassium current equations have been developed to include the rapid (IKr) and slow (IKs) components of the delayed rectifier current based on the data of of Heath and Terrar, along with data from Sanguinetti and Jurkiewicz. Length- and tension-dependent changes in mechanical and electrophysiological processes have been incorporated as described recently by Kohl et al. Drug receptor interactions have started to be developed, using the sodium channel as the first target. The new model has been tested against experimental data on action potential clamp, and on force-interval and duration-interval relations; it has been found to reliably reproduce experimental observations.
</para>

<para>
The original paper reference is cited below:
</para>

<para>
Improved guinea-pig ventricular cell model incorporating a diadic space, IKr and IKs, and length- and tension-dependent processes, Denis Noble, Anthony Varghese, Peter Kohl and Penelope Noble, 1998, <emphasis>Can J Cardiol</emphasis>, 14, 123-134.  <ulink url="http://www.ncbi.nlm.nih.gov/pubmed/9487284">PubMed ID: 9487284</ulink>
</para>

<informalfigure float="0" id="fig_cell_diagram">
<mediaobject>
  <imageobject>
    <objectinfo>
      <title>cell diagram of Noble'98 model showing ionic currents, pumps and exchangers within the sarcolemma and the sarcoplasmic reticulum</title>
    </objectinfo>
    <imagedata fileref="noble_1998c.png"/>
  </imageobject>
</mediaobject>
<caption>A schematic diagram describing the current flows across the cell membrane that are captured in the Noble'98 model.</caption>
</informalfigure>

</sect1>
</article>
</documentation>
   <units name="micrometre">
      <unit units="metre" prefix="micro"/>
   </units>
   <units name="per_micrometre">
      <unit units="micrometre" exponent="-1"/>
   </units>
   <units name="per_second">
      <unit units="second" exponent="-1"/>
   </units>
   <units name="millivolt">
      <unit units="volt" prefix="milli"/>
   </units>
   <units name="per_millivolt">
      <unit units="volt" prefix="milli" exponent="-1"/>
   </units>
   <units name="per_millivolt_second">
      <unit units="millivolt" exponent="-1"/>
      <unit units="second" exponent="-1"/>
   </units>
   <units name="microS">
      <unit units="siemens" prefix="micro"/>
   </units>
   <units name="microF">
      <unit units="farad" prefix="micro"/>
   </units>
   <units name="nanoA">
      <unit units="ampere" prefix="nano"/>
   </units>
   <units name="mA_nA">
      <unit units="ampere" prefix="milli"/>
      <unit units="nanoA" exponent="-1"/>
   </units>
   <units name="nanoA_per_millimolar">
      <unit units="ampere" prefix="nano"/>
      <unit units="millimolar" exponent="-1"/>
   </units>
   <units name="millimolar">
      <unit units="mole" prefix="milli"/>
      <unit units="litre" exponent="-1"/>
   </units>
   <units name="millimolar4">
      <unit units="millimolar" exponent="4"/>
   </units>
   <units name="millimolar_per_second">
      <unit units="millimolar"/>
      <unit units="second" exponent="-1"/>
   </units>
   <units name="per_millimolar_second">
      <unit units="millimolar" exponent="-1"/>
      <unit units="second" exponent="-1"/>
   </units>
   <units name="joule_per_kilomole_kelvin">
      <unit units="joule"/>
      <unit units="mole" prefix="kilo" exponent="-1"/>
      <unit units="kelvin" exponent="-1"/>
   </units>
   <units name="coulomb_per_mole">
      <unit units="coulomb"/>
      <unit units="mole" exponent="-1"/>
   </units>
   <units name="micrometre3">
      <unit units="metre" prefix="micro" exponent="3"/>
   </units>
   <units name="litre_micrometre3">
      <unit units="litre"/>
      <unit units="micrometre3" exponent="-1"/>
   </units>
   <component name="environment">
      <variable units="second" public_interface="out" name="time"/>
   </component>
   <component name="membrane">
      <variable units="millivolt" public_interface="out" name="V" initial_value="-92.849333"/>
      <variable units="joule_per_kilomole_kelvin" public_interface="out" name="R" initial_value="8314.472"/>
      <variable units="kelvin" public_interface="out" name="T" initial_value="310"/>
      <variable units="coulomb_per_mole" public_interface="out" name="F" initial_value="96485.3415"/>
      <variable units="microF" name="Cm" initial_value="9.5e-5"/>
      <variable units="second" public_interface="in" name="time"/>
      <variable units="nanoA" public_interface="in" name="i_K1"/>
      <variable units="nanoA" public_interface="in" name="i_to"/>
      <variable units="nanoA" public_interface="in" name="i_Kr"/>
      <variable units="nanoA" public_interface="in" name="i_Ks"/>
      <variable units="nanoA" public_interface="in" name="i_K_ATP"/>
      <variable units="nanoA" public_interface="in" name="i_K_ACh"/>
      <variable units="nanoA" public_interface="in" name="i_Ca_L_K_cyt"/>
      <variable units="nanoA" public_interface="in" name="i_Ca_L_K_ds"/>
      <variable units="nanoA" public_interface="in" name="i_NaK"/>
      <variable units="nanoA" public_interface="in" name="i_Na"/>
      <variable units="nanoA" public_interface="in" name="i_b_Na"/>
      <variable units="nanoA" public_interface="in" name="i_p_Na"/>
      <variable units="nanoA" public_interface="in" name="i_Ca_L_Na_cyt"/>
      <variable units="nanoA" public_interface="in" name="i_Ca_L_Na_ds"/>
      <variable units="nanoA" public_interface="in" name="i_NaCa_cyt"/>
      <variable units="nanoA" public_interface="in" name="i_NaCa_ds"/>
      <variable units="nanoA" public_interface="in" name="i_Ca_L_Ca_cyt"/>
      <variable units="nanoA" public_interface="in" name="i_Ca_L_Ca_ds"/>
      <variable units="nanoA" public_interface="in" name="i_b_Ca"/>
      <variable units="nanoA" public_interface="in" name="i_stretch"/>
      <variable units="nanoA" name="i_Stim"/>
      <variable units="second" name="stim_start" initial_value="0.1"/>
      <variable units="second" name="stim_end" initial_value="10000"/>
      <variable units="second" name="stim_period" initial_value="1"/>
      <variable units="second" name="stim_duration" initial_value="0.003"/>
      <variable units="nanoA" name="stim_amplitude" initial_value="-3"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_Stim</ci>
            <piecewise>
               <piece>
                  <ci>stim_amplitude</ci>
                  <apply>
                     <and/>
                     <apply>
                        <geq/>
                        <ci>time</ci>
                        <ci>stim_start</ci>
                     </apply>
                     <apply>
                        <leq/>
                        <ci>time</ci>
                        <ci>stim_end</ci>
                     </apply>
                     <apply>
                        <leq/>
                        <apply>
                           <minus/>
                           <apply>
                              <minus/>
                              <ci>time</ci>
                              <ci>stim_start</ci>
                           </apply>
                           <apply>
                              <times/>
                              <apply>
                                 <floor/>
                                 <apply>
                                    <divide/>
                                    <apply>
                                       <minus/>
                                       <ci>time</ci>
                                       <ci>stim_start</ci>
                                    </apply>
                                    <ci>stim_period</ci>
                                 </apply>
                              </apply>
                              <ci>stim_period</ci>
                           </apply>
                        </apply>
                        <ci>stim_duration</ci>
                     </apply>
                  </apply>
               </piece>
               <otherwise>
                  <cn cellml:units="nanoA">0</cn>
               </otherwise>
            </piecewise>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>V</ci>
            </apply>
            <apply>
               <times/>
               <apply>
                  <divide/>
                  <apply>
                     <minus/>
                     <cn cellml:units="dimensionless">1</cn>
                  </apply>
                  <ci>Cm</ci>
               </apply>
               <apply>
                  <plus/>
                  <ci>i_Stim</ci>
                  <ci>i_K1</ci>
                  <ci>i_to</ci>
                  <ci>i_Kr</ci>
                  <ci>i_Ks</ci>
                  <ci>i_K_ATP</ci>
                  <ci>i_K_ACh</ci>
                  <ci>i_NaK</ci>
                  <ci>i_Na</ci>
                  <ci>i_b_Na</ci>
                  <ci>i_p_Na</ci>
                  <ci>i_Ca_L_Na_cyt</ci>
                  <ci>i_Ca_L_Na_ds</ci>
                  <ci>i_NaCa_cyt</ci>
                  <ci>i_NaCa_ds</ci>
                  <ci>i_Ca_L_Ca_cyt</ci>
                  <ci>i_Ca_L_Ca_ds</ci>
                  <ci>i_Ca_L_K_cyt</ci>
                  <ci>i_Ca_L_K_ds</ci>
                  <ci>i_b_Ca</ci>
                  <ci>i_stretch</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="reversal_potentials">
      <variable units="millivolt" public_interface="out" name="E_Na"/>
      <variable units="millivolt" public_interface="out" name="E_K"/>
      <variable units="millivolt" public_interface="out" name="E_Ks"/>
      <variable units="millivolt" public_interface="out" name="E_Ca"/>
      <variable units="millivolt" public_interface="out" name="E_mh"/>
      <variable units="dimensionless" name="P_kna" initial_value="0.03"/>
      <variable units="millimolar" public_interface="in" name="K_o"/>
      <variable units="millimolar" public_interface="in" name="Na_o"/>
      <variable units="millimolar" public_interface="in" name="K_i"/>
      <variable units="millimolar" public_interface="in" name="Na_i"/>
      <variable units="millimolar" public_interface="in" name="Ca_o"/>
      <variable units="millimolar" public_interface="in" name="Ca_i"/>
      <variable units="joule_per_kilomole_kelvin" public_interface="in" name="R"/>
      <variable units="coulomb_per_mole" public_interface="in" name="F"/>
      <variable units="kelvin" public_interface="in" name="T"/>
      <variable units="second" public_interface="in" name="time"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <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/>
            <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_Ks</ci>
            <apply>
               <times/>
               <apply>
                  <divide/>
                  <apply>
                     <times/>
                     <ci>R</ci>
                     <ci>T</ci>
                  </apply>
                  <ci>F</ci>
               </apply>
               <apply>
                  <ln/>
                  <apply>
                     <divide/>
                     <apply>
                        <plus/>
                        <ci>K_o</ci>
                        <apply>
                           <times/>
                           <ci>P_kna</ci>
                           <ci>Na_o</ci>
                        </apply>
                     </apply>
                     <apply>
                        <plus/>
                        <ci>K_i</ci>
                        <apply>
                           <times/>
                           <ci>P_kna</ci>
                           <ci>Na_i</ci>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>E_Ca</ci>
            <apply>
               <times/>
               <apply>
                  <divide/>
                  <apply>
                     <times/>
                     <cn cellml:units="dimensionless">0.5</cn>
                     <ci>R</ci>
                     <ci>T</ci>
                  </apply>
                  <ci>F</ci>
               </apply>
               <apply>
                  <ln/>
                  <apply>
                     <divide/>
                     <ci>Ca_o</ci>
                     <ci>Ca_i</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>E_mh</ci>
            <apply>
               <times/>
               <apply>
                  <divide/>
                  <apply>
                     <times/>
                     <ci>R</ci>
                     <ci>T</ci>
                  </apply>
                  <ci>F</ci>
               </apply>
               <apply>
                  <ln/>
                  <apply>
                     <divide/>
                     <apply>
                        <plus/>
                        <ci>Na_o</ci>
                        <apply>
                           <times/>
                           <cn cellml:units="dimensionless">0.12</cn>
                           <ci>K_o</ci>
                        </apply>
                     </apply>
                     <apply>
                        <plus/>
                        <ci>Na_i</ci>
                        <apply>
                           <times/>
                           <cn cellml:units="dimensionless">0.12</cn>
                           <ci>K_i</ci>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="time_independent_potassium_current">
      <variable units="nanoA" public_interface="out" name="i_K1"/>
      <variable units="millimolar" public_interface="out" name="K_mk1" initial_value="10"/>
      <variable units="microS" name="g_K1" initial_value="0.5"/>
      <variable units="millimolar" public_interface="in" name="K_o"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="millivolt" public_interface="in" name="E_K"/>
      <variable units="joule_per_kilomole_kelvin" public_interface="in" name="R"/>
      <variable units="coulomb_per_mole" public_interface="in" name="F"/>
      <variable units="kelvin" public_interface="in" name="T"/>
      <variable units="second" public_interface="in" name="time"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_K1</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <ci>g_K1</ci>
                        <ci>K_o</ci>
                     </apply>
                     <apply>
                        <plus/>
                        <ci>K_o</ci>
                        <ci>K_mk1</ci>
                     </apply>
                  </apply>
                  <apply>
                     <minus/>
                     <ci>V</ci>
                     <ci>E_K</ci>
                  </apply>
               </apply>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <times/>
                           <apply>
                              <minus/>
                              <apply>
                                 <minus/>
                                 <ci>V</ci>
                                 <ci>E_K</ci>
                              </apply>
                              <cn cellml:units="millivolt">10</cn>
                           </apply>
                           <ci>F</ci>
                           <cn cellml:units="dimensionless">1.25</cn>
                        </apply>
                        <apply>
                           <times/>
                           <ci>R</ci>
                           <ci>T</ci>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="rapid_delayed_rectifier_potassium_current">
      <variable units="nanoA" public_interface="out" name="i_Kr"/>
      <variable units="microS" name="g_Kr1" initial_value="0.0021"/>
      <variable units="microS" name="g_Kr2" initial_value="0.0013"/>
      <variable units="millivolt" public_interface="in" name="E_K"/>
      <variable units="second" public_interface="in" private_interface="out" name="time"/>
      <variable units="millivolt" public_interface="in" private_interface="out" name="V"/>
      <variable units="dimensionless" private_interface="in" name="xr1"/>
      <variable units="dimensionless" private_interface="in" name="xr2"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_Kr</ci>
            <apply>
               <times/>
               <apply>
                  <divide/>
                  <apply>
                     <times/>
                     <apply>
                        <plus/>
                        <apply>
                           <times/>
                           <ci>g_Kr1</ci>
                           <ci>xr1</ci>
                        </apply>
                        <apply>
                           <times/>
                           <ci>g_Kr2</ci>
                           <ci>xr2</ci>
                        </apply>
                     </apply>
                     <cn cellml:units="dimensionless">1</cn>
                  </apply>
                  <apply>
                     <plus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <plus/>
                              <ci>V</ci>
                              <cn cellml:units="millivolt">9</cn>
                           </apply>
                           <cn cellml:units="millivolt">22.4</cn>
                        </apply>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>E_K</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="rapid_delayed_rectifier_potassium_current_xr1_gate">
      <variable units="dimensionless" public_interface="out" name="xr1" initial_value="1.03e-5"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="per_second" name="alpha_xr1"/>
      <variable units="per_second" name="beta_xr1"/>
      <variable units="second" public_interface="in" name="time"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>alpha_xr1</ci>
            <apply>
               <divide/>
               <cn cellml:units="per_second">50</cn>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <minus/>
                           <apply>
                              <minus/>
                              <ci>V</ci>
                              <cn cellml:units="millivolt">5</cn>
                           </apply>
                        </apply>
                        <cn cellml:units="millivolt">9</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>beta_xr1</ci>
            <apply>
               <times/>
               <cn cellml:units="per_second">0.05</cn>
               <apply>
                  <exp/>
                  <apply>
                     <divide/>
                     <apply>
                        <minus/>
                        <apply>
                           <minus/>
                           <ci>V</ci>
                           <cn cellml:units="millivolt">20</cn>
                        </apply>
                     </apply>
                     <cn cellml:units="millivolt">15</cn>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>xr1</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>alpha_xr1</ci>
                  <apply>
                     <minus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <ci>xr1</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>beta_xr1</ci>
                  <ci>xr1</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="rapid_delayed_rectifier_potassium_current_xr2_gate">
      <variable units="dimensionless" public_interface="out" name="xr2" initial_value="2e-7"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="per_second" name="alpha_xr2"/>
      <variable units="per_second" name="beta_xr2"/>
      <variable units="second" public_interface="in" name="time"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>alpha_xr2</ci>
            <apply>
               <divide/>
               <cn cellml:units="per_second">50</cn>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <minus/>
                           <apply>
                              <minus/>
                              <ci>V</ci>
                              <cn cellml:units="millivolt">5</cn>
                           </apply>
                        </apply>
                        <cn cellml:units="millivolt">9</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>beta_xr2</ci>
            <apply>
               <times/>
               <cn cellml:units="per_second">0.4</cn>
               <apply>
                  <exp/>
                  <apply>
                     <minus/>
                     <apply>
                        <power/>
                        <apply>
                           <divide/>
                           <apply>
                              <plus/>
                              <ci>V</ci>
                              <cn cellml:units="millivolt">30</cn>
                           </apply>
                           <cn cellml:units="millivolt">30</cn>
                        </apply>
                        <cn cellml:units="dimensionless">3</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>xr2</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>alpha_xr2</ci>
                  <apply>
                     <minus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <ci>xr2</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>beta_xr2</ci>
                  <ci>xr2</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="slow_delayed_rectifier_potassium_current">
      <variable units="nanoA" public_interface="out" name="i_Ks"/>
      <variable units="microS" name="g_Ks" initial_value="0.0026"/>
      <variable units="millivolt" public_interface="in" name="E_Ks"/>
      <variable units="second" public_interface="in" private_interface="out" name="time"/>
      <variable units="millivolt" public_interface="in" private_interface="out" name="V"/>
      <variable units="dimensionless" private_interface="in" name="xs"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_Ks</ci>
            <apply>
               <times/>
               <ci>g_Ks</ci>
               <apply>
                  <power/>
                  <ci>xs</ci>
                  <cn cellml:units="dimensionless">2</cn>
               </apply>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>E_Ks</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="slow_delayed_rectifier_potassium_current_xs_gate">
      <variable units="dimensionless" public_interface="out" name="xs" initial_value="0.001302"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="per_second" name="alpha_xs"/>
      <variable units="per_second" name="beta_xs"/>
      <variable units="second" public_interface="in" name="time"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>alpha_xs</ci>
            <apply>
               <divide/>
               <cn cellml:units="per_second">14</cn>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <minus/>
                           <apply>
                              <minus/>
                              <ci>V</ci>
                              <cn cellml:units="millivolt">40</cn>
                           </apply>
                        </apply>
                        <cn cellml:units="millivolt">9</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>beta_xs</ci>
            <apply>
               <times/>
               <cn cellml:units="per_second">1</cn>
               <apply>
                  <exp/>
                  <apply>
                     <divide/>
                     <apply>
                        <minus/>
                        <ci>V</ci>
                     </apply>
                     <cn cellml:units="millivolt">45</cn>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>xs</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>alpha_xs</ci>
                  <apply>
                     <minus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <ci>xs</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>beta_xs</ci>
                  <ci>xs</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="ATP_dependent_potassium_current">
      <variable units="nanoA" public_interface="out" name="i_K_ATP"/>
      <variable units="microS" name="g_K_ATP" initial_value="0"/>
      <variable units="millimolar" name="K_ATP" initial_value="0.1"/>
      <variable units="millimolar" name="ATP" initial_value="0"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="second" public_interface="in" name="time"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_K_ATP</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <ci>g_K_ATP</ci>
                  <apply>
                     <plus/>
                     <ci>V</ci>
                     <cn cellml:units="millivolt">80</cn>
                  </apply>
               </apply>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <power/>
                     <apply>
                        <divide/>
                        <ci>ATP</ci>
                        <ci>K_ATP</ci>
                     </apply>
                     <cn cellml:units="dimensionless">2</cn>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="sodium_activated_potassium_current">
      <variable units="nanoA" public_interface="out" name="i_KNa"/>
      <variable units="microS" name="g_K_Na" initial_value="0"/>
      <variable units="millimolar" name="K_kna" initial_value="20"/>
      <variable units="millimolar" public_interface="in" name="Na_i"/>
      <variable units="millivolt" public_interface="in" name="E_K"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="second" public_interface="in" name="time"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_KNa</ci>
            <apply>
               <times/>
               <apply>
                  <divide/>
                  <apply>
                     <times/>
                     <ci>g_K_Na</ci>
                     <ci>Na_i</ci>
                  </apply>
                  <apply>
                     <plus/>
                     <ci>Na_i</ci>
                     <ci>K_kna</ci>
                  </apply>
               </apply>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>E_K</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="fast_sodium_current">
      <variable units="nanoA" public_interface="out" name="i_Na"/>
      <variable units="microS" name="g_Na" initial_value="2.5"/>
      <variable units="millivolt" public_interface="in" name="E_mh"/>
      <variable units="second" public_interface="in" private_interface="out" name="time"/>
      <variable units="millivolt" public_interface="in" private_interface="out" name="V"/>
      <variable units="dimensionless" private_interface="in" name="m"/>
      <variable units="dimensionless" private_interface="in" name="h"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_Na</ci>
            <apply>
               <times/>
               <ci>g_Na</ci>
               <apply>
                  <power/>
                  <ci>m</ci>
                  <cn cellml:units="dimensionless">3</cn>
               </apply>
               <ci>h</ci>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>E_mh</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="fast_sodium_current_m_gate">
      <variable units="dimensionless" public_interface="out" name="m" initial_value="0.0016203"/>
      <variable units="per_second" name="alpha_m"/>
      <variable units="per_second" name="beta_m"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="second" public_interface="in" name="time"/>
      <variable units="millivolt" name="delta_m" initial_value="1e-5"/>
      <variable units="millivolt" name="E0_m"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>E0_m</ci>
            <apply>
               <plus/>
               <ci>V</ci>
               <cn cellml:units="millivolt">41</cn>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>alpha_m</ci>
            <piecewise>
               <piece>
                  <cn cellml:units="per_second">2000</cn>
                  <apply>
                     <lt/>
                     <apply>
                        <abs/>
                        <ci>E0_m</ci>
                     </apply>
                     <ci>delta_m</ci>
                  </apply>
               </piece>
               <otherwise>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <cn cellml:units="per_millivolt_second">200</cn>
                        <ci>E0_m</ci>
                     </apply>
                     <apply>
                        <minus/>
                        <cn cellml:units="dimensionless">1</cn>
                        <apply>
                           <exp/>
                           <apply>
                              <times/>
                              <apply>
                                 <minus/>
                                 <cn cellml:units="per_millivolt">0.1</cn>
                              </apply>
                              <ci>E0_m</ci>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </otherwise>
            </piecewise>
         </apply>
         <apply>
            <eq/>
            <ci>beta_m</ci>
            <apply>
               <times/>
               <cn cellml:units="per_second">8000</cn>
               <apply>
                  <exp/>
                  <apply>
                     <times/>
                     <apply>
                        <minus/>
                        <cn cellml:units="per_millivolt">0.056</cn>
                     </apply>
                     <apply>
                        <plus/>
                        <ci>V</ci>
                        <cn cellml:units="millivolt">66</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>m</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>alpha_m</ci>
                  <apply>
                     <minus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <ci>m</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>beta_m</ci>
                  <ci>m</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="fast_sodium_current_h_gate">
      <variable units="dimensionless" public_interface="out" name="h" initial_value="0.9944036"/>
      <variable units="per_second" name="alpha_h"/>
      <variable units="per_second" name="beta_h"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="second" public_interface="in" name="time"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>alpha_h</ci>
            <apply>
               <times/>
               <cn cellml:units="per_second">20</cn>
               <apply>
                  <exp/>
                  <apply>
                     <times/>
                     <apply>
                        <minus/>
                        <cn cellml:units="per_millivolt">0.125</cn>
                     </apply>
                     <apply>
                        <plus/>
                        <ci>V</ci>
                        <cn cellml:units="millivolt">75</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>beta_h</ci>
            <apply>
               <divide/>
               <cn cellml:units="per_second">2000</cn>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <times/>
                     <cn cellml:units="dimensionless">320</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <times/>
                           <apply>
                              <minus/>
                              <cn cellml:units="per_millivolt">0.1</cn>
                           </apply>
                           <apply>
                              <plus/>
                              <ci>V</ci>
                              <cn cellml:units="millivolt">75</cn>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>h</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>alpha_h</ci>
                  <apply>
                     <minus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <ci>h</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>beta_h</ci>
                  <ci>h</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="persistent_sodium_current">
      <variable units="nanoA" public_interface="out" name="i_p_Na"/>
      <variable units="microS" name="g_pna" initial_value="0.004"/>
      <variable units="millivolt" public_interface="in" name="E_Na"/>
      <variable units="second" public_interface="in" name="time"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_p_Na</ci>
            <apply>
               <times/>
               <apply>
                  <divide/>
                  <apply>
                     <times/>
                     <ci>g_pna</ci>
                     <cn cellml:units="dimensionless">1</cn>
                  </apply>
                  <apply>
                     <plus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <minus/>
                              <apply>
                                 <plus/>
                                 <ci>V</ci>
                                 <cn cellml:units="millivolt">52</cn>
                              </apply>
                           </apply>
                           <cn cellml:units="millivolt">8</cn>
                        </apply>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>E_Na</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="sodium_background_current">
      <variable units="nanoA" public_interface="out" name="i_b_Na"/>
      <variable units="microS" name="g_bna" initial_value="0.0006"/>
      <variable units="millimolar" public_interface="in" name="Na_i"/>
      <variable units="millivolt" public_interface="in" name="E_Na"/>
      <variable units="second" public_interface="in" name="time"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_b_Na</ci>
            <apply>
               <times/>
               <ci>g_bna</ci>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>E_Na</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="L_type_Ca_channel">
      <variable units="nanoA" public_interface="out" name="i_Ca_L_Ca_cyt"/>
      <variable units="nanoA" public_interface="out" name="i_Ca_L_K_cyt"/>
      <variable units="nanoA" public_interface="out" name="i_Ca_L_Na_cyt"/>
      <variable units="nanoA" public_interface="out" name="i_Ca_L_Ca_ds"/>
      <variable units="nanoA" public_interface="out" name="i_Ca_L_K_ds"/>
      <variable units="nanoA" public_interface="out" name="i_Ca_L_Na_ds"/>
      <variable units="nanoA" public_interface="out" name="i_Ca_L"/>
      <variable units="nanoA_per_millimolar" name="P_Ca_L" initial_value="0.1"/>
      <variable units="dimensionless" name="P_CaK" initial_value="0.002"/>
      <variable units="dimensionless" name="P_CaNa" initial_value="0.01"/>
      <variable units="millimolar" public_interface="in" name="Ca_o"/>
      <variable units="millimolar" public_interface="in" private_interface="out" name="Ca_i"/>
      <variable units="millimolar" public_interface="in" private_interface="out" name="Ca_ds"/>
      <variable units="millimolar" public_interface="in" name="K_o"/>
      <variable units="millimolar" public_interface="in" name="K_i"/>
      <variable units="millimolar" public_interface="in" name="Na_o"/>
      <variable units="millimolar" public_interface="in" name="Na_i"/>
      <variable units="joule_per_kilomole_kelvin" public_interface="in" name="R"/>
      <variable units="coulomb_per_mole" public_interface="in" name="F"/>
      <variable units="kelvin" public_interface="in" name="T"/>
      <variable units="second" public_interface="in" private_interface="out" name="time"/>
      <variable units="millivolt" public_interface="in" private_interface="out" name="V"/>
      <variable units="dimensionless" private_interface="in" name="d"/>
      <variable units="dimensionless" private_interface="in" name="f"/>
      <variable units="dimensionless" private_interface="in" name="f2"/>
      <variable units="dimensionless" private_interface="in" name="f2ds"/>
      <variable units="millimolar" private_interface="out" name="Km_f2" initial_value="100000"/>
      <variable units="millimolar" private_interface="out" name="Km_f2ds" initial_value="0.001"/>
      <variable units="per_second" public_interface="out" private_interface="out" name="R_decay" initial_value="20"/>
      <variable units="dimensionless" name="FrICa" initial_value="1"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_Ca_L_Ca_cyt</ci>
            <apply>
               <times/>
               <apply>
                  <divide/>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <apply>
                           <minus/>
                           <cn cellml:units="dimensionless">1</cn>
                           <ci>FrICa</ci>
                        </apply>
                        <cn cellml:units="dimensionless">4</cn>
                        <ci>P_Ca_L</ci>
                        <ci>d</ci>
                        <ci>f</ci>
                        <ci>f2</ci>
                        <apply>
                           <minus/>
                           <ci>V</ci>
                           <cn cellml:units="millivolt">50</cn>
                        </apply>
                        <ci>F</ci>
                     </apply>
                     <apply>
                        <times/>
                        <ci>R</ci>
                        <ci>T</ci>
                     </apply>
                  </apply>
                  <apply>
                     <minus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <times/>
                              <apply>
                                 <minus/>
                                 <apply>
                                    <minus/>
                                    <ci>V</ci>
                                    <cn cellml:units="millivolt">50</cn>
                                 </apply>
                              </apply>
                              <ci>F</ci>
                              <cn cellml:units="dimensionless">2</cn>
                           </apply>
                           <apply>
                              <times/>
                              <ci>R</ci>
                              <ci>T</ci>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <minus/>
                  <apply>
                     <times/>
                     <ci>Ca_i</ci>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <times/>
                              <cn cellml:units="millivolt">100</cn>
                              <ci>F</ci>
                           </apply>
                           <apply>
                              <times/>
                              <ci>R</ci>
                              <ci>T</ci>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
                  <apply>
                     <times/>
                     <ci>Ca_o</ci>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <times/>
                              <apply>
                                 <minus/>
                                 <apply>
                                    <minus/>
                                    <ci>V</ci>
                                    <cn cellml:units="millivolt">50</cn>
                                 </apply>
                              </apply>
                              <ci>F</ci>
                              <cn cellml:units="dimensionless">2</cn>
                           </apply>
                           <apply>
                              <times/>
                              <ci>R</ci>
                              <ci>T</ci>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>i_Ca_L_K_cyt</ci>
            <apply>
               <times/>
               <apply>
                  <divide/>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <apply>
                           <minus/>
                           <cn cellml:units="dimensionless">1</cn>
                           <ci>FrICa</ci>
                        </apply>
                        <ci>P_CaK</ci>
                        <ci>P_Ca_L</ci>
                        <ci>d</ci>
                        <ci>f</ci>
                        <ci>f2</ci>
                        <apply>
                           <minus/>
                           <ci>V</ci>
                           <cn cellml:units="millivolt">50</cn>
                        </apply>
                        <ci>F</ci>
                     </apply>
                     <apply>
                        <times/>
                        <ci>R</ci>
                        <ci>T</ci>
                     </apply>
                  </apply>
                  <apply>
                     <minus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <times/>
                              <apply>
                                 <minus/>
                                 <apply>
                                    <minus/>
                                    <ci>V</ci>
                                    <cn cellml:units="millivolt">50</cn>
                                 </apply>
                              </apply>
                              <ci>F</ci>
                           </apply>
                           <apply>
                              <times/>
                              <ci>R</ci>
                              <ci>T</ci>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <minus/>
                  <apply>
                     <times/>
                     <ci>K_i</ci>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <times/>
                              <cn cellml:units="millivolt">50</cn>
                              <ci>F</ci>
                           </apply>
                           <apply>
                              <times/>
                              <ci>R</ci>
                              <ci>T</ci>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
                  <apply>
                     <times/>
                     <ci>K_o</ci>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <times/>
                              <apply>
                                 <minus/>
                                 <apply>
                                    <minus/>
                                    <ci>V</ci>
                                    <cn cellml:units="millivolt">50</cn>
                                 </apply>
                              </apply>
                              <ci>F</ci>
                           </apply>
                           <apply>
                              <times/>
                              <ci>R</ci>
                              <ci>T</ci>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>i_Ca_L_Na_cyt</ci>
            <apply>
               <times/>
               <apply>
                  <divide/>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <apply>
                           <minus/>
                           <cn cellml:units="dimensionless">1</cn>
                           <ci>FrICa</ci>
                        </apply>
                        <ci>P_CaNa</ci>
                        <ci>P_Ca_L</ci>
                        <ci>d</ci>
                        <ci>f</ci>
                        <ci>f2</ci>
                        <apply>
                           <minus/>
                           <ci>V</ci>
                           <cn cellml:units="millivolt">50</cn>
                        </apply>
                        <ci>F</ci>
                     </apply>
                     <apply>
                        <times/>
                        <ci>R</ci>
                        <ci>T</ci>
                     </apply>
                  </apply>
                  <apply>
                     <minus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <times/>
                              <apply>
                                 <minus/>
                                 <apply>
                                    <minus/>
                                    <ci>V</ci>
                                    <cn cellml:units="millivolt">50</cn>
                                 </apply>
                              </apply>
                              <ci>F</ci>
                           </apply>
                           <apply>
                              <times/>
                              <ci>R</ci>
                              <ci>T</ci>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <minus/>
                  <apply>
                     <times/>
                     <ci>Na_i</ci>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <times/>
                              <cn cellml:units="millivolt">50</cn>
                              <ci>F</ci>
                           </apply>
                           <apply>
                              <times/>
                              <ci>R</ci>
                              <ci>T</ci>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
                  <apply>
                     <times/>
                     <ci>Na_o</ci>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <times/>
                              <apply>
                                 <minus/>
                                 <apply>
                                    <minus/>
                                    <ci>V</ci>
                                    <cn cellml:units="millivolt">50</cn>
                                 </apply>
                              </apply>
                              <ci>F</ci>
                           </apply>
                           <apply>
                              <times/>
                              <ci>R</ci>
                              <ci>T</ci>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>i_Ca_L_Ca_ds</ci>
            <apply>
               <times/>
               <apply>
                  <divide/>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <ci>FrICa</ci>
                        <cn cellml:units="dimensionless">4</cn>
                        <ci>P_Ca_L</ci>
                        <ci>d</ci>
                        <ci>f</ci>
                        <ci>f2ds</ci>
                        <apply>
                           <minus/>
                           <ci>V</ci>
                           <cn cellml:units="millivolt">50</cn>
                        </apply>
                        <ci>F</ci>
                     </apply>
                     <apply>
                        <times/>
                        <ci>R</ci>
                        <ci>T</ci>
                     </apply>
                  </apply>
                  <apply>
                     <minus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <times/>
                              <apply>
                                 <minus/>
                                 <apply>
                                    <minus/>
                                    <ci>V</ci>
                                    <cn cellml:units="millivolt">50</cn>
                                 </apply>
                              </apply>
                              <ci>F</ci>
                              <cn cellml:units="dimensionless">2</cn>
                           </apply>
                           <apply>
                              <times/>
                              <ci>R</ci>
                              <ci>T</ci>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <minus/>
                  <apply>
                     <times/>
                     <ci>Ca_i</ci>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <times/>
                              <cn cellml:units="millivolt">100</cn>
                              <ci>F</ci>
                           </apply>
                           <apply>
                              <times/>
                              <ci>R</ci>
                              <ci>T</ci>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
                  <apply>
                     <times/>
                     <ci>Ca_o</ci>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <times/>
                              <apply>
                                 <minus/>
                                 <apply>
                                    <minus/>
                                    <ci>V</ci>
                                    <cn cellml:units="millivolt">50</cn>
                                 </apply>
                              </apply>
                              <ci>F</ci>
                              <cn cellml:units="dimensionless">2</cn>
                           </apply>
                           <apply>
                              <times/>
                              <ci>R</ci>
                              <ci>T</ci>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>i_Ca_L_K_ds</ci>
            <apply>
               <times/>
               <apply>
                  <divide/>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <ci>FrICa</ci>
                        <ci>P_CaK</ci>
                        <ci>P_Ca_L</ci>
                        <ci>d</ci>
                        <ci>f</ci>
                        <ci>f2ds</ci>
                        <apply>
                           <minus/>
                           <ci>V</ci>
                           <cn cellml:units="millivolt">50</cn>
                        </apply>
                        <ci>F</ci>
                     </apply>
                     <apply>
                        <times/>
                        <ci>R</ci>
                        <ci>T</ci>
                     </apply>
                  </apply>
                  <apply>
                     <minus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <times/>
                              <apply>
                                 <minus/>
                                 <apply>
                                    <minus/>
                                    <ci>V</ci>
                                    <cn cellml:units="millivolt">50</cn>
                                 </apply>
                              </apply>
                              <ci>F</ci>
                           </apply>
                           <apply>
                              <times/>
                              <ci>R</ci>
                              <ci>T</ci>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <minus/>
                  <apply>
                     <times/>
                     <ci>K_i</ci>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <times/>
                              <cn cellml:units="millivolt">50</cn>
                              <ci>F</ci>
                           </apply>
                           <apply>
                              <times/>
                              <ci>R</ci>
                              <ci>T</ci>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
                  <apply>
                     <times/>
                     <ci>K_o</ci>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <times/>
                              <apply>
                                 <minus/>
                                 <apply>
                                    <minus/>
                                    <ci>V</ci>
                                    <cn cellml:units="millivolt">50</cn>
                                 </apply>
                              </apply>
                              <ci>F</ci>
                           </apply>
                           <apply>
                              <times/>
                              <ci>R</ci>
                              <ci>T</ci>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>i_Ca_L_Na_ds</ci>
            <apply>
               <times/>
               <apply>
                  <divide/>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <ci>FrICa</ci>
                        <ci>P_CaNa</ci>
                        <ci>P_Ca_L</ci>
                        <ci>d</ci>
                        <ci>f</ci>
                        <ci>f2ds</ci>
                        <apply>
                           <minus/>
                           <ci>V</ci>
                           <cn cellml:units="millivolt">50</cn>
                        </apply>
                        <ci>F</ci>
                     </apply>
                     <apply>
                        <times/>
                        <ci>R</ci>
                        <ci>T</ci>
                     </apply>
                  </apply>
                  <apply>
                     <minus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <times/>
                              <apply>
                                 <minus/>
                                 <apply>
                                    <minus/>
                                    <ci>V</ci>
                                    <cn cellml:units="millivolt">50</cn>
                                 </apply>
                              </apply>
                              <ci>F</ci>
                           </apply>
                           <apply>
                              <times/>
                              <ci>R</ci>
                              <ci>T</ci>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <minus/>
                  <apply>
                     <times/>
                     <ci>Na_i</ci>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <times/>
                              <cn cellml:units="millivolt">50</cn>
                              <ci>F</ci>
                           </apply>
                           <apply>
                              <times/>
                              <ci>R</ci>
                              <ci>T</ci>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
                  <apply>
                     <times/>
                     <ci>Na_o</ci>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <times/>
                              <apply>
                                 <minus/>
                                 <apply>
                                    <minus/>
                                    <ci>V</ci>
                                    <cn cellml:units="millivolt">50</cn>
                                 </apply>
                              </apply>
                              <ci>F</ci>
                           </apply>
                           <apply>
                              <times/>
                              <ci>R</ci>
                              <ci>T</ci>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>i_Ca_L</ci>
            <apply>
               <plus/>
               <ci>i_Ca_L_Ca_cyt</ci>
               <ci>i_Ca_L_K_cyt</ci>
               <ci>i_Ca_L_Na_cyt</ci>
               <ci>i_Ca_L_Ca_ds</ci>
               <ci>i_Ca_L_K_ds</ci>
               <ci>i_Ca_L_Na_ds</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="L_type_Ca_channel_d_gate">
      <variable units="dimensionless" public_interface="out" name="d" initial_value="0"/>
      <variable units="per_second" name="alpha_d"/>
      <variable units="per_second" name="beta_d"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="second" public_interface="in" name="time"/>
      <variable units="millivolt" name="E0_d"/>
      <variable units="dimensionless" name="speed_d" initial_value="3"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>E0_d</ci>
            <apply>
               <minus/>
               <apply>
                  <plus/>
                  <ci>V</ci>
                  <cn cellml:units="millivolt">24</cn>
               </apply>
               <cn cellml:units="millivolt">5</cn>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>alpha_d</ci>
            <piecewise>
               <piece>
                  <cn cellml:units="per_second">120</cn>
                  <apply>
                     <lt/>
                     <apply>
                        <abs/>
                        <ci>E0_d</ci>
                     </apply>
                     <cn cellml:units="millivolt">0.0001</cn>
                  </apply>
               </piece>
               <otherwise>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <cn cellml:units="per_millivolt_second">30</cn>
                        <ci>E0_d</ci>
                     </apply>
                     <apply>
                        <minus/>
                        <cn cellml:units="dimensionless">1</cn>
                        <apply>
                           <exp/>
                           <apply>
                              <divide/>
                              <apply>
                                 <minus/>
                                 <ci>E0_d</ci>
                              </apply>
                              <cn cellml:units="millivolt">4</cn>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </otherwise>
            </piecewise>
         </apply>
         <apply>
            <eq/>
            <ci>beta_d</ci>
            <piecewise>
               <piece>
                  <cn cellml:units="per_second">120</cn>
                  <apply>
                     <lt/>
                     <apply>
                        <abs/>
                        <ci>E0_d</ci>
                     </apply>
                     <cn cellml:units="millivolt">0.0001</cn>
                  </apply>
               </piece>
               <otherwise>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <cn cellml:units="per_millivolt_second">12</cn>
                        <ci>E0_d</ci>
                     </apply>
                     <apply>
                        <minus/>
                        <apply>
                           <exp/>
                           <apply>
                              <divide/>
                              <ci>E0_d</ci>
                              <cn cellml:units="millivolt">10</cn>
                           </apply>
                        </apply>
                        <cn cellml:units="dimensionless">1</cn>
                     </apply>
                  </apply>
               </otherwise>
            </piecewise>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>d</ci>
            </apply>
            <apply>
               <times/>
               <ci>speed_d</ci>
               <apply>
                  <minus/>
                  <apply>
                     <times/>
                     <ci>alpha_d</ci>
                     <apply>
                        <minus/>
                        <cn cellml:units="dimensionless">1</cn>
                        <ci>d</ci>
                     </apply>
                  </apply>
                  <apply>
                     <times/>
                     <ci>beta_d</ci>
                     <ci>d</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="L_type_Ca_channel_f_gate">
      <variable units="dimensionless" public_interface="out" name="f" initial_value="1"/>
      <variable units="per_second" name="alpha_f"/>
      <variable units="per_second" name="beta_f"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="second" public_interface="in" name="time"/>
      <variable units="dimensionless" name="speed_f" initial_value="0.3"/>
      <variable units="millivolt" name="delta_f" initial_value="0.0001"/>
      <variable units="millivolt" name="E0_f"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>E0_f</ci>
            <apply>
               <plus/>
               <ci>V</ci>
               <cn cellml:units="millivolt">34</cn>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>alpha_f</ci>
            <piecewise>
               <piece>
                  <cn cellml:units="per_second">25</cn>
                  <apply>
                     <lt/>
                     <apply>
                        <abs/>
                        <ci>E0_f</ci>
                     </apply>
                     <ci>delta_f</ci>
                  </apply>
               </piece>
               <otherwise>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <cn cellml:units="per_millivolt_second">6.25</cn>
                        <ci>E0_f</ci>
                     </apply>
                     <apply>
                        <minus/>
                        <apply>
                           <exp/>
                           <apply>
                              <divide/>
                              <ci>E0_f</ci>
                              <cn cellml:units="millivolt">4</cn>
                           </apply>
                        </apply>
                        <cn cellml:units="dimensionless">1</cn>
                     </apply>
                  </apply>
               </otherwise>
            </piecewise>
         </apply>
         <apply>
            <eq/>
            <ci>beta_f</ci>
            <apply>
               <divide/>
               <cn cellml:units="per_second">12</cn>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <times/>
                           <apply>
                              <minus/>
                              <cn cellml:units="dimensionless">1</cn>
                           </apply>
                           <apply>
                              <plus/>
                              <ci>V</ci>
                              <cn cellml:units="millivolt">34</cn>
                           </apply>
                        </apply>
                        <cn cellml:units="millivolt">4</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>f</ci>
            </apply>
            <apply>
               <times/>
               <ci>speed_f</ci>
               <apply>
                  <minus/>
                  <apply>
                     <times/>
                     <ci>alpha_f</ci>
                     <apply>
                        <minus/>
                        <cn cellml:units="dimensionless">1</cn>
                        <ci>f</ci>
                     </apply>
                  </apply>
                  <apply>
                     <times/>
                     <ci>beta_f</ci>
                     <ci>f</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="L_type_Ca_channel_f2_gate">
      <variable units="dimensionless" public_interface="out" name="f2" initial_value="0.9349197"/>
      <variable units="millimolar" public_interface="in" name="Km_f2"/>
      <variable units="millimolar" public_interface="in" name="Ca_i"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="second" public_interface="in" name="time"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>f2</ci>
            </apply>
            <apply>
               <minus/>
               <cn cellml:units="per_second">1</cn>
               <apply>
                  <times/>
                  <cn cellml:units="per_second">1</cn>
                  <apply>
                     <plus/>
                     <apply>
                        <divide/>
                        <ci>Ca_i</ci>
                        <apply>
                           <plus/>
                           <ci>Km_f2</ci>
                           <ci>Ca_i</ci>
                        </apply>
                     </apply>
                     <ci>f2</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="L_type_Ca_channel_f2ds_gate">
      <variable units="dimensionless" public_interface="out" name="f2ds" initial_value="0.9651958"/>
      <variable units="millimolar" public_interface="in" name="Km_f2ds"/>
      <variable units="per_second" public_interface="in" name="R_decay"/>
      <variable units="millimolar" public_interface="in" name="Ca_ds"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="second" public_interface="in" name="time"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>f2ds</ci>
            </apply>
            <apply>
               <times/>
               <ci>R_decay</ci>
               <apply>
                  <minus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <plus/>
                     <apply>
                        <divide/>
                        <ci>Ca_ds</ci>
                        <apply>
                           <plus/>
                           <ci>Km_f2ds</ci>
                           <ci>Ca_ds</ci>
                        </apply>
                     </apply>
                     <ci>f2ds</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="calcium_background_current">
      <variable units="nanoA" public_interface="out" name="i_b_Ca"/>
      <variable units="microS" name="g_bca" initial_value="0.00025"/>
      <variable units="millivolt" public_interface="in" name="E_Ca"/>
      <variable units="second" public_interface="in" name="time"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_b_Ca</ci>
            <apply>
               <times/>
               <ci>g_bca</ci>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>E_Ca</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="transient_outward_current">
      <variable units="nanoA" public_interface="out" name="i_to"/>
      <variable units="microS" name="g_to" initial_value="0.005"/>
      <variable units="dimensionless" name="g_tos" initial_value="0"/>
      <variable units="millivolt" public_interface="in" name="E_K"/>
      <variable units="millivolt" public_interface="in" private_interface="out" name="V"/>
      <variable units="second" public_interface="in" private_interface="out" name="time"/>
      <variable units="dimensionless" private_interface="in" name="s"/>
      <variable units="dimensionless" private_interface="in" name="r"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_to</ci>
            <apply>
               <times/>
               <ci>g_to</ci>
               <apply>
                  <plus/>
                  <ci>g_tos</ci>
                  <apply>
                     <times/>
                     <ci>s</ci>
                     <apply>
                        <minus/>
                        <cn cellml:units="dimensionless">1</cn>
                        <ci>g_tos</ci>
                     </apply>
                  </apply>
               </apply>
               <ci>r</ci>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>E_K</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="transient_outward_current_s_gate">
      <variable units="dimensionless" public_interface="out" name="s" initial_value="0.9948645"/>
      <variable units="per_second" name="alpha_s"/>
      <variable units="per_second" name="beta_s"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="second" public_interface="in" name="time"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>alpha_s</ci>
            <apply>
               <times/>
               <cn cellml:units="per_second">0.033</cn>
               <apply>
                  <exp/>
                  <apply>
                     <divide/>
                     <apply>
                        <minus/>
                        <ci>V</ci>
                     </apply>
                     <cn cellml:units="millivolt">17</cn>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>beta_s</ci>
            <apply>
               <divide/>
               <cn cellml:units="per_second">33</cn>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <times/>
                        <apply>
                           <minus/>
                           <cn cellml:units="per_millivolt">0.125</cn>
                        </apply>
                        <apply>
                           <plus/>
                           <ci>V</ci>
                           <cn cellml:units="millivolt">10</cn>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>s</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>alpha_s</ci>
                  <apply>
                     <minus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <ci>s</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>beta_s</ci>
                  <ci>s</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="transient_outward_current_r_gate">
      <variable units="dimensionless" public_interface="out" name="r" initial_value="0"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="second" public_interface="in" name="time"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>r</ci>
            </apply>
            <apply>
               <times/>
               <cn cellml:units="per_second">333</cn>
               <apply>
                  <minus/>
                  <apply>
                     <divide/>
                     <cn cellml:units="dimensionless">1</cn>
                     <apply>
                        <plus/>
                        <cn cellml:units="dimensionless">1</cn>
                        <apply>
                           <exp/>
                           <apply>
                              <divide/>
                              <apply>
                                 <minus/>
                                 <apply>
                                    <plus/>
                                    <ci>V</ci>
                                    <cn cellml:units="millivolt">4</cn>
                                 </apply>
                              </apply>
                              <cn cellml:units="millivolt">5</cn>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
                  <ci>r</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="ACh_dependent_potassium_current">
      <variable units="nanoA" public_interface="out" name="i_K_ACh"/>
      <variable units="microS" name="g_KACh" initial_value="0"/>
      <variable units="millimolar" name="ACh" initial_value="5"/>
      <variable units="millimolar" name="K_D" initial_value="0.00013"/>
      <variable units="millivolt" public_interface="in" name="E_K"/>
      <variable units="millimolar" public_interface="in" name="K_o"/>
      <variable units="millimolar" public_interface="in" name="K_mk1"/>
      <variable units="joule_per_kilomole_kelvin" public_interface="in" name="R"/>
      <variable units="coulomb_per_mole" public_interface="in" name="F"/>
      <variable units="kelvin" public_interface="in" name="T"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="second" public_interface="in" private_interface="out" name="time"/>
      <variable units="dimensionless" private_interface="in" name="x_ACh"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_K_ACh</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <ci>g_KACh</ci>
                  <apply>
                     <divide/>
                     <ci>K_o</ci>
                     <apply>
                        <plus/>
                        <ci>K_o</ci>
                        <ci>K_mk1</ci>
                     </apply>
                  </apply>
                  <ci>x_ACh</ci>
                  <apply>
                     <divide/>
                     <cn cellml:units="dimensionless">1</cn>
                     <apply>
                        <plus/>
                        <cn cellml:units="dimensionless">1</cn>
                        <apply>
                           <power/>
                           <apply>
                              <divide/>
                              <ci>K_D</ci>
                              <ci>ACh</ci>
                           </apply>
                           <cn cellml:units="dimensionless">2</cn>
                        </apply>
                     </apply>
                  </apply>
                  <apply>
                     <minus/>
                     <ci>V</ci>
                     <ci>E_K</ci>
                  </apply>
               </apply>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <times/>
                           <cn cellml:units="dimensionless">2</cn>
                           <ci>F</ci>
                           <apply>
                              <minus/>
                              <ci>V</ci>
                              <apply>
                                 <plus/>
                                 <ci>E_K</ci>
                                 <cn cellml:units="millivolt">10</cn>
                              </apply>
                           </apply>
                        </apply>
                        <apply>
                           <times/>
                           <ci>R</ci>
                           <ci>T</ci>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="ACh_dependent_potassium_current_xACh_gate">
      <variable units="dimensionless" public_interface="out" name="x_ACh" initial_value="0"/>
      <variable units="per_second" name="alpha_ACh" initial_value="0.5"/>
      <variable units="per_second" name="beta_ACh" initial_value="0.5"/>
      <variable units="second" public_interface="in" name="time"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>x_ACh</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>alpha_ACh</ci>
                  <apply>
                     <minus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <ci>x_ACh</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>beta_ACh</ci>
                  <ci>x_ACh</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="sodium_potassium_pump">
      <variable units="nanoA" public_interface="out" name="i_NaK"/>
      <variable units="nanoA" name="i_NaK_max" initial_value="0.7"/>
      <variable units="millimolar" name="K_mK" initial_value="1"/>
      <variable units="millimolar" name="K_mNa" initial_value="40"/>
      <variable units="millimolar" public_interface="in" name="K_o"/>
      <variable units="millimolar" public_interface="in" name="Na_i"/>
      <variable units="second" public_interface="in" name="time"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_NaK</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <ci>i_NaK_max</ci>
                        <ci>K_o</ci>
                     </apply>
                     <apply>
                        <plus/>
                        <ci>K_mK</ci>
                        <ci>K_o</ci>
                     </apply>
                  </apply>
                  <ci>Na_i</ci>
               </apply>
               <apply>
                  <plus/>
                  <ci>K_mNa</ci>
                  <ci>Na_i</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="sodium_calcium_exchanger">
      <variable units="nanoA" public_interface="out" name="i_NaCa"/>
      <variable units="nanoA" public_interface="out" name="i_NaCa_cyt"/>
      <variable units="nanoA" public_interface="out" name="i_NaCa_ds"/>
      <variable units="nanoA" name="k_NaCa" initial_value="0.0005"/>
      <variable units="dimensionless" name="n_NaCa" initial_value="3"/>
      <variable units="dimensionless" name="d_NaCa" initial_value="0"/>
      <variable units="dimensionless" name="gamma" initial_value="0.5"/>
      <variable units="dimensionless" name="FRiNaCa" initial_value="0.001"/>
      <variable units="millimolar" public_interface="in" name="Na_i"/>
      <variable units="millimolar" public_interface="in" name="Na_o"/>
      <variable units="millimolar" public_interface="in" name="Ca_i"/>
      <variable units="millimolar" public_interface="in" name="Ca_ds"/>
      <variable units="millimolar" public_interface="in" name="Ca_o"/>
      <variable units="joule_per_kilomole_kelvin" public_interface="in" name="R"/>
      <variable units="coulomb_per_mole" public_interface="in" name="F"/>
      <variable units="kelvin" public_interface="in" name="T"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="second" public_interface="in" name="time"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_NaCa_cyt</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <apply>
                     <minus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <ci>FRiNaCa</ci>
                  </apply>
                  <ci>k_NaCa</ci>
                  <apply>
                     <minus/>
                     <apply>
                        <times/>
                        <apply>
                           <exp/>
                           <apply>
                              <divide/>
                              <apply>
                                 <times/>
                                 <ci>gamma</ci>
                                 <apply>
                                    <minus/>
                                    <ci>n_NaCa</ci>
                                    <cn cellml:units="dimensionless">2</cn>
                                 </apply>
                                 <ci>V</ci>
                                 <ci>F</ci>
                              </apply>
                              <apply>
                                 <times/>
                                 <ci>R</ci>
                                 <ci>T</ci>
                              </apply>
                           </apply>
                        </apply>
                        <apply>
                           <power/>
                           <ci>Na_i</ci>
                           <ci>n_NaCa</ci>
                        </apply>
                        <ci>Ca_o</ci>
                     </apply>
                     <apply>
                        <times/>
                        <apply>
                           <exp/>
                           <apply>
                              <divide/>
                              <apply>
                                 <times/>
                                 <apply>
                                    <minus/>
                                    <ci>gamma</ci>
                                    <cn cellml:units="dimensionless">1</cn>
                                 </apply>
                                 <apply>
                                    <minus/>
                                    <ci>n_NaCa</ci>
                                    <cn cellml:units="dimensionless">2</cn>
                                 </apply>
                                 <ci>V</ci>
                                 <ci>F</ci>
                              </apply>
                              <apply>
                                 <times/>
                                 <ci>R</ci>
                                 <ci>T</ci>
                              </apply>
                           </apply>
                        </apply>
                        <apply>
                           <power/>
                           <ci>Na_o</ci>
                           <ci>n_NaCa</ci>
                        </apply>
                        <ci>Ca_i</ci>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <apply>
                     <plus/>
                     <cn cellml:units="millimolar4">1</cn>
                     <apply>
                        <times/>
                        <ci>d_NaCa</ci>
                        <apply>
                           <plus/>
                           <apply>
                              <times/>
                              <ci>Ca_i</ci>
                              <apply>
                                 <power/>
                                 <ci>Na_o</ci>
                                 <ci>n_NaCa</ci>
                              </apply>
                           </apply>
                           <apply>
                              <times/>
                              <ci>Ca_o</ci>
                              <apply>
                                 <power/>
                                 <ci>Na_i</ci>
                                 <ci>n_NaCa</ci>
                              </apply>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
                  <apply>
                     <plus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <apply>
                        <divide/>
                        <ci>Ca_i</ci>
                        <cn cellml:units="millimolar">0.0069</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>i_NaCa_ds</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <ci>FRiNaCa</ci>
                  <ci>k_NaCa</ci>
                  <apply>
                     <minus/>
                     <apply>
                        <times/>
                        <apply>
                           <exp/>
                           <apply>
                              <divide/>
                              <apply>
                                 <times/>
                                 <ci>gamma</ci>
                                 <apply>
                                    <minus/>
                                    <ci>n_NaCa</ci>
                                    <cn cellml:units="dimensionless">2</cn>
                                 </apply>
                                 <ci>V</ci>
                                 <ci>F</ci>
                              </apply>
                              <apply>
                                 <times/>
                                 <ci>R</ci>
                                 <ci>T</ci>
                              </apply>
                           </apply>
                        </apply>
                        <apply>
                           <power/>
                           <ci>Na_i</ci>
                           <ci>n_NaCa</ci>
                        </apply>
                        <ci>Ca_o</ci>
                     </apply>
                     <apply>
                        <times/>
                        <apply>
                           <exp/>
                           <apply>
                              <divide/>
                              <apply>
                                 <times/>
                                 <apply>
                                    <minus/>
                                    <ci>gamma</ci>
                                    <cn cellml:units="dimensionless">1</cn>
                                 </apply>
                                 <apply>
                                    <minus/>
                                    <ci>n_NaCa</ci>
                                    <cn cellml:units="dimensionless">2</cn>
                                 </apply>
                                 <ci>V</ci>
                                 <ci>F</ci>
                              </apply>
                              <apply>
                                 <times/>
                                 <ci>R</ci>
                                 <ci>T</ci>
                              </apply>
                           </apply>
                        </apply>
                        <apply>
                           <power/>
                           <ci>Na_o</ci>
                           <ci>n_NaCa</ci>
                        </apply>
                        <ci>Ca_ds</ci>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <apply>
                     <plus/>
                     <cn cellml:units="millimolar4">1</cn>
                     <apply>
                        <times/>
                        <ci>d_NaCa</ci>
                        <apply>
                           <plus/>
                           <apply>
                              <times/>
                              <ci>Ca_ds</ci>
                              <apply>
                                 <power/>
                                 <ci>Na_o</ci>
                                 <ci>n_NaCa</ci>
                              </apply>
                           </apply>
                           <apply>
                              <times/>
                              <ci>Ca_o</ci>
                              <apply>
                                 <power/>
                                 <ci>Na_i</ci>
                                 <ci>n_NaCa</ci>
                              </apply>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
                  <apply>
                     <plus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <apply>
                        <divide/>
                        <ci>Ca_ds</ci>
                        <cn cellml:units="millimolar">0.0069</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>i_NaCa</ci>
            <apply>
               <plus/>
               <ci>i_NaCa_cyt</ci>
               <ci>i_NaCa_ds</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="sarcoplasmic_reticulum_calcium_pump">
      <variable units="millimolar_per_second" public_interface="out" name="i_up"/>
      <variable units="dimensionless" name="K_1"/>
      <variable units="millimolar" name="K_2"/>
      <variable units="millimolar" name="K_cyca" initial_value="0.0003"/>
      <variable units="dimensionless" name="K_xcs" initial_value="0.4"/>
      <variable units="millimolar" name="K_srca" initial_value="0.5"/>
      <variable units="millimolar_per_second" name="alpha_up" initial_value="0.4"/>
      <variable units="millimolar_per_second" name="beta_up" initial_value="0.03"/>
      <variable units="millimolar" public_interface="in" name="Ca_i"/>
      <variable units="millimolar" public_interface="in" name="Ca_up"/>
      <variable units="second" public_interface="in" name="time"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>K_1</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <ci>K_cyca</ci>
                  <ci>K_xcs</ci>
               </apply>
               <ci>K_srca</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>K_2</ci>
            <apply>
               <plus/>
               <ci>Ca_i</ci>
               <apply>
                  <times/>
                  <ci>Ca_up</ci>
                  <ci>K_1</ci>
               </apply>
               <apply>
                  <times/>
                  <ci>K_cyca</ci>
                  <ci>K_xcs</ci>
               </apply>
               <ci>K_cyca</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>i_up</ci>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <apply>
                     <divide/>
                     <ci>Ca_i</ci>
                     <ci>K_2</ci>
                  </apply>
                  <ci>alpha_up</ci>
               </apply>
               <apply>
                  <times/>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <ci>Ca_up</ci>
                        <ci>K_1</ci>
                     </apply>
                     <ci>K_2</ci>
                  </apply>
                  <ci>beta_up</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="calcium_translocation">
      <variable units="millimolar_per_second" public_interface="out" name="i_trans"/>
      <variable units="millimolar" public_interface="in" name="Ca_rel"/>
      <variable units="millimolar" public_interface="in" name="Ca_up"/>
      <variable units="second" public_interface="in" name="time"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_trans</ci>
            <apply>
               <times/>
               <cn cellml:units="per_second">50</cn>
               <apply>
                  <minus/>
                  <ci>Ca_up</ci>
                  <ci>Ca_rel</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="calcium_release">
      <variable units="millimolar_per_second" public_interface="out" name="i_rel"/>
      <variable units="dimensionless" name="VoltDep"/>
      <variable units="dimensionless" name="RegBindSite"/>
      <variable units="dimensionless" name="CaiReg"/>
      <variable units="dimensionless" name="CadsReg"/>
      <variable units="per_second" name="ActRate"/>
      <variable units="per_second" name="InactRate"/>
      <variable units="per_second" name="K_leak_rate" initial_value="0.005"/>
      <variable units="per_second" name="SRLeak"/>
      <variable units="per_second" name="K_m_rel" initial_value="250"/>
      <variable units="millimolar" name="K_m_Ca_cyt" initial_value="0.0005"/>
      <variable units="millimolar" name="K_m_Ca_ds" initial_value="0.01"/>
      <variable units="dimensionless" name="PrecFrac"/>
      <variable units="dimensionless" name="ActFrac" initial_value="0.001914"/>
      <variable units="dimensionless" name="ProdFrac" initial_value="0.2854569"/>
      <variable units="dimensionless" name="SpeedRel"/>
      <variable units="per_micrometre" name="gama_SR_SL" initial_value="2.5"/>
      <variable units="dimensionless" name="gama_SR_IT" initial_value="2.5"/>
      <variable units="second" public_interface="in" name="time"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="millimolar" public_interface="in" name="Ca_i"/>
      <variable units="millimolar" public_interface="in" name="Ca_ds"/>
      <variable units="millimolar" public_interface="in" name="Ca_rel"/>
      <variable units="dimensionless" public_interface="in" name="isometric_tension"/>
      <variable units="micrometre" public_interface="in" name="sarcomere_length"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>VoltDep</ci>
            <apply>
               <exp/>
               <apply>
                  <times/>
                  <cn cellml:units="per_millivolt">0.08</cn>
                  <apply>
                     <minus/>
                     <ci>V</ci>
                     <cn cellml:units="millivolt">40</cn>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>CaiReg</ci>
            <apply>
               <divide/>
               <ci>Ca_i</ci>
               <apply>
                  <plus/>
                  <ci>Ca_i</ci>
                  <ci>K_m_Ca_cyt</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>CadsReg</ci>
            <apply>
               <divide/>
               <ci>Ca_ds</ci>
               <apply>
                  <plus/>
                  <ci>Ca_ds</ci>
                  <ci>K_m_Ca_ds</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>RegBindSite</ci>
            <apply>
               <plus/>
               <ci>CaiReg</ci>
               <apply>
                  <times/>
                  <apply>
                     <minus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <ci>CaiReg</ci>
                  </apply>
                  <ci>CadsReg</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>ActRate</ci>
            <apply>
               <plus/>
               <apply>
                  <times/>
                  <cn cellml:units="per_second">0</cn>
                  <ci>VoltDep</ci>
               </apply>
               <apply>
                  <times/>
                  <cn cellml:units="per_second">500</cn>
                  <apply>
                     <power/>
                     <ci>RegBindSite</ci>
                     <cn cellml:units="dimensionless">2</cn>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>InactRate</ci>
            <apply>
               <plus/>
               <cn cellml:units="per_second">60</cn>
               <apply>
                  <times/>
                  <cn cellml:units="per_second">500</cn>
                  <apply>
                     <power/>
                     <ci>RegBindSite</ci>
                     <cn cellml:units="dimensionless">2</cn>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>SpeedRel</ci>
            <piecewise>
               <piece>
                  <cn cellml:units="dimensionless">5</cn>
                  <apply>
                     <lt/>
                     <ci>V</ci>
                     <apply>
                        <minus/>
                        <cn cellml:units="millivolt">50</cn>
                     </apply>
                  </apply>
               </piece>
               <otherwise>
                  <cn cellml:units="dimensionless">1</cn>
               </otherwise>
            </piecewise>
         </apply>
         <apply>
            <eq/>
            <ci>SRLeak</ci>
            <piecewise>
               <piece>
                  <apply>
                     <times/>
                     <ci>K_leak_rate</ci>
                     <apply>
                        <exp/>
                        <apply>
                           <times/>
                           <ci>gama_SR_IT</ci>
                           <ci>isometric_tension</ci>
                        </apply>
                     </apply>
                  </apply>
                  <apply>
                     <gt/>
                     <ci>isometric_tension</ci>
                     <cn cellml:units="dimensionless">0</cn>
                  </apply>
               </piece>
               <otherwise>
                  <apply>
                     <times/>
                     <ci>K_leak_rate</ci>
                     <apply>
                        <exp/>
                        <apply>
                           <times/>
                           <ci>gama_SR_SL</ci>
                           <ci>sarcomere_length</ci>
                        </apply>
                     </apply>
                  </apply>
               </otherwise>
            </piecewise>
         </apply>
         <apply>
            <eq/>
            <ci>PrecFrac</ci>
            <apply>
               <minus/>
               <apply>
                  <minus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <ci>ActFrac</ci>
               </apply>
               <ci>ProdFrac</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>ActFrac</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>PrecFrac</ci>
                  <ci>SpeedRel</ci>
                  <ci>ActRate</ci>
               </apply>
               <apply>
                  <times/>
                  <ci>ActFrac</ci>
                  <ci>SpeedRel</ci>
                  <ci>InactRate</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>ProdFrac</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>ActFrac</ci>
                  <ci>SpeedRel</ci>
                  <ci>InactRate</ci>
               </apply>
               <apply>
                  <times/>
                  <ci>SpeedRel</ci>
                  <cn cellml:units="per_second">1</cn>
                  <ci>ProdFrac</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>i_rel</ci>
            <apply>
               <times/>
               <apply>
                  <plus/>
                  <apply>
                     <times/>
                     <apply>
                        <power/>
                        <apply>
                           <divide/>
                           <ci>ActFrac</ci>
                           <apply>
                              <plus/>
                              <ci>ActFrac</ci>
                              <cn cellml:units="dimensionless">0.25</cn>
                           </apply>
                        </apply>
                        <cn cellml:units="dimensionless">2</cn>
                     </apply>
                     <ci>K_m_rel</ci>
                  </apply>
                  <ci>SRLeak</ci>
               </apply>
               <ci>Ca_rel</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="extracellular_sodium_concentration">
      <variable units="millimolar" public_interface="out" name="Na_o" initial_value="140"/>
   </component>
   <component name="intracellular_sodium_concentration">
      <variable units="millimolar" public_interface="out" name="Na_i" initial_value="7.3321223"/>
      <variable units="micrometre3" public_interface="in" name="V_i"/>
      <variable units="second" public_interface="in" name="time"/>
      <variable units="coulomb_per_mole" public_interface="in" name="F"/>
      <variable units="nanoA" public_interface="in" name="i_Na"/>
      <variable units="nanoA" public_interface="in" name="i_b_Na"/>
      <variable units="nanoA" public_interface="in" name="i_p_Na"/>
      <variable units="nanoA" public_interface="in" name="i_Ca_L_Na_cyt"/>
      <variable units="nanoA" public_interface="in" name="i_Ca_L_Na_ds"/>
      <variable units="nanoA" public_interface="in" name="i_NaK"/>
      <variable units="nanoA" public_interface="in" name="i_NaCa_cyt"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>Na_i</ci>
            </apply>
            <apply>
               <times/>
               <apply>
                  <divide/>
                  <apply>
                     <minus/>
                     <cn cellml:units="mA_nA">1</cn>
                  </apply>
                  <apply>
                     <times/>
                     <cn cellml:units="litre_micrometre3">1</cn>
                     <ci>V_i</ci>
                     <ci>F</ci>
                  </apply>
               </apply>
               <apply>
                  <plus/>
                  <ci>i_Na</ci>
                  <ci>i_p_Na</ci>
                  <ci>i_b_Na</ci>
                  <apply>
                     <times/>
                     <cn cellml:units="dimensionless">3</cn>
                     <ci>i_NaK</ci>
                  </apply>
                  <apply>
                     <times/>
                     <cn cellml:units="dimensionless">3</cn>
                     <ci>i_NaCa_cyt</ci>
                  </apply>
                  <ci>i_Ca_L_Na_cyt</ci>
                  <ci>i_Ca_L_Na_ds</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="extracellular_calcium_concentration">
      <variable units="millimolar" public_interface="out" name="Ca_o" initial_value="2"/>
   </component>
   <component name="extracellular_potassium_concentration">
      <variable units="millimolar" public_interface="out" name="K_o" initial_value="4"/>
   </component>
   <component name="intracellular_potassium_concentration">
      <variable units="millimolar" public_interface="out" name="K_i" initial_value="136.5644281"/>
      <variable units="micrometre3" public_interface="in" name="V_i"/>
      <variable units="nanoA" public_interface="in" name="i_K1"/>
      <variable units="nanoA" public_interface="in" name="i_Kr"/>
      <variable units="nanoA" public_interface="in" name="i_Ks"/>
      <variable units="nanoA" public_interface="in" name="i_Ca_L_K_cyt"/>
      <variable units="nanoA" public_interface="in" name="i_Ca_L_K_ds"/>
      <variable units="nanoA" public_interface="in" name="i_to"/>
      <variable units="nanoA" public_interface="in" name="i_NaK"/>
      <variable units="second" public_interface="in" name="time"/>
      <variable units="coulomb_per_mole" public_interface="in" name="F"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>K_i</ci>
            </apply>
            <apply>
               <times/>
               <apply>
                  <divide/>
                  <apply>
                     <minus/>
                     <cn cellml:units="mA_nA">1</cn>
                  </apply>
                  <apply>
                     <times/>
                     <cn cellml:units="litre_micrometre3">1</cn>
                     <ci>V_i</ci>
                     <ci>F</ci>
                  </apply>
               </apply>
               <apply>
                  <minus/>
                  <apply>
                     <plus/>
                     <ci>i_K1</ci>
                     <ci>i_Kr</ci>
                     <ci>i_Ks</ci>
                     <ci>i_Ca_L_K_cyt</ci>
                     <ci>i_Ca_L_K_ds</ci>
                     <ci>i_to</ci>
                  </apply>
                  <apply>
                     <times/>
                     <cn cellml:units="dimensionless">2</cn>
                     <ci>i_NaK</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="intracellular_calcium_concentration">
      <variable units="millimolar" public_interface="out" name="Ca_i" initial_value="1.4e-5"/>
      <variable units="millimolar" public_interface="out" name="Ca_ds" initial_value="1.88e-5"/>
      <variable units="millimolar" public_interface="out" name="Ca_up" initial_value="0.4531889"/>
      <variable units="millimolar" public_interface="out" name="Ca_rel" initial_value="0.4481927"/>
      <variable units="millimolar" public_interface="out" name="Ca_Calmod" initial_value="0.0005555"/>
      <variable units="millimolar" public_interface="out" name="Ca_Trop" initial_value="0.0002"/>
      <variable units="millimolar" public_interface="out" name="Calmod" initial_value="0.02"/>
      <variable units="millimolar" public_interface="out" name="Trop" initial_value="0.05"/>
      <variable units="per_millimolar_second" name="alpha_Calmod" initial_value="100000"/>
      <variable units="per_second" name="beta_Calmod" initial_value="50"/>
      <variable units="per_millimolar_second" name="alpha_Trop"/>
      <variable units="per_second" name="beta_Trop" initial_value="200"/>
      <variable units="per_micrometre" name="gama_Trop_SL" initial_value="1.5"/>
      <variable units="per_millimolar_second" name="KTrop" initial_value="5000"/>
      <variable units="micrometre" name="radius" initial_value="12"/>
      <variable units="micrometre" name="length" initial_value="74"/>
      <variable units="micrometre3" public_interface="out" name="V_i"/>
      <variable units="micrometre3" name="V_Cell"/>
      <variable units="dimensionless" name="V_i_ratio"/>
      <variable units="dimensionless" name="V_ds_ratio" initial_value="0.1"/>
      <variable units="dimensionless" name="V_rel_ratio" initial_value="0.1"/>
      <variable units="dimensionless" name="V_e_ratio" initial_value="0.4"/>
      <variable units="dimensionless" name="V_up_ratio" initial_value="0.01"/>
      <variable units="per_second" name="Kdecay" initial_value="10"/>
      <variable units="millimolar_per_second" public_interface="in" name="i_up"/>
      <variable units="millimolar_per_second" public_interface="in" name="i_trans"/>
      <variable units="millimolar_per_second" public_interface="in" name="i_rel"/>
      <variable units="nanoA" public_interface="in" name="i_NaCa_cyt"/>
      <variable units="nanoA" public_interface="in" name="i_NaCa_ds"/>
      <variable units="nanoA" public_interface="in" name="i_Ca_L_Ca_cyt"/>
      <variable units="nanoA" public_interface="in" name="i_Ca_L_Ca_ds"/>
      <variable units="nanoA" public_interface="in" name="i_b_Ca"/>
      <variable units="coulomb_per_mole" public_interface="in" name="F"/>
      <variable units="dimensionless" public_interface="in" name="isometric_tension"/>
      <variable units="micrometre" public_interface="in" name="sarcomere_length"/>
      <variable units="second" public_interface="in" name="time"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>V_Cell</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <cn cellml:units="dimensionless">3.141592654</cn>
                  <apply>
                     <power/>
                     <apply>
                        <divide/>
                        <ci>radius</ci>
                        <cn cellml:units="dimensionless">1000</cn>
                     </apply>
                     <cn cellml:units="dimensionless">2</cn>
                  </apply>
                  <ci>length</ci>
               </apply>
               <cn cellml:units="dimensionless">1000</cn>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>V_i_ratio</ci>
            <apply>
               <minus/>
               <apply>
                  <minus/>
                  <apply>
                     <minus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <ci>V_e_ratio</ci>
                  </apply>
                  <ci>V_up_ratio</ci>
               </apply>
               <ci>V_rel_ratio</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>V_i</ci>
            <apply>
               <times/>
               <ci>V_Cell</ci>
               <ci>V_i_ratio</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>Ca_i</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <minus/>
                  <apply>
                     <minus/>
                     <apply>
                        <plus/>
                        <apply>
                           <times/>
                           <apply>
                              <divide/>
                              <apply>
                                 <minus/>
                                 <cn cellml:units="mA_nA">1</cn>
                              </apply>
                              <apply>
                                 <times/>
                                 <cn cellml:units="dimensionless">2</cn>
                                 <cn cellml:units="litre_micrometre3">1</cn>
                                 <ci>V_i</ci>
                                 <ci>F</ci>
                              </apply>
                           </apply>
                           <apply>
                              <minus/>
                              <apply>
                                 <minus/>
                                 <apply>
                                    <plus/>
                                    <ci>i_Ca_L_Ca_cyt</ci>
                                    <ci>i_b_Ca</ci>
                                 </apply>
                                 <apply>
                                    <times/>
                                    <cn cellml:units="dimensionless">2</cn>
                                    <ci>i_NaCa_cyt</ci>
                                 </apply>
                              </apply>
                              <apply>
                                 <times/>
                                 <cn cellml:units="dimensionless">2</cn>
                                 <ci>i_NaCa_ds</ci>
                              </apply>
                           </apply>
                        </apply>
                        <apply>
                           <times/>
                           <ci>Ca_ds</ci>
                           <ci>V_ds_ratio</ci>
                           <ci>Kdecay</ci>
                        </apply>
                        <apply>
                           <divide/>
                           <apply>
                              <times/>
                              <ci>i_rel</ci>
                              <ci>V_rel_ratio</ci>
                           </apply>
                           <ci>V_i_ratio</ci>
                        </apply>
                     </apply>
                     <apply>
                        <diff/>
                        <bvar>
                           <ci>time</ci>
                        </bvar>
                        <ci>Ca_Calmod</ci>
                     </apply>
                  </apply>
                  <apply>
                     <diff/>
                     <bvar>
                        <ci>time</ci>
                     </bvar>
                     <ci>Ca_Trop</ci>
                  </apply>
               </apply>
               <ci>i_up</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>Ca_ds</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <divide/>
                  <apply>
                     <times/>
                     <apply>
                        <minus/>
                        <cn cellml:units="mA_nA">1</cn>
                     </apply>
                     <ci>i_Ca_L_Ca_ds</ci>
                  </apply>
                  <apply>
                     <times/>
                     <cn cellml:units="dimensionless">2</cn>
                     <cn cellml:units="litre_micrometre3">1</cn>
                     <ci>V_ds_ratio</ci>
                     <ci>V_i</ci>
                     <ci>F</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>Ca_ds</ci>
                  <ci>Kdecay</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>Ca_up</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <apply>
                     <divide/>
                     <ci>V_i_ratio</ci>
                     <ci>V_up_ratio</ci>
                  </apply>
                  <ci>i_up</ci>
               </apply>
               <ci>i_trans</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>Ca_rel</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <apply>
                     <divide/>
                     <ci>V_up_ratio</ci>
                     <ci>V_rel_ratio</ci>
                  </apply>
                  <ci>i_trans</ci>
               </apply>
               <ci>i_rel</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>Ca_Calmod</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>alpha_Calmod</ci>
                  <ci>Ca_i</ci>
                  <apply>
                     <minus/>
                     <ci>Calmod</ci>
                     <ci>Ca_Calmod</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>beta_Calmod</ci>
                  <ci>Ca_Calmod</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>alpha_Trop</ci>
            <apply>
               <times/>
               <ci>KTrop</ci>
               <apply>
                  <exp/>
                  <apply>
                     <times/>
                     <ci>gama_Trop_SL</ci>
                     <ci>sarcomere_length</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>Ca_Trop</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>alpha_Trop</ci>
                  <ci>Ca_i</ci>
                  <apply>
                     <minus/>
                     <ci>Trop</ci>
                     <ci>Ca_Trop</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>beta_Trop</ci>
                  <ci>Ca_Trop</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="contraction">
      <variable units="dimensionless" name="light_chain" initial_value="3.32e-5"/>
      <variable units="dimensionless" name="cross_bridge" initial_value="8.09e-5"/>
      <variable units="per_second" name="KCont1" initial_value="12000"/>
      <variable units="dimensionless" name="XCont2"/>
      <variable units="dimensionless" name="XCont1"/>
      <variable units="per_second" name="KCont2" initial_value="100"/>
      <variable units="per_second" name="KCont3" initial_value="60"/>
      <variable units="per_second" name="KCont4" initial_value="25"/>
      <variable units="micrometre" public_interface="out" name="sarcomere_length" initial_value="2"/>
      <variable units="per_micrometre" name="cross_bridge_density" initial_value="0.05"/>
      <variable units="dimensionless" name="tension_rest"/>
      <variable units="dimensionless" name="tension_active"/>
      <variable units="micrometre" name="overlap"/>
      <variable units="dimensionless" name="cross_bridge_availability"/>
      <variable units="dimensionless" public_interface="out" name="isometric_tension"/>
      <variable units="millimolar" public_interface="in" name="Ca_Calmod"/>
      <variable units="millimolar" public_interface="in" name="Ca_Trop"/>
      <variable units="millimolar" public_interface="in" name="Calmod"/>
      <variable units="millimolar" public_interface="in" name="Trop"/>
      <variable units="second" public_interface="in" name="time"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>XCont2</ci>
            <apply>
               <divide/>
               <ci>Ca_Calmod</ci>
               <ci>Calmod</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>XCont1</ci>
            <apply>
               <divide/>
               <ci>Ca_Trop</ci>
               <ci>Trop</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>light_chain</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>KCont1</ci>
                  <apply>
                     <power/>
                     <ci>XCont1</ci>
                     <cn cellml:units="dimensionless">2</cn>
                  </apply>
                  <ci>XCont2</ci>
                  <apply>
                     <minus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <ci>light_chain</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>KCont2</ci>
                  <ci>light_chain</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>cross_bridge</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>KCont3</ci>
                  <ci>light_chain</ci>
                  <apply>
                     <minus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <ci>cross_bridge</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>KCont4</ci>
                  <ci>cross_bridge</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>tension_rest</ci>
            <apply>
               <times/>
               <cn cellml:units="dimensionless">2e-4</cn>
               <apply>
                  <exp/>
                  <apply>
                     <times/>
                     <cn cellml:units="per_micrometre">2</cn>
                     <ci>sarcomere_length</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>tension_active</ci>
            <piecewise>
               <piece>
                  <apply>
                     <minus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <times/>
                           <apply>
                              <minus/>
                              <cn cellml:units="per_micrometre">3</cn>
                           </apply>
                           <apply>
                              <minus/>
                              <ci>sarcomere_length</ci>
                              <cn cellml:units="micrometre">1</cn>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
                  <apply>
                     <gt/>
                     <ci>sarcomere_length</ci>
                     <cn cellml:units="micrometre">1</cn>
                  </apply>
               </piece>
               <otherwise>
                  <cn cellml:units="dimensionless">0</cn>
               </otherwise>
            </piecewise>
         </apply>
         <apply>
            <eq/>
            <ci>overlap</ci>
            <piecewise>
               <piece>
                  <apply>
                     <minus/>
                     <cn cellml:units="micrometre">1</cn>
                     <apply>
                        <times/>
                        <cn cellml:units="dimensionless">0.625</cn>
                        <apply>
                           <minus/>
                           <ci>sarcomere_length</ci>
                           <cn cellml:units="micrometre">2</cn>
                        </apply>
                     </apply>
                  </apply>
                  <apply>
                     <gt/>
                     <ci>sarcomere_length</ci>
                     <cn cellml:units="micrometre">2</cn>
                  </apply>
               </piece>
               <otherwise>
                  <cn cellml:units="micrometre">1</cn>
               </otherwise>
            </piecewise>
         </apply>
         <apply>
            <eq/>
            <ci>cross_bridge_availability</ci>
            <apply>
               <times/>
               <ci>tension_active</ci>
               <ci>overlap</ci>
               <ci>cross_bridge_density</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>isometric_tension</ci>
            <apply>
               <plus/>
               <apply>
                  <times/>
                  <ci>cross_bridge</ci>
                  <ci>cross_bridge_availability</ci>
               </apply>
               <ci>tension_rest</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="stretch_current">
      <variable units="nanoA" public_interface="out" name="i_stretch"/>
      <variable units="nanoA" name="i_Ca_stretch"/>
      <variable units="nanoA" name="i_K_stretch"/>
      <variable units="nanoA" name="i_Na_stretch"/>
      <variable units="nanoA" name="i_An_stretch"/>
      <variable units="per_micrometre" name="gama_SAC_SL" initial_value="2.5"/>
      <variable units="dimensionless" name="gama_SAC_IT" initial_value="2.5"/>
      <variable units="micrometre" name="SLHST" initial_value="2"/>
      <variable units="dimensionless" name="ITHST" initial_value="1"/>
      <variable units="microS" name="g_Ca_stretch" initial_value="0.01"/>
      <variable units="microS" name="g_K_stretch" initial_value="0.01"/>
      <variable units="microS" name="g_Na_stretch" initial_value="0.01"/>
      <variable units="microS" name="g_An_stretch" initial_value="0.01"/>
      <variable units="millivolt" name="E_An_stretch" initial_value="-20"/>
      <variable units="dimensionless" name="f_stretch"/>
      <variable units="dimensionless" public_interface="in" name="isometric_tension"/>
      <variable units="micrometre" public_interface="in" name="sarcomere_length"/>
      <variable units="millivolt" public_interface="in" name="E_Ca"/>
      <variable units="millivolt" public_interface="in" name="E_K"/>
      <variable units="millivolt" public_interface="in" name="E_Na"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="second" public_interface="in" name="time"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>f_stretch</ci>
            <piecewise>
               <piece>
                  <apply>
                     <divide/>
                     <cn cellml:units="dimensionless">1</cn>
                     <apply>
                        <plus/>
                        <cn cellml:units="dimensionless">1</cn>
                        <apply>
                           <exp/>
                           <apply>
                              <times/>
                              <apply>
                                 <minus/>
                                 <cn cellml:units="dimensionless">2</cn>
                              </apply>
                              <ci>gama_SAC_IT</ci>
                              <apply>
                                 <minus/>
                                 <ci>isometric_tension</ci>
                                 <ci>ITHST</ci>
                              </apply>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
                  <apply>
                     <gt/>
                     <ci>isometric_tension</ci>
                     <cn cellml:units="dimensionless">0</cn>
                  </apply>
               </piece>
               <otherwise>
                  <apply>
                     <divide/>
                     <cn cellml:units="dimensionless">1</cn>
                     <apply>
                        <plus/>
                        <cn cellml:units="dimensionless">1</cn>
                        <apply>
                           <exp/>
                           <apply>
                              <times/>
                              <apply>
                                 <minus/>
                                 <cn cellml:units="dimensionless">2</cn>
                              </apply>
                              <ci>gama_SAC_SL</ci>
                              <apply>
                                 <minus/>
                                 <ci>sarcomere_length</ci>
                                 <ci>SLHST</ci>
                              </apply>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </otherwise>
            </piecewise>
         </apply>
         <apply>
            <eq/>
            <ci>i_Ca_stretch</ci>
            <apply>
               <times/>
               <ci>g_Ca_stretch</ci>
               <ci>f_stretch</ci>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>E_Ca</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>i_Na_stretch</ci>
            <apply>
               <times/>
               <ci>g_Na_stretch</ci>
               <ci>f_stretch</ci>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>E_Na</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>i_K_stretch</ci>
            <apply>
               <times/>
               <ci>g_K_stretch</ci>
               <ci>f_stretch</ci>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>E_K</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>i_An_stretch</ci>
            <apply>
               <times/>
               <ci>g_An_stretch</ci>
               <ci>f_stretch</ci>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>E_An_stretch</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>i_stretch</ci>
            <apply>
               <plus/>
               <ci>i_Ca_stretch</ci>
               <ci>i_Na_stretch</ci>
               <ci>i_K_stretch</ci>
               <ci>i_An_stretch</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="fibroblast">
      <variable units="microS" name="g_fibro_junct" initial_value="2.9e-4"/>
      <variable units="microS" name="g_fibro" initial_value="2e-4"/>
      <variable units="microF" name="c_fibro" initial_value="1e-5"/>
      <variable units="microS" name="g_fibro_stretch" initial_value="0"/>
      <variable units="millivolt" name="E_fibro_stretch" initial_value="0"/>
      <variable units="millivolt" name="V_fibro" initial_value="-20"/>
      <variable units="nanoA" name="i_fibro"/>
      <variable units="nanoA" name="i_fibro_junct"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="second" public_interface="in" name="time"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_fibro</ci>
            <apply>
               <plus/>
               <apply>
                  <times/>
                  <ci>g_fibro</ci>
                  <apply>
                     <plus/>
                     <ci>V_fibro</ci>
                     <cn cellml:units="millivolt">20</cn>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>g_fibro_stretch</ci>
                  <apply>
                     <minus/>
                     <ci>V_fibro</ci>
                     <ci>E_fibro_stretch</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>i_fibro_junct</ci>
            <apply>
               <times/>
               <apply>
                  <minus/>
                  <ci>g_fibro_junct</ci>
               </apply>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>V_fibro</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>V_fibro</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <apply>
                     <plus/>
                     <ci>i_fibro</ci>
                     <ci>i_fibro_junct</ci>
                  </apply>
               </apply>
               <ci>c_fibro</ci>
            </apply>
         </apply>
      </math>
   </component>
   <group>
      <relationship_ref relationship="containment"/>
      <component_ref component="membrane">
         <component_ref component="fast_sodium_current">
            <component_ref component="fast_sodium_current_m_gate"/>
            <component_ref component="fast_sodium_current_h_gate"/>
         </component_ref>
         <component_ref component="L_type_Ca_channel">
            <component_ref component="L_type_Ca_channel_d_gate"/>
            <component_ref component="L_type_Ca_channel_f_gate"/>
            <component_ref component="L_type_Ca_channel_f2_gate"/>
            <component_ref component="L_type_Ca_channel_f2ds_gate"/>
         </component_ref>
         <component_ref component="rapid_delayed_rectifier_potassium_current">
            <component_ref component="rapid_delayed_rectifier_potassium_current_xr1_gate"/>
            <component_ref component="rapid_delayed_rectifier_potassium_current_xr2_gate"/>
         </component_ref>
         <component_ref component="slow_delayed_rectifier_potassium_current">
            <component_ref component="slow_delayed_rectifier_potassium_current_xs_gate"/>
         </component_ref>
         <component_ref component="time_independent_potassium_current"/>
         <component_ref component="ATP_dependent_potassium_current"/>
         <component_ref component="sodium_activated_potassium_current"/>
         <component_ref component="persistent_sodium_current"/>
         <component_ref component="sodium_background_current"/>
         <component_ref component="calcium_background_current"/>
         <component_ref component="transient_outward_current">
            <component_ref component="transient_outward_current_s_gate"/>
            <component_ref component="transient_outward_current_r_gate"/>
         </component_ref>
         <component_ref component="ACh_dependent_potassium_current">
            <component_ref component="ACh_dependent_potassium_current_xACh_gate"/>
         </component_ref>
         <component_ref component="sodium_potassium_pump"/>
         <component_ref component="sodium_calcium_exchanger"/>
         <component_ref component="sarcoplasmic_reticulum_calcium_pump"/>
         <component_ref component="calcium_translocation"/>
         <component_ref component="calcium_release"/>
         <component_ref component="extracellular_sodium_concentration"/>
         <component_ref component="intracellular_sodium_concentration"/>
         <component_ref component="extracellular_potassium_concentration"/>
         <component_ref component="intracellular_potassium_concentration"/>
         <component_ref component="extracellular_calcium_concentration"/>
         <component_ref component="intracellular_calcium_concentration"/>
         <component_ref component="reversal_potentials"/>
         <component_ref component="contraction"/>
         <component_ref component="stretch_current"/>
         <component_ref component="fibroblast"/>
      </component_ref>
   </group>
   <group>
      <relationship_ref relationship="encapsulation"/>
      <component_ref component="fast_sodium_current">
         <component_ref component="fast_sodium_current_m_gate"/>
         <component_ref component="fast_sodium_current_h_gate"/>
      </component_ref>
      <component_ref component="L_type_Ca_channel">
         <component_ref component="L_type_Ca_channel_d_gate"/>
         <component_ref component="L_type_Ca_channel_f_gate"/>
         <component_ref component="L_type_Ca_channel_f2_gate"/>
         <component_ref component="L_type_Ca_channel_f2ds_gate"/>
      </component_ref>
      <component_ref component="transient_outward_current">
         <component_ref component="transient_outward_current_s_gate"/>
         <component_ref component="transient_outward_current_r_gate"/>
      </component_ref>
      <component_ref component="rapid_delayed_rectifier_potassium_current">
         <component_ref component="rapid_delayed_rectifier_potassium_current_xr1_gate"/>
         <component_ref component="rapid_delayed_rectifier_potassium_current_xr2_gate"/>
      </component_ref>
      <component_ref component="slow_delayed_rectifier_potassium_current">
         <component_ref component="slow_delayed_rectifier_potassium_current_xs_gate"/>
      </component_ref>
      <component_ref component="ACh_dependent_potassium_current">
         <component_ref component="ACh_dependent_potassium_current_xACh_gate"/>
      </component_ref>
   </group>
   <connection>
      <map_components component_2="environment" component_1="membrane"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="environment" component_1="reversal_potentials"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="environment" component_1="fast_sodium_current"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="environment" component_1="L_type_Ca_channel"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="environment" component_1="rapid_delayed_rectifier_potassium_current"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="environment" component_1="slow_delayed_rectifier_potassium_current"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="environment" component_1="ATP_dependent_potassium_current"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="environment" component_1="ACh_dependent_potassium_current"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="environment" component_1="sodium_activated_potassium_current"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="environment" component_1="time_independent_potassium_current"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="environment" component_1="persistent_sodium_current"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="environment" component_1="sodium_background_current"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="environment" component_1="calcium_background_current"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="environment" component_1="transient_outward_current"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="environment" component_1="sodium_potassium_pump"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="environment" component_1="sodium_calcium_exchanger"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="environment" component_1="intracellular_sodium_concentration"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="environment" component_1="intracellular_potassium_concentration"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="environment" component_1="intracellular_calcium_concentration"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="environment" component_1="sarcoplasmic_reticulum_calcium_pump"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="environment" component_1="calcium_translocation"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="environment" component_1="calcium_release"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="environment" component_1="contraction"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="environment" component_1="stretch_current"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="environment" component_1="fibroblast"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="membrane" component_1="calcium_release"/>
      <map_variables variable_2="V" variable_1="V"/>
   </connection>
   <connection>
      <map_components component_2="reversal_potentials" component_1="membrane"/>
      <map_variables variable_2="R" variable_1="R"/>
      <map_variables variable_2="T" variable_1="T"/>
      <map_variables variable_2="F" variable_1="F"/>
   </connection>
   <connection>
      <map_components component_2="fast_sodium_current" component_1="membrane"/>
      <map_variables variable_2="V" variable_1="V"/>
      <map_variables variable_2="i_Na" variable_1="i_Na"/>
   </connection>
   <connection>
      <map_components component_2="L_type_Ca_channel" component_1="membrane"/>
      <map_variables variable_2="V" variable_1="V"/>
      <map_variables variable_2="i_Ca_L_Ca_cyt" variable_1="i_Ca_L_Ca_cyt"/>
      <map_variables variable_2="i_Ca_L_Ca_ds" variable_1="i_Ca_L_Ca_ds"/>
      <map_variables variable_2="i_Ca_L_Na_cyt" variable_1="i_Ca_L_Na_cyt"/>
      <map_variables variable_2="i_Ca_L_Na_ds" variable_1="i_Ca_L_Na_ds"/>
      <map_variables variable_2="i_Ca_L_K_cyt" variable_1="i_Ca_L_K_cyt"/>
      <map_variables variable_2="i_Ca_L_K_ds" variable_1="i_Ca_L_K_ds"/>
      <map_variables variable_2="R" variable_1="R"/>
      <map_variables variable_2="T" variable_1="T"/>
      <map_variables variable_2="F" variable_1="F"/>
   </connection>
   <connection>
      <map_components component_2="rapid_delayed_rectifier_potassium_current" component_1="membrane"/>
      <map_variables variable_2="V" variable_1="V"/>
      <map_variables variable_2="i_Kr" variable_1="i_Kr"/>
   </connection>
   <connection>
      <map_components component_2="slow_delayed_rectifier_potassium_current" component_1="membrane"/>
      <map_variables variable_2="V" variable_1="V"/>
      <map_variables variable_2="i_Ks" variable_1="i_Ks"/>
   </connection>
   <connection>
      <map_components component_2="transient_outward_current" component_1="membrane"/>
      <map_variables variable_2="V" variable_1="V"/>
      <map_variables variable_2="i_to" variable_1="i_to"/>
   </connection>
   <connection>
      <map_components component_2="time_independent_potassium_current" component_1="membrane"/>
      <map_variables variable_2="V" variable_1="V"/>
      <map_variables variable_2="i_K1" variable_1="i_K1"/>
      <map_variables variable_2="R" variable_1="R"/>
      <map_variables variable_2="T" variable_1="T"/>
      <map_variables variable_2="F" variable_1="F"/>
   </connection>
   <connection>
      <map_components component_2="sodium_activated_potassium_current" component_1="membrane"/>
      <map_variables variable_2="V" variable_1="V"/>
   </connection>
   <connection>
      <map_components component_2="ATP_dependent_potassium_current" component_1="membrane"/>
      <map_variables variable_2="V" variable_1="V"/>
      <map_variables variable_2="i_K_ATP" variable_1="i_K_ATP"/>
   </connection>
   <connection>
      <map_components component_2="ACh_dependent_potassium_current" component_1="membrane"/>
      <map_variables variable_2="V" variable_1="V"/>
      <map_variables variable_2="i_K_ACh" variable_1="i_K_ACh"/>
      <map_variables variable_2="R" variable_1="R"/>
      <map_variables variable_2="T" variable_1="T"/>
      <map_variables variable_2="F" variable_1="F"/>
   </connection>
   <connection>
      <map_components component_2="persistent_sodium_current" component_1="membrane"/>
      <map_variables variable_2="V" variable_1="V"/>
      <map_variables variable_2="i_p_Na" variable_1="i_p_Na"/>
   </connection>
   <connection>
      <map_components component_2="sodium_background_current" component_1="membrane"/>
      <map_variables variable_2="V" variable_1="V"/>
      <map_variables variable_2="i_b_Na" variable_1="i_b_Na"/>
   </connection>
   <connection>
      <map_components component_2="calcium_background_current" component_1="membrane"/>
      <map_variables variable_2="V" variable_1="V"/>
      <map_variables variable_2="i_b_Ca" variable_1="i_b_Ca"/>
   </connection>
   <connection>
      <map_components component_2="sodium_potassium_pump" component_1="membrane"/>
      <map_variables variable_2="i_NaK" variable_1="i_NaK"/>
   </connection>
   <connection>
      <map_components component_2="sodium_calcium_exchanger" component_1="membrane"/>
      <map_variables variable_2="V" variable_1="V"/>
      <map_variables variable_2="i_NaCa_cyt" variable_1="i_NaCa_cyt"/>
      <map_variables variable_2="i_NaCa_ds" variable_1="i_NaCa_ds"/>
      <map_variables variable_2="R" variable_1="R"/>
      <map_variables variable_2="T" variable_1="T"/>
      <map_variables variable_2="F" variable_1="F"/>
   </connection>
   <connection>
      <map_components component_2="stretch_current" component_1="membrane"/>
      <map_variables variable_2="V" variable_1="V"/>
      <map_variables variable_2="i_stretch" variable_1="i_stretch"/>
   </connection>
   <connection>
      <map_components component_2="fibroblast" component_1="membrane"/>
      <map_variables variable_2="V" variable_1="V"/>
   </connection>
   <connection>
      <map_components component_2="intracellular_sodium_concentration" component_1="membrane"/>
      <map_variables variable_2="F" variable_1="F"/>
   </connection>
   <connection>
      <map_components component_2="intracellular_potassium_concentration" component_1="membrane"/>
      <map_variables variable_2="F" variable_1="F"/>
   </connection>
   <connection>
      <map_components component_2="intracellular_calcium_concentration" component_1="membrane"/>
      <map_variables variable_2="F" variable_1="F"/>
   </connection>
   <connection>
      <map_components component_2="sodium_background_current" component_1="reversal_potentials"/>
      <map_variables variable_2="E_Na" variable_1="E_Na"/>
   </connection>
   <connection>
      <map_components component_2="persistent_sodium_current" component_1="reversal_potentials"/>
      <map_variables variable_2="E_Na" variable_1="E_Na"/>
   </connection>
   <connection>
      <map_components component_2="time_independent_potassium_current" component_1="reversal_potentials"/>
      <map_variables variable_2="E_K" variable_1="E_K"/>
   </connection>
   <connection>
      <map_components component_2="rapid_delayed_rectifier_potassium_current" component_1="reversal_potentials"/>
      <map_variables variable_2="E_K" variable_1="E_K"/>
   </connection>
   <connection>
      <map_components component_2="slow_delayed_rectifier_potassium_current" component_1="reversal_potentials"/>
      <map_variables variable_2="E_Ks" variable_1="E_Ks"/>
   </connection>
   <connection>
      <map_components component_2="sodium_activated_potassium_current" component_1="reversal_potentials"/>
      <map_variables variable_2="E_K" variable_1="E_K"/>
   </connection>
   <connection>
      <map_components component_2="transient_outward_current" component_1="reversal_potentials"/>
      <map_variables variable_2="E_K" variable_1="E_K"/>
   </connection>
   <connection>
      <map_components component_2="ACh_dependent_potassium_current" component_1="reversal_potentials"/>
      <map_variables variable_2="E_K" variable_1="E_K"/>
   </connection>
   <connection>
      <map_components component_2="calcium_background_current" component_1="reversal_potentials"/>
      <map_variables variable_2="E_Ca" variable_1="E_Ca"/>
   </connection>
   <connection>
      <map_components component_2="fast_sodium_current" component_1="reversal_potentials"/>
      <map_variables variable_2="E_mh" variable_1="E_mh"/>
   </connection>
   <connection>
      <map_components component_2="extracellular_potassium_concentration" component_1="reversal_potentials"/>
      <map_variables variable_2="K_o" variable_1="K_o"/>
   </connection>
   <connection>
      <map_components component_2="extracellular_sodium_concentration" component_1="reversal_potentials"/>
      <map_variables variable_2="Na_o" variable_1="Na_o"/>
   </connection>
   <connection>
      <map_components component_2="intracellular_potassium_concentration" component_1="reversal_potentials"/>
      <map_variables variable_2="K_i" variable_1="K_i"/>
   </connection>
   <connection>
      <map_components component_2="intracellular_sodium_concentration" component_1="reversal_potentials"/>
      <map_variables variable_2="Na_i" variable_1="Na_i"/>
   </connection>
   <connection>
      <map_components component_2="extracellular_calcium_concentration" component_1="reversal_potentials"/>
      <map_variables variable_2="Ca_o" variable_1="Ca_o"/>
   </connection>
   <connection>
      <map_components component_2="intracellular_calcium_concentration" component_1="reversal_potentials"/>
      <map_variables variable_2="Ca_i" variable_1="Ca_i"/>
   </connection>
   <connection>
      <map_components component_2="stretch_current" component_1="reversal_potentials"/>
      <map_variables variable_2="E_Ca" variable_1="E_Ca"/>
      <map_variables variable_2="E_K" variable_1="E_K"/>
      <map_variables variable_2="E_Na" variable_1="E_Na"/>
   </connection>
   <connection>
      <map_components component_2="extracellular_potassium_concentration" component_1="time_independent_potassium_current"/>
      <map_variables variable_2="K_o" variable_1="K_o"/>
   </connection>
   <connection>
      <map_components component_2="extracellular_potassium_concentration" component_1="ACh_dependent_potassium_current"/>
      <map_variables variable_2="K_o" variable_1="K_o"/>
   </connection>
   <connection>
      <map_components component_2="intracellular_potassium_concentration" component_1="time_independent_potassium_current"/>
      <map_variables variable_2="i_K1" variable_1="i_K1"/>
   </connection>
   <connection>
      <map_components component_2="intracellular_potassium_concentration" component_1="rapid_delayed_rectifier_potassium_current"/>
      <map_variables variable_2="i_Kr" variable_1="i_Kr"/>
   </connection>
   <connection>
      <map_components component_2="intracellular_potassium_concentration" component_1="slow_delayed_rectifier_potassium_current"/>
      <map_variables variable_2="i_Ks" variable_1="i_Ks"/>
   </connection>
   <connection>
      <map_components component_2="intracellular_sodium_concentration" component_1="fast_sodium_current"/>
      <map_variables variable_2="i_Na" variable_1="i_Na"/>
   </connection>
   <connection>
      <map_components component_2="intracellular_sodium_concentration" component_1="sodium_background_current"/>
      <map_variables variable_2="i_b_Na" variable_1="i_b_Na"/>
      <map_variables variable_2="Na_i" variable_1="Na_i"/>
   </connection>
   <connection>
      <map_components component_2="intracellular_sodium_concentration" component_1="persistent_sodium_current"/>
      <map_variables variable_2="i_p_Na" variable_1="i_p_Na"/>
   </connection>
   <connection>
      <map_components component_2="intracellular_sodium_concentration" component_1="L_type_Ca_channel"/>
      <map_variables variable_2="i_Ca_L_Na_cyt" variable_1="i_Ca_L_Na_cyt"/>
      <map_variables variable_2="i_Ca_L_Na_ds" variable_1="i_Ca_L_Na_ds"/>
      <map_variables variable_2="Na_i" variable_1="Na_i"/>
   </connection>
   <connection>
      <map_components component_2="intracellular_calcium_concentration" component_1="L_type_Ca_channel"/>
      <map_variables variable_2="i_Ca_L_Ca_cyt" variable_1="i_Ca_L_Ca_cyt"/>
      <map_variables variable_2="i_Ca_L_Ca_ds" variable_1="i_Ca_L_Ca_ds"/>
      <map_variables variable_2="Ca_i" variable_1="Ca_i"/>
      <map_variables variable_2="Ca_ds" variable_1="Ca_ds"/>
   </connection>
   <connection>
      <map_components component_2="intracellular_potassium_concentration" component_1="L_type_Ca_channel"/>
      <map_variables variable_2="i_Ca_L_K_cyt" variable_1="i_Ca_L_K_cyt"/>
      <map_variables variable_2="i_Ca_L_K_ds" variable_1="i_Ca_L_K_ds"/>
      <map_variables variable_2="K_i" variable_1="K_i"/>
   </connection>
   <connection>
      <map_components component_2="extracellular_potassium_concentration" component_1="L_type_Ca_channel"/>
      <map_variables variable_2="K_o" variable_1="K_o"/>
   </connection>
   <connection>
      <map_components component_2="extracellular_sodium_concentration" component_1="L_type_Ca_channel"/>
      <map_variables variable_2="Na_o" variable_1="Na_o"/>
   </connection>
   <connection>
      <map_components component_2="extracellular_calcium_concentration" component_1="L_type_Ca_channel"/>
      <map_variables variable_2="Ca_o" variable_1="Ca_o"/>
   </connection>
   <connection>
      <map_components component_2="intracellular_calcium_concentration" component_1="calcium_background_current"/>
      <map_variables variable_2="i_b_Ca" variable_1="i_b_Ca"/>
   </connection>
   <connection>
      <map_components component_2="intracellular_sodium_concentration" component_1="sodium_activated_potassium_current"/>
      <map_variables variable_2="Na_i" variable_1="Na_i"/>
   </connection>
   <connection>
      <map_components component_2="intracellular_sodium_concentration" component_1="sodium_potassium_pump"/>
      <map_variables variable_2="i_NaK" variable_1="i_NaK"/>
      <map_variables variable_2="Na_i" variable_1="Na_i"/>
   </connection>
   <connection>
      <map_components component_2="extracellular_potassium_concentration" component_1="sodium_potassium_pump"/>
      <map_variables variable_2="K_o" variable_1="K_o"/>
   </connection>
   <connection>
      <map_components component_2="intracellular_potassium_concentration" component_1="sodium_potassium_pump"/>
      <map_variables variable_2="i_NaK" variable_1="i_NaK"/>
   </connection>
   <connection>
      <map_components component_2="intracellular_potassium_concentration" component_1="transient_outward_current"/>
      <map_variables variable_2="i_to" variable_1="i_to"/>
   </connection>
   <connection>
      <map_components component_2="intracellular_calcium_concentration" component_1="sodium_calcium_exchanger"/>
      <map_variables variable_2="i_NaCa_cyt" variable_1="i_NaCa_cyt"/>
      <map_variables variable_2="i_NaCa_ds" variable_1="i_NaCa_ds"/>
      <map_variables variable_2="Ca_i" variable_1="Ca_i"/>
      <map_variables variable_2="Ca_ds" variable_1="Ca_ds"/>
   </connection>
   <connection>
      <map_components component_2="intracellular_sodium_concentration" component_1="sodium_calcium_exchanger"/>
      <map_variables variable_2="i_NaCa_cyt" variable_1="i_NaCa_cyt"/>
      <map_variables variable_2="Na_i" variable_1="Na_i"/>
   </connection>
   <connection>
      <map_components component_2="extracellular_sodium_concentration" component_1="sodium_calcium_exchanger"/>
      <map_variables variable_2="Na_o" variable_1="Na_o"/>
   </connection>
   <connection>
      <map_components component_2="extracellular_calcium_concentration" component_1="sodium_calcium_exchanger"/>
      <map_variables variable_2="Ca_o" variable_1="Ca_o"/>
   </connection>
   <connection>
      <map_components component_2="intracellular_calcium_concentration" component_1="sarcoplasmic_reticulum_calcium_pump"/>
      <map_variables variable_2="i_up" variable_1="i_up"/>
      <map_variables variable_2="Ca_i" variable_1="Ca_i"/>
      <map_variables variable_2="Ca_up" variable_1="Ca_up"/>
   </connection>
   <connection>
      <map_components component_2="intracellular_calcium_concentration" component_1="calcium_translocation"/>
      <map_variables variable_2="i_trans" variable_1="i_trans"/>
      <map_variables variable_2="Ca_rel" variable_1="Ca_rel"/>
      <map_variables variable_2="Ca_up" variable_1="Ca_up"/>
   </connection>
   <connection>
      <map_components component_2="intracellular_calcium_concentration" component_1="calcium_release"/>
      <map_variables variable_2="i_rel" variable_1="i_rel"/>
      <map_variables variable_2="Ca_rel" variable_1="Ca_rel"/>
      <map_variables variable_2="Ca_i" variable_1="Ca_i"/>
      <map_variables variable_2="Ca_ds" variable_1="Ca_ds"/>
   </connection>
   <connection>
      <map_components component_2="contraction" component_1="intracellular_calcium_concentration"/>
      <map_variables variable_2="Ca_Calmod" variable_1="Ca_Calmod"/>
      <map_variables variable_2="Ca_Trop" variable_1="Ca_Trop"/>
      <map_variables variable_2="Calmod" variable_1="Calmod"/>
      <map_variables variable_2="Trop" variable_1="Trop"/>
      <map_variables variable_2="isometric_tension" variable_1="isometric_tension"/>
      <map_variables variable_2="sarcomere_length" variable_1="sarcomere_length"/>
   </connection>
   <connection>
      <map_components component_2="contraction" component_1="calcium_release"/>
      <map_variables variable_2="isometric_tension" variable_1="isometric_tension"/>
      <map_variables variable_2="sarcomere_length" variable_1="sarcomere_length"/>
   </connection>
   <connection>
      <map_components component_2="contraction" component_1="stretch_current"/>
      <map_variables variable_2="isometric_tension" variable_1="isometric_tension"/>
      <map_variables variable_2="sarcomere_length" variable_1="sarcomere_length"/>
   </connection>
   <connection>
      <map_components component_2="rapid_delayed_rectifier_potassium_current_xr1_gate" component_1="rapid_delayed_rectifier_potassium_current"/>
      <map_variables variable_2="xr1" variable_1="xr1"/>
      <map_variables variable_2="time" variable_1="time"/>
      <map_variables variable_2="V" variable_1="V"/>
   </connection>
   <connection>
      <map_components component_2="rapid_delayed_rectifier_potassium_current_xr2_gate" component_1="rapid_delayed_rectifier_potassium_current"/>
      <map_variables variable_2="xr2" variable_1="xr2"/>
      <map_variables variable_2="time" variable_1="time"/>
      <map_variables variable_2="V" variable_1="V"/>
   </connection>
   <connection>
      <map_components component_2="slow_delayed_rectifier_potassium_current_xs_gate" component_1="slow_delayed_rectifier_potassium_current"/>
      <map_variables variable_2="xs" variable_1="xs"/>
      <map_variables variable_2="time" variable_1="time"/>
      <map_variables variable_2="V" variable_1="V"/>
   </connection>
   <connection>
      <map_components component_2="ACh_dependent_potassium_current_xACh_gate" component_1="ACh_dependent_potassium_current"/>
      <map_variables variable_2="x_ACh" variable_1="x_ACh"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="fast_sodium_current_m_gate" component_1="fast_sodium_current"/>
      <map_variables variable_2="m" variable_1="m"/>
      <map_variables variable_2="time" variable_1="time"/>
      <map_variables variable_2="V" variable_1="V"/>
   </connection>
   <connection>
      <map_components component_2="fast_sodium_current_h_gate" component_1="fast_sodium_current"/>
      <map_variables variable_2="h" variable_1="h"/>
      <map_variables variable_2="time" variable_1="time"/>
      <map_variables variable_2="V" variable_1="V"/>
   </connection>
   <connection>
      <map_components component_2="L_type_Ca_channel_d_gate" component_1="L_type_Ca_channel"/>
      <map_variables variable_2="d" variable_1="d"/>
      <map_variables variable_2="time" variable_1="time"/>
      <map_variables variable_2="V" variable_1="V"/>
   </connection>
   <connection>
      <map_components component_2="L_type_Ca_channel_f_gate" component_1="L_type_Ca_channel"/>
      <map_variables variable_2="f" variable_1="f"/>
      <map_variables variable_2="time" variable_1="time"/>
      <map_variables variable_2="V" variable_1="V"/>
   </connection>
   <connection>
      <map_components component_2="L_type_Ca_channel_f2_gate" component_1="L_type_Ca_channel"/>
      <map_variables variable_2="f2" variable_1="f2"/>
      <map_variables variable_2="time" variable_1="time"/>
      <map_variables variable_2="V" variable_1="V"/>
      <map_variables variable_2="Ca_i" variable_1="Ca_i"/>
      <map_variables variable_2="Km_f2" variable_1="Km_f2"/>
   </connection>
   <connection>
      <map_components component_2="L_type_Ca_channel_f2ds_gate" component_1="L_type_Ca_channel"/>
      <map_variables variable_2="f2ds" variable_1="f2ds"/>
      <map_variables variable_2="time" variable_1="time"/>
      <map_variables variable_2="V" variable_1="V"/>
      <map_variables variable_2="Ca_ds" variable_1="Ca_ds"/>
      <map_variables variable_2="Km_f2ds" variable_1="Km_f2ds"/>
      <map_variables variable_2="R_decay" variable_1="R_decay"/>
   </connection>
   <connection>
      <map_components component_2="transient_outward_current_s_gate" component_1="transient_outward_current"/>
      <map_variables variable_2="s" variable_1="s"/>
      <map_variables variable_2="time" variable_1="time"/>
      <map_variables variable_2="V" variable_1="V"/>
   </connection>
   <connection>
      <map_components component_2="transient_outward_current_r_gate" component_1="transient_outward_current"/>
      <map_variables variable_2="r" variable_1="r"/>
      <map_variables variable_2="time" variable_1="time"/>
      <map_variables variable_2="V" variable_1="V"/>
   </connection>
   <connection>
      <map_components component_2="intracellular_calcium_concentration" component_1="intracellular_sodium_concentration"/>
      <map_variables variable_2="V_i" variable_1="V_i"/>
   </connection>
   <connection>
      <map_components component_2="intracellular_calcium_concentration" component_1="intracellular_potassium_concentration"/>
      <map_variables variable_2="V_i" variable_1="V_i"/>
   </connection>
   <connection>
      <map_components component_2="time_independent_potassium_current" component_1="ACh_dependent_potassium_current"/>
      <map_variables variable_2="K_mk1" variable_1="K_mk1"/>
   </connection>




<rdf:RDF>
  <rdf:Seq rdf:about="rdf:#0ce2c5c0-6262-4a24-917e-1fc336b0cabb">
    <rdf:li rdf:resource="rdf:#492a8d76-a4ab-4919-8506-fd55bb48c9e6"/>
    <rdf:li rdf:resource="rdf:#8bb3cc2e-3939-4d46-b384-50f8f57d72d5"/>
    <rdf:li rdf:resource="rdf:#8bda7cea-82ea-48a2-a242-f068440dd53b"/>
    <rdf:li rdf:resource="rdf:#1c28a87b-ecfc-4a17-9e63-1069ce345cd1"/>
  </rdf:Seq>
  <rdf:Description rdf:about="rdf:#72025205-c2ee-46ba-945e-6ad6e24afc85">
    <vCard:Given>D</vCard:Given>
    <vCard:Family>Noble</vCard:Family>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#$hm4db">
    <ns7:boundIntervals rdf:resource="rdf:#$im4db"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#db5ec880-ef58-4c15-8e0d-36bf8e7e497f">
    <vCard:Given>Penny</vCard:Given>
    <vCard:Family>Noble</vCard:Family>
    <vCard:Other>J</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#9c09dc78-7f49-4457-93ed-c7b9955a38c8">
    <vCard:Given>A</vCard:Given>
    <vCard:Family>Varghese</vCard:Family>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#a5f2a5b3-2edd-47c1-8582-1eab7e8f752b">
    <vCard:Given>James</vCard:Given>
    <vCard:Family>Lawson</vCard:Family>
    <vCard:Other>Richard</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#d7beb4e1-8910-4c28-b560-61c8d1d4a814">
    <dcterms:W3CDTF>2007-12-04T15:11:55+13:00</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#2ecb01cc-a373-40fa-9501-4aa93ccd5d8a">
    <dcterms:modified rdf:resource="rdf:#a96b2bdc-cbd5-4502-ac7c-203b452d742b"/>
    <rdf:value>Penny Noble was informed by a community member that INaCads was not included in the dCads/dt ODE and hence the system was out of balance. This has been fixed in version 07. An alteration was also made to the model to remove two differential expressions in the intracellular_calcium concentration component to allow PCEnv 0.2 to run the model.</rdf:value>
    <cmeta:modifier rdf:resource="rdf:#dec66365-7c06-4d35-8dfb-28989da96412"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#6b5feb3f-e3a7-44d9-ab1c-7874f76bb698">
    <vCard:FN>James Lawson</vCard:FN>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#f7a88ef0-b25b-4f1b-8c84-925bdd5b9e1c">
    <vCard:N rdf:resource="rdf:#16333cab-2375-4624-a3e1-16842ffa5ddc"/>
  </rdf:Description>
  <rdf:Description rdf:about="#noble_model_1998_stretch">
    <ns7:simulation rdf:resource="rdf:#$hm4db"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#8bda7cea-82ea-48a2-a242-f068440dd53b">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#7342d718-2ba4-41c7-8a1d-ff647241975f"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#9f557c82-27e9-4db6-93a8-a9bd42afe2a6">
    <vCard:ORG rdf:resource="rdf:#f663fe1d-2bd6-497b-94d0-eef64ccedde3"/>
    <vCard:EMAIL rdf:resource="rdf:#8c5aca23-4409-46cb-bbca-4ca7826d19e6"/>
    <vCard:N rdf:resource="rdf:#db5ec880-ef58-4c15-8e0d-36bf8e7e497f"/>
  </rdf:Description>
  <rdf:Description rdf:about="">
    <dc:publisher>Department of Physiology, Anatomy &amp; Genetics, University of Oxford</dc:publisher>
    <cmeta:comment rdf:resource="rdf:#1b1008fc-cda4-422b-a815-3a82b5e9c32a"/>
    <dcterms:created rdf:resource="rdf:#1f496d47-a63b-480e-941e-53eb648706ec"/>
    <dc:creator rdf:resource="rdf:#9f557c82-27e9-4db6-93a8-a9bd42afe2a6"/>
    <cmeta:modification rdf:resource="rdf:#6761e204-596b-44f6-b8cd-51afc006050e"/>
    <cmeta:modification rdf:resource="rdf:#2daea026-24aa-4913-8020-2620ff205c30"/>
    <cmeta:modification rdf:resource="rdf:#2ecb01cc-a373-40fa-9501-4aa93ccd5d8a"/>
    <cmeta:modification rdf:resource="rdf:#43217e3a-8f02-4e5e-bf24-4d9a005b8777"/>
    <cmeta:modification rdf:resource="rdf:#4fa87bf8-9c73-415d-ac3a-b97e8e3887d6"/>
    <cmeta:modification rdf:resource="rdf:#bb802211-58b6-4242-94b0-8eec067e1049"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#1f496d47-a63b-480e-941e-53eb648706ec">
    <dcterms:W3CDTF>2006-05-03T00:00:00+00:00</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#f663fe1d-2bd6-497b-94d0-eef64ccedde3">
    <vCard:Orgname>Oxford University</vCard:Orgname>
    <vCard:Orgunit>Department of Physiology, Anatomy &amp; Genetics</vCard:Orgunit>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#16333cab-2375-4624-a3e1-16842ffa5ddc">
    <vCard:Given>James</vCard:Given>
    <vCard:Family>Lawson</vCard:Family>
    <vCard:Other>Richard</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#1bbdcfe9-2109-4062-86de-7fc13439a418">
    <dcterms:W3CDTF>2007-07-18T12:59:52+12:00</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#6761e204-596b-44f6-b8cd-51afc006050e">
    <dcterms:modified rdf:resource="rdf:#1bca4043-64c7-44e2-97ec-092c2d8c44cb"/>
    <rdf:value>Note that the missing i_NCX_ds component was added to the dCai/dt rather than the dCads/dt equation in order to achieve balance in the model.</rdf:value>
    <cmeta:modifier rdf:resource="rdf:#e3a54851-287a-46f2-9a61-0f77ae0c1458"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#$im4db">
    <rdf:rest rdf:resource="http://www.w3.org/1999/02/22-rdf-syntax-ns#nil"/>
    <rdf:first rdf:resource="rdf:#$jm4db"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#4fa87bf8-9c73-415d-ac3a-b97e8e3887d6">
    <dcterms:modified rdf:resource="rdf:#86d71735-52ed-4a55-a0a2-1cc257a8217f"/>
    <rdf:value>This model has been curated by Penny Noble</rdf:value>
    <cmeta:modifier rdf:resource="rdf:#e6b753c8-3d08-4984-9f32-e577c2597f13"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#a96b2bdc-cbd5-4502-ac7c-203b452d742b">
    <dcterms:W3CDTF>2007-09-06T13:01:21+12:00</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#2daea026-24aa-4913-8020-2620ff205c30">
    <dcterms:modified rdf:resource="rdf:#d7beb4e1-8910-4c28-b560-61c8d1d4a814"/>
    <rdf:value>Units checked, curated.</rdf:value>
    <cmeta:modifier rdf:resource="rdf:#19aeec40-a962-4cbe-9ab2-65f9b1438afd"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#1b1008fc-cda4-422b-a815-3a82b5e9c32a">
    <dc:creator rdf:resource="rdf:#6b5feb3f-e3a7-44d9-ab1c-7874f76bb698"/>
    <rdf:value>This version was created by Penny Noble and is known to read into COR. James Lawson fixed an error (19/04/07) that prevented it reading in PCEnv (unsupported predefined operator diff in component intracellular_calcium_concentration. This variant of version 06 is the STRETCH model, which includes acetylcholine dependent modulation and the electromechanical processes simulated by Variant 02, and also includes simulation of stretch and fibroblasts. This model is parameterised the same as the BASIC model, with all additional variables set to zero.&#13;
To simulate multiple action potentials with this model the following integrator parameters must be used: max step size must be 0.001 or less, max point density must be 10,000 or more.</rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#43217e3a-8f02-4e5e-bf24-4d9a005b8777">
    <dcterms:modified rdf:resource="rdf:#1bbdcfe9-2109-4062-86de-7fc13439a418"/>
    <rdf:value>Changed stim_end in stimulus protocol from 10.1 to 1,000,000 seconds</rdf:value>
    <cmeta:modifier rdf:resource="rdf:#e2a1ec0e-f041-4b1c-af55-110096409042"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#8bb3cc2e-3939-4d46-b384-50f8f57d72d5">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#9c09dc78-7f49-4457-93ed-c7b9955a38c8"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#7f2d6d6b-344c-4a7d-a15b-a4747f3898ca">
    <dc:title>Canadian Journal of Cardiology</dc:title>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#cfbad6be-cc46-4275-a900-a78cb6fd8bd4">
    <dcterms:W3CDTF>2007-06-22T13:59:26+12:00</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#8c52a9a5-590f-4c2e-b319-a344696e33f1">
    <vCard:Given>Penny</vCard:Given>
    <vCard:Family>Noble</vCard:Family>
    <vCard:Other>J</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#$6wmB13">
    <ns7:nonstandard-pointDensity>100000</ns7:nonstandard-pointDensity>
    <ns7:maximumStepSize>0.001</ns7:maximumStepSize>
    <ns7:nonstandard-algorithm>bdf15</ns7:nonstandard-algorithm>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#e3a54851-287a-46f2-9a61-0f77ae0c1458">
    <vCard:N rdf:resource="rdf:#feb02716-3dd9-4ca6-8307-90e863c534d1"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#86d71735-52ed-4a55-a0a2-1cc257a8217f">
    <dcterms:W3CDTF>2007-02-06T13:00:00+13:00</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#dec66365-7c06-4d35-8dfb-28989da96412">
    <vCard:N rdf:resource="rdf:#b3346d0b-66a3-4213-a4e0-321b7948aa9d"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#bb802211-58b6-4242-94b0-8eec067e1049">
    <dcterms:modified rdf:resource="rdf:#cfbad6be-cc46-4275-a900-a78cb6fd8bd4"/>
    <rdf:value>Several variables have been given cmeta:id's to allow creation of a PCEnv session file.</rdf:value>
    <cmeta:modifier rdf:resource="rdf:#f7a88ef0-b25b-4f1b-8c84-925bdd5b9e1c"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#492a8d76-a4ab-4919-8506-fd55bb48c9e6">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#72025205-c2ee-46ba-945e-6ad6e24afc85"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#feb02716-3dd9-4ca6-8307-90e863c534d1">
    <vCard:Given>Penny</vCard:Given>
    <vCard:Family>Noble</vCard:Family>
    <vCard:Other/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#34f1c61b-8e98-4ed5-a473-3ba951668467">
    <vCard:Given>P</vCard:Given>
    <vCard:Family>Noble</vCard:Family>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#1c28a87b-ecfc-4a17-9e63-1069ce345cd1">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#34f1c61b-8e98-4ed5-a473-3ba951668467"/>
  </rdf:Description>
  <rdf:Description rdf:about="#noble_varghese_kohl_noble_1998_stretch">
    <ns7:simulation rdf:resource="rdf:#$4wmB13"/>
    <bqs:reference rdf:resource="rdf:#89bb8cae-ce4c-427e-b18f-fe308d30dd9e"/>
<bqs:reference rdf:parseType="Resource">
  <dc:subject rdf:parseType="Resource">
    <bqs:subject_type>keyword</bqs:subject_type>
    <rdf:value>
      <rdf:Bag>
        <rdf:li>diadic space</rdf:li>
        <rdf:li>guinea-pig</rdf:li>
        <rdf:li>excitation-contraction coupling</rdf:li>
        <rdf:li>ventricular myocyte</rdf:li>
        <rdf:li>calcium dynamics</rdf:li>
        <rdf:li>electrophysiology</rdf:li>
        <rdf:li>cardiac</rdf:li>
      </rdf:Bag>
    </rdf:value>
  </dc:subject>
</bqs:reference>
    <cmeta:comment rdf:resource="rdf:#56512c4d-6eaf-4d72-9205-e876c9ea5e70"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#56512c4d-6eaf-4d72-9205-e876c9ea5e70">
    <dc:creator rdf:resource="rdf:#5ba8104a-bf61-4827-8ecf-ec222459ce5b"/>
    <rdf:value>This file is a CellML description of Noble, Carghese, Kohl and Noble's 1998 extension of their 1991 guinea pig ventricular cell model. This model incorporates a diadic space, IKr and IKs currents, and length and tension dependent processes. Variant 3 is parameterised for the 'stretch' model.</rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#e6b753c8-3d08-4984-9f32-e577c2597f13">
    <vCard:N rdf:resource="rdf:#8c52a9a5-590f-4c2e-b319-a344696e33f1"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#7342d718-2ba4-41c7-8a1d-ff647241975f">
    <vCard:Given>P</vCard:Given>
    <vCard:Family>Kohl</vCard:Family>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#89bb8cae-ce4c-427e-b18f-fe308d30dd9e">
    <bqs:Pubmed_id>9487284</bqs:Pubmed_id>
    <bqs:JournalArticle rdf:resource="rdf:#7412123a-3acd-4f20-aea4-11f746133c40"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#07a73155-cbd2-4f61-848e-2c11a2608cf0">
    <dcterms:W3CDTF>1998-01-00 00:00</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#8c5aca23-4409-46cb-bbca-4ca7826d19e6">
    <rdf:type rdf:resource="http://imc.org/vCard/3.0#internet"/>
    <rdf:value>penny.noble@dpag.ox.ac.uk</rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#$jm4db">
    <ns7:nonstandard-pointDensity>50000</ns7:nonstandard-pointDensity>
    <ns7:maximumStepSize>0.001</ns7:maximumStepSize>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#1bca4043-64c7-44e2-97ec-092c2d8c44cb">
    <dcterms:W3CDTF>2008-05-14T03:16:52+12:00</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#c261f321-6315-485c-b701-bff0ae199be0">
    <vCard:Given>Penny</vCard:Given>
    <vCard:Family>Noble</vCard:Family>
    <vCard:Other>J</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#e2a1ec0e-f041-4b1c-af55-110096409042">
    <vCard:N rdf:resource="rdf:#a5f2a5b3-2edd-47c1-8582-1eab7e8f752b"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#$4wmB13">
    <ns7:boundIntervals rdf:resource="rdf:#$5wmB13"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#5ba8104a-bf61-4827-8ecf-ec222459ce5b">
    <vCard:FN>James Lawson</vCard:FN>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#19aeec40-a962-4cbe-9ab2-65f9b1438afd">
    <vCard:N rdf:resource="rdf:#c261f321-6315-485c-b701-bff0ae199be0"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#7412123a-3acd-4f20-aea4-11f746133c40">
    <dc:creator rdf:resource="rdf:#0ce2c5c0-6262-4a24-917e-1fc336b0cabb"/>
    <dc:title>Improved guinea-pig ventricular cell model incorporating a diadic space, IKr and IKs, and length- and tension-dependent processes</dc:title>
    <bqs:volume>14(1)</bqs:volume>
    <bqs:first_page>123</bqs:first_page>
    <bqs:Journal rdf:resource="rdf:#7f2d6d6b-344c-4a7d-a15b-a4747f3898ca"/>
    <dcterms:issued rdf:resource="rdf:#07a73155-cbd2-4f61-848e-2c11a2608cf0"/>
    <bqs:last_page>134</bqs:last_page>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#$5wmB13">
    <rdf:rest rdf:resource="http://www.w3.org/1999/02/22-rdf-syntax-ns#nil"/>
    <rdf:first rdf:resource="rdf:#$6wmB13"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#b3346d0b-66a3-4213-a4e0-321b7948aa9d">
    <vCard:Given>James</vCard:Given>
    <vCard:Family>Lawson</vCard:Family>
    <vCard:Other>Richard</vCard:Other>
  </rdf:Description>
</rdf:RDF>
</model>