Location: Faber, Rudy, 2000 @ 124c7736bff0 / faber_rudy_modified_version_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_modified_version_2000.cellml

<?xml version='1.0' encoding='utf-8'?>
<!--
This CellML file was generated on 31/07/2008 at 3:02:14 at p.m. 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#" name="LR_Dynamic_model_2000" cmeta:id="LR_Dynamic_model_2000">
<documentation xmlns="http://cellml.org/tmp-documentation">
<article>
  <articleinfo>
  <title>The Most Recent Verison of the Luo-Rudy Mammalian Ventricular Model II</title>
  <author>
    <firstname>Catherine</firstname>
          <surname>Lloyd</surname>
    <affiliation>
      <shortaffil>Auckland Bioengineering Institute</shortaffil>
    </affiliation>
  </author>
</articleinfo>
  <section id="sec_status">
    <title>Model Status</title>
    <para>
            This model has been curated and is known to run in PCEnv and COR to recreate the published results.  The units have been checked and are known to be consistent.  The CellML model has been based on a translation of the LRd07 Matlab code which can be downloaded from the Rudy lab website.  It is basically the Faber and Rudy 2000 model with the addition of a transient potassium current (Ito) and an updated Irel current.
          </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.  Please note that 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 Yoram Rudy (To view this website, follow this link: <ulink url="http://rudylab.wustl.edu/research/cell/methodology/cellmodels/LRd/history.htm"> http://rudylab.wustl.edu/research/cell/methodology/cellmodels/LRd/history.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>

<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 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_cell_diagram_LRd">
<mediaobject>
  <imageobject>
    <objectinfo>
      <title>cell diagram</title>
    </objectinfo>
    <imagedata fileref="faber_2000_new.png"/>
  </imageobject>
</mediaobject>
<caption>In this most recent version of the 1994 Luo-Rudy Dynamic Ventricular Model, the model has been extended to include a transient outward potassium current (I<subscript>to</subscript>) and also the eqautions for the I<subscript>rel</subscript> current have been updated.</caption>
</informalfigure>

<informalfigure float="0" id="state_diagram">
<mediaobject>
  <imageobject>
    <objectinfo>
      <title>state diagram</title>
    </objectinfo>
    <imagedata fileref="livshitz_rudy_2007.png"/>
  </imageobject>
</mediaobject>
<caption>In this model, Livshitz and Rudy have merged the Luo-Rudy dynamic (LRd) and Hund-Rudy dynamic
(HRd) models of mammalian ventricular AP.  This is a schematic diagram of the currents and the ion exchanges described by this merged model.</caption>
</informalfigure>

<para>
All the original paper references are cited below.  Note that the CellML model described here is not based on a single publication per se.  Rather, it is mainly based upon the 2000 Faber and Rudy model, with the addition of a transient outward potassium current (I<subscript>to</subscript>) and also the eqautions for the I<subscript>rel</subscript> current have been updated.  These modifications were based on the LRd07 Matlab code, which can be downloaded from the <ulink url="http://rudylab.wustl.edu/research/cell/methodology/cellmodels/HRd/HRd07.html">Rudy Research Lab website</ulink>
</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>

<para>
<ulink url="http://ajpheart.physiology.org/cgi/content/abstract/292/6/H2854">Regulation of Ca2+ and electrical alternans in cardiac myocytes: role of CAMKII and repolarizing currents</ulink>, Leonid M. Livshitz and Yoram Rudy, 2007, <ulink url="http://ajpheart.physiology.org/">
            <emphasis>American Journal of Physiology: Heart and Circulatory Physiology</emphasis>
          </ulink>, 292, H2854-H2866. (<ulink url="http://ajpheart.physiology.org/cgi/content/full/292/6/H2854"> Full text</ulink> and <ulink url="http://ajpheart.physiology.org/cgi/reprint/292/6/H2854">PDF</ulink> versions of this article are available to subscribers on the <emphasis>American Journal of Physiology: Heart and Circulatory Physiology</emphasis> website).  <ulink url="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?cmd=Retrieve&amp;db=PubMed&amp;list_uids=17277017&amp;dopt=Abstract">PubMed ID: 17277017</ulink>
</para>

</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="mM_per_millivolt_ms">
      <unit units="mole" prefix="milli"/>
      <unit units="litre" exponent="-1"/>
      <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_per_mV">
      <unit units="mole" prefix="milli"/>
      <unit units="litre" exponent="-1"/>
      <unit units="volt" prefix="milli" exponent="-1"/>
   </units>
   <units name="millimolar">
      <unit units="mole" prefix="milli"/>
      <unit units="litre" exponent="-1"/>
   </units>
   <units name="millimolar2">
      <unit units="millimolar" exponent="2"/>
   </units>
   <units name="millimolar3">
      <unit units="millimolar" exponent="3"/>
   </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_NaCa"/>
      <variable units="microA_per_microF" public_interface="in" name="i_K1"/>
      <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_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="10"/>
      <variable units="ms" name="stim_end" initial_value="9e6"/>
      <variable units="ms" name="stim_period" initial_value="300"/>
      <variable units="ms" name="stim_duration" initial_value="0.5"/>
      <variable units="microA_per_microF" name="stim_amplitude" initial_value="-80"/>
      <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_K1</ci>
                  <ci>i_K_ATP</ci>
                  <ci>i_to</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="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>
            <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>
         <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>
            <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>
         </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">32</cn>
                           </apply>
                           <cn cellml:units="millivolt">8</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="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="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="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>
                  <divide/>
                  <apply>
                     <times/>
                     <ci>I_NaK</ci>
                     <ci>f_NaK</ci>
                  </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>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <divide/>
                     <ci>K_mKo</ci>
                     <ci>Ko</ci>
                  </apply>
               </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="in" 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="ms" name="tau_tr" initial_value="120"/>
      <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="CaT" initial_value="0.01544711"/>
      <variable units="millimolar" public_interface="out" name="Cao" initial_value="1.8"/>
      <variable units="millimolar" public_interface="out" name="Ca_JSR"/>
      <variable units="millimolar" public_interface="out" name="JSR" initial_value="8.290468"/>
      <variable units="millimolar" public_interface="out" name="NSR" initial_value="1.516756041281"/>
      <variable units="millimolar" public_interface="out" name="Cai"/>
      <variable units="millimolar" public_interface="out" name="bjsr"/>
      <variable units="millimolar2" public_interface="out" name="cjsr"/>
      <variable units="millimolar" public_interface="out" name="bmyo"/>
      <variable units="millimolar2" public_interface="out" name="cmyo"/>
      <variable units="millimolar3" public_interface="out" name="dmyo"/>
      <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_mTRPN" initial_value="0.0005"/>
      <variable units="millimolar" name="K_mCMDN" initial_value="0.00238"/>
      <variable units="millimolar" name="TRPN_max" initial_value="0.07"/>
      <variable units="millimolar" name="CMDN_max" initial_value="0.05"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>bjsr</ci>
            <apply>
               <plus/>
               <ci>CSQN_max</ci>
               <ci>K_mCSQN</ci>
               <apply>
                  <minus/>
                  <ci>JSR</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>cjsr</ci>
            <apply>
               <times/>
               <ci>K_mCSQN</ci>
               <ci>JSR</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>Ca_JSR</ci>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <apply>
                     <root/>
                     <apply>
                        <plus/>
                        <apply>
                           <power/>
                           <ci>bjsr</ci>
                           <cn cellml:units="dimensionless">2</cn>
                        </apply>
                        <apply>
                           <times/>
                           <cn cellml:units="dimensionless">4</cn>
                           <ci>cjsr</ci>
                        </apply>
                     </apply>
                  </apply>
                  <ci>bjsr</ci>
               </apply>
               <cn cellml:units="dimensionless">2</cn>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>bmyo</ci>
            <apply>
               <plus/>
               <ci>CMDN_max</ci>
               <ci>TRPN_max</ci>
               <ci>K_mCMDN</ci>
               <ci>K_mTRPN</ci>
               <apply>
                  <minus/>
                  <ci>CaT</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>cmyo</ci>
            <apply>
               <plus/>
               <apply>
                  <times/>
                  <ci>K_mCMDN</ci>
                  <ci>K_mTRPN</ci>
               </apply>
               <apply>
                  <times/>
                  <ci>TRPN_max</ci>
                  <ci>K_mCMDN</ci>
               </apply>
               <apply>
                  <times/>
                  <ci>CMDN_max</ci>
                  <ci>K_mTRPN</ci>
               </apply>
               <apply>
                  <times/>
                  <apply>
                     <minus/>
                     <ci>CaT</ci>
                  </apply>
                  <apply>
                     <plus/>
                     <ci>K_mTRPN</ci>
                     <ci>K_mCMDN</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>dmyo</ci>
            <apply>
               <times/>
               <apply>
                  <minus/>
                  <ci>K_mTRPN</ci>
               </apply>
               <ci>K_mCMDN</ci>
               <ci>CaT</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>Cai</ci>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <apply>
                     <divide/>
                     <cn cellml:units="dimensionless">2</cn>
                     <cn cellml:units="dimensionless">3</cn>
                  </apply>
                  <apply>
                     <root/>
                     <apply>
                        <minus/>
                        <apply>
                           <power/>
                           <ci>bmyo</ci>
                           <cn cellml:units="dimensionless">2</cn>
                        </apply>
                        <apply>
                           <times/>
                           <cn cellml:units="dimensionless">3</cn>
                           <ci>cmyo</ci>
                        </apply>
                     </apply>
                  </apply>
                  <apply>
                     <cos/>
                     <apply>
                        <divide/>
                        <apply>
                           <arccos/>
                           <apply>
                              <divide/>
                              <apply>
                                 <minus/>
                                 <apply>
                                    <times/>
                                    <cn cellml:units="dimensionless">9</cn>
                                    <ci>bmyo</ci>
                                    <ci>cmyo</ci>
                                 </apply>
                                 <apply>
                                    <plus/>
                                    <apply>
                                       <times/>
                                       <cn cellml:units="dimensionless">2</cn>
                                       <apply>
                                          <power/>
                                          <ci>bmyo</ci>
                                          <cn cellml:units="dimensionless">3</cn>
                                       </apply>
                                    </apply>
                                    <apply>
                                       <times/>
                                       <cn cellml:units="dimensionless">27</cn>
                                       <ci>dmyo</ci>
                                    </apply>
                                 </apply>
                              </apply>
                              <apply>
                                 <times/>
                                 <cn cellml:units="dimensionless">2</cn>
                                 <apply>
                                    <root/>
                                    <apply>
                                       <power/>
                                       <apply>
                                          <minus/>
                                          <apply>
                                             <power/>
                                             <ci>bmyo</ci>
                                             <cn cellml:units="dimensionless">2</cn>
                                          </apply>
                                          <apply>
                                             <times/>
                                             <cn cellml:units="dimensionless">3</cn>
                                             <ci>cmyo</ci>
                                          </apply>
                                       </apply>
                                       <cn cellml:units="dimensionless">3</cn>
                                    </apply>
                                 </apply>
                              </apply>
                           </apply>
                        </apply>
                        <cn cellml:units="dimensionless">3</cn>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <divide/>
                  <ci>bmyo</ci>
                  <cn cellml:units="dimensionless">3</cn>
               </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>NSR</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>i_tr</ci>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>NSR</ci>
                  <ci>Ca_JSR</ci>
               </apply>
               <ci>tau_tr</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>JSR</ci>
            </apply>
            <apply>
               <minus/>
               <ci>i_tr</ci>
               <ci>i_rel</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>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>CaT</ci>
            </apply>
            <apply>
               <plus/>
               <apply>
                  <divide/>
                  <apply>
                     <times/>
                     <apply>
                        <minus/>
                        <cn cellml:units="microF_per_cm2">1</cn>
                     </apply>
                     <ci>A_cap</ci>
                     <apply>
                        <plus/>
                        <ci>i_CaCa</ci>
                        <ci>i_Ca_T</ci>
                        <ci>i_p_Ca</ci>
                        <ci>i_Ca_b</ci>
                        <apply>
                           <times/>
                           <apply>
                              <minus/>
                              <cn cellml:units="dimensionless">2</cn>
                           </apply>
                           <ci>i_NaCa</ci>
                        </apply>
                     </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>
      </math>
   </component>
   <component name="CICR_current">
      <variable units="mM_per_millivolt_ms" name="kappa" initial_value="0.125"/>
      <variable units="ms" name="tau" initial_value="4.75"/>
      <variable units="millimolar" name="K_relss" initial_value="1"/>
      <variable units="dimensionless" name="qn" initial_value="9"/>
      <variable units="millimolar_per_mV" public_interface="out" name="alpha_rel"/>
      <variable units="ms" public_interface="out" name="tau_rel"/>
      <variable units="millimolar_per_ms" public_interface="out" name="I_relss"/>
      <variable units="millimolar_per_ms" public_interface="out" name="i_rel" initial_value="-0.000000000000000000474565"/>
      <variable units="microA_per_microF" public_interface="in" name="i_CaCa"/>
      <variable units="millimolar" public_interface="in" name="Ca_JSR"/>
      <variable units="ms" public_interface="in" name="time"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>alpha_rel</ci>
            <apply>
               <times/>
               <ci>kappa</ci>
               <ci>tau</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>I_relss</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <ci>i_CaCa</ci>
                  <ci>alpha_rel</ci>
               </apply>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <power/>
                     <apply>
                        <divide/>
                        <ci>K_relss</ci>
                        <ci>Ca_JSR</ci>
                     </apply>
                     <ci>qn</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>tau_rel</ci>
            <apply>
               <divide/>
               <ci>tau</ci>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <divide/>
                     <cn cellml:units="millimolar">0.0123</cn>
                     <ci>Ca_JSR</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>i_rel</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <apply>
                     <plus/>
                     <ci>I_relss</ci>
                     <ci>i_rel</ci>
                  </apply>
               </apply>
               <ci>tau_rel</ci>
            </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_ns_K"/>
      <variable units="microA_per_microF" public_interface="in" name="I_st"/>
      <variable units="microA_per_microF" public_interface="in" name="i_K_ATP"/>
      <variable units="microA_per_microF" public_interface="in" name="i_to"/>
	  <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_K_ATP</ci>
                     <ci>i_to</ci>
					 <ci>i_Kp</ci>
                     <ci>i_K_Na</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="transient_outward_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="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="environment" component_1="CICR_current"/>
      <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="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="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="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"/>
   </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="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="transient_outward_current"/>
      <map_variables variable_2="i_to" variable_1="i_to"/>
   </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="transient_outward_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="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>
   <connection>
      <map_components component_2="calcium_dynamics" component_1="CICR_current"/>
      <map_variables variable_2="i_rel" variable_1="i_rel"/>
      <map_variables variable_2="Ca_JSR" variable_1="Ca_JSR"/>
   </connection>
   <connection>
      <map_components component_2="L_type_Ca_channel" component_1="CICR_current"/>
      <map_variables variable_2="i_CaCa" variable_1="i_CaCa"/>
   </connection>
<rdf:RDF>
  <rdf:Seq rdf:about="rdf:#f0c2ed0f-950b-457a-87f5-16f24a7b5a90">
    <rdf:li rdf:resource="rdf:#cd2bb52f-284f-43de-b024-9fe10f80ab2e"/>
    <rdf:li rdf:resource="rdf:#b743931e-2d65-4bab-96f0-5ec40bdc4366"/>
  </rdf:Seq>
  <rdf:Description rdf:about="rdf:#069cbe50-344f-438d-9321-954372ce4388">
    <vCard:FN/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#324dc3d7-c193-43e3-81d5-dfae211d50f8">
    <vCard:Orgname>The University of Auckland</vCard:Orgname>
    <vCard:Orgunit>Auckland Bioengineering Institute</vCard:Orgunit>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#bda92893-302a-4ffd-86d4-d78d1a5d07dd">
    <vCard:FN>Catherine Lloyd</vCard:FN>
  </rdf:Description>
  <rdf:Description rdf:about="">
    <dc:publisher>Auckland Bioengineering Institute, The University of Auckland</dc:publisher>
    <cmeta:comment rdf:resource="rdf:#7578c32c-b0c8-442a-aa5a-7a4d90184272"/>
    <dcterms:created rdf:resource="rdf:#32666255-7f00-47c5-a4c5-eae91eab1c97"/>
    <dc:creator rdf:resource="rdf:#4ccd9d28-15f4-4041-ab5e-e815fc7bace8"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#7578c32c-b0c8-442a-aa5a-7a4d90184272">
    <dc:creator rdf:resource="rdf:#bda92893-302a-4ffd-86d4-d78d1a5d07dd"/>
    <rdf:value>This model has been curated and is known to run in PCEnv and COR to recreate the published results.  The units have been checked and are known to be consistent.&#13;
&#13;
The CellML model has been based on a translation of the LRd07 Matlab code which can be downloaded from the Rudy lab website.  It is basically the Faber and Rudy 2000 model with the addition of a transient potassium current (Ito) and an updated Irel current.</rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#914db032-065c-4620-a03c-e67de44fc090">
    <rdf:type rdf:resource="http://imc.org/vCard/3.0#internet"/>
    <rdf:value>j.ashton@auckland.ac.nz</rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="#LR_Dynamic_model_2000">
    <bqs:reference rdf:resource="rdf:#a8b27314-2526-4722-af9b-ecaf4cb46854"/>
<bqs:reference rdf:parseType="Resource">
  <dc:subject rdf:parseType="Resource">
    <bqs:subject_type>keyword</bqs:subject_type>
    <rdf:value>
      <rdf:Bag>
        <rdf:li>excitation-contraction coupling</rdf:li>
        <rdf:li>Ventricular Myocyte</rdf:li>
        <rdf:li>calcium</rdf:li>
        <rdf:li>electrophysiology</rdf:li>
        <rdf:li>cardiac</rdf:li>
        <rdf:li>pathology</rdf:li>
        <rdf:li>ventricular myocyte</rdf:li>
      </rdf:Bag>
    </rdf:value>
  </dc:subject>
</bqs:reference>
    <cmeta:comment rdf:resource="rdf:#8fde73b5-5891-4fc8-bfef-b3c0eb7a750a"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#cd2bb52f-284f-43de-b024-9fe10f80ab2e">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#aeccb117-b108-4dbf-b3a7-f82bd484e701"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#84df8fca-47ca-4db4-b632-f3f72250b527">
    <vCard:Given>Y</vCard:Given>
    <vCard:Family>Rudy</vCard:Family>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#9da9f3b0-68a4-45b9-a70f-173c3f89c5da">
    <vCard:Given>Jesse</vCard:Given>
    <vCard:Family>Ashton</vCard:Family>
    <vCard:Other/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#8fde73b5-5891-4fc8-bfef-b3c0eb7a750a">
    <dc:creator rdf:resource="rdf:#069cbe50-344f-438d-9321-954372ce4388"/>
    <rdf:value/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#c0125ce0-0f7e-44b8-a0b7-e6471cc6a42d">
    <dc:creator rdf:resource="rdf:#f0c2ed0f-950b-457a-87f5-16f24a7b5a90"/>
    <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:#71c978ef-b44c-4ca7-8ff1-369b1524a64a"/>
    <dcterms:issued rdf:resource="rdf:#e27e0947-cc2c-443f-86a3-6a6ea4b3b563"/>
    <bqs:last_page>2404</bqs:last_page>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#e27e0947-cc2c-443f-86a3-6a6ea4b3b563">
    <dcterms:W3CDTF>2000-05-00 00:00</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#71c978ef-b44c-4ca7-8ff1-369b1524a64a">
    <dc:title>Biophysical Journal</dc:title>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#32666255-7f00-47c5-a4c5-eae91eab1c97">
    <dcterms:W3CDTF>2008-08-01T00:00:00+00:00</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#a8b27314-2526-4722-af9b-ecaf4cb46854">
    <bqs:Pubmed_id>10777735</bqs:Pubmed_id>
    <bqs:JournalArticle rdf:resource="rdf:#c0125ce0-0f7e-44b8-a0b7-e6471cc6a42d"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#b743931e-2d65-4bab-96f0-5ec40bdc4366">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#84df8fca-47ca-4db4-b632-f3f72250b527"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#aeccb117-b108-4dbf-b3a7-f82bd484e701">
    <vCard:Given>G</vCard:Given>
    <vCard:Family>Faber</vCard:Family>
    <vCard:Other>M</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#4ccd9d28-15f4-4041-ab5e-e815fc7bace8">
    <vCard:ORG rdf:resource="rdf:#324dc3d7-c193-43e3-81d5-dfae211d50f8"/>
    <vCard:EMAIL rdf:resource="rdf:#914db032-065c-4620-a03c-e67de44fc090"/>
    <vCard:N rdf:resource="rdf:#9da9f3b0-68a4-45b9-a70f-173c3f89c5da"/>
  </rdf:Description>
</rdf:RDF>
</model>