Location: Faber, Rudy, 2000 @ 124c7736bff0 / faber_rudy_2000.cellml

Author:
Catherine Lloyd <c.lloyd@auckland.ac.nz>
Date:
2009-10-07 13:34:08+13:00
Desc:
Added the slightly modified version of the Faber and Rudy as a subexposure of the model - not as the main exposure.
Permanent Source URI:
https://models.cellml.org/workspace/faber_rudy_2000/rawfile/124c7736bff0d65eb458994e1cd4b08a1e74f089/faber_rudy_2000.cellml

<?xml version='1.0' encoding='utf-8'?>
<!--
This CellML file was generated on 21/04/2008 at 14:12:36 using:

COR (0.9.31.915)
Copyright 2002-2008 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="faber_rudy_2000_version5" name="faber_rudy_2000_version5">
<documentation xmlns="http://cellml.org/tmp-documentation">
<article>
  <articleinfo>
  <title>Updated Luo-Rudy Mammalian Ventricular Model II</title>
  <author>
    <firstname>Penny</firstname>
          <surname>Noble</surname>
    <affiliation>
      <shortaffil>Oxford University</shortaffil>
    </affiliation>
  </author>
</articleinfo>
  <section id="sec_status">
    <title>Model Status</title>
    <para>
            This model has been curated and validated by Penny Noble of Oxford University and is known to run in PCEnv and COR.
          </para>
  </section>
  <sect1 id="sec_structure">
<title>Model Structure</title>

<para>A subsequent pair of papers by Ching-hsing Luo and Yoram Rudy (1994a; 1994b) further developed the Luo-Rudy I model by addressing some of the issues which were not investigated in their original model.  In particular, the LR-II model incorporates a more thorough description of the processes which regulate intracellular calcium ion concentration and the movement of calcium ions through the cell and to and from the sarcoplasmic reticulum (see figure 1 below).
</para>

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

<para>
<ulink url="http://circres.ahajournals.org/cgi/content/abstract/74/6/1071">A Dynamic Model of the Cardiac Ventricular Action Potential - Simulations of Ionic Currents and Concentration Changes</ulink>, Ching-hsing Luo and Yoram Rudy, 1994, <ulink url="http://circres.ahajournals.org/">
            <emphasis>Circulation Research</emphasis>
          </ulink>, 74, 1071-1097.
</para>

<para>
Since the model was published in 1994, it has been updated several times, incorporating new experimental data taken from papers written by other authors.  These modifications have been considered in the model presented here.  To illustrate how the Luo-Rudy II dynamic model has been developed, see the sequence of diagrams.  The 1997 Shaw and Rudy model diagram is identical to the 1999 Viswanathan model diagram since the only difference between the two models is that the Viswanathan model introduces a second activation gate (X<subscript>s2</subscript>) into the formulation of I<subscript>Ks</subscript>.  In addition, the development of the Luo-Rudy II model has been documented by Case Western Reserve University (To view this website, follow this link: <ulink url="http://www.cwru.edu/med/CBRTC/LRdOnline/development.htm"> http://www.cwru.edu/med/CBRTC/LRdOnline/development.htm</ulink>  
</para>

<informalfigure float="0" id="fig_cell_diagram_LR_II">
<mediaobject>
  <imageobject>
    <objectinfo>
      <title>cell diagram of the LR-II model showing ionic currents, pumps and exchangers within the sarcolemma and the sarcoplasmic reticulum</title>
    </objectinfo>
    <imagedata fileref="luo_1994.png"/>
  </imageobject>
</mediaobject>
<caption>A schematic diagram describing the current flows across the cell membrane that are captured in the original, 1994 Luo-Rudy II, dynamic model.</caption>
</informalfigure>

<informalfigure float="0" id="fig_cell_diagram_zeng">
<mediaobject>
  <imageobject>
    <objectinfo>
      <title>cell diagram of the zeng model showing ionic currents, pumps and exchangers within the sarcolemma and the sarcoplasmic reticulum</title>
    </objectinfo>
    <imagedata fileref="zeng_1995.png"/>
  </imageobject>
</mediaobject>
<caption>A schematic diagram showing the updated Luo-Rudy II, dynamic model.  In their 1995 ventricular model Zeng <emphasis>et al</emphasis> have included the slow and fast components of the delayed rectifier K<superscript>+</superscript> current, I<subscript>Kr</subscript> and I<subscript>Ks</subscript>.  A T-type Calcium channel (I<subscript>Ca(T)</subscript>) is also added.</caption>
</informalfigure>

<informalfigure float="0" id="fig_cell_diagram_shaw">
<mediaobject>
  <imageobject>
    <objectinfo>
      <title>cell diagram of the shaw and viswanathan models showing ionic currents, pumps and exchangers within the sarcolemma and the sarcoplasmic reticulum</title>
    </objectinfo>
    <imagedata fileref="shaw_1997.png"/>
  </imageobject>
</mediaobject>
<caption>
In their 1997 model Shaw and Rudy have built upon the 1995 Zeng <emphasis>et al</emphasis> version of the 1994 Luo-Rudy Dynamic Ventricular Model.  In addition to the slow and fast components of the delayed rectifier K<superscript>+</superscript> current, (I<subscript>Kr</subscript> and I<subscript>Ks</subscript>) and the  T-type Calcium channel (I<subscript>Ca(T)</subscript>), Shaw and Rudy have included an ATP-dependent potassium current (I<subscript>K(ATP)</subscript>) in their model.  In turn, Viswanathan <emphasis>et al</emphasis> develop this model to include two activation gates (X<subscript>s1</subscript> and X<subscript>s2</subscript>) into the formulation of I<subscript>Ks</subscript>.</caption>
</informalfigure>

<para>
All the original paper references are cited below:
</para>

<para>
<ulink url="http://circres.ahajournals.org/cgi/content/abstract/77/1/140">Two Components of the Delayed Rectifier K<superscript>+</superscript> Current in Ventricular Myocytes of the Guinea Pig Type</ulink>, Jinglin Zeng, Kenneth R. Laurita, David S. Rosenbaum and Yoram Rudy, 1995, <ulink url="http://circres.ahajournals.org/">
            <emphasis>Circulation Research</emphasis>
          </ulink>, 77, 140-152.  (The <ulink url="http://circres.ahajournals.org/cgi/content/full/77/1/140">full text</ulink> version of the article is available for Journal Members on the Circulation Research website.)  <ulink url="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?cmd=Retrieve&amp;db=PubMed&amp;list_uids=7788872&amp;dopt=Abstract">PubMed ID: 7788872</ulink>  
</para>

<para>
<ulink url="http://www.sciencedirect.com/science?_ob=ArticleURL&amp;_udi=B6T14-3R37WSS-17&amp;_user=140507&amp;_coverDate=08%2F31%2F1997&amp;_alid=185463716&amp;_rdoc=1&amp;_fmt=summary&amp;_orig=search&amp;_qd=1&amp;_cdi=4880&amp;_sort=d&amp;_docanchor=&amp;view=c&amp;_acct=C000011498&amp;_version=1&amp;_urlVersion=0&amp;_userid=140507&amp;md5=ad03ef8fddd6cc6d11f8484e421991f2">Electrophysiological effects of acute myocardial ischemia: a theoretical study of altered cell excitability and action potential duration</ulink>, Robin M. Shaw and Yoram Rudy, 1997, <ulink url="http://www.sciencedirect.com/science?_ob=JournalURL&amp;_cdi=4880&amp;_auth=y&amp;_acct=C000011498&amp;_version=1&amp;_urlVersion=0&amp;_userid=140507&amp;md5=845774525765eaf31dec5d2fc2a614f6">
            <emphasis>Cardiovascular Research</emphasis>
          </ulink>, 35, 256-272.  <ulink url="http://www.sciencedirect.com/science?_ob=ArticleURL&amp;_udi=B6T14-3R37WSS-17&amp;_coverDate=08%2F31%2F1997&amp;_alid=185463716&amp;_rdoc=1&amp;_fmt=&amp;_orig=search&amp;_qd=1&amp;_cdi=4880&amp;_sort=d&amp;view=c&amp;_acct=C000011498&amp;_version=1&amp;_urlVersion=0&amp;_userid=140507&amp;md5=2b061b7dec7050a8c72e2eba382d5e97">Full text</ulink> and PDF versions of the article are available to subscribers on the Science Direct website.  <ulink url="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?cmd=Retrieve&amp;db=PubMed&amp;list_uids=9349389&amp;dopt=Abstract">PubMed ID: 9349389</ulink>    
</para>

<para>
<ulink url="http://circ.ahajournals.org/cgi/content/abstract/99/18/2466">Effects of I<subscript>Kr</subscript> and I<subscript>Ks</subscript> Heterogeneity on Action Potential Duration and Its Rate Dependence: A Simulation Study</ulink>, Prakash C. Viswanathan, BE; Robin M. Shaw, PhD; and Yoram Rudy, PhD, 1999, <ulink url="http://circ.ahajournals.org/">
            <emphasis>Circulation</emphasis>
          </ulink>, 99, 2466-2474.  (<ulink url="http://www.circ.ahajournals.org/cgi/content/full/99/18/2466">Full text</ulink> and <ulink url="http://www.circ.ahajournals.org/cgi/reprint/99/18/2466.pdf">PDF</ulink> versions of the article are available for Journal Members on the Circulation website.)  <ulink url="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?cmd=Retrieve&amp;db=PubMed&amp;list_uids=10318671&amp;dopt=Abstract">PubMed ID: 10318671</ulink>      
</para>

<para>
<ulink url="http://www.biophysj.org/cgi/content/abstract/78/5/2392">Action Potential and Contractility in [Na<superscript>+</superscript>]<subscript>i</subscript> Overloaded Cardiac Myocytes</ulink>, Gregory M. Faber and Yoram Rudy, 2000, <ulink url="http://www.biophysj.org/">
            <emphasis>Biophysical Journal</emphasis>
          </ulink>, 78, 2392-2404.  (<ulink url="http://www.biophysj.org/cgi/content/full/78/5/2392">Full text</ulink> and <ulink url="http://www.biophysj.org/cgi/reprint/78/5/2392.pdf">PDF</ulink> versions of the article are available for Journal Members on the Biophysical Journal website.)  <ulink url="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?cmd=Retrieve&amp;db=PubMed&amp;list_uids=10777735&amp;dopt=Abstract">PubMed ID: 10777735</ulink>   
</para>

<informalfigure float="0" id="fig_cell_diagram_faber">
<mediaobject>
  <imageobject>
    <objectinfo>
      <title>cell diagram of the shaw and viswanathan models showing ionic currents, pumps and exchangers within the sarcolemma and the sarcoplasmic reticulum</title>
    </objectinfo>
    <imagedata fileref="faber_2000.png"/>
  </imageobject>
</mediaobject>
<caption>In this most recent version of the 1994 Luo-Rudy Dynamic Ventricular Model, Faber and Rudy have built upon the 1999 Viswanathan model and they have added a sodium-dependent potassium current (I<subscript>K(Na)</subscript>).</caption>
</informalfigure>

<informalfigure float="0" id="fig_cellml_rendering">
<mediaobject>
  <imageobject>
    <objectinfo>
      <title>the cellml rendering of the Luo-Rudy II model</title>
    </objectinfo>
    <imagedata fileref="cellml_rendering.gif"/>
  </imageobject>
</mediaobject>
<caption>The network defined in the CellML description of the most recent version of the Luo-Rudy II model. A key describing the significance of the shapes of the components and the colours of the connections between them is in the <ulink url="${HTML_EXMPL_GRAPHICAL_NOTATION}">notation guide</ulink>.  For simplicity, not all the variables are shown.</caption>
</informalfigure>


</sect1>
</article>
</documentation>
   <units name="ms">
      <unit units="second" prefix="milli"/>
   </units>
   <units name="per_ms">
      <unit units="ms" 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="ms_per_millivolt">
      <unit units="ms"/>
      <unit units="volt" prefix="milli" exponent="-1"/>
   </units>
   <units name="per_millivolt_ms">
      <unit units="millivolt" exponent="-1"/>
      <unit units="ms" exponent="-1"/>
   </units>
   <units name="millivolt_per_ms">
      <unit units="volt" prefix="milli"/>
      <unit units="ms" exponent="-1"/>
   </units>
   <units name="milliS_per_microF">
      <unit units="siemens" prefix="milli"/>
      <unit units="farad" prefix="micro" exponent="-1"/>
   </units>
   <units name="cm2">
      <unit units="metre" prefix="centi" exponent="2"/>
   </units>
   <units name="milliS_per_cm2">
      <unit units="siemens" prefix="milli"/>
      <unit units="metre" prefix="centi" exponent="-2"/>
   </units>
   <units name="nanoS_per_cm2">
      <unit units="siemens" prefix="nano"/>
      <unit units="metre" prefix="centi" exponent="-2"/>
   </units>
   <units name="microF_per_cm2">
      <unit units="microF"/>
      <unit units="metre" prefix="centi" exponent="-2"/>
   </units>
   <units name="microF">
      <unit units="farad" prefix="micro"/>
   </units>
   <units name="microA_per_microF">
      <unit units="ampere" prefix="micro"/>
      <unit units="farad" prefix="micro" exponent="-1"/>
   </units>
   <units name="millimolar_per_ms">
      <unit units="mole" prefix="milli"/>
      <unit units="litre" exponent="-1"/>
      <unit units="ms" 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="micromolar">
      <unit units="mole" prefix="micro"/>
      <unit units="litre" 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="cm_per_second">
      <unit units="metre" prefix="centi"/>
      <unit units="second" exponent="-1"/>
   </units>
   <units name="mm">
      <unit units="metre" prefix="milli"/>
   </units>
   <units name="mm2">
      <unit units="metre" prefix="milli" exponent="2"/>
   </units>
   <units name="micro_litre">
      <unit units="litre" prefix="micro"/>
   </units>
   <units name="litre_per_farad_millisecond">
      <unit units="litre"/>
      <unit units="farad" exponent="-1"/>
      <unit units="ms" exponent="-1"/>
   </units>
   <component name="environment">
      <variable units="ms" public_interface="out" name="time"/>
   </component>
   <component name="cell">
      <variable units="millivolt" public_interface="out" name="V" initial_value="-85.126196430406"/>
      <variable units="joule_per_kilomole_kelvin" public_interface="out" name="R" initial_value="8314"/>
      <variable units="kelvin" public_interface="out" name="T" initial_value="310"/>
      <variable units="coulomb_per_mole" public_interface="out" name="F" initial_value="96485"/>
      <variable units="microA_per_microF" public_interface="out" name="I_st"/>
      <variable units="ms" public_interface="in" name="time"/>
      <variable units="microA_per_microF" public_interface="in" name="i_Na"/>
      <variable units="microA_per_microF" public_interface="in" name="i_Ca_L"/>
      <variable units="microA_per_microF" public_interface="in" name="i_Ca_T"/>
      <variable units="microA_per_microF" public_interface="in" name="i_Kr"/>
      <variable units="microA_per_microF" public_interface="in" name="i_Ks"/>
      <variable units="microA_per_microF" public_interface="in" name="i_K_Na"/>
      <variable units="microA_per_microF" public_interface="in" name="i_K_ATP"/>
      <variable units="microA_per_microF" public_interface="in" name="i_to"/>
      <variable units="microA_per_microF" public_interface="in" name="i_NaCa"/>
      <variable units="microA_per_microF" public_interface="in" name="i_K1"/>
      <variable units="microA_per_microF" public_interface="in" name="i_Kp"/>
      <variable units="microA_per_microF" public_interface="in" name="i_p_Ca"/>
      <variable units="microA_per_microF" public_interface="in" name="i_Na_b"/>
      <variable units="microA_per_microF" public_interface="in" name="i_Ca_b"/>
      <variable units="microA_per_microF" public_interface="in" name="i_NaK"/>
      <variable units="microA_per_microF" public_interface="in" name="i_ns_Ca"/>
      <variable units="microA_per_microF" public_interface="out" name="dVdt"/>
      <variable units="ms" name="stim_start" initial_value="100"/>
      <variable units="ms" name="stim_end" initial_value="9e6"/>
      <variable units="ms" name="stim_period" initial_value="1000"/>
      <variable units="ms" name="stim_duration" initial_value="2"/>
      <variable units="microA_per_microF" name="stim_amplitude" initial_value="-25.5"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>I_st</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="microA_per_microF">0</cn>
               </otherwise>
            </piecewise>
         </apply>
         <apply>
            <eq/>
            <ci>dVdt</ci>
            <apply>
               <minus/>
               <apply>
                  <plus/>
                  <ci>i_Na</ci>
                  <ci>i_Ca_L</ci>
                  <ci>i_Ca_T</ci>
                  <ci>i_Kr</ci>
                  <ci>i_Ks</ci>
                  <ci>i_K_Na</ci>
                  <ci>i_K_ATP</ci>
                  <ci>i_to</ci>
                  <ci>i_K1</ci>
                  <ci>i_Kp</ci>
                  <ci>i_NaCa</ci>
                  <ci>i_p_Ca</ci>
                  <ci>i_Na_b</ci>
                  <ci>i_Ca_b</ci>
                  <ci>i_NaK</ci>
                  <ci>i_ns_Ca</ci>
                  <ci>I_st</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>V</ci>
            </apply>
            <ci>dVdt</ci>
         </apply>
      </math>
   </component>
   <component name="fast_sodium_current">
      <variable units="microA_per_microF" public_interface="out" name="i_Na"/>
      <variable units="millivolt" public_interface="out" name="E_Na"/>
      <variable units="milliS_per_microF" name="g_Na" initial_value="16"/>
      <variable units="ms" public_interface="in" private_interface="out" name="time"/>
      <variable units="millivolt" public_interface="in" private_interface="out" name="V"/>
      <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="millimolar" public_interface="in" name="Nai"/>
      <variable units="millimolar" public_interface="in" name="Nao"/>
      <variable units="dimensionless" private_interface="in" name="m"/>
      <variable units="dimensionless" private_interface="in" name="h"/>
      <variable units="dimensionless" private_interface="in" name="j"/>
      <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>Nao</ci>
                     <ci>Nai</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <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>
               <ci>j</ci>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>E_Na</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="fast_sodium_current_m_gate">
      <variable units="dimensionless" public_interface="out" name="m" initial_value="0.001513296212"/>
      <variable units="per_ms" name="alpha_m"/>
      <variable units="per_ms" name="beta_m"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="millivolt" name="delta_m" initial_value="1e-5"/>
      <variable units="millivolt" name="E0_m"/>
      <variable units="ms" public_interface="in" name="time"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>E0_m</ci>
            <apply>
               <plus/>
               <ci>V</ci>
               <cn cellml:units="millivolt">47.13</cn>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>alpha_m</ci>
            <piecewise>
               <piece>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <cn cellml:units="per_millivolt_ms">0.32</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>
                  <apply>
                     <geq/>
                     <apply>
                        <abs/>
                        <ci>E0_m</ci>
                     </apply>
                     <ci>delta_m</ci>
                  </apply>
               </piece>
               <otherwise>
                  <cn cellml:units="per_ms">3.2</cn>
               </otherwise>
            </piecewise>
         </apply>
         <apply>
            <eq/>
            <ci>beta_m</ci>
            <apply>
               <times/>
               <cn cellml:units="per_ms">0.08</cn>
               <apply>
                  <exp/>
                  <apply>
                     <divide/>
                     <apply>
                        <minus/>
                        <ci>V</ci>
                     </apply>
                     <cn cellml:units="millivolt">11</cn>
                  </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.985318710339"/>
      <variable units="per_ms" name="alpha_h"/>
      <variable units="per_ms" name="beta_h"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="ms" public_interface="in" name="time"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>alpha_h</ci>
            <piecewise>
               <piece>
                  <apply>
                     <times/>
                     <cn cellml:units="per_ms">0.135</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <plus/>
                              <cn cellml:units="millivolt">80</cn>
                              <ci>V</ci>
                           </apply>
                           <apply>
                              <minus/>
                              <cn cellml:units="millivolt">6.8</cn>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
                  <apply>
                     <lt/>
                     <ci>V</ci>
                     <apply>
                        <minus/>
                        <cn cellml:units="millivolt">40</cn>
                     </apply>
                  </apply>
               </piece>
               <otherwise>
                  <cn cellml:units="per_ms">0</cn>
               </otherwise>
            </piecewise>
         </apply>
         <apply>
            <eq/>
            <ci>beta_h</ci>
            <piecewise>
               <piece>
                  <apply>
                     <plus/>
                     <apply>
                        <times/>
                        <cn cellml:units="per_ms">3.56</cn>
                        <apply>
                           <exp/>
                           <apply>
                              <times/>
                              <cn cellml:units="per_millivolt">0.079</cn>
                              <ci>V</ci>
                           </apply>
                        </apply>
                     </apply>
                     <apply>
                        <times/>
                        <cn cellml:units="per_ms">310000</cn>
                        <apply>
                           <exp/>
                           <apply>
                              <times/>
                              <cn cellml:units="per_millivolt">0.35</cn>
                              <ci>V</ci>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
                  <apply>
                     <lt/>
                     <ci>V</ci>
                     <apply>
                        <minus/>
                        <cn cellml:units="millivolt">40</cn>
                     </apply>
                  </apply>
               </piece>
               <otherwise>
                  <apply>
                     <divide/>
                     <cn cellml:units="per_ms">1</cn>
                     <apply>
                        <times/>
                        <cn cellml:units="dimensionless">0.13</cn>
                        <apply>
                           <plus/>
                           <cn cellml:units="dimensionless">1</cn>
                           <apply>
                              <exp/>
                              <apply>
                                 <divide/>
                                 <apply>
                                    <plus/>
                                    <ci>V</ci>
                                    <cn cellml:units="millivolt">10.66</cn>
                                 </apply>
                                 <apply>
                                    <minus/>
                                    <cn cellml:units="millivolt">11.1</cn>
                                 </apply>
                              </apply>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </otherwise>
            </piecewise>
         </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="fast_sodium_current_j_gate">
      <variable units="dimensionless" public_interface="out" name="j" initial_value="0.990731711845"/>
      <variable units="per_ms" name="alpha_j"/>
      <variable units="per_ms" name="beta_j"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="ms" public_interface="in" name="time"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>alpha_j</ci>
            <piecewise>
               <piece>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <apply>
                           <minus/>
                           <apply>
                              <plus/>
                              <apply>
                                 <times/>
                                 <cn cellml:units="per_millivolt_ms">127140</cn>
                                 <apply>
                                    <exp/>
                                    <apply>
                                       <times/>
                                       <cn cellml:units="per_millivolt">0.2444</cn>
                                       <ci>V</ci>
                                    </apply>
                                 </apply>
                              </apply>
                              <apply>
                                 <times/>
                                 <cn cellml:units="per_millivolt_ms">3.474e-5</cn>
                                 <apply>
                                    <exp/>
                                    <apply>
                                       <times/>
                                       <apply>
                                          <minus/>
                                          <cn cellml:units="per_millivolt">0.04391</cn>
                                       </apply>
                                       <ci>V</ci>
                                    </apply>
                                 </apply>
                              </apply>
                           </apply>
                        </apply>
                        <apply>
                           <plus/>
                           <ci>V</ci>
                           <cn cellml:units="millivolt">37.78</cn>
                        </apply>
                     </apply>
                     <apply>
                        <plus/>
                        <cn cellml:units="dimensionless">1</cn>
                        <apply>
                           <exp/>
                           <apply>
                              <times/>
                              <cn cellml:units="per_millivolt">0.311</cn>
                              <apply>
                                 <plus/>
                                 <ci>V</ci>
                                 <cn cellml:units="millivolt">79.23</cn>
                              </apply>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
                  <apply>
                     <lt/>
                     <ci>V</ci>
                     <apply>
                        <minus/>
                        <cn cellml:units="millivolt">40</cn>
                     </apply>
                  </apply>
               </piece>
               <otherwise>
                  <cn cellml:units="per_ms">0</cn>
               </otherwise>
            </piecewise>
         </apply>
         <apply>
            <eq/>
            <ci>beta_j</ci>
            <piecewise>
               <piece>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <cn cellml:units="per_ms">0.1212</cn>
                        <apply>
                           <exp/>
                           <apply>
                              <times/>
                              <apply>
                                 <minus/>
                                 <cn cellml:units="per_millivolt">0.01052</cn>
                              </apply>
                              <ci>V</ci>
                           </apply>
                        </apply>
                     </apply>
                     <apply>
                        <plus/>
                        <cn cellml:units="dimensionless">1</cn>
                        <apply>
                           <exp/>
                           <apply>
                              <times/>
                              <apply>
                                 <minus/>
                                 <cn cellml:units="per_millivolt">0.1378</cn>
                              </apply>
                              <apply>
                                 <plus/>
                                 <ci>V</ci>
                                 <cn cellml:units="millivolt">40.14</cn>
                              </apply>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
                  <apply>
                     <lt/>
                     <ci>V</ci>
                     <apply>
                        <minus/>
                        <cn cellml:units="millivolt">40</cn>
                     </apply>
                  </apply>
               </piece>
               <otherwise>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <cn cellml:units="per_ms">0.3</cn>
                        <apply>
                           <exp/>
                           <apply>
                              <times/>
                              <apply>
                                 <minus/>
                                 <cn cellml:units="per_millivolt">2.535e-7</cn>
                              </apply>
                              <ci>V</ci>
                           </apply>
                        </apply>
                     </apply>
                     <apply>
                        <plus/>
                        <cn cellml:units="dimensionless">1</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">32</cn>
                              </apply>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </otherwise>
            </piecewise>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>j</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>alpha_j</ci>
                  <apply>
                     <minus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <ci>j</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>beta_j</ci>
                  <ci>j</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="L_type_Ca_channel">
      <variable units="microA_per_microF" public_interface="out" name="i_Ca_L"/>
      <variable units="microA_per_microF" public_interface="out" name="i_CaCa"/>
      <variable units="microA_per_microF" public_interface="out" name="i_CaK"/>
      <variable units="microA_per_microF" public_interface="out" name="i_CaNa"/>
      <variable units="dimensionless" public_interface="out" name="gamma_Nai" initial_value="0.75"/>
      <variable units="dimensionless" public_interface="out" name="gamma_Nao" initial_value="0.75"/>
      <variable units="dimensionless" public_interface="out" name="gamma_Ki" initial_value="0.75"/>
      <variable units="dimensionless" public_interface="out" name="gamma_Ko" initial_value="0.75"/>
      <variable units="dimensionless" name="gamma_Cai" initial_value="1"/>
      <variable units="dimensionless" name="gamma_Cao" initial_value="0.341"/>
      <variable units="microA_per_microF" name="I_CaCa"/>
      <variable units="microA_per_microF" name="I_CaK"/>
      <variable units="microA_per_microF" name="I_CaNa"/>
      <variable units="litre_per_farad_millisecond" name="P_Ca" initial_value="0.00054"/>
      <variable units="litre_per_farad_millisecond" name="P_Na" initial_value="6.75e-7"/>
      <variable units="litre_per_farad_millisecond" name="P_K" initial_value="1.93e-7"/>
      <variable units="ms" public_interface="in" private_interface="out" name="time"/>
      <variable units="millivolt" public_interface="in" private_interface="out" name="V"/>
      <variable units="millimolar" public_interface="in" private_interface="out" name="Cai"/>
      <variable units="joule_per_kilomole_kelvin" public_interface="in" name="R"/>
      <variable units="kelvin" public_interface="in" name="T"/>
      <variable units="coulomb_per_mole" public_interface="in" name="F"/>
      <variable units="millimolar" public_interface="in" name="Cao"/>
      <variable units="millimolar" public_interface="in" name="Nao"/>
      <variable units="millimolar" public_interface="in" name="Ko"/>
      <variable units="millimolar" public_interface="in" name="Nai"/>
      <variable units="millimolar" public_interface="in" name="Ki"/>
      <variable units="dimensionless" private_interface="in" name="d"/>
      <variable units="dimensionless" private_interface="in" name="f"/>
      <variable units="dimensionless" private_interface="in" name="f_Ca"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>I_CaCa</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <ci>P_Ca</ci>
                        <apply>
                           <power/>
                           <cn cellml:units="dimensionless">2</cn>
                           <cn cellml:units="dimensionless">2</cn>
                        </apply>
                        <ci>V</ci>
                        <apply>
                           <power/>
                           <ci>F</ci>
                           <cn cellml:units="dimensionless">2</cn>
                        </apply>
                     </apply>
                     <apply>
                        <times/>
                        <ci>R</ci>
                        <ci>T</ci>
                     </apply>
                  </apply>
                  <apply>
                     <minus/>
                     <apply>
                        <times/>
                        <ci>gamma_Cai</ci>
                        <ci>Cai</ci>
                        <apply>
                           <exp/>
                           <apply>
                              <divide/>
                              <apply>
                                 <times/>
                                 <cn cellml:units="dimensionless">2</cn>
                                 <ci>V</ci>
                                 <ci>F</ci>
                              </apply>
                              <apply>
                                 <times/>
                                 <ci>R</ci>
                                 <ci>T</ci>
                              </apply>
                           </apply>
                        </apply>
                     </apply>
                     <apply>
                        <times/>
                        <ci>gamma_Cao</ci>
                        <ci>Cao</ci>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <minus/>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <times/>
                           <cn cellml:units="dimensionless">2</cn>
                           <ci>V</ci>
                           <ci>F</ci>
                        </apply>
                        <apply>
                           <times/>
                           <ci>R</ci>
                           <ci>T</ci>
                        </apply>
                     </apply>
                  </apply>
                  <cn cellml:units="dimensionless">1</cn>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>I_CaNa</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <ci>P_Na</ci>
                        <apply>
                           <power/>
                           <cn cellml:units="dimensionless">1</cn>
                           <cn cellml:units="dimensionless">2</cn>
                        </apply>
                        <ci>V</ci>
                        <apply>
                           <power/>
                           <ci>F</ci>
                           <cn cellml:units="dimensionless">2</cn>
                        </apply>
                     </apply>
                     <apply>
                        <times/>
                        <ci>R</ci>
                        <ci>T</ci>
                     </apply>
                  </apply>
                  <apply>
                     <minus/>
                     <apply>
                        <times/>
                        <ci>gamma_Nai</ci>
                        <ci>Nai</ci>
                        <apply>
                           <exp/>
                           <apply>
                              <divide/>
                              <apply>
                                 <times/>
                                 <cn cellml:units="dimensionless">1</cn>
                                 <ci>V</ci>
                                 <ci>F</ci>
                              </apply>
                              <apply>
                                 <times/>
                                 <ci>R</ci>
                                 <ci>T</ci>
                              </apply>
                           </apply>
                        </apply>
                     </apply>
                     <apply>
                        <times/>
                        <ci>gamma_Nao</ci>
                        <ci>Nao</ci>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <minus/>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <times/>
                           <cn cellml:units="dimensionless">1</cn>
                           <ci>V</ci>
                           <ci>F</ci>
                        </apply>
                        <apply>
                           <times/>
                           <ci>R</ci>
                           <ci>T</ci>
                        </apply>
                     </apply>
                  </apply>
                  <cn cellml:units="dimensionless">1</cn>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>I_CaK</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <ci>P_K</ci>
                        <apply>
                           <power/>
                           <cn cellml:units="dimensionless">1</cn>
                           <cn cellml:units="dimensionless">2</cn>
                        </apply>
                        <ci>V</ci>
                        <apply>
                           <power/>
                           <ci>F</ci>
                           <cn cellml:units="dimensionless">2</cn>
                        </apply>
                     </apply>
                     <apply>
                        <times/>
                        <ci>R</ci>
                        <ci>T</ci>
                     </apply>
                  </apply>
                  <apply>
                     <minus/>
                     <apply>
                        <times/>
                        <ci>gamma_Ki</ci>
                        <ci>Ki</ci>
                        <apply>
                           <exp/>
                           <apply>
                              <divide/>
                              <apply>
                                 <times/>
                                 <cn cellml:units="dimensionless">1</cn>
                                 <ci>V</ci>
                                 <ci>F</ci>
                              </apply>
                              <apply>
                                 <times/>
                                 <ci>R</ci>
                                 <ci>T</ci>
                              </apply>
                           </apply>
                        </apply>
                     </apply>
                     <apply>
                        <times/>
                        <ci>gamma_Ko</ci>
                        <ci>Ko</ci>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <minus/>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <times/>
                           <cn cellml:units="dimensionless">1</cn>
                           <ci>V</ci>
                           <ci>F</ci>
                        </apply>
                        <apply>
                           <times/>
                           <ci>R</ci>
                           <ci>T</ci>
                        </apply>
                     </apply>
                  </apply>
                  <cn cellml:units="dimensionless">1</cn>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>i_CaCa</ci>
            <apply>
               <times/>
               <ci>d</ci>
               <ci>f</ci>
               <ci>f_Ca</ci>
               <ci>I_CaCa</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>i_CaNa</ci>
            <apply>
               <times/>
               <ci>d</ci>
               <ci>f</ci>
               <ci>f_Ca</ci>
               <ci>I_CaNa</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>i_CaK</ci>
            <apply>
               <times/>
               <ci>d</ci>
               <ci>f</ci>
               <ci>f_Ca</ci>
               <ci>I_CaK</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>i_Ca_L</ci>
            <apply>
               <plus/>
               <ci>i_CaCa</ci>
               <ci>i_CaK</ci>
               <ci>i_CaNa</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="L_type_Ca_channel_d_gate">
      <variable units="dimensionless" public_interface="out" name="d" initial_value="0.000005906564"/>
      <variable units="per_ms" name="alpha_d"/>
      <variable units="per_ms" name="beta_d"/>
      <variable units="dimensionless" name="d_infinity"/>
      <variable units="ms" name="tau_d"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="ms" public_interface="in" name="time"/>
      <variable units="millivolt" name="E0_d"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>E0_d</ci>
            <apply>
               <plus/>
               <ci>V</ci>
               <cn cellml:units="millivolt">10</cn>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>d_infinity</ci>
            <apply>
               <divide/>
               <cn cellml:units="dimensionless">1</cn>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <minus/>
                           <ci>E0_d</ci>
                        </apply>
                        <cn cellml:units="millivolt">6.24</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>tau_d</ci>
            <piecewise>
               <piece>
                  <apply>
                     <divide/>
                     <cn cellml:units="ms">1</cn>
                     <apply>
                        <times/>
                        <cn cellml:units="dimensionless">0.035</cn>
                        <cn cellml:units="dimensionless">6.24</cn>
                        <cn cellml:units="dimensionless">2</cn>
                     </apply>
                  </apply>
                  <apply>
                     <lt/>
                     <apply>
                        <abs/>
                        <ci>E0_d</ci>
                     </apply>
                     <cn cellml:units="millivolt">1e-5</cn>
                  </apply>
               </piece>
               <otherwise>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <cn cellml:units="ms">1</cn>
                        <ci>d_infinity</ci>
                        <apply>
                           <minus/>
                           <cn cellml:units="dimensionless">1</cn>
                           <apply>
                              <exp/>
                              <apply>
                                 <divide/>
                                 <apply>
                                    <minus/>
                                    <ci>E0_d</ci>
                                 </apply>
                                 <cn cellml:units="millivolt">6.24</cn>
                              </apply>
                           </apply>
                        </apply>
                     </apply>
                     <apply>
                        <times/>
                        <cn cellml:units="per_millivolt">0.035</cn>
                        <ci>E0_d</ci>
                     </apply>
                  </apply>
               </otherwise>
            </piecewise>
         </apply>
         <apply>
            <eq/>
            <ci>alpha_d</ci>
            <apply>
               <divide/>
               <ci>d_infinity</ci>
               <ci>tau_d</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>beta_d</ci>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <ci>d_infinity</ci>
               </apply>
               <ci>tau_d</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>d</ci>
            </apply>
            <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>
      </math>
   </component>
   <component name="L_type_Ca_channel_f_gate">
      <variable units="dimensionless" public_interface="out" name="f" initial_value="0.999390880784"/>
      <variable units="per_ms" name="alpha_f"/>
      <variable units="per_ms" name="beta_f"/>
      <variable units="dimensionless" name="f_infinity"/>
      <variable units="ms" name="tau_f"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="ms" public_interface="in" name="time"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>f_infinity</ci>
            <apply>
               <plus/>
               <apply>
                  <divide/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <plus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <plus/>
                              <ci>V</ci>
                              <cn cellml:units="millivolt">35.06</cn>
                           </apply>
                           <cn cellml:units="millivolt">8.6</cn>
                        </apply>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <divide/>
                  <cn cellml:units="dimensionless">0.6</cn>
                  <apply>
                     <plus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <minus/>
                              <cn cellml:units="millivolt">50</cn>
                              <ci>V</ci>
                           </apply>
                           <cn cellml:units="millivolt">20</cn>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>tau_f</ci>
            <apply>
               <divide/>
               <cn cellml:units="ms">1</cn>
               <apply>
                  <plus/>
                  <apply>
                     <times/>
                     <cn cellml:units="dimensionless">0.0197</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <minus/>
                           <apply>
                              <power/>
                              <apply>
                                 <times/>
                                 <cn cellml:units="per_millivolt">0.0337</cn>
                                 <apply>
                                    <plus/>
                                    <ci>V</ci>
                                    <cn cellml:units="millivolt">10</cn>
                                 </apply>
                              </apply>
                              <cn cellml:units="dimensionless">2</cn>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
                  <cn cellml:units="dimensionless">0.02</cn>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>alpha_f</ci>
            <apply>
               <divide/>
               <ci>f_infinity</ci>
               <ci>tau_f</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>beta_f</ci>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <ci>f_infinity</ci>
               </apply>
               <ci>tau_f</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>f</ci>
            </apply>
            <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>
      </math>
   </component>
   <component name="L_type_Ca_channel_f_Ca_gate">
      <variable units="dimensionless" public_interface="out" name="f_Ca"/>
      <variable units="millimolar" name="Km_Ca" initial_value="0.0006"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="millimolar" public_interface="in" name="Cai"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>f_Ca</ci>
            <apply>
               <divide/>
               <cn cellml:units="dimensionless">1</cn>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <divide/>
                     <ci>Cai</ci>
                     <ci>Km_Ca</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="T_type_Ca_channel">
      <variable units="microA_per_microF" public_interface="out" name="i_Ca_T"/>
      <variable units="milliS_per_microF" name="g_CaT" initial_value="0.05"/>
      <variable units="millivolt" public_interface="in" name="E_Ca"/>
      <variable units="ms" 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="b"/>
      <variable units="dimensionless" private_interface="in" name="g"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_Ca_T</ci>
            <apply>
               <times/>
               <ci>g_CaT</ci>
               <ci>b</ci>
               <ci>b</ci>
               <ci>g</ci>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>E_Ca</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="T_type_Ca_channel_b_gate">
      <variable units="dimensionless" public_interface="out" name="b" initial_value="0.001378275288"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="ms" public_interface="in" name="time"/>
      <variable units="dimensionless" name="b_inf"/>
      <variable units="ms" name="tau_b"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>b_inf</ci>
            <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">14</cn>
                           </apply>
                        </apply>
                        <cn cellml:units="millivolt">10.8</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>tau_b</ci>
            <apply>
               <plus/>
               <cn cellml:units="ms">3.7</cn>
               <apply>
                  <divide/>
                  <cn cellml:units="ms">6.1</cn>
                  <apply>
                     <plus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <plus/>
                              <ci>V</ci>
                              <cn cellml:units="millivolt">25</cn>
                           </apply>
                           <cn cellml:units="millivolt">4.5</cn>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>b</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>b_inf</ci>
                  <ci>b</ci>
               </apply>
               <ci>tau_b</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="T_type_Ca_channel_g_gate">
      <variable units="dimensionless" public_interface="out" name="g" initial_value="0.988597502434"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="ms" public_interface="in" name="time"/>
      <variable units="dimensionless" name="g_inf"/>
      <variable units="ms" name="tau_g"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>g_inf</ci>
            <apply>
               <divide/>
               <cn cellml:units="dimensionless">1</cn>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <plus/>
                           <ci>V</ci>
                           <cn cellml:units="millivolt">60</cn>
                        </apply>
                        <cn cellml:units="millivolt">5.6</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>tau_g</ci>
            <piecewise>
               <piece>
                  <apply>
                     <plus/>
                     <apply>
                        <times/>
                        <apply>
                           <minus/>
                           <cn cellml:units="ms_per_millivolt">0.875</cn>
                        </apply>
                        <ci>V</ci>
                     </apply>
                     <cn cellml:units="ms">12</cn>
                  </apply>
                  <apply>
                     <leq/>
                     <ci>V</ci>
                     <cn cellml:units="millivolt">0</cn>
                  </apply>
               </piece>
               <otherwise>
                  <cn cellml:units="ms">12</cn>
               </otherwise>
            </piecewise>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>g</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>g_inf</ci>
                  <ci>g</ci>
               </apply>
               <ci>tau_g</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="rapid_delayed_rectifier_potassium_current">
      <variable units="microA_per_microF" public_interface="out" name="i_Kr"/>
      <variable units="milliS_per_microF" name="g_Kr"/>
      <variable units="milliS_per_microF" name="G_Kr" initial_value="0.02614"/>
      <variable units="dimensionless" name="Rect"/>
      <variable units="millivolt" public_interface="in" name="E_K"/>
      <variable units="ms" public_interface="in" private_interface="out" name="time"/>
      <variable units="millivolt" public_interface="in" private_interface="out" name="V"/>
      <variable units="millimolar" public_interface="in" name="Ko"/>
      <variable units="dimensionless" private_interface="in" name="xr"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>g_Kr</ci>
            <apply>
               <times/>
               <ci>G_Kr</ci>
               <apply>
                  <root/>
                  <apply>
                     <divide/>
                     <ci>Ko</ci>
                     <cn cellml:units="millimolar">5.4</cn>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>Rect</ci>
            <apply>
               <divide/>
               <cn cellml:units="dimensionless">1</cn>
               <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>
         <apply>
            <eq/>
            <ci>i_Kr</ci>
            <apply>
               <times/>
               <ci>g_Kr</ci>
               <ci>xr</ci>
               <ci>Rect</ci>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>E_K</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="rapid_delayed_rectifier_potassium_current_xr_gate">
      <variable units="dimensionless" public_interface="out" name="xr" initial_value="0.000207067204"/>
      <variable units="dimensionless" name="xr_infinity"/>
      <variable units="ms" name="tau_xr"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="ms" public_interface="in" name="time"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>xr_infinity</ci>
            <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">21.5</cn>
                           </apply>
                        </apply>
                        <cn cellml:units="millivolt">7.5</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>tau_xr</ci>
            <apply>
               <divide/>
               <cn cellml:units="ms">1</cn>
               <apply>
                  <plus/>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <cn cellml:units="per_millivolt">0.00138</cn>
                        <apply>
                           <plus/>
                           <ci>V</ci>
                           <cn cellml:units="millivolt">14.2</cn>
                        </apply>
                     </apply>
                     <apply>
                        <minus/>
                        <cn cellml:units="dimensionless">1</cn>
                        <apply>
                           <exp/>
                           <apply>
                              <times/>
                              <apply>
                                 <minus/>
                                 <cn cellml:units="per_millivolt">0.123</cn>
                              </apply>
                              <apply>
                                 <plus/>
                                 <ci>V</ci>
                                 <cn cellml:units="millivolt">14.2</cn>
                              </apply>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <cn cellml:units="per_millivolt">0.00061</cn>
                        <apply>
                           <plus/>
                           <ci>V</ci>
                           <cn cellml:units="millivolt">38.9</cn>
                        </apply>
                     </apply>
                     <apply>
                        <minus/>
                        <apply>
                           <exp/>
                           <apply>
                              <times/>
                              <cn cellml:units="per_millivolt">0.145</cn>
                              <apply>
                                 <plus/>
                                 <ci>V</ci>
                                 <cn cellml:units="millivolt">38.9</cn>
                              </apply>
                           </apply>
                        </apply>
                        <cn cellml:units="dimensionless">1</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>xr</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>xr_infinity</ci>
                  <ci>xr</ci>
               </apply>
               <ci>tau_xr</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="slow_delayed_rectifier_potassium_current">
      <variable units="microA_per_microF" public_interface="out" name="i_Ks"/>
      <variable units="milliS_per_microF" name="g_Ks"/>
      <variable units="milliS_per_microF" name="G_Ks" initial_value="0.433"/>
      <variable units="millivolt" name="E_Ks"/>
      <variable units="dimensionless" name="PNaK" initial_value="0.01833"/>
      <variable units="ms" public_interface="in" private_interface="out" name="time"/>
      <variable units="millivolt" public_interface="in" private_interface="out" name="V"/>
      <variable units="joule_per_kilomole_kelvin" public_interface="in" name="R"/>
      <variable units="kelvin" public_interface="in" name="T"/>
      <variable units="coulomb_per_mole" public_interface="in" name="F"/>
      <variable units="millimolar" public_interface="in" name="Ko"/>
      <variable units="millimolar" public_interface="in" name="Ki"/>
      <variable units="millimolar" public_interface="in" name="Nao"/>
      <variable units="millimolar" public_interface="in" name="Nai"/>
      <variable units="millimolar" public_interface="in" name="Cai"/>
      <variable units="dimensionless" private_interface="in" name="xs1"/>
      <variable units="dimensionless" private_interface="in" name="xs2"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <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>Ko</ci>
                        <apply>
                           <times/>
                           <ci>PNaK</ci>
                           <ci>Nao</ci>
                        </apply>
                     </apply>
                     <apply>
                        <plus/>
                        <ci>Ki</ci>
                        <apply>
                           <times/>
                           <ci>PNaK</ci>
                           <ci>Nai</ci>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>g_Ks</ci>
            <apply>
               <times/>
               <ci>G_Ks</ci>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <divide/>
                     <cn cellml:units="dimensionless">0.6</cn>
                     <apply>
                        <plus/>
                        <cn cellml:units="dimensionless">1</cn>
                        <apply>
                           <power/>
                           <apply>
                              <divide/>
                              <cn cellml:units="millimolar">3.8e-5</cn>
                              <ci>Cai</ci>
                           </apply>
                           <cn cellml:units="dimensionless">1.4</cn>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>i_Ks</ci>
            <apply>
               <times/>
               <ci>g_Ks</ci>
               <ci>xs1</ci>
               <ci>xs2</ci>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>E_Ks</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="slow_delayed_rectifier_potassium_current_xs1_gate">
      <variable units="dimensionless" public_interface="out" name="xs1" initial_value="0.007136102382"/>
      <variable units="dimensionless" name="xs1_infinity"/>
      <variable units="ms" name="tau_xs1"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="ms" public_interface="in" name="time"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>xs1_infinity</ci>
            <apply>
               <divide/>
               <cn cellml:units="dimensionless">1</cn>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <minus/>
                           <apply>
                              <minus/>
                              <ci>V</ci>
                              <cn cellml:units="millivolt">1.5</cn>
                           </apply>
                        </apply>
                        <cn cellml:units="millivolt">16.7</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>tau_xs1</ci>
            <apply>
               <divide/>
               <cn cellml:units="ms">1</cn>
               <apply>
                  <plus/>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <cn cellml:units="per_millivolt">7.19e-5</cn>
                        <apply>
                           <plus/>
                           <ci>V</ci>
                           <cn cellml:units="millivolt">30</cn>
                        </apply>
                     </apply>
                     <apply>
                        <minus/>
                        <cn cellml:units="dimensionless">1</cn>
                        <apply>
                           <exp/>
                           <apply>
                              <times/>
                              <apply>
                                 <minus/>
                                 <cn cellml:units="per_millivolt">0.148</cn>
                              </apply>
                              <apply>
                                 <plus/>
                                 <ci>V</ci>
                                 <cn cellml:units="millivolt">30</cn>
                              </apply>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <cn cellml:units="per_millivolt">0.000131</cn>
                        <apply>
                           <plus/>
                           <ci>V</ci>
                           <cn cellml:units="millivolt">30</cn>
                        </apply>
                     </apply>
                     <apply>
                        <minus/>
                        <apply>
                           <exp/>
                           <apply>
                              <times/>
                              <cn cellml:units="per_millivolt">0.0687</cn>
                              <apply>
                                 <plus/>
                                 <ci>V</ci>
                                 <cn cellml:units="millivolt">30</cn>
                              </apply>
                           </apply>
                        </apply>
                        <cn cellml:units="dimensionless">1</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>xs1</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>xs1_infinity</ci>
                  <ci>xs1</ci>
               </apply>
               <ci>tau_xs1</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="slow_delayed_rectifier_potassium_current_xs2_gate">
      <variable units="dimensionless" public_interface="out" name="xs2" initial_value="0.039518996812"/>
      <variable units="dimensionless" name="xs2_infinity"/>
      <variable units="ms" name="tau_xs2"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="ms" public_interface="in" name="time"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>xs2_infinity</ci>
            <apply>
               <divide/>
               <cn cellml:units="dimensionless">1</cn>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <minus/>
                           <apply>
                              <minus/>
                              <ci>V</ci>
                              <cn cellml:units="millivolt">1.5</cn>
                           </apply>
                        </apply>
                        <cn cellml:units="millivolt">16.7</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>tau_xs2</ci>
            <apply>
               <divide/>
               <cn cellml:units="ms">4</cn>
               <apply>
                  <plus/>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <cn cellml:units="per_millivolt">7.19e-5</cn>
                        <apply>
                           <plus/>
                           <ci>V</ci>
                           <cn cellml:units="millivolt">30</cn>
                        </apply>
                     </apply>
                     <apply>
                        <minus/>
                        <cn cellml:units="dimensionless">1</cn>
                        <apply>
                           <exp/>
                           <apply>
                              <times/>
                              <apply>
                                 <minus/>
                                 <cn cellml:units="per_millivolt">0.148</cn>
                              </apply>
                              <apply>
                                 <plus/>
                                 <ci>V</ci>
                                 <cn cellml:units="millivolt">30</cn>
                              </apply>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <cn cellml:units="per_millivolt">0.000131</cn>
                        <apply>
                           <plus/>
                           <ci>V</ci>
                           <cn cellml:units="millivolt">30</cn>
                        </apply>
                     </apply>
                     <apply>
                        <minus/>
                        <apply>
                           <exp/>
                           <apply>
                              <times/>
                              <cn cellml:units="per_millivolt">0.0687</cn>
                              <apply>
                                 <plus/>
                                 <ci>V</ci>
                                 <cn cellml:units="millivolt">30</cn>
                              </apply>
                           </apply>
                        </apply>
                        <cn cellml:units="dimensionless">1</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>xs2</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>xs2_infinity</ci>
                  <ci>xs2</ci>
               </apply>
               <ci>tau_xs2</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="time_independent_potassium_current">
      <variable units="microA_per_microF" public_interface="out" name="i_K1"/>
      <variable units="millivolt" public_interface="out" private_interface="out" name="E_K"/>
      <variable units="milliS_per_microF" name="g_K1"/>
      <variable units="milliS_per_microF" name="G_K1" initial_value="0.75"/>
      <variable units="millivolt" public_interface="in" private_interface="out" name="V"/>
      <variable units="millimolar" public_interface="in" name="Ko"/>
      <variable units="millimolar" public_interface="in" name="Ki"/>
      <variable units="joule_per_kilomole_kelvin" public_interface="in" name="R"/>
      <variable units="kelvin" public_interface="in" name="T"/>
      <variable units="coulomb_per_mole" public_interface="in" name="F"/>
      <variable units="dimensionless" private_interface="in" name="K1_infinity"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>g_K1</ci>
            <apply>
               <times/>
               <ci>G_K1</ci>
               <apply>
                  <root/>
                  <apply>
                     <divide/>
                     <ci>Ko</ci>
                     <cn cellml:units="millimolar">5.4</cn>
                  </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>Ko</ci>
                     <ci>Ki</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>i_K1</ci>
            <apply>
               <times/>
               <ci>g_K1</ci>
               <ci>K1_infinity</ci>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>E_K</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="time_independent_potassium_current_K1_gate">
      <variable units="dimensionless" public_interface="out" name="K1_infinity"/>
      <variable units="per_ms" name="alpha_K1"/>
      <variable units="per_ms" name="beta_K1"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="millivolt" public_interface="in" name="E_K"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>alpha_K1</ci>
            <apply>
               <divide/>
               <cn cellml:units="per_ms">1.02</cn>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <times/>
                        <cn cellml:units="per_millivolt">0.2385</cn>
                        <apply>
                           <minus/>
                           <apply>
                              <minus/>
                              <ci>V</ci>
                              <ci>E_K</ci>
                           </apply>
                           <cn cellml:units="millivolt">59.215</cn>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>beta_K1</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <cn cellml:units="per_ms">1</cn>
                  <apply>
                     <plus/>
                     <apply>
                        <times/>
                        <cn cellml:units="dimensionless">0.49124</cn>
                        <apply>
                           <exp/>
                           <apply>
                              <times/>
                              <cn cellml:units="per_millivolt">0.08032</cn>
                              <apply>
                                 <plus/>
                                 <apply>
                                    <minus/>
                                    <ci>V</ci>
                                    <ci>E_K</ci>
                                 </apply>
                                 <cn cellml:units="millivolt">5.476</cn>
                              </apply>
                           </apply>
                        </apply>
                     </apply>
                     <apply>
                        <exp/>
                        <apply>
                           <times/>
                           <cn cellml:units="per_millivolt">0.06175</cn>
                           <apply>
                              <minus/>
                              <apply>
                                 <minus/>
                                 <ci>V</ci>
                                 <ci>E_K</ci>
                              </apply>
                              <cn cellml:units="millivolt">594.31</cn>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <times/>
                        <apply>
                           <minus/>
                           <cn cellml:units="per_millivolt">0.5143</cn>
                        </apply>
                        <apply>
                           <plus/>
                           <apply>
                              <minus/>
                              <ci>V</ci>
                              <ci>E_K</ci>
                           </apply>
                           <cn cellml:units="millivolt">4.753</cn>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>K1_infinity</ci>
            <apply>
               <divide/>
               <ci>alpha_K1</ci>
               <apply>
                  <plus/>
                  <ci>alpha_K1</ci>
                  <ci>beta_K1</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="plateau_potassium_current">
      <variable units="microA_per_microF" public_interface="out" name="i_Kp"/>
      <variable units="milliS_per_microF" name="g_Kp" initial_value="0.00552"/>
      <variable units="dimensionless" name="Kp"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="millivolt" public_interface="in" name="E_K"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>Kp</ci>
            <apply>
               <divide/>
               <cn cellml:units="dimensionless">1</cn>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <minus/>
                           <cn cellml:units="millivolt">7.488</cn>
                           <ci>V</ci>
                        </apply>
                        <cn cellml:units="millivolt">5.98</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>i_Kp</ci>
            <apply>
               <times/>
               <ci>g_Kp</ci>
               <ci>Kp</ci>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>E_K</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="sodium_activated_potassium_current">
      <variable units="microA_per_microF" public_interface="out" name="i_K_Na"/>
      <variable units="milliS_per_microF" name="g_K_Na" initial_value="0.12848"/>
      <variable units="dimensionless" name="nKNa" initial_value="2.8"/>
      <variable units="dimensionless" name="pona"/>
      <variable units="dimensionless" name="pov"/>
      <variable units="millimolar" name="kdKNa" initial_value="66"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="millimolar" public_interface="in" name="Nai"/>
      <variable units="millivolt" public_interface="in" name="E_K"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>pona</ci>
            <apply>
               <divide/>
               <cn cellml:units="dimensionless">0.85</cn>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <power/>
                     <apply>
                        <divide/>
                        <ci>kdKNa</ci>
                        <ci>Nai</ci>
                     </apply>
                     <ci>nKNa</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>pov</ci>
            <apply>
               <minus/>
               <cn cellml:units="dimensionless">0.8</cn>
               <apply>
                  <divide/>
                  <cn cellml:units="dimensionless">0.65</cn>
                  <apply>
                     <plus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <plus/>
                              <ci>V</ci>
                              <cn cellml:units="millivolt">125</cn>
                           </apply>
                           <cn cellml:units="millivolt">15</cn>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>i_K_Na</ci>
            <apply>
               <times/>
               <ci>g_K_Na</ci>
               <ci>pona</ci>
               <ci>pov</ci>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>E_K</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="ATP_sensitive_potassium_current">
      <variable units="microA_per_microF" public_interface="out" name="i_K_ATP"/>
      <variable units="milliS_per_microF" name="g_K_ATP"/>
      <variable units="dimensionless" name="i_K_ATP_on" initial_value="1"/>
      <variable units="dimensionless" name="nATP" initial_value="0.24"/>
      <variable units="dimensionless" name="nicholsarea" initial_value="5e-5"/>
      <variable units="millimolar" name="ATPi" initial_value="3"/>
      <variable units="dimensionless" name="hATP" initial_value="2"/>
      <variable units="millimolar" name="kATP" initial_value="0.00025"/>
      <variable units="dimensionless" name="pATP"/>
      <variable units="milliS_per_microF" name="GKbaraATP"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="millimolar" public_interface="in" name="Ko"/>
      <variable units="millivolt" public_interface="in" name="E_K"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>g_K_ATP</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <ci>i_K_ATP_on</ci>
                  <cn cellml:units="milliS_per_microF">0.000193</cn>
               </apply>
               <ci>nicholsarea</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>pATP</ci>
            <apply>
               <divide/>
               <cn cellml:units="dimensionless">1</cn>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <power/>
                     <apply>
                        <divide/>
                        <ci>ATPi</ci>
                        <ci>kATP</ci>
                     </apply>
                     <ci>hATP</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>GKbaraATP</ci>
            <apply>
               <times/>
               <ci>g_K_ATP</ci>
               <ci>pATP</ci>
               <apply>
                  <power/>
                  <apply>
                     <divide/>
                     <ci>Ko</ci>
                     <cn cellml:units="millimolar">4</cn>
                  </apply>
                  <ci>nATP</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>i_K_ATP</ci>
            <apply>
               <times/>
               <ci>GKbaraATP</ci>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>E_K</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="transient_outward_current">
      <variable units="microA_per_microF" public_interface="out" name="i_to"/>
      <variable units="milliS_per_microF" name="g_to"/>
      <variable units="millivolt" public_interface="in" name="E_K"/>
      <variable units="millivolt" public_interface="in" private_interface="out" name="V"/>
      <variable units="ms" public_interface="in" private_interface="out" name="time"/>
      <variable units="dimensionless" name="rvdv"/>
      <variable units="dimensionless" private_interface="in" name="zdv"/>
      <variable units="dimensionless" private_interface="in" name="ydv"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>g_to</ci>
            <apply>
               <times/>
               <cn cellml:units="dimensionless">0</cn>
               <cn cellml:units="milliS_per_microF">0.5</cn>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>rvdv</ci>
            <apply>
               <exp/>
               <apply>
                  <divide/>
                  <ci>V</ci>
                  <cn cellml:units="millivolt">100</cn>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>i_to</ci>
            <apply>
               <times/>
               <ci>g_to</ci>
               <apply>
                  <power/>
                  <ci>zdv</ci>
                  <cn cellml:units="dimensionless">3</cn>
               </apply>
               <ci>ydv</ci>
               <ci>rvdv</ci>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>E_K</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="transient_outward_current_zdv_gate">
      <variable units="dimensionless" public_interface="out" name="zdv" initial_value="0.014537782303"/>
      <variable units="per_ms" name="alpha_zdv"/>
      <variable units="per_ms" name="beta_zdv"/>
      <variable units="ms" name="tau_zdv"/>
      <variable units="dimensionless" name="zdv_ss"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="ms" public_interface="in" name="time"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>alpha_zdv</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <cn cellml:units="per_ms">10</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <minus/>
                           <ci>V</ci>
                           <cn cellml:units="millivolt">40</cn>
                        </apply>
                        <cn cellml:units="millivolt">25</cn>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <minus/>
                           <ci>V</ci>
                           <cn cellml:units="millivolt">40</cn>
                        </apply>
                        <cn cellml:units="millivolt">25</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>beta_zdv</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <cn cellml:units="per_ms">10</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <minus/>
                           <apply>
                              <plus/>
                              <ci>V</ci>
                              <cn cellml:units="millivolt">90</cn>
                           </apply>
                        </apply>
                        <cn cellml:units="millivolt">25</cn>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <minus/>
                           <apply>
                              <plus/>
                              <ci>V</ci>
                              <cn cellml:units="millivolt">90</cn>
                           </apply>
                        </apply>
                        <cn cellml:units="millivolt">25</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>tau_zdv</ci>
            <apply>
               <divide/>
               <cn cellml:units="dimensionless">1</cn>
               <apply>
                  <plus/>
                  <ci>alpha_zdv</ci>
                  <ci>beta_zdv</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>zdv_ss</ci>
            <apply>
               <divide/>
               <ci>alpha_zdv</ci>
               <apply>
                  <plus/>
                  <ci>alpha_zdv</ci>
                  <ci>beta_zdv</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>zdv</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>zdv_ss</ci>
                  <ci>zdv</ci>
               </apply>
               <ci>tau_zdv</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="transient_outward_current_ydv_gate">
      <variable units="dimensionless" public_interface="out" name="ydv" initial_value="0.99993940527"/>
      <variable units="per_ms" name="alpha_ydv"/>
      <variable units="per_ms" name="beta_ydv"/>
      <variable units="ms" name="tau_ydv"/>
      <variable units="dimensionless" name="ydv_ss"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="ms" public_interface="in" name="time"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>alpha_ydv</ci>
            <apply>
               <divide/>
               <cn cellml:units="per_ms">0.015</cn>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <plus/>
                           <ci>V</ci>
                           <cn cellml:units="millivolt">60</cn>
                        </apply>
                        <cn cellml:units="millivolt">5</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>beta_ydv</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <cn cellml:units="per_ms">0.1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <plus/>
                           <ci>V</ci>
                           <cn cellml:units="millivolt">25</cn>
                        </apply>
                        <cn cellml:units="millivolt">5</cn>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <plus/>
                           <ci>V</ci>
                           <cn cellml:units="millivolt">25</cn>
                        </apply>
                        <cn cellml:units="millivolt">5</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>tau_ydv</ci>
            <apply>
               <divide/>
               <cn cellml:units="dimensionless">1</cn>
               <apply>
                  <plus/>
                  <ci>alpha_ydv</ci>
                  <ci>beta_ydv</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>ydv_ss</ci>
            <apply>
               <divide/>
               <ci>alpha_ydv</ci>
               <apply>
                  <plus/>
                  <ci>alpha_ydv</ci>
                  <ci>beta_ydv</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>ydv</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>ydv_ss</ci>
                  <ci>ydv</ci>
               </apply>
               <ci>tau_ydv</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="sarcolemmal_calcium_pump">
      <variable units="microA_per_microF" public_interface="out" name="i_p_Ca"/>
      <variable units="millimolar" name="K_mpCa" initial_value="0.0005"/>
      <variable units="microA_per_microF" name="I_pCa" initial_value="1.15"/>
      <variable units="millimolar" public_interface="in" name="Cai"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_p_Ca</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <ci>I_pCa</ci>
                  <ci>Cai</ci>
               </apply>
               <apply>
                  <plus/>
                  <ci>K_mpCa</ci>
                  <ci>Cai</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="sodium_background_current">
      <variable units="microA_per_microF" public_interface="out" name="i_Na_b"/>
      <variable units="milliS_per_microF" name="g_Nab" initial_value="0.004"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="millivolt" public_interface="in" name="E_Na"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_Na_b</ci>
            <apply>
               <times/>
               <ci>g_Nab</ci>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>E_Na</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="calcium_background_current">
      <variable units="microA_per_microF" public_interface="out" name="i_Ca_b"/>
      <variable units="milliS_per_microF" name="g_Cab" initial_value="0.003016"/>
      <variable units="millivolt" public_interface="out" name="E_Ca"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="joule_per_kilomole_kelvin" public_interface="in" name="R"/>
      <variable units="kelvin" public_interface="in" name="T"/>
      <variable units="coulomb_per_mole" public_interface="in" name="F"/>
      <variable units="millimolar" public_interface="in" name="Cai"/>
      <variable units="millimolar" public_interface="in" name="Cao"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>E_Ca</ci>
            <apply>
               <times/>
               <apply>
                  <divide/>
                  <apply>
                     <times/>
                     <ci>R</ci>
                     <ci>T</ci>
                  </apply>
                  <apply>
                     <times/>
                     <cn cellml:units="dimensionless">2</cn>
                     <ci>F</ci>
                  </apply>
               </apply>
               <apply>
                  <ln/>
                  <apply>
                     <divide/>
                     <ci>Cao</ci>
                     <ci>Cai</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>i_Ca_b</ci>
            <apply>
               <times/>
               <ci>g_Cab</ci>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>E_Ca</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="sodium_potassium_pump">
      <variable units="microA_per_microF" public_interface="out" name="i_NaK"/>
      <variable units="microA_per_microF" name="I_NaK" initial_value="2.25"/>
      <variable units="dimensionless" name="f_NaK"/>
      <variable units="millimolar" name="K_mNai" initial_value="10"/>
      <variable units="millimolar" name="K_mKo" initial_value="1.5"/>
      <variable units="dimensionless" name="sigma"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="joule_per_kilomole_kelvin" public_interface="in" name="R"/>
      <variable units="kelvin" public_interface="in" name="T"/>
      <variable units="coulomb_per_mole" public_interface="in" name="F"/>
      <variable units="millimolar" public_interface="in" name="Nai"/>
      <variable units="millimolar" public_interface="in" name="Nao"/>
      <variable units="millimolar" public_interface="in" name="Ko"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>sigma</ci>
            <apply>
               <times/>
               <apply>
                  <divide/>
                  <cn cellml:units="dimensionless">1</cn>
                  <cn cellml:units="dimensionless">7</cn>
               </apply>
               <apply>
                  <minus/>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <ci>Nao</ci>
                        <cn cellml:units="millimolar">67.3</cn>
                     </apply>
                  </apply>
                  <cn cellml:units="dimensionless">1</cn>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>f_NaK</ci>
            <apply>
               <divide/>
               <cn cellml:units="dimensionless">1</cn>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <times/>
                     <cn cellml:units="dimensionless">0.1245</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <times/>
                              <apply>
                                 <minus/>
                                 <cn cellml:units="dimensionless">0.1</cn>
                              </apply>
                              <ci>V</ci>
                              <ci>F</ci>
                           </apply>
                           <apply>
                              <times/>
                              <ci>R</ci>
                              <ci>T</ci>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
                  <apply>
                     <times/>
                     <cn cellml:units="dimensionless">0.0365</cn>
                     <ci>sigma</ci>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <times/>
                              <apply>
                                 <minus/>
                                 <ci>V</ci>
                              </apply>
                              <ci>F</ci>
                           </apply>
                           <apply>
                              <times/>
                              <ci>R</ci>
                              <ci>T</ci>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>i_NaK</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <ci>I_NaK</ci>
                        <ci>f_NaK</ci>
                        <cn cellml:units="dimensionless">1</cn>
                     </apply>
                     <apply>
                        <plus/>
                        <cn cellml:units="dimensionless">1</cn>
                        <apply>
                           <power/>
                           <apply>
                              <divide/>
                              <ci>K_mNai</ci>
                              <ci>Nai</ci>
                           </apply>
                           <cn cellml:units="dimensionless">2</cn>
                        </apply>
                     </apply>
                  </apply>
                  <ci>Ko</ci>
               </apply>
               <apply>
                  <plus/>
                  <ci>Ko</ci>
                  <ci>K_mKo</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="non_specific_calcium_activated_current">
      <variable units="microA_per_microF" public_interface="out" name="i_ns_Ca"/>
      <variable units="microA_per_microF" public_interface="out" name="i_ns_Na"/>
      <variable units="microA_per_microF" public_interface="out" name="i_ns_K"/>
      <variable units="litre_per_farad_millisecond" name="P_ns_Ca"/>
      <variable units="dimensionless" public_interface="in" name="gamma_Nai"/>
      <variable units="dimensionless" public_interface="in" name="gamma_Nao"/>
      <variable units="dimensionless" public_interface="in" name="gamma_Ki"/>
      <variable units="dimensionless" public_interface="in" name="gamma_Ko"/>
      <variable units="joule_per_kilomole_kelvin" public_interface="in" name="R"/>
      <variable units="kelvin" public_interface="in" name="T"/>
      <variable units="coulomb_per_mole" public_interface="in" name="F"/>
      <variable units="millimolar" public_interface="in" name="Cao"/>
      <variable units="millimolar" public_interface="in" name="Nao"/>
      <variable units="millimolar" public_interface="in" name="Ko"/>
      <variable units="millimolar" public_interface="in" name="Nai"/>
      <variable units="millimolar" public_interface="in" name="Ki"/>
      <variable units="microA_per_microF" name="I_ns_Na"/>
      <variable units="microA_per_microF" name="I_ns_K"/>
      <variable units="millimolar" name="K_m_ns_Ca" initial_value="0.0012"/>
      <variable units="ms" public_interface="in" name="time"/>
      <variable units="millimolar" public_interface="in" name="Cai"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>P_ns_Ca</ci>
            <cn cellml:units="litre_per_farad_millisecond">1.75e-7</cn>
         </apply>
         <apply>
            <eq/>
            <ci>I_ns_Na</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <ci>P_ns_Ca</ci>
                        <apply>
                           <power/>
                           <cn cellml:units="dimensionless">1</cn>
                           <cn cellml:units="dimensionless">2</cn>
                        </apply>
                        <ci>V</ci>
                        <apply>
                           <power/>
                           <ci>F</ci>
                           <cn cellml:units="dimensionless">2</cn>
                        </apply>
                     </apply>
                     <apply>
                        <times/>
                        <ci>R</ci>
                        <ci>T</ci>
                     </apply>
                  </apply>
                  <apply>
                     <minus/>
                     <apply>
                        <times/>
                        <ci>gamma_Nai</ci>
                        <ci>Nai</ci>
                        <apply>
                           <exp/>
                           <apply>
                              <divide/>
                              <apply>
                                 <times/>
                                 <cn cellml:units="dimensionless">1</cn>
                                 <ci>V</ci>
                                 <ci>F</ci>
                              </apply>
                              <apply>
                                 <times/>
                                 <ci>R</ci>
                                 <ci>T</ci>
                              </apply>
                           </apply>
                        </apply>
                     </apply>
                     <apply>
                        <times/>
                        <ci>gamma_Nao</ci>
                        <ci>Nao</ci>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <minus/>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <times/>
                           <cn cellml:units="dimensionless">1</cn>
                           <ci>V</ci>
                           <ci>F</ci>
                        </apply>
                        <apply>
                           <times/>
                           <ci>R</ci>
                           <ci>T</ci>
                        </apply>
                     </apply>
                  </apply>
                  <cn cellml:units="dimensionless">1</cn>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>I_ns_K</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <ci>P_ns_Ca</ci>
                        <apply>
                           <power/>
                           <cn cellml:units="dimensionless">1</cn>
                           <cn cellml:units="dimensionless">2</cn>
                        </apply>
                        <ci>V</ci>
                        <apply>
                           <power/>
                           <ci>F</ci>
                           <cn cellml:units="dimensionless">2</cn>
                        </apply>
                     </apply>
                     <apply>
                        <times/>
                        <ci>R</ci>
                        <ci>T</ci>
                     </apply>
                  </apply>
                  <apply>
                     <minus/>
                     <apply>
                        <times/>
                        <ci>gamma_Ki</ci>
                        <ci>Ki</ci>
                        <apply>
                           <exp/>
                           <apply>
                              <divide/>
                              <apply>
                                 <times/>
                                 <cn cellml:units="dimensionless">1</cn>
                                 <ci>V</ci>
                                 <ci>F</ci>
                              </apply>
                              <apply>
                                 <times/>
                                 <ci>R</ci>
                                 <ci>T</ci>
                              </apply>
                           </apply>
                        </apply>
                     </apply>
                     <apply>
                        <times/>
                        <ci>gamma_Ko</ci>
                        <ci>Ko</ci>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <minus/>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <times/>
                           <cn cellml:units="dimensionless">1</cn>
                           <ci>V</ci>
                           <ci>F</ci>
                        </apply>
                        <apply>
                           <times/>
                           <ci>R</ci>
                           <ci>T</ci>
                        </apply>
                     </apply>
                  </apply>
                  <cn cellml:units="dimensionless">1</cn>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>i_ns_Na</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <ci>I_ns_Na</ci>
                  <cn cellml:units="dimensionless">1</cn>
               </apply>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <power/>
                     <apply>
                        <divide/>
                        <ci>K_m_ns_Ca</ci>
                        <ci>Cai</ci>
                     </apply>
                     <cn cellml:units="dimensionless">3</cn>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>i_ns_K</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <ci>I_ns_K</ci>
                  <cn cellml:units="dimensionless">1</cn>
               </apply>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <power/>
                     <apply>
                        <divide/>
                        <ci>K_m_ns_Ca</ci>
                        <ci>Cai</ci>
                     </apply>
                     <cn cellml:units="dimensionless">3</cn>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>i_ns_Ca</ci>
            <apply>
               <plus/>
               <ci>i_ns_Na</ci>
               <ci>i_ns_K</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="Na_Ca_exchanger">
      <variable units="microA_per_microF" public_interface="out" name="i_NaCa"/>
      <variable units="microA_per_microF" name="c1" initial_value="0.00025"/>
      <variable units="dimensionless" name="c2" initial_value="0.0001"/>
      <variable units="dimensionless" name="gamma" initial_value="0.15"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="joule_per_kilomole_kelvin" public_interface="in" name="R"/>
      <variable units="kelvin" public_interface="in" name="T"/>
      <variable units="coulomb_per_mole" public_interface="in" name="F"/>
      <variable units="millimolar" public_interface="in" name="Nai"/>
      <variable units="millimolar" public_interface="in" name="Nao"/>
      <variable units="millimolar" public_interface="in" name="Cai"/>
      <variable units="millimolar" public_interface="in" name="Cao"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_NaCa</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <ci>c1</ci>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <times/>
                           <apply>
                              <minus/>
                              <ci>gamma</ci>
                              <cn cellml:units="dimensionless">1</cn>
                           </apply>
                           <ci>V</ci>
                           <ci>F</ci>
                        </apply>
                        <apply>
                           <times/>
                           <ci>R</ci>
                           <ci>T</ci>
                        </apply>
                     </apply>
                  </apply>
                  <apply>
                     <minus/>
                     <apply>
                        <times/>
                        <apply>
                           <exp/>
                           <apply>
                              <divide/>
                              <apply>
                                 <times/>
                                 <ci>V</ci>
                                 <ci>F</ci>
                              </apply>
                              <apply>
                                 <times/>
                                 <ci>R</ci>
                                 <ci>T</ci>
                              </apply>
                           </apply>
                        </apply>
                        <apply>
                           <power/>
                           <ci>Nai</ci>
                           <cn cellml:units="dimensionless">3</cn>
                        </apply>
                        <ci>Cao</ci>
                     </apply>
                     <apply>
                        <times/>
                        <apply>
                           <power/>
                           <ci>Nao</ci>
                           <cn cellml:units="dimensionless">3</cn>
                        </apply>
                        <ci>Cai</ci>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <plus/>
                  <cn cellml:units="millimolar4">1</cn>
                  <apply>
                     <times/>
                     <ci>c2</ci>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <times/>
                              <apply>
                                 <minus/>
                                 <ci>gamma</ci>
                                 <cn cellml:units="dimensionless">1</cn>
                              </apply>
                              <ci>V</ci>
                              <ci>F</ci>
                           </apply>
                           <apply>
                              <times/>
                              <ci>R</ci>
                              <ci>T</ci>
                           </apply>
                        </apply>
                     </apply>
                     <apply>
                        <plus/>
                        <apply>
                           <times/>
                           <apply>
                              <exp/>
                              <apply>
                                 <divide/>
                                 <apply>
                                    <times/>
                                    <ci>V</ci>
                                    <ci>F</ci>
                                 </apply>
                                 <apply>
                                    <times/>
                                    <ci>R</ci>
                                    <ci>T</ci>
                                 </apply>
                              </apply>
                           </apply>
                           <apply>
                              <power/>
                              <ci>Nai</ci>
                              <cn cellml:units="dimensionless">3</cn>
                           </apply>
                           <ci>Cao</ci>
                        </apply>
                        <apply>
                           <times/>
                           <apply>
                              <power/>
                              <ci>Nao</ci>
                              <cn cellml:units="dimensionless">3</cn>
                           </apply>
                           <ci>Cai</ci>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="calcium_dynamics">
      <variable units="millimolar_per_ms" public_interface="out" name="i_rel"/>
      <variable units="millimolar_per_ms" public_interface="out" name="i_up"/>
      <variable units="millimolar_per_ms" public_interface="out" name="i_leak"/>
      <variable units="millimolar_per_ms" public_interface="out" name="i_tr"/>
      <variable units="per_ms" name="G_rel"/>
      <variable units="per_ms" name="G_rel_Viswanathan"/>
      <variable units="per_ms" name="G_rel_max" initial_value="60"/>
      <variable units="per_ms" name="G_rel_overload" initial_value="4"/>
      <variable units="ms" name="tau_tr" initial_value="180"/>
      <variable units="millimolar" name="K_mrel" initial_value="0.0008"/>
      <variable units="millimolar" name="delta_Ca_ith" initial_value="0.00018"/>
      <variable units="millimolar" name="CSQN_max" initial_value="10"/>
      <variable units="millimolar" name="K_mCSQN" initial_value="0.8"/>
      <variable units="millimolar" name="K_mup" initial_value="0.00092"/>
      <variable units="per_ms" name="K_leak"/>
      <variable units="millimolar_per_ms" name="I_up" initial_value="0.00875"/>
      <variable units="millimolar" name="Ca_NSR_max" initial_value="15"/>
      <variable units="ms" public_interface="in" name="time"/>
      <variable units="coulomb_per_mole" public_interface="in" name="F"/>
      <variable units="millimolar" public_interface="out" name="Cai" initial_value="0.000091526844"/>
      <variable units="millimolar" public_interface="out" name="Cao" initial_value="1.8"/>
      <variable units="millimolar" public_interface="out" name="Ca_JSR" initial_value="1.024584190405"/>
      <variable units="millimolar" public_interface="out" name="Ca_NSR" initial_value="1.516756041281"/>
      <variable units="micro_litre" public_interface="in" name="V_myo"/>
      <variable units="cm2" public_interface="in" name="A_cap"/>
      <variable units="micro_litre" public_interface="in" name="V_JSR"/>
      <variable units="micro_litre" public_interface="in" name="V_NSR"/>
      <variable units="microA_per_microF" public_interface="in" name="i_NaCa"/>
      <variable units="microA_per_microF" public_interface="in" name="i_CaCa"/>
      <variable units="microA_per_microF" public_interface="in" name="i_p_Ca"/>
      <variable units="microA_per_microF" public_interface="in" name="i_Ca_b"/>
      <variable units="microA_per_microF" public_interface="in" name="i_Ca_T"/>
      <variable units="millimolar" name="K_mTn" initial_value="0.0005"/>
      <variable units="millimolar" name="K_mCMDN" initial_value="0.00238"/>
      <variable units="millimolar" name="Tn_max" initial_value="0.07"/>
      <variable units="millimolar" name="CMDN_max" initial_value="0.05"/>
      <variable units="microA_per_microF" public_interface="in" name="dVdt"/>
      <variable units="dimensionless" name="APtrack" initial_value="0"/>
      <variable units="dimensionless" name="APtrack2" initial_value="0"/>
      <variable units="dimensionless" name="APtrack3" initial_value="0.00012947996"/>
      <variable units="millimolar" name="Cainfluxtrack" initial_value="-0"/>
      <variable units="dimensionless" name="OVRLDtrack" initial_value="0"/>
      <variable units="dimensionless" name="OVRLDtrack2" initial_value="0"/>
      <variable units="dimensionless" name="OVRLDtrack3" initial_value="0"/>
      <variable units="dimensionless" name="CSQNthresh" initial_value="0.7"/>
      <variable units="dimensionless" name="Logicthresh" initial_value="0.98"/>
      <variable units="dimensionless" name="RyRopen"/>
      <variable units="dimensionless" name="RyRclose"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>APtrack</ci>
            </apply>
            <piecewise>
               <piece>
                  <apply>
                     <minus/>
                     <apply>
                        <times/>
                        <cn cellml:units="per_ms">100</cn>
                        <apply>
                           <minus/>
                           <cn cellml:units="dimensionless">1</cn>
                           <ci>APtrack</ci>
                        </apply>
                     </apply>
                     <apply>
                        <times/>
                        <cn cellml:units="per_ms">0.5</cn>
                        <ci>APtrack</ci>
                     </apply>
                  </apply>
                  <apply>
                     <gt/>
                     <ci>dVdt</ci>
                     <cn cellml:units="millivolt_per_ms">150</cn>
                  </apply>
               </piece>
               <otherwise>
                  <apply>
                     <times/>
                     <apply>
                        <minus/>
                        <cn cellml:units="per_ms">0.5</cn>
                     </apply>
                     <ci>APtrack</ci>
                  </apply>
               </otherwise>
            </piecewise>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>APtrack2</ci>
            </apply>
            <piecewise>
               <piece>
                  <apply>
                     <minus/>
                     <apply>
                        <times/>
                        <cn cellml:units="per_ms">100</cn>
                        <apply>
                           <minus/>
                           <cn cellml:units="dimensionless">1</cn>
                           <ci>APtrack2</ci>
                        </apply>
                     </apply>
                     <apply>
                        <times/>
                        <cn cellml:units="per_ms">0.5</cn>
                        <ci>APtrack2</ci>
                     </apply>
                  </apply>
                  <apply>
                     <and/>
                     <apply>
                        <lt/>
                        <ci>APtrack</ci>
                        <cn cellml:units="dimensionless">0.2</cn>
                     </apply>
                     <apply>
                        <gt/>
                        <ci>APtrack</ci>
                        <cn cellml:units="dimensionless">0.18</cn>
                     </apply>
                  </apply>
               </piece>
               <otherwise>
                  <apply>
                     <times/>
                     <apply>
                        <minus/>
                        <cn cellml:units="per_ms">0.5</cn>
                     </apply>
                     <ci>APtrack2</ci>
                  </apply>
               </otherwise>
            </piecewise>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>APtrack3</ci>
            </apply>
            <piecewise>
               <piece>
                  <apply>
                     <minus/>
                     <apply>
                        <times/>
                        <cn cellml:units="per_ms">100</cn>
                        <apply>
                           <minus/>
                           <cn cellml:units="dimensionless">1</cn>
                           <ci>APtrack3</ci>
                        </apply>
                     </apply>
                     <apply>
                        <times/>
                        <cn cellml:units="per_ms">0.5</cn>
                        <ci>APtrack3</ci>
                     </apply>
                  </apply>
                  <apply>
                     <and/>
                     <apply>
                        <lt/>
                        <ci>APtrack</ci>
                        <cn cellml:units="dimensionless">0.2</cn>
                     </apply>
                     <apply>
                        <gt/>
                        <ci>APtrack</ci>
                        <cn cellml:units="dimensionless">0.18</cn>
                     </apply>
                  </apply>
               </piece>
               <otherwise>
                  <apply>
                     <times/>
                     <apply>
                        <minus/>
                        <cn cellml:units="per_ms">0.01</cn>
                     </apply>
                     <ci>APtrack3</ci>
                  </apply>
               </otherwise>
            </piecewise>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>Cainfluxtrack</ci>
            </apply>
            <piecewise>
               <piece>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <apply>
                           <minus/>
                           <cn cellml:units="microF_per_cm2">1</cn>
                        </apply>
                        <ci>A_cap</ci>
                        <apply>
                           <plus/>
                           <apply>
                              <minus/>
                              <apply>
                                 <plus/>
                                 <ci>i_CaCa</ci>
                                 <ci>i_Ca_T</ci>
                              </apply>
                              <apply>
                                 <times/>
                                 <cn cellml:units="dimensionless">2</cn>
                                 <ci>i_NaCa</ci>
                              </apply>
                           </apply>
                           <ci>i_p_Ca</ci>
                           <ci>i_Ca_b</ci>
                        </apply>
                     </apply>
                     <apply>
                        <times/>
                        <cn cellml:units="dimensionless">2</cn>
                        <ci>V_myo</ci>
                        <ci>F</ci>
                     </apply>
                  </apply>
                  <apply>
                     <gt/>
                     <ci>APtrack</ci>
                     <cn cellml:units="dimensionless">0.2</cn>
                  </apply>
               </piece>
               <piece>
                  <cn cellml:units="millimolar_per_ms">0</cn>
                  <apply>
                     <and/>
                     <apply>
                        <gt/>
                        <ci>APtrack2</ci>
                        <cn cellml:units="dimensionless">0.01</cn>
                     </apply>
                     <apply>
                        <leq/>
                        <ci>APtrack</ci>
                        <cn cellml:units="dimensionless">0.2</cn>
                     </apply>
                  </apply>
               </piece>
               <otherwise>
                  <apply>
                     <times/>
                     <apply>
                        <minus/>
                        <cn cellml:units="per_ms">0.5</cn>
                     </apply>
                     <ci>Cainfluxtrack</ci>
                  </apply>
               </otherwise>
            </piecewise>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>OVRLDtrack</ci>
            </apply>
            <piecewise>
               <piece>
                  <apply>
                     <times/>
                     <cn cellml:units="per_ms">50</cn>
                     <apply>
                        <minus/>
                        <cn cellml:units="dimensionless">1</cn>
                        <ci>OVRLDtrack</ci>
                     </apply>
                  </apply>
                  <apply>
                     <and/>
                     <apply>
                        <gt/>
                        <apply>
                           <divide/>
                           <cn cellml:units="dimensionless">1</cn>
                           <apply>
                              <plus/>
                              <cn cellml:units="dimensionless">1</cn>
                              <apply>
                                 <divide/>
                                 <ci>K_mCSQN</ci>
                                 <ci>Ca_JSR</ci>
                              </apply>
                           </apply>
                        </apply>
                        <ci>CSQNthresh</ci>
                     </apply>
                     <apply>
                        <lt/>
                        <ci>OVRLDtrack3</ci>
                        <cn cellml:units="dimensionless">0.37</cn>
                     </apply>
                     <apply>
                        <lt/>
                        <ci>APtrack3</ci>
                        <cn cellml:units="dimensionless">0.37</cn>
                     </apply>
                  </apply>
               </piece>
               <otherwise>
                  <apply>
                     <times/>
                     <apply>
                        <minus/>
                        <cn cellml:units="per_ms">0.5</cn>
                     </apply>
                     <ci>OVRLDtrack</ci>
                  </apply>
               </otherwise>
            </piecewise>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>OVRLDtrack2</ci>
            </apply>
            <piecewise>
               <piece>
                  <apply>
                     <times/>
                     <cn cellml:units="per_ms">50</cn>
                     <apply>
                        <minus/>
                        <cn cellml:units="dimensionless">1</cn>
                        <ci>OVRLDtrack2</ci>
                     </apply>
                  </apply>
                  <apply>
                     <and/>
                     <apply>
                        <gt/>
                        <ci>OVRLDtrack</ci>
                        <ci>Logicthresh</ci>
                     </apply>
                     <apply>
                        <lt/>
                        <ci>OVRLDtrack2</ci>
                        <ci>Logicthresh</ci>
                     </apply>
                  </apply>
               </piece>
               <otherwise>
                  <apply>
                     <times/>
                     <apply>
                        <minus/>
                        <cn cellml:units="per_ms">0.5</cn>
                     </apply>
                     <ci>OVRLDtrack2</ci>
                  </apply>
               </otherwise>
            </piecewise>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>OVRLDtrack3</ci>
            </apply>
            <piecewise>
               <piece>
                  <apply>
                     <times/>
                     <cn cellml:units="per_ms">50</cn>
                     <apply>
                        <minus/>
                        <cn cellml:units="dimensionless">1</cn>
                        <ci>OVRLDtrack3</ci>
                     </apply>
                  </apply>
                  <apply>
                     <and/>
                     <apply>
                        <gt/>
                        <ci>OVRLDtrack</ci>
                        <ci>Logicthresh</ci>
                     </apply>
                     <apply>
                        <lt/>
                        <ci>OVRLDtrack3</ci>
                        <ci>Logicthresh</ci>
                     </apply>
                  </apply>
               </piece>
               <otherwise>
                  <apply>
                     <times/>
                     <apply>
                        <minus/>
                        <cn cellml:units="per_ms">0.01</cn>
                     </apply>
                     <ci>OVRLDtrack3</ci>
                  </apply>
               </otherwise>
            </piecewise>
         </apply>
         <apply>
            <eq/>
            <ci>G_rel_Viswanathan</ci>
            <piecewise>
               <piece>
                  <apply>
                     <times/>
                     <apply>
                        <divide/>
                        <apply>
                           <times/>
                           <ci>G_rel_max</ci>
                           <apply>
                              <minus/>
                              <ci>Cainfluxtrack</ci>
                              <ci>delta_Ca_ith</ci>
                           </apply>
                        </apply>
                        <apply>
                           <minus/>
                           <apply>
                              <plus/>
                              <ci>K_mrel</ci>
                              <ci>Cainfluxtrack</ci>
                           </apply>
                           <ci>delta_Ca_ith</ci>
                        </apply>
                     </apply>
                     <apply>
                        <minus/>
                        <cn cellml:units="dimensionless">1</cn>
                        <ci>APtrack2</ci>
                     </apply>
                     <ci>APtrack2</ci>
                  </apply>
                  <apply>
                     <gt/>
                     <ci>Cainfluxtrack</ci>
                     <ci>delta_Ca_ith</ci>
                  </apply>
               </piece>
               <piece>
                  <apply>
                     <times/>
                     <ci>G_rel_overload</ci>
                     <apply>
                        <minus/>
                        <cn cellml:units="dimensionless">1</cn>
                        <ci>OVRLDtrack2</ci>
                     </apply>
                     <ci>OVRLDtrack2</ci>
                  </apply>
                  <apply>
                     <and/>
                     <apply>
                        <leq/>
                        <ci>Cainfluxtrack</ci>
                        <ci>delta_Ca_ith</ci>
                     </apply>
                     <apply>
                        <gt/>
                        <ci>OVRLDtrack2</ci>
                        <cn cellml:units="dimensionless">0</cn>
                     </apply>
                  </apply>
               </piece>
               <otherwise>
                  <cn cellml:units="per_ms">0</cn>
               </otherwise>
            </piecewise>
         </apply>
         <apply>
            <eq/>
            <ci>RyRopen</ci>
            <apply>
               <divide/>
               <cn cellml:units="dimensionless">1</cn>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <times/>
                     <apply>
                        <exp/>
                        <cn cellml:units="dimensionless">2</cn>
                     </apply>
                     <ci>APtrack2</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>RyRclose</ci>
            <apply>
               <minus/>
               <cn cellml:units="dimensionless">1</cn>
               <ci>RyRopen</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>G_rel</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <ci>RyRopen</ci>
                  <ci>RyRclose</ci>
                  <cn cellml:units="per_ms">150</cn>
               </apply>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <plus/>
                           <apply>
                              <minus/>
                              <apply>
                                 <plus/>
                                 <ci>i_CaCa</ci>
                                 <ci>i_Ca_T</ci>
                              </apply>
                              <apply>
                                 <times/>
                                 <cn cellml:units="dimensionless">2</cn>
                                 <ci>i_NaCa</ci>
                              </apply>
                           </apply>
                           <ci>i_p_Ca</ci>
                           <ci>i_Ca_b</ci>
                           <cn cellml:units="microA_per_microF">5</cn>
                        </apply>
                        <cn cellml:units="microA_per_microF">0.9</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>i_rel</ci>
            <apply>
               <times/>
               <ci>G_rel</ci>
               <apply>
                  <minus/>
                  <ci>Ca_JSR</ci>
                  <ci>Cai</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>i_up</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <ci>I_up</ci>
                  <ci>Cai</ci>
               </apply>
               <apply>
                  <plus/>
                  <ci>Cai</ci>
                  <ci>K_mup</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>K_leak</ci>
            <apply>
               <divide/>
               <ci>I_up</ci>
               <ci>Ca_NSR_max</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>i_leak</ci>
            <apply>
               <times/>
               <ci>K_leak</ci>
               <ci>Ca_NSR</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>i_tr</ci>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>Ca_NSR</ci>
                  <ci>Ca_JSR</ci>
               </apply>
               <ci>tau_tr</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>Ca_JSR</ci>
            </apply>
            <apply>
               <times/>
               <apply>
                  <divide/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <plus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <apply>
                        <divide/>
                        <apply>
                           <times/>
                           <ci>CSQN_max</ci>
                           <ci>K_mCSQN</ci>
                        </apply>
                        <apply>
                           <power/>
                           <apply>
                              <plus/>
                              <ci>K_mCSQN</ci>
                              <ci>Ca_JSR</ci>
                           </apply>
                           <cn cellml:units="dimensionless">2</cn>
                        </apply>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <minus/>
                  <ci>i_tr</ci>
                  <ci>i_rel</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>Ca_NSR</ci>
            </apply>
            <apply>
               <plus/>
               <apply>
                  <minus/>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <apply>
                           <minus/>
                           <ci>i_tr</ci>
                        </apply>
                        <ci>V_JSR</ci>
                     </apply>
                     <ci>V_NSR</ci>
                  </apply>
                  <ci>i_leak</ci>
               </apply>
               <ci>i_up</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>Cai</ci>
            </apply>
            <apply>
               <times/>
               <apply>
                  <divide/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <plus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <apply>
                        <divide/>
                        <apply>
                           <times/>
                           <ci>CMDN_max</ci>
                           <ci>K_mCMDN</ci>
                        </apply>
                        <apply>
                           <power/>
                           <apply>
                              <plus/>
                              <ci>K_mCMDN</ci>
                              <ci>Cai</ci>
                           </apply>
                           <cn cellml:units="dimensionless">2</cn>
                        </apply>
                     </apply>
                     <apply>
                        <divide/>
                        <apply>
                           <times/>
                           <ci>Tn_max</ci>
                           <ci>K_mTn</ci>
                        </apply>
                        <apply>
                           <power/>
                           <apply>
                              <plus/>
                              <ci>K_mTn</ci>
                              <ci>Cai</ci>
                           </apply>
                           <cn cellml:units="dimensionless">2</cn>
                        </apply>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <plus/>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <apply>
                           <minus/>
                           <cn cellml:units="microF_per_cm2">1</cn>
                        </apply>
                        <ci>A_cap</ci>
                        <apply>
                           <plus/>
                           <apply>
                              <minus/>
                              <apply>
                                 <plus/>
                                 <ci>i_CaCa</ci>
                                 <ci>i_Ca_T</ci>
                              </apply>
                              <apply>
                                 <times/>
                                 <cn cellml:units="dimensionless">2</cn>
                                 <ci>i_NaCa</ci>
                              </apply>
                           </apply>
                           <ci>i_p_Ca</ci>
                           <ci>i_Ca_b</ci>
                        </apply>
                     </apply>
                     <apply>
                        <times/>
                        <cn cellml:units="dimensionless">2</cn>
                        <ci>V_myo</ci>
                        <ci>F</ci>
                     </apply>
                  </apply>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <ci>i_rel</ci>
                        <ci>V_JSR</ci>
                     </apply>
                     <ci>V_myo</ci>
                  </apply>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <apply>
                           <minus/>
                           <ci>i_leak</ci>
                           <ci>i_up</ci>
                        </apply>
                        <ci>V_NSR</ci>
                     </apply>
                     <ci>V_myo</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="ionic_concentrations">
      <variable units="millimolar" public_interface="out" name="Nai" initial_value="10"/>
      <variable units="millimolar" public_interface="out" name="Nao" initial_value="140"/>
      <variable units="millimolar" public_interface="out" name="Ki" initial_value="144.473230653346"/>
      <variable units="millimolar" public_interface="out" name="Ko" initial_value="5.4"/>
      <variable units="cm2" public_interface="in" name="A_cap"/>
      <variable units="micro_litre" public_interface="in" name="V_myo"/>
      <variable units="ms" public_interface="in" name="time"/>
      <variable units="coulomb_per_mole" public_interface="in" name="F"/>
      <variable units="microA_per_microF" public_interface="in" name="i_Na"/>
      <variable units="microA_per_microF" public_interface="in" name="i_CaNa"/>
      <variable units="microA_per_microF" public_interface="in" name="i_Na_b"/>
      <variable units="microA_per_microF" public_interface="in" name="i_ns_Na"/>
      <variable units="microA_per_microF" public_interface="in" name="i_NaCa"/>
      <variable units="microA_per_microF" public_interface="in" name="i_NaK"/>
      <variable units="microA_per_microF" public_interface="in" name="i_CaK"/>
      <variable units="microA_per_microF" public_interface="in" name="i_Kr"/>
      <variable units="microA_per_microF" public_interface="in" name="i_Ks"/>
      <variable units="microA_per_microF" public_interface="in" name="i_K1"/>
      <variable units="microA_per_microF" public_interface="in" name="i_Kp"/>
      <variable units="microA_per_microF" public_interface="in" name="i_K_Na"/>
      <variable units="microA_per_microF" public_interface="in" name="i_K_ATP"/>
      <variable units="microA_per_microF" public_interface="in" name="i_ns_K"/>
      <variable units="microA_per_microF" public_interface="in" name="i_to"/>
      <variable units="microA_per_microF" public_interface="in" name="I_st"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>Nai</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <apply>
                     <minus/>
                     <cn cellml:units="microF_per_cm2">1</cn>
                  </apply>
                  <apply>
                     <plus/>
                     <ci>i_Na</ci>
                     <ci>i_CaNa</ci>
                     <ci>i_Na_b</ci>
                     <ci>i_ns_Na</ci>
                     <apply>
                        <times/>
                        <ci>i_NaCa</ci>
                        <cn cellml:units="dimensionless">3</cn>
                     </apply>
                     <apply>
                        <times/>
                        <ci>i_NaK</ci>
                        <cn cellml:units="dimensionless">3</cn>
                     </apply>
                  </apply>
                  <ci>A_cap</ci>
               </apply>
               <apply>
                  <times/>
                  <ci>V_myo</ci>
                  <ci>F</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>Ki</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <apply>
                     <minus/>
                     <cn cellml:units="microF_per_cm2">1</cn>
                  </apply>
                  <apply>
                     <plus/>
                     <ci>I_st</ci>
                     <ci>i_CaK</ci>
                     <ci>i_Kr</ci>
                     <ci>i_Ks</ci>
                     <ci>i_K1</ci>
                     <ci>i_Kp</ci>
                     <ci>i_K_Na</ci>
                     <ci>i_K_ATP</ci>
                     <ci>i_to</ci>
                     <ci>i_ns_K</ci>
                     <apply>
                        <times/>
                        <apply>
                           <minus/>
                           <ci>i_NaK</ci>
                        </apply>
                        <cn cellml:units="dimensionless">2</cn>
                     </apply>
                  </apply>
                  <ci>A_cap</ci>
               </apply>
               <apply>
                  <times/>
                  <ci>V_myo</ci>
                  <ci>F</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="geometry">
      <variable units="mm" name="preplength" initial_value="0.1"/>
      <variable units="mm" name="radius" initial_value="0.011"/>
      <variable units="micro_litre" name="volume"/>
      <variable units="micro_litre" public_interface="out" name="V_myo"/>
      <variable units="micro_litre" public_interface="out" name="V_JSR"/>
      <variable units="micro_litre" public_interface="out" name="V_NSR"/>
      <variable units="cm2" public_interface="out" name="A_cap" initial_value="1.534e-4"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>volume</ci>
            <apply>
               <times/>
               <pi/>
               <ci>preplength</ci>
               <apply>
                  <power/>
                  <ci>radius</ci>
                  <cn cellml:units="dimensionless">2</cn>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>V_myo</ci>
            <apply>
               <times/>
               <cn cellml:units="dimensionless">0.68</cn>
               <ci>volume</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>V_JSR</ci>
            <apply>
               <times/>
               <cn cellml:units="dimensionless">0.0048</cn>
               <ci>volume</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>V_NSR</ci>
            <apply>
               <times/>
               <cn cellml:units="dimensionless">0.0552</cn>
               <ci>volume</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="transmembrane_currents"/>
   <component name="intracellular_currents"/>
   <component name="Ions_n_reversal_potentials"/>
   <group>
      <relationship_ref relationship="containment"/>
      <component_ref component="cell">
         <component_ref component="geometry"/>
         <component_ref component="transmembrane_currents">
            <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="fast_sodium_current_j_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_f_Ca_gate"/>
            </component_ref>
            <component_ref component="T_type_Ca_channel">
               <component_ref component="T_type_Ca_channel_b_gate"/>
               <component_ref component="T_type_Ca_channel_g_gate"/>
            </component_ref>
            <component_ref component="rapid_delayed_rectifier_potassium_current">
               <component_ref component="rapid_delayed_rectifier_potassium_current_xr_gate"/>
            </component_ref>
            <component_ref component="slow_delayed_rectifier_potassium_current">
               <component_ref component="slow_delayed_rectifier_potassium_current_xs1_gate"/>
               <component_ref component="slow_delayed_rectifier_potassium_current_xs2_gate"/>
            </component_ref>
            <component_ref component="time_independent_potassium_current">
               <component_ref component="time_independent_potassium_current_K1_gate"/>
            </component_ref>
            <component_ref component="transient_outward_current">
               <component_ref component="transient_outward_current_zdv_gate"/>
               <component_ref component="transient_outward_current_ydv_gate"/>
            </component_ref>
            <component_ref component="Na_Ca_exchanger"/>
            <component_ref component="plateau_potassium_current"/>
            <component_ref component="sodium_activated_potassium_current"/>
            <component_ref component="ATP_sensitive_potassium_current"/>
            <component_ref component="sodium_background_current"/>
            <component_ref component="calcium_background_current"/>
            <component_ref component="sodium_potassium_pump"/>
            <component_ref component="non_specific_calcium_activated_current"/>
         </component_ref>
         <component_ref component="intracellular_currents">
            <component_ref component="sarcolemmal_calcium_pump"/>
         </component_ref>
         <component_ref component="Ions_n_reversal_potentials">
            <component_ref component="ionic_concentrations"/>
            <component_ref component="calcium_dynamics"/>
         </component_ref>
      </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="fast_sodium_current_j_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_f_Ca_gate"/>
      </component_ref>
      <component_ref component="T_type_Ca_channel">
         <component_ref component="T_type_Ca_channel_b_gate"/>
         <component_ref component="T_type_Ca_channel_g_gate"/>
      </component_ref>
      <component_ref component="rapid_delayed_rectifier_potassium_current">
         <component_ref component="rapid_delayed_rectifier_potassium_current_xr_gate"/>
      </component_ref>
      <component_ref component="slow_delayed_rectifier_potassium_current">
         <component_ref component="slow_delayed_rectifier_potassium_current_xs1_gate"/>
         <component_ref component="slow_delayed_rectifier_potassium_current_xs2_gate"/>
      </component_ref>
      <component_ref component="time_independent_potassium_current">
         <component_ref component="time_independent_potassium_current_K1_gate"/>
      </component_ref>
      <component_ref component="transient_outward_current">
         <component_ref component="transient_outward_current_zdv_gate"/>
         <component_ref component="transient_outward_current_ydv_gate"/>
      </component_ref>
   </group>
   <connection>
      <map_components component_2="environment" component_1="cell"/>
      <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="T_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="transient_outward_current"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="environment" component_1="non_specific_calcium_activated_current"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="environment" component_1="ionic_concentrations"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="environment" component_1="calcium_dynamics"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="fast_sodium_current" component_1="cell"/>
      <map_variables variable_2="V" variable_1="V"/>
      <map_variables variable_2="i_Na" variable_1="i_Na"/>
      <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="L_type_Ca_channel" component_1="cell"/>
      <map_variables variable_2="V" variable_1="V"/>
      <map_variables variable_2="i_Ca_L" variable_1="i_Ca_L"/>
      <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="T_type_Ca_channel" component_1="cell"/>
      <map_variables variable_2="V" variable_1="V"/>
      <map_variables variable_2="i_Ca_T" variable_1="i_Ca_T"/>
   </connection>
   <connection>
      <map_components component_2="rapid_delayed_rectifier_potassium_current" component_1="cell"/>
      <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="cell"/>
      <map_variables variable_2="V" variable_1="V"/>
      <map_variables variable_2="i_Ks" variable_1="i_Ks"/>
      <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="cell"/>
      <map_variables variable_2="V" variable_1="V"/>
      <map_variables variable_2="i_K_Na" variable_1="i_K_Na"/>
   </connection>
   <connection>
      <map_components component_2="ATP_sensitive_potassium_current" component_1="cell"/>
      <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="transient_outward_current" component_1="cell"/>
      <map_variables variable_2="V" variable_1="V"/>
      <map_variables variable_2="i_to" variable_1="i_to"/>
   </connection>
   <connection>
      <map_components component_2="Na_Ca_exchanger" component_1="cell"/>
      <map_variables variable_2="V" variable_1="V"/>
      <map_variables variable_2="i_NaCa" variable_1="i_NaCa"/>
      <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="time_independent_potassium_current" component_1="cell"/>
      <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="plateau_potassium_current" component_1="cell"/>
      <map_variables variable_2="V" variable_1="V"/>
      <map_variables variable_2="i_Kp" variable_1="i_Kp"/>
   </connection>
   <connection>
      <map_components component_2="sarcolemmal_calcium_pump" component_1="cell"/>
      <map_variables variable_2="i_p_Ca" variable_1="i_p_Ca"/>
   </connection>
   <connection>
      <map_components component_2="sodium_background_current" component_1="cell"/>
      <map_variables variable_2="V" variable_1="V"/>
      <map_variables variable_2="i_Na_b" variable_1="i_Na_b"/>
   </connection>
   <connection>
      <map_components component_2="calcium_background_current" component_1="cell"/>
      <map_variables variable_2="V" variable_1="V"/>
      <map_variables variable_2="i_Ca_b" variable_1="i_Ca_b"/>
      <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_potassium_pump" component_1="cell"/>
      <map_variables variable_2="V" variable_1="V"/>
      <map_variables variable_2="i_NaK" variable_1="i_NaK"/>
      <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="non_specific_calcium_activated_current" component_1="cell"/>
      <map_variables variable_2="V" variable_1="V"/>
      <map_variables variable_2="R" variable_1="R"/>
      <map_variables variable_2="T" variable_1="T"/>
      <map_variables variable_2="F" variable_1="F"/>
      <map_variables variable_2="i_ns_Ca" variable_1="i_ns_Ca"/>
   </connection>
   <connection>
      <map_components component_2="ionic_concentrations" component_1="cell"/>
      <map_variables variable_2="F" variable_1="F"/>
      <map_variables variable_2="I_st" variable_1="I_st"/>
   </connection>
   <connection>
      <map_components component_2="calcium_dynamics" component_1="cell"/>
      <map_variables variable_2="F" variable_1="F"/>
      <map_variables variable_2="dVdt" variable_1="dVdt"/>
   </connection>
   <connection>
      <map_components component_2="ionic_concentrations" component_1="fast_sodium_current"/>
      <map_variables variable_2="i_Na" variable_1="i_Na"/>
      <map_variables variable_2="Nao" variable_1="Nao"/>
      <map_variables variable_2="Nai" variable_1="Nai"/>
   </connection>
   <connection>
      <map_components component_2="sodium_background_current" component_1="fast_sodium_current"/>
      <map_variables variable_2="E_Na" variable_1="E_Na"/>
   </connection>
   <connection>
      <map_components component_2="ionic_concentrations" component_1="L_type_Ca_channel"/>
      <map_variables variable_2="Nao" variable_1="Nao"/>
      <map_variables variable_2="Nai" variable_1="Nai"/>
      <map_variables variable_2="Ko" variable_1="Ko"/>
      <map_variables variable_2="Ki" variable_1="Ki"/>
      <map_variables variable_2="i_CaNa" variable_1="i_CaNa"/>
      <map_variables variable_2="i_CaK" variable_1="i_CaK"/>
   </connection>
   <connection>
      <map_components component_2="calcium_dynamics" component_1="L_type_Ca_channel"/>
      <map_variables variable_2="Cai" variable_1="Cai"/>
      <map_variables variable_2="Cao" variable_1="Cao"/>
      <map_variables variable_2="i_CaCa" variable_1="i_CaCa"/>
   </connection>
   <connection>
      <map_components component_2="calcium_dynamics" component_1="T_type_Ca_channel"/>
      <map_variables variable_2="i_Ca_T" variable_1="i_Ca_T"/>
   </connection>
   <connection>
      <map_components component_2="ionic_concentrations" component_1="slow_delayed_rectifier_potassium_current"/>
      <map_variables variable_2="Ki" variable_1="Ki"/>
      <map_variables variable_2="Ko" variable_1="Ko"/>
      <map_variables variable_2="Nao" variable_1="Nao"/>
      <map_variables variable_2="Nai" variable_1="Nai"/>
      <map_variables variable_2="i_Ks" variable_1="i_Ks"/>
   </connection>
   <connection>
      <map_components component_2="calcium_dynamics" component_1="slow_delayed_rectifier_potassium_current"/>
      <map_variables variable_2="Cai" variable_1="Cai"/>
   </connection>
   <connection>
      <map_components component_2="ionic_concentrations" component_1="rapid_delayed_rectifier_potassium_current"/>
      <map_variables variable_2="Ko" variable_1="Ko"/>
      <map_variables variable_2="i_Kr" variable_1="i_Kr"/>
   </connection>
   <connection>
      <map_components component_2="ionic_concentrations" component_1="transient_outward_current"/>
      <map_variables variable_2="i_to" variable_1="i_to"/>
   </connection>
   <connection>
      <map_components component_2="ionic_concentrations" component_1="Na_Ca_exchanger"/>
      <map_variables variable_2="Nai" variable_1="Nai"/>
      <map_variables variable_2="Nao" variable_1="Nao"/>
      <map_variables variable_2="i_NaCa" variable_1="i_NaCa"/>
   </connection>
   <connection>
      <map_components component_2="calcium_dynamics" component_1="Na_Ca_exchanger"/>
      <map_variables variable_2="Cai" variable_1="Cai"/>
      <map_variables variable_2="Cao" variable_1="Cao"/>
      <map_variables variable_2="i_NaCa" variable_1="i_NaCa"/>
   </connection>
   <connection>
      <map_components component_2="ionic_concentrations" component_1="time_independent_potassium_current"/>
      <map_variables variable_2="Ki" variable_1="Ki"/>
      <map_variables variable_2="Ko" variable_1="Ko"/>
      <map_variables variable_2="i_K1" variable_1="i_K1"/>
   </connection>
   <connection>
      <map_components component_2="ionic_concentrations" component_1="sodium_activated_potassium_current"/>
      <map_variables variable_2="Nai" variable_1="Nai"/>
      <map_variables variable_2="i_K_Na" variable_1="i_K_Na"/>
   </connection>
   <connection>
      <map_components component_2="ionic_concentrations" component_1="ATP_sensitive_potassium_current"/>
      <map_variables variable_2="Ko" variable_1="Ko"/>
      <map_variables variable_2="i_K_ATP" variable_1="i_K_ATP"/>
   </connection>
   <connection>
      <map_components component_2="time_independent_potassium_current" component_1="plateau_potassium_current"/>
      <map_variables variable_2="E_K" variable_1="E_K"/>
   </connection>
   <connection>
      <map_components component_2="time_independent_potassium_current" component_1="rapid_delayed_rectifier_potassium_current"/>
      <map_variables variable_2="E_K" variable_1="E_K"/>
   </connection>
   <connection>
      <map_components component_2="time_independent_potassium_current" component_1="sodium_activated_potassium_current"/>
      <map_variables variable_2="E_K" variable_1="E_K"/>
   </connection>
   <connection>
      <map_components component_2="time_independent_potassium_current" component_1="ATP_sensitive_potassium_current"/>
      <map_variables variable_2="E_K" variable_1="E_K"/>
   </connection>
   <connection>
      <map_components component_2="time_independent_potassium_current" component_1="transient_outward_current"/>
      <map_variables variable_2="E_K" variable_1="E_K"/>
   </connection>
   <connection>
      <map_components component_2="ionic_concentrations" component_1="plateau_potassium_current"/>
      <map_variables variable_2="i_Kp" variable_1="i_Kp"/>
   </connection>
   <connection>
      <map_components component_2="calcium_dynamics" component_1="sarcolemmal_calcium_pump"/>
      <map_variables variable_2="i_p_Ca" variable_1="i_p_Ca"/>
      <map_variables variable_2="Cai" variable_1="Cai"/>
   </connection>
   <connection>
      <map_components component_2="calcium_background_current" component_1="T_type_Ca_channel"/>
      <map_variables variable_2="E_Ca" variable_1="E_Ca"/>
   </connection>
   <connection>
      <map_components component_2="ionic_concentrations" component_1="sodium_background_current"/>
      <map_variables variable_2="i_Na_b" variable_1="i_Na_b"/>
   </connection>
   <connection>
      <map_components component_2="calcium_dynamics" component_1="calcium_background_current"/>
      <map_variables variable_2="Cai" variable_1="Cai"/>
      <map_variables variable_2="Cao" variable_1="Cao"/>
      <map_variables variable_2="i_Ca_b" variable_1="i_Ca_b"/>
   </connection>
   <connection>
      <map_components component_2="ionic_concentrations" component_1="sodium_potassium_pump"/>
      <map_variables variable_2="Nai" variable_1="Nai"/>
      <map_variables variable_2="Nao" variable_1="Nao"/>
      <map_variables variable_2="Ko" variable_1="Ko"/>
      <map_variables variable_2="i_NaK" variable_1="i_NaK"/>
   </connection>
   <connection>
      <map_components component_2="ionic_concentrations" component_1="non_specific_calcium_activated_current"/>
      <map_variables variable_2="Nao" variable_1="Nao"/>
      <map_variables variable_2="Nai" variable_1="Nai"/>
      <map_variables variable_2="Ko" variable_1="Ko"/>
      <map_variables variable_2="Ki" variable_1="Ki"/>
      <map_variables variable_2="i_ns_Na" variable_1="i_ns_Na"/>
      <map_variables variable_2="i_ns_K" variable_1="i_ns_K"/>
   </connection>
   <connection>
      <map_components component_2="calcium_dynamics" component_1="non_specific_calcium_activated_current"/>
      <map_variables variable_2="Cai" variable_1="Cai"/>
      <map_variables variable_2="Cao" variable_1="Cao"/>
   </connection>
   <connection>
      <map_components component_2="non_specific_calcium_activated_current" component_1="L_type_Ca_channel"/>
      <map_variables variable_2="gamma_Nao" variable_1="gamma_Nao"/>
      <map_variables variable_2="gamma_Nai" variable_1="gamma_Nai"/>
      <map_variables variable_2="gamma_Ko" variable_1="gamma_Ko"/>
      <map_variables variable_2="gamma_Ki" variable_1="gamma_Ki"/>
   </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="fast_sodium_current_j_gate" component_1="fast_sodium_current"/>
      <map_variables variable_2="j" variable_1="j"/>
      <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_f_Ca_gate" component_1="L_type_Ca_channel"/>
      <map_variables variable_2="f_Ca" variable_1="f_Ca"/>
      <map_variables variable_2="Cai" variable_1="Cai"/>
      <map_variables variable_2="V" variable_1="V"/>
   </connection>
   <connection>
      <map_components component_2="T_type_Ca_channel_b_gate" component_1="T_type_Ca_channel"/>
      <map_variables variable_2="b" variable_1="b"/>
      <map_variables variable_2="time" variable_1="time"/>
      <map_variables variable_2="V" variable_1="V"/>
   </connection>
   <connection>
      <map_components component_2="T_type_Ca_channel_g_gate" component_1="T_type_Ca_channel"/>
      <map_variables variable_2="g" variable_1="g"/>
      <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_xr_gate" component_1="rapid_delayed_rectifier_potassium_current"/>
      <map_variables variable_2="xr" variable_1="xr"/>
      <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_xs1_gate" component_1="slow_delayed_rectifier_potassium_current"/>
      <map_variables variable_2="xs1" variable_1="xs1"/>
      <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_xs2_gate" component_1="slow_delayed_rectifier_potassium_current"/>
      <map_variables variable_2="xs2" variable_1="xs2"/>
      <map_variables variable_2="time" variable_1="time"/>
      <map_variables variable_2="V" variable_1="V"/>
   </connection>
   <connection>
      <map_components component_2="time_independent_potassium_current_K1_gate" component_1="time_independent_potassium_current"/>
      <map_variables variable_2="K1_infinity" variable_1="K1_infinity"/>
      <map_variables variable_2="E_K" variable_1="E_K"/>
      <map_variables variable_2="V" variable_1="V"/>
   </connection>
   <connection>
      <map_components component_2="transient_outward_current_ydv_gate" component_1="transient_outward_current"/>
      <map_variables variable_2="ydv" variable_1="ydv"/>
      <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_zdv_gate" component_1="transient_outward_current"/>
      <map_variables variable_2="zdv" variable_1="zdv"/>
      <map_variables variable_2="time" variable_1="time"/>
      <map_variables variable_2="V" variable_1="V"/>
   </connection>
   <connection>
      <map_components component_2="ionic_concentrations" component_1="geometry"/>
      <map_variables variable_2="V_myo" variable_1="V_myo"/>
      <map_variables variable_2="A_cap" variable_1="A_cap"/>
   </connection>
   <connection>
      <map_components component_2="geometry" component_1="calcium_dynamics"/>
      <map_variables variable_2="V_myo" variable_1="V_myo"/>
      <map_variables variable_2="A_cap" variable_1="A_cap"/>
      <map_variables variable_2="V_JSR" variable_1="V_JSR"/>
      <map_variables variable_2="V_NSR" variable_1="V_NSR"/>
   </connection>
<rdf:RDF>
  <rdf:Seq rdf:about="rdf:#9b6b12f4-7c11-4a19-a5d0-b72a173411ff">
    <rdf:li rdf:resource="rdf:#dba32688-0867-4a2a-acfa-18cbf0261435"/>
    <rdf:li rdf:resource="rdf:#e7b77551-d534-4391-8f1a-73f244372afd"/>
  </rdf:Seq>
  <rdf:Description rdf:about="rdf:#52b20981-fb3e-431a-be2c-d0abd9b365f6">
    <dc:creator rdf:resource="rdf:#05739223-6287-4bc4-b98b-4191b106c2ae"/>
    <rdf:value>Since the model was published in 1994, it has been updated several times, incorporating new experimental data taken from papers written by other authors.  These modifications have been considered in the model presented here.  To illustrate how the Luo-Rudy II dynamic model has been developed, see the sequence of diagrams.  The 1997 Shaw and Rudy model diagram is identical to the 1999 Viswanathan model diagram since the only difference between the two models is that the Viswanathan model introduces a second activation gate (X&lt;subscript&gt;s2&lt;/subscript&gt;) into the formulation of I&lt;subscript&gt;Ks&lt;/subscript&gt;.</rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#b2535049-2e20-4454-9c79-7b61b3fb2d79">
    <vCard:N rdf:resource="rdf:#250e831d-159c-4178-8bfd-ab2851fbaab4"/>
  </rdf:Description>
  <rdf:Description rdf:about="">
    <dc:publisher/>
    <cmeta:comment rdf:resource="rdf:#52b20981-fb3e-431a-be2c-d0abd9b365f6"/>
    <dcterms:created rdf:resource="rdf:#b0fa4fcf-9058-49b1-ac1e-51d5308f83df"/>
    <dc:creator rdf:resource="rdf:#2535490c-7843-4eb9-a433-ee7226b78fbe"/>
    <cmeta:modification rdf:resource="rdf:#c5b9a0b9-6e00-4220-97fc-89679f4b3e98"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#ef395efd-5ff2-458e-a30e-e592db663345">
    <vCard:FN/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#1f9ded7d-a130-4df8-afbf-2b38b71f8826">
    <dcterms:W3CDTF>2008-05-14T01:30:59+12:00</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#05739223-6287-4bc4-b98b-4191b106c2ae">
    <vCard:FN>Penny Noble</vCard:FN>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#a1a005fa-0d09-4031-b9f0-91e01567551e">
    <dc:creator rdf:resource="rdf:#ef395efd-5ff2-458e-a30e-e592db663345"/>
    <rdf:value/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#250e831d-159c-4178-8bfd-ab2851fbaab4">
    <vCard:Given>Penny</vCard:Given>
    <vCard:Family>Noble</vCard:Family>
    <vCard:Other/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#4a50ef2c-c23a-4166-8930-916a286d79e5">
    <vCard:Given>Y</vCard:Given>
    <vCard:Family>Rudy</vCard:Family>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#2535490c-7843-4eb9-a433-ee7226b78fbe">
    <vCard:ORG rdf:resource="rdf:#d9a861c5-da9e-48cf-bf39-fed2ac4822c1"/>
    <vCard:EMAIL rdf:resource="rdf:#1a939bbc-bfbc-4823-a844-d4f33cb394ca"/>
    <vCard:N rdf:resource="rdf:#63f4146b-264b-482d-bff4-38ab536b8c57"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#4fbe3ba6-328e-4c4d-9249-8a2a965ecaea">
    <bqs:Pubmed_id>10777735</bqs:Pubmed_id>
    <bqs:JournalArticle rdf:resource="rdf:#dd1796f4-c286-4c30-9e0d-ba79d2f0e2d2"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#e7b77551-d534-4391-8f1a-73f244372afd">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#4a50ef2c-c23a-4166-8930-916a286d79e5"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#dd1796f4-c286-4c30-9e0d-ba79d2f0e2d2">
    <dc:creator rdf:resource="rdf:#9b6b12f4-7c11-4a19-a5d0-b72a173411ff"/>
    <dc:title>Action potential and contractility changes in [Na(+)](i) overloaded cardiac myocytes: a simulation study.</dc:title>
    <bqs:volume>78</bqs:volume>
    <bqs:first_page>2392</bqs:first_page>
    <bqs:Journal rdf:resource="rdf:#7bbc1652-a7fd-41eb-a3f0-7900d6d7ef93"/>
    <dcterms:issued rdf:resource="rdf:#6080d8cb-d870-4d50-90b8-9ae31af5b294"/>
    <bqs:last_page>2404</bqs:last_page>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#63f4146b-264b-482d-bff4-38ab536b8c57">
    <vCard:Given>Penny</vCard:Given>
    <vCard:Family>Noble</vCard:Family>
    <vCard:Other/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#b0fa4fcf-9058-49b1-ac1e-51d5308f83df">
    <dcterms:W3CDTF>2008-05-14T00:00:00+00:00</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#6080d8cb-d870-4d50-90b8-9ae31af5b294">
    <dcterms:W3CDTF>2000-05-00 00:00</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#d9a861c5-da9e-48cf-bf39-fed2ac4822c1">
    <vCard:Orgname>University of Oxford</vCard:Orgname>
    <vCard:Orgunit/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#1a939bbc-bfbc-4823-a844-d4f33cb394ca">
    <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:#7bbc1652-a7fd-41eb-a3f0-7900d6d7ef93">
    <dc:title>Biophysics Journal</dc:title>
  </rdf:Description>
  <rdf:Description rdf:about="#faber_rudy_2000_version5">
    <bqs:reference rdf:resource="rdf:#4fbe3ba6-328e-4c4d-9249-8a2a965ecaea"/>
    <cmeta:comment rdf:resource="rdf:#a1a005fa-0d09-4031-b9f0-91e01567551e"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#7a861dea-0ca9-408e-b817-b6e4f5722634">
    <vCard:Given>GM</vCard:Given>
    <vCard:Family>Faber</vCard:Family>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#c5b9a0b9-6e00-4220-97fc-89679f4b3e98">
    <dcterms:modified rdf:resource="rdf:#1f9ded7d-a130-4df8-afbf-2b38b71f8826"/>
    <rdf:value>Curation of units and correction of some code. Still getting results for ICaL that differ from the publication although equations are correct.</rdf:value>
    <cmeta:modifier rdf:resource="rdf:#b2535049-2e20-4454-9c79-7b61b3fb2d79"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#dba32688-0867-4a2a-acfa-18cbf0261435">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#7a861dea-0ca9-408e-b817-b6e4f5722634"/>
  </rdf:Description>
</rdf:RDF>
</model>