Location: Rice, Jafri, Winslow, 2000 @ 4533913da7a3 / rice_jafri_winslow_2000.cellml

Author:
pmr2.import <nobody@models.cellml.org>
Date:
2009-06-17 15:43:04+12:00
Desc:
committing version02 of rice_jafri_winslow_2000
Permanent Source URI:
https://models.cellml.org/workspace/rice_jafri_winslow_2000/rawfile/4533913da7a3d2c499b34891dc5df5d240f99ab1/rice_jafri_winslow_2000.cellml

<?xml version='1.0' encoding='utf-8'?>
<!--  FILE :  rice_model2_2000_raw.xml

CREATED :  24th June 2002

LAST MODIFIED : 9th April 2003

AUTHOR :  Catherine Lloyd
          The Bioengineering Institute
          The University of Auckland
          
MODEL STATUS :  This model conforms to the CellML 1.0 Specification released on
                10th August 2001, and the 16/01/2002 CellML Metadata 1.0
                Specification.

DESCRIPTION :  This file contains a CellML description of Rice et al's 2000
               detailed model of ventricular cardiac myocytes.  It simulates
               action potentials, Ca2+-handling mechanisms and isometric force
               generation by the myofilaments.

CHANGES:  
  18/07/2002 - CML - Added more metadata.
  09/04/2003 - AAC - Added publication date information.  
--><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:pathway_editor="http://www.physiome.com/pathway_editor/1.0#" xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" pathway_editor:rendering_config_file="rice_model2_2000_CellMLrender.xml" cmeta:id="rice_jafri_winslow_2000_version02" name="rice_jafri_winslow_2000_version02">
<documentation xmlns="http://cellml.org/tmp-documentation">
<article>
  <articleinfo>
  <title>Modelling Interval-Force Relations in Cardiac Muscle</title>
  <author>
    <firstname>Catherine</firstname>
          <surname>Lloyd</surname>
    <affiliation>
      <shortaffil>Bioengineering Institute, University of Auckland</shortaffil>
    </affiliation>
  </author>
</articleinfo>
  <section id="sec_status">
    <title>Model Status</title>
    <para>
            This is the original unchecked version of the model imported from the previous
            CellML model repository, 24-Jan-2006.
          </para>
  </section>
  <sect1 id="sec_structure">
<title>Model Structure</title>

<para>
Short-term interval-force (I-F) relations describe the dependence of cardiac muscle contraction strength for short interbeat intervals.  Two phenomena that characterise these short-term I-F relations are restitution and postextrasystolic potentiation.  In their 2000 paper, J. Jeremy Rice, M. Saleet Jafri and Raimond L. Winslow use two modelling approaches to study short-term interval-force relations in cardiac muscle.  Their first approach employs a relatively simple (or minimal), discrete-time model of excitation-contraction coupling.  This model explains restitution and postextrasystolic potentiation in terms of the total amount of Ca<superscript>2+</superscript> stored in the sarcoplasmic reticulum (SR) and a feature in which Ca<superscript>2+</superscript> stored in the SR slowly becomes available for release into the sarcoplasm (see <xref linkend="fig_simple_diagram"/> below).
</para>

<informalfigure float="0" id="fig_simple_diagram">
          <mediaobject>
            <imageobject>
              <objectinfo>
                <title>simple_model</title>
              </objectinfo>
              <imagedata fileref="rice_2000a.png"/>
            </imageobject>
            <caption>
Schematic diagram of a simple, discrete-time model of Ca<superscript>2+</superscript>-handling in restitution and postextrasystolic potentiation. In<subscript>SC</subscript> and Out<subscript>SC</subscript> represent Ca<superscript>2+</superscript> influx and efflux across the sarcolemma, and In<subscript>SR</subscript> and Out<subscript>SR</subscript> represent sarcoplasmic reticulum Ca<superscript>2+</superscript> release and uptake. 
</caption>
          </mediaobject>
        </informalfigure>

<para>
In their second modelling approach, J. Jeremy Rice, M. Saleet Jafri and Raimond L. Winslow develop a more detailed ventricular cell model that is able to simulate action potentials, Ca<superscript>2+</superscript>-handling mechanisms, and isometric force generation by the myofilaments.  This single comprehensive (or maximal) model combines two previously developed models (see <ulink url="${HTML_EXMPL_JRW_MODEL}"> The Jafri-Rice-Winslow Ventricular Model, 1998</ulink> and <ulink url="${HTML_EXMPL_RICE_MODEL}"> Cooperative Mechanisms in Cardiac Muscle, Rice <emphasis>et al.</emphasis>, 1999</ulink>).  The first ventricular cell model describes membrane currents and Ca<superscript>2+</superscript>-handling, and the second model of myofilaments provides the force generation of muscle contraction (see <xref linkend="fig_detailed_diagram"/> below).  The ventricular cell model provides the Ca<superscript>2+</superscript> needed to drive the myofilament contraction.  This is a feedforward pathway in which Ca<superscript>2+</superscript> binds to troponin and initiates force generation by triggering myofilament contraction.  A feedback pathway also links the two models, as the affinity of troponin for Ca<superscript>2+</superscript> is a function of developed force. 
</para>

<para>
<ulink url="http://ajpheart.physiology.org/cgi/content/abstract/278/3/H913">Modeling short-term interval-force relations in cardiac muscle</ulink> J. Jeremy Rice, M. Saleet Jafri and Raimond L. Winslow, 2000, <ulink url="http://ajpheart.physiology.org/">
            <emphasis>American Journal of Physiology</emphasis>
          </ulink>, 278, H913-H931.  (<ulink url="http://ajpheart.physiology.org/cgi/content/full/278/3/H913">Full text</ulink> and <ulink url="http://ajpheart.physiology.org/cgi/reprint/278/3/H913.pdf">PDF</ulink> versions of the article are available for Journal Members on the American Journal of Physiology website.)  <ulink url="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?cmd=Retrieve&amp;db=PubMed&amp;list_uids=10710361&amp;dopt=Abstract">PubMed ID: 10710361</ulink> 
</para>

<informalfigure float="0" id="fig_detailed_diagram">
<mediaobject>
  <imageobject>
    <objectinfo>
      <title>detailed_model</title>
    </objectinfo>
    <imagedata fileref="rice_2000b.png"/>
  </imageobject>
</mediaobject>
<caption>Schematic diagram of a detailed, single ventricular cell model which is able to simulate action potentials, Ca<superscript>2+</superscript>-handling and isometric force generation by myofilaments.</caption>
</informalfigure>

<para>
Simulation results generated from the model support experimental data.  They suggest that short-term interval-force relations result mainly from the interplay between ryanodine receptor adaptation and SR Ca<superscript>2+</superscript> loading with additional contributions from membrane currents and myofilament activation.
</para>

</sect1>
</article>
</documentation>
  
  
  <!--
    Below, we define some additional units for association with variables and
    constants within the model. The identifiers are fairly self-explanatory.
  -->
  <units name="millisecond">
    <unit units="second" prefix="milli"/>
  </units>
  
  <units name="per_millisecond">
    <unit units="second" prefix="milli" exponent="-1"/>
  </units>
  
  <units name="per_second">
    <unit units="second" exponent="-1"/>
  </units>
  
  <units name="micro_litre">
    <unit units="litre" prefix="micro"/>
  </units>
  
  <units name="micrometre">
    <unit units="metre" prefix="micro"/>
  </units>
  
  <units name="cm2">
    <unit units="metre" prefix="centi" exponent="2"/>
  </units>
  
  <units name="cm_per_second">
    <unit units="metre" prefix="centi"/>
    <unit units="second" exponent="-1"/>
  </units>
  
  <units name="cm_per_millisecond">
    <unit units="metre" prefix="centi"/>
    <unit units="millisecond" exponent="-1"/>
  </units>
  
  <units name="millivolt">
    <unit units="volt" prefix="milli"/>
  </units>
  
  <units name="per_millivolt">
    <unit units="millivolt" exponent="-1"/>
  </units>
  
  <units name="per_millivolt_millisecond">
    <unit units="millivolt" exponent="-1"/>
    <unit units="millisecond" exponent="-1"/>
  </units>
  
  <units name="millimolar">
    <unit units="mole" prefix="milli"/>
    <unit units="litre" exponent="-1"/>
  </units>
  
  <units name="micromolar">
    <unit units="mole" prefix="micro"/>
    <unit units="litre" exponent="-1"/>
  </units>
  
  <units name="micromolar_per_second">
    <unit units="micromolar"/>
    <unit units="second" exponent="-1"/>
  </units>
  
  <units name="micromolar_per_millisecond">
    <unit units="micromolar"/>
    <unit units="millisecond" exponent="-1"/>
  </units>
  
  <units name="per_micromolar_per_second">
    <unit units="micromolar" exponent="-1"/>
    <unit units="second" exponent="-1"/>
  </units>
  
  <units name="per_micromolar3_per_millisecond">
    <unit units="micromolar" exponent="-3"/>
    <unit units="millisecond" exponent="-1"/>
  </units>
  
  <units name="per_micromolar4_per_millisecond">
    <unit units="micromolar" exponent="-4"/>
    <unit units="millisecond" exponent="-1"/>
  </units>
  
  <units name="microF_per_cm2">
    <unit units="farad" prefix="micro"/>
    <unit units="metre" prefix="centi" exponent="-2"/>
  </units>

  <units name="microA_per_microF">
    <unit units="ampere" prefix="micro"/>
    <unit units="farad" prefix="micro" exponent="-1"/>
  </units>
  
  <units name="milliS_per_microF">
    <unit units="siemens" prefix="milli"/>
    <unit units="farad" prefix="micro" exponent="-1"/>
  </units>
  
  <units name="joule_per_mole_kelvin">
    <unit units="joule"/>
    <unit units="mole" exponent="-1"/>
    <unit units="kelvin" exponent="-1"/>
  </units>
  
  <units name="coulomb_per_millimole">
    <unit units="coulomb"/>
    <unit units="mole" prefix="milli" exponent="-1"/>
  </units>
  
  <units name="N_per_mm2">
    <unit units="newton"/>
    <unit units="metre" prefix="milli" exponent="-2"/>
  </units>
  
  <!--
    The "environment" component is used to declare variables that are used by
    all or most of the other components, in this case just "time".
  -->
  
  <component name="environment">
    <variable units="millisecond" public_interface="out" name="time"/>
  </component>
  
  <component name="cell_geometry">
    <variable units="cm2" public_interface="out" name="A_cap" initial_value="1.534E-4"/>
    <variable units="micro_litre" public_interface="out" name="V_myo" initial_value="25.84E-6"/>
    <variable units="micro_litre" public_interface="out" name="V_JSR" initial_value="0.16E-6"/>
    <variable units="micro_litre" public_interface="out" name="V_NSR" initial_value="1.40E-6"/>
    <variable units="micro_litre" public_interface="out" name="V_SS" initial_value="0.495E-6"/>
  </component>
  
  <!--
    The "membrane" component is really the `root' node of our model.
    It defines the action potential variable "V" among other things.
  -->
  <component name="membrane">
    <variable units="millivolt" public_interface="out" name="V" initial_value="-84.1638"/>
    <variable units="joule_per_mole_kelvin" public_interface="out" name="R" initial_value="8.314"/>
    <variable units="kelvin" public_interface="out" name="T" initial_value="310.0"/>
    <variable units="coulomb_per_millimole" public_interface="out" name="F" initial_value="96.5"/>
   
    <variable units="microF_per_cm2" name="Cm" initial_value="1.0"/>    
   
    <variable units="millisecond" 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"/>
    <variable units="microA_per_microF" public_interface="in" name="i_Ca_K"/>
    <variable units="microA_per_microF" public_interface="in" name="i_K"/>
    <variable units="microA_per_microF" public_interface="in" name="i_K1"/>
    <variable units="microA_per_microF" public_interface="in" name="i_NaCa"/>
    <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"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="membrane_voltage_diff_eq">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> V </ci>
        </apply>
        <apply>
          <times/>
          <apply>
            <divide/>
            <cn cellml:units="dimensionless"> -1.0 </cn>
            <ci> Cm </ci>
          </apply>
          <apply>
            <plus/>
            <ci> i_Na </ci>
            <ci> i_Ca </ci>
            <ci> i_Ca_K </ci>
            <ci> i_K </ci>
            <ci> i_K1 </ci>
            <ci> i_Kp </ci>
            <ci> i_NaCa </ci>
            <ci> i_NaK </ci>
            <ci> i_ns_Ca </ci>
            <ci> i_p_Ca </ci>
            <ci> i_Ca_b </ci>
            <ci> i_Na_b </ci>
          </apply>
        </apply>
      </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="12.8"/>
    
    <variable units="millisecond" public_interface="in" private_interface="out" name="time"/>
    <variable units="millivolt" public_interface="in" private_interface="out" name="V"/>
    <variable units="joule_per_mole_kelvin" public_interface="in" name="R"/>
    <variable units="kelvin" public_interface="in" name="T"/>
    <variable units="coulomb_per_millimole" public_interface="in" name="F"/>
    <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 id="i_Na_calculation">
        <eq/>
        <ci> i_Na </ci>
        <apply>
          <times/>
          <ci> g_Na </ci>
          <apply>
            <power/>
            <ci> m </ci>
            <cn cellml:units="dimensionless"> 3.0 </cn>
          </apply>
          <ci> h </ci>
          <ci> j </ci>
          <apply>
            <minus/>
            <ci> V </ci>
            <ci> E_Na </ci>
          </apply>
        </apply>
      </apply> 
      
      <apply id="E_Na_calculation">
        <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>
    </math>
  </component>
  
  <component name="fast_sodium_current_m_gate">
    <variable units="dimensionless" public_interface="out" name="m" initial_value="0.0326174"/>
        
    <variable units="per_millisecond" name="alpha_m"/>
    <variable units="per_millisecond" name="beta_m"/>

    <variable units="millivolt" public_interface="in" name="V"/>
    <variable units="millisecond" public_interface="in" name="time"/>
   
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="alpha_m_calculation">
        <eq/>
        <ci> alpha_m </ci>
        <apply>
          <times/>
          <cn cellml:units="per_millivolt_millisecond"> 0.32 </cn>
          <apply>
            <divide/>
            <apply>
              <plus/>
              <ci> V </ci>
              <cn cellml:units="millivolt"> 47.13 </cn>
            </apply>
            <apply>
              <minus/>
              <cn cellml:units="dimensionless"> 1.0 </cn>
              <apply>
                <exp/>
                <apply>
                  <times/>
                  <cn cellml:units="per_millivolt"> -0.1 </cn>
                  <apply>
                    <plus/>
                    <ci> V </ci>
                    <cn cellml:units="millivolt"> 47.13 </cn>
                  </apply>
                </apply>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply id="beta_m_calculation">
        <eq/>
        <ci> beta_m </ci>
        <apply>
          <times/>
          <cn cellml:units="per_millisecond"> 0.08 </cn>
          <apply>
            <exp/>
            <apply>
              <divide/>
              <apply>
                <minus/>
                <ci> V </ci>
              </apply>
              <cn cellml:units="millivolt"> 11.0 </cn>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply id="dm_dt">
        <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="per_millivolt"> 1.0 </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.962369"/>
        
    <variable units="per_millisecond" name="alpha_h"/>
    <variable units="per_millisecond" name="beta_h"/>

    <variable units="millivolt" public_interface="in" name="V"/>
    <variable units="millisecond" public_interface="in" name="time"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="alpha_h_calculation">
        <eq/>
        <ci> alpha_h </ci>
        <piecewise>
          <piece>
            <apply>
              <times/>
              <cn cellml:units="per_millisecond"> 0.135 </cn>
              <apply>
                <exp/>
                <apply>
                  <divide/>
                  <apply>
                    <plus/>
                    <cn cellml:units="millivolt"> 80.0 </cn>
                    <ci> V </ci>
                  </apply>
                  <cn cellml:units="millivolt"> -6.8 </cn>
                </apply>
              </apply>
            </apply>
            <apply>
              <lt/>
              <ci> V </ci> 
              <cn cellml:units="millivolt"> -40.0 </cn>
            </apply>
          </piece>
          <otherwise>
            <cn cellml:units="per_millisecond"> 0.0 </cn>
          </otherwise>
        </piecewise>
      </apply>
       
       
      <apply id="beta_h_calculation">
        <eq/>
        <ci> beta_h </ci>
        <piecewise>
          <piece>
            <apply>
              <plus/>
              <apply>
                <times/>
                <cn cellml:units="per_millisecond"> 3.56 </cn>
                <apply>
                  <exp/>
                  <apply>
                    <times/>
                    <cn cellml:units="millivolt"> 0.079 </cn>
                    <ci> V </ci>
                  </apply>
                </apply>
              </apply>
              <apply>
                <times/>
                <cn cellml:units="per_millisecond"> 310000.0 </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>
              <cn cellml:units="millivolt"> -40.0 </cn>
            </apply>
          </piece>
          <otherwise>
            <apply>
              <divide/>
              <cn cellml:units="dimensionless"> 1.0 </cn>
              <apply>
                <times/>
                <cn cellml:units="millisecond"> 0.13 </cn>
                <apply>
                  <plus/>
                  <cn cellml:units="dimensionless"> 1.0 </cn>
                  <apply>
                    <exp/>
                    <apply>
                      <divide/>
                      <apply>
                        <minus/>
                        <apply>
                          <plus/>
                          <ci> V </ci>
                          <cn cellml:units="millivolt"> 10.66 </cn>
                        </apply>
                      </apply>
                      <cn cellml:units="millivolt"> 11.1 </cn>
                    </apply>
                  </apply>
                </apply>
              </apply>
            </apply>
          </otherwise>
        </piecewise>
      </apply>
      
      <apply id="dh_dt">
        <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.0 </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.975658"/>
        
    <variable units="per_millisecond" name="alpha_j"/>
    <variable units="per_millisecond" name="beta_j"/>
    
    <variable units="millivolt" public_interface="in" name="V"/>
    <variable units="millisecond" public_interface="in" name="time"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="alpha_j_calculation">
        <eq/>
        <ci> alpha_j </ci>
        <piecewise>
          <piece>
            <apply>
              <times/>
              <apply>
                <minus/>
                <apply>
                  <times/>
                  <cn cellml:units="per_millivolt_millisecond"> -127140.0
</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_millisecond"> 0.00003474 </cn>
                  <apply>
                    <exp/>
                    <apply>
                      <times/>
                      <cn cellml:units="per_millivolt"> -0.04391 </cn>
                      <ci> V </ci>
                    </apply>
                  </apply>
                </apply>
              </apply>
              <apply>
                <divide/>
                <apply>
                  <plus/>
                  <ci> V </ci>
                  <cn cellml:units="millivolt"> 37.78 </cn>
                </apply>
                <apply>
                  <plus/>
                  <cn cellml:units="dimensionless"> 1.0 </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>
            <apply>
              <lt/>
              <ci> V </ci>
              <cn cellml:units="millivolt"> -40.0 </cn>
            </apply>
          </piece>
          <otherwise>
            <cn cellml:units="per_millisecond"> 0.0 </cn>
          </otherwise>
        </piecewise>
      </apply>
       
      <apply id="beta_j_calculation">
        <eq/>
        <ci> beta_j </ci>
        <piecewise>
          <piece>
            <apply>
              <divide/>
              <apply>
                <times/>
                <cn cellml:units="per_millisecond"> 0.1212 </cn>
                <apply>
                  <exp/>
                  <apply>
                    <times/>
                    <cn cellml:units="per_millivolt"> -0.01052 </cn>
                    <ci> V </ci>
                  </apply>
                </apply>
              </apply>
              <apply>
                <plus/>
                <cn cellml:units="dimensionless"> 1.0 </cn>
                <apply>
                  <exp/>
                  <apply>
                    <times/>
                    <cn cellml:units="per_millivolt"> -0.1378 </cn>
                    <apply>
                      <plus/>       
                      <ci> V </ci>
                      <cn cellml:units="millivolt"> 40.14 </cn>
                    </apply>
                  </apply>
                </apply>
              </apply>
            </apply>
            <apply>
              <lt/>
              <ci> V </ci>
              <cn cellml:units="millivolt"> -40.0 </cn>
            </apply>
          </piece>
          <otherwise>
            <apply>
              <divide/>
              <apply>
                <times/>
                <cn cellml:units="per_millisecond"> 0.3 </cn>
                <apply>
                  <exp/>
                  <apply>
                    <times/>
                    <cn cellml:units="per_millivolt"> -0.0000002535 </cn>
                    <ci> V </ci>                                                
                  </apply>
                </apply>
              </apply>
              <apply>
                <plus/>
                <cn cellml:units="dimensionless"> 1.0 </cn>
                <apply>
                  <exp/>
                  <apply>
                    <times/>
                    <cn cellml:units="per_millivolt"> -0.1 </cn>
                    <apply>
                      <plus/>
                      <ci> V </ci>
                      <cn cellml:units="millivolt"> 32.0 </cn>
                    </apply>
                  </apply>
                </apply>
              </apply>
            </apply>                                                            
          </otherwise> 
        </piecewise>
      </apply>
      
      <apply id="dj_dt">
        <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.0 </cn>
              <ci> j </ci>
            </apply>
          </apply>
          <apply>
            <times/>
            <ci> beta_j </ci>
            <ci> j </ci>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="time_dependent_potassium_current">
    <variable units="microA_per_microF" public_interface="out" name="i_K"/>
    
    <variable units="milliS_per_microF" name="g_K"/>
    <variable units="millivolt" name="E_K"/>
    <variable units="cm_per_second" name="PNa_K" initial_value="0.01833"/>
    
    <variable units="millisecond" 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="millimolar" public_interface="in" name="Ki"/>
    <variable units="millimolar" public_interface="in" name="Nao"/>
    <variable units="millimolar" public_interface="in" name="Nai"/>
    <variable units="joule_per_mole_kelvin" public_interface="in" name="R"/>
    <variable units="kelvin" public_interface="in" name="T"/>
    <variable units="coulomb_per_millimole" public_interface="in" name="F"/>
    
    <variable units="dimensionless" private_interface="in" name="X"/>
    <variable units="dimensionless" private_interface="in" name="Xi"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="g_K_calculation">
        <eq/>
        <ci> g_K </ci>
        <apply>
          <times/>
          <cn cellml:units="milliS_per_microF"> 0.1128 </cn>
          <apply>
            <root/>
            <apply>
              <divide/>
              <ci> Ko </ci>
              <cn cellml:units="millimolar"> 5.4 </cn>
            </apply>
          </apply>                                                              
        </apply>
      </apply>
      
      <apply id="E_K_calculation">
        <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/>
              <apply>
                <plus/>
                <ci> Ko </ci>
                <apply>
                  <times/>
                  <ci> PNa_K </ci>
                  <ci> Nao </ci>
                </apply>
              </apply>
              <apply>
                <plus/>
                <ci> Ki </ci>
                <apply>
                  <times/>
                  <ci> PNa_K </ci>
                  <ci> Nai </ci>
                </apply>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply id="i_K_calculation">
        <eq/>
        <ci> i_K </ci>
        <apply>
          <times/>
          <ci> g_K </ci>
          <ci> Xi </ci>
          <apply>
            <power/>
            <ci> X </ci>
            <cn cellml:units="dimensionless"> 2.0 </cn>
          </apply>
          <apply>
            <minus/>
            <ci> V </ci>
            <ci> E_K </ci>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="time_dependent_potassium_current_X_gate">
    <variable units="dimensionless" public_interface="out" name="X" initial_value="0.000337175"/>
    
    <variable units="per_millisecond" name="alpha_X"/>
    <variable units="per_millisecond" name="beta_X"/>
    
    <variable units="millivolt" public_interface="in" name="V"/>
    <variable units="millisecond" public_interface="in" name="time"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="alpha_X_calculation">
        <eq/>
        <ci> alpha_X </ci>
        <apply>
          <divide/>
          <apply>
            <times/>
            <cn cellml:units="per_millivolt_millisecond"> 0.0000719 </cn>
            <apply>
              <plus/>
              <ci> V </ci>
              <cn cellml:units="millivolt"> 30.0 </cn>
            </apply>
          </apply>
          <apply>
            <minus/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
              <exp/>
              <apply>
                <times/>
                <cn cellml:units="per_millivolt"> -0.148 </cn>
                <apply>
                  <plus/>
                  <ci> V </ci>
                  <cn cellml:units="millivolt"> 30.0 </cn>
                </apply>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply id="beta_X_calculation">
        <eq/>
        <ci> beta_X </ci>
        <apply>
          <divide/>
          <apply>
            <times/>
            <cn cellml:units="per_millivolt_millisecond"> 0.000131 </cn>
            <apply>
              <plus/>
              <ci> V </ci>
              <cn cellml:units="millivolt"> 30.0 </cn>
            </apply>
          </apply>
          <apply>
            <plus/>
            <cn cellml:units="dimensionless"> -1.0 </cn>
            <apply>
              <exp/>
              <apply>
                <times/>
                <cn cellml:units="per_millivolt"> 0.0687 </cn>
                <apply>
                  <plus/>
                  <ci> V </ci>
                  <cn cellml:units="millivolt"> 30.0 </cn>
                </apply>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply id="dX_dt">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> X </ci>
        </apply>
        <apply>
          <minus/>
          <apply>
            <times/>
            <ci> alpha_X </ci>
            <apply>
              <minus/>
              <cn cellml:units="dimensionless"> 1.0 </cn>
              <ci> X </ci>
            </apply>
          </apply>
          <apply>
            <times/>
            <ci> beta_X </ci>
            <ci> X </ci>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
      
  <component name="time_dependent_potassium_current_Xi_gate">
    <variable units="dimensionless" public_interface="out" name="Xi"/>

    <variable units="millivolt" public_interface="in" name="V"/>
    <variable units="millisecond" public_interface="in" name="time"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="Xi_calculation">
        <eq/>
        <ci> Xi </ci>
        <apply>
          <divide/>
          <cn cellml:units="dimensionless"> 1.0 </cn>
          <apply>
            <plus/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
              <exp/>
              <apply>
                <divide/>
                <apply>
                  <minus/>
                  <ci> V </ci>
                  <cn cellml:units="millivolt"> 40.0 </cn>
                </apply>
                <cn cellml:units="millivolt"> 40.0 </cn>
              </apply>
            </apply>
          </apply>
        </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_K1"/>
    
    <variable units="milliS_per_microF" name="g_K1"/>
    
    <variable units="millisecond" 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="millimolar" public_interface="in" name="Ki"/>
    <variable units="joule_per_mole_kelvin" public_interface="in" name="R"/>
    <variable units="kelvin" public_interface="in" name="T"/>
    <variable units="coulomb_per_millimole" public_interface="in" name="F"/> 
    
    <variable units="dimensionless" private_interface="in" name="K1_infinity"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="g_K1_calculation">
        <eq/>
        <ci> g_K1 </ci>
        <apply>
          <times/>
          <cn cellml:units="milliS_per_microF"> 0.75 </cn>
          <apply>
            <root/>
            <apply>
              <divide/>
              <ci> Ko </ci>
              <cn cellml:units="millimolar"> 5.4 </cn>
            </apply>
          </apply>
        </apply>
      </apply>
        
      <apply id="E_K1_calculation">
        <eq/>
        <ci> E_K1 </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 id="i_K1_calculation">
        <eq/>
        <ci> i_K1 </ci>
        <apply>
          <times/>
          <ci> g_K1 </ci>
          <ci> K1_infinity </ci>
          <apply>
            <minus/>     
            <ci> V </ci>
            <ci> E_K1 </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_millisecond" name="alpha_K1"/>
    <variable units="per_millisecond" name="beta_K1"/>
    
    <variable units="millivolt" public_interface="in" name="V"/>
    <variable units="millisecond" public_interface="in" name="time"/>
    <variable units="millivolt" public_interface="in" name="E_K1"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="alpha_K1_calculation">
        <eq/>
        <ci> alpha_K1 </ci>
        <apply>
          <divide/>
          <cn cellml:units="per_millisecond"> 1.02 </cn>
          <apply>
            <plus/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
              <exp/>
              <apply>
                <times/>
                <cn cellml:units="per_millivolt"> 0.2385 </cn>
                <apply>
                  <minus/>
                  <ci> V </ci>
                  <apply>
                    <plus/>
                    <ci> E_K1 </ci>
                    <cn cellml:units="millivolt"> 59.215 </cn>
                  </apply>
                </apply>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
          
      <apply id="beta_K1_calculation">
        <eq/>
        <ci> beta_K1 </ci>
        <apply>
          <divide/>
          <apply>
            <plus/>
            <apply>
              <times/>
              <cn cellml:units="per_millisecond"> 0.49124 </cn>
              <apply>
                <exp/>
                <apply>
                  <times/>
                  <cn cellml:units="dimensionless"> 0.08032 </cn>
                  <apply>
                    <plus/>
                    <apply>
                      <minus/>
                      <ci> V </ci>
                      <ci> E_K1 </ci>
                    </apply>
                    <cn cellml:units="millivolt"> 5.476 </cn>
                  </apply>
                </apply>
              </apply>
            </apply>
            <apply>
              <exp/>
              <apply>
                <times/>
                <cn cellml:units="dimensionless"> 0.06175 </cn>
                <apply>
                  <minus/>
                  <apply>
                    <minus/>
                    <ci> V </ci>
                    <ci> E_K1 </ci>
                  </apply>
                  <cn cellml:units="millivolt"> 594.31 </cn>
                </apply>
              </apply>
            </apply>
          </apply>
          <apply>
            <plus/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
              <exp/>
              <apply>
                <times/>
                <cn cellml:units="per_millivolt"> -0.5143 </cn>
                <apply>
                  <plus/>
                  <apply>
                    <minus/>
                    <ci> V </ci>
                    <ci> E_K1 </ci>
                  </apply>
                  <cn cellml:units="millivolt"> 4.753 </cn>
                </apply>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
          
      <apply id="K1_infinity_calculation">
        <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="millivolt" name="E_Kp"/>
    <variable units="milliS_per_microF" name="g_Kp" initial_value="0.00828"/>
    <variable units="dimensionless" name="Kp"/>
    
    <variable units="millisecond" public_interface="in" name="time"/>
    <variable units="millivolt" public_interface="in" name="V"/>
    <variable units="millivolt" public_interface="in" name="E_K1"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="E_Kp_calculation">
        <eq/>
        <ci> E_Kp </ci>
        <ci> E_K1 </ci>
      </apply>
      
      <apply id="Kp_calculation">
        <eq/>
        <ci> Kp </ci>
        <apply>
          <divide/>
          <cn cellml:units="dimensionless"> 1.0 </cn>
          <apply>
            <plus/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
              <exp/>
              <apply>
                <divide/>
                <apply>
                  <minus/>
                  <cn cellml:units="millivolt"> 7.488 </cn>
                  <ci> V </ci>
                </apply>
                <cn cellml:units="millivolt"> 5.98 </cn>
              </apply>
            </apply> 
          </apply>
        </apply>
      </apply>
              
      <apply id="i_Kp_calculation">
        <eq/>
        <ci> i_Kp </ci>
        <apply>
          <times/>
          <ci> g_Kp </ci>
          <ci> Kp </ci>
          <apply>
            <minus/>
            <ci> V </ci>
            <ci> E_Kp </ci>
          </apply>
        </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="k_NaCa" initial_value="9000.0"/>
    <variable units="millimolar" name="K_mNa" initial_value="87.5"/>
    <variable units="millimolar" name="K_mCa" initial_value="1.38"/>
    <variable units="dimensionless" name="k_sat" initial_value="0.1"/>
    <variable units="dimensionless" name="eta" initial_value="0.35"/>
    
    <variable units="millisecond" public_interface="in" name="time"/>
    <variable units="millivolt" public_interface="in" name="V"/>
    <variable units="joule_per_mole_kelvin" public_interface="in" name="R"/>
    <variable units="kelvin" public_interface="in" name="T"/>
    <variable units="coulomb_per_millimole" public_interface="in" name="F"/>
    <variable units="millimolar" public_interface="in" name="Nai"/>
    <variable units="millimolar" public_interface="in" name="Nao"/>
    <variable units="micromolar" public_interface="in" name="Cai"/>
    <variable units="millimolar" public_interface="in" name="Cao"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="Na_Ca_exchanger">
        <eq/>
        <ci> i_NaCa </ci>
        <apply>
          <times/> 
          <ci> k_NaCa </ci>
          <apply>
            <divide/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
              <plus/>
              <apply>
                <power/>
                <ci> K_mNa </ci>
                <cn cellml:units="dimensionless"> 3.0 </cn>
              </apply>
              <apply>
                <power/>
                <ci> Nao </ci>
                <cn cellml:units="dimensionless"> 3.0 </cn>
              </apply>
            </apply>
          </apply>
          <apply>
            <divide/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
              <plus/>
              <ci> K_mCa </ci>
              <ci> Cao </ci>
            </apply>
          </apply>
          <apply>
            <divide/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
              <plus/>
              <cn cellml:units="dimensionless"> 1.0 </cn>
              <apply>
                <times/>
                <ci> k_sat </ci>
                <apply>
                  <exp/>
                  <apply>
                    <times/>
                    <apply>
                      <minus/>
                      <ci> eta </ci>
                      <cn cellml:units="dimensionless"> 1.0 </cn>
                    </apply>
                    <apply>
                      <divide/>
                      <apply>
                        <times/>
                        <ci> V </ci>
                        <ci> F </ci>
                      </apply>
                      <apply>
                        <times/>
                        <ci> R </ci>
                        <ci> T </ci>
                      </apply>
                    </apply>
                  </apply>
                </apply>
              </apply>
            </apply>
          </apply>
          <apply>
            <minus/>
            <apply>
              <times/>
              <apply>
                <exp/>
                <apply>
                  <times/>
                  <ci> eta </ci>
                  <apply>
                    <divide/>
                    <apply>
                      <times/>
                      <ci> V </ci>
                      <ci> F </ci>
                    </apply>
                    <apply>
                      <times/>
                      <ci> R </ci>
                      <ci> T </ci>
                    </apply>
                  </apply>
                </apply>
              </apply>
              <apply>
                <power/>
                <ci> Nai </ci>
                <cn cellml:units="dimensionless"> 3.0 </cn>
              </apply>
              <ci> Cao </ci>
            </apply>
            <apply>
              <times/>
              <apply>
                <exp/>
                <apply>
                  <times/>
                  <apply>
                    <minus/>
                    <ci> eta </ci>
                    <cn cellml:units="dimensionless"> 1.0 </cn>
                  </apply>
                  <apply>
                    <divide/>
                    <apply>
                      <times/>
                      <ci> V </ci>
                      <ci> F </ci>
                    </apply>
                    <apply>
                      <times/>
                      <ci> R </ci>
                      <ci> T </ci>
                    </apply>
                  </apply>
                </apply>
              </apply>
              <apply>
                <power/>
                <ci> Nao </ci>
                <cn cellml:units="dimensionless"> 3.0 </cn>
              </apply>
              <ci> Cai </ci>
            </apply>
          </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.47"/>
    <variable units="dimensionless" name="f_NaK"/>
    <variable units="millimolar" name="K_mNai" initial_value="10.0"/>
    <variable units="millimolar" name="K_mKo" initial_value="1.5"/>
    <variable units="dimensionless" name="sigma"/>
    
    <variable units="millisecond" public_interface="in" name="time"/>
    <variable units="millivolt" public_interface="in" name="V"/>
    <variable units="joule_per_mole_kelvin" public_interface="in" name="R"/>
    <variable units="kelvin" public_interface="in" name="T"/>
    <variable units="coulomb_per_millimole" 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 id="f_NaK_calculation">
        <eq/>
        <ci> f_NaK </ci>
        <apply>
          <divide/>
          <cn cellml:units="dimensionless"> 1.0 </cn>
          <apply>
            <plus/>
            <apply>
              <plus/>
              <cn cellml:units="dimensionless"> 1.0 </cn>
              <apply>
                <times/>
                <cn cellml:units="dimensionless"> 0.1245 </cn>
                <apply>
                  <exp/>
                  <apply>
                    <times/>
                    <cn cellml:units="dimensionless"> -0.1 </cn>
                    <apply>
                      <divide/>
                      <apply>
                        <times/>
                        <ci> V </ci>
                        <ci> F </ci>
                      </apply>
                      <apply>
                        <times/>
                        <ci> R </ci>
                        <ci> T </ci>
                      </apply>
                    </apply>
                  </apply>
                </apply>
              </apply>
            </apply>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> 0.0365 </cn>
              <ci> sigma </ci>
              <apply>
                <exp/>
                <apply>
                  <minus/>
                  <apply>
                    <divide/>
                    <apply>
                      <times/>
                      <ci> V </ci>
                      <ci> F </ci>
                    </apply>
                    <apply>
                      <times/>
                      <ci> R </ci>
                      <ci> T </ci>
                    </apply>
                  </apply>
                </apply>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply id="sigma_calculation">
        <eq/>
        <ci> sigma </ci>
        <apply>
          <times/>
          <apply>
            <divide/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <cn cellml:units="dimensionless"> 7.0 </cn>
          </apply>
          <apply>
            <minus/>
            <apply>
              <exp/>
              <apply>
                <divide/>
                <ci> Nao </ci>
                <cn cellml:units="dimensionless"> 67.3 </cn>
              </apply>
            </apply>
            <cn cellml:units="dimensionless"> 1.0 </cn>
          </apply>
        </apply>
      </apply>
      
      <apply id="i_NaK_calculation">
        <eq/>
        <ci> i_NaK </ci>
        <apply>
          <times/>
          <ci> I_NaK </ci>
          <ci> f_NaK </ci>
          <apply>
            <divide/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
              <plus/>
              <cn cellml:units="dimensionless"> 1.0 </cn>
              <apply>
                <power/>
                <apply>
                  <divide/>
                  <ci> K_mNai </ci>
                  <ci> Nai </ci>
                </apply>
                <cn cellml:units="dimensionless"> 1.5 </cn>
              </apply>
            </apply>
          </apply>
          <apply>
            <divide/>
            <ci> Ko </ci>
            <apply>
              <plus/>
              <ci> Ko </ci>
              <ci> K_mKo </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="microA_per_microF" name="I_ns_Ca"/>
    <variable units="microA_per_microF" name="I_ns_Na"/>
    <variable units="microA_per_microF" name="I_ns_K"/>
    <variable units="micromolar" name="K_m_ns_Ca" initial_value="1.2"/>
    <variable units="cm_per_second" name="Pns_Na" initial_value="0.0"/>
    <variable units="cm_per_second" name="Pns_K" initial_value="0.0"/>
    
    <variable units="millisecond" public_interface="in" name="time"/>
    <variable units="micromolar" public_interface="in" name="Cai"/>
    <variable units="micromolar" public_interface="in" name="Nai"/>
    <variable units="micromolar" public_interface="in" name="Ki"/>
    <variable units="micromolar" public_interface="in" name="Nao"/>
    <variable units="micromolar" public_interface="in" name="Ko"/>
    <variable units="millivolt" public_interface="in" name="V"/>
    <variable units="joule_per_mole_kelvin" public_interface="in" name="R"/>
    <variable units="kelvin" public_interface="in" name="T"/>
    <variable units="coulomb_per_millimole" public_interface="in" name="F"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="i_ns_Na_calculation">
        <eq/>
        <ci> i_ns_Na </ci>
        <apply>
          <times/>
          <ci> I_ns_Na </ci>
          <apply>
            <divide/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
              <plus/>
              <cn cellml:units="dimensionless"> 1.0 </cn>
              <apply>
                <power/>
                <apply>
                  <divide/>
                  <ci> K_m_ns_Ca </ci>
                  <ci> Cai </ci>
                </apply>
                <cn cellml:units="dimensionless"> 3.0 </cn>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply id="i_ns_K_calculation">
        <eq/>
        <ci> i_ns_K </ci>
        <apply>
          <times/>
          <ci> I_ns_K </ci>
          <apply>
            <divide/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
              <plus/>
              <cn cellml:units="dimensionless"> 1.0 </cn>
              <apply>
                <power/>
                <apply>
                  <divide/>
                  <ci> K_m_ns_Ca </ci>
                  <ci> Cai </ci>
                </apply>
                <cn cellml:units="dimensionless"> 3.0 </cn>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply id="i_ns_Ca_calculation">
        <eq/>
        <ci> i_ns_Ca </ci>
        <apply>
          <plus/>
          <ci> i_ns_Na </ci>
          <ci> i_ns_K </ci>
        </apply>
      </apply>
      
      <apply id="I_ns_Na_calculation">
        <eq/>
        <ci> I_ns_Na </ci>
        <apply>
          <times/>
          <ci> Pns_Na </ci>
          <apply>
            <divide/>
            <apply>
              <times/>
              <ci> V </ci>
              <apply>
                <power/>
                <ci> F </ci>
                <cn cellml:units="dimensionless"> 2.0 </cn>
              </apply>
            </apply>
            <apply>
              <times/>
              <ci> R </ci>
              <ci> T </ci>
            </apply>
          </apply>
          <apply>
            <divide/>
            <apply>
              <minus/>
              <apply>
                <times/>
                <cn cellml:units="dimensionless"> 0.75 </cn>
                <ci> Nai </ci>
                <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>
              <apply>
                <times/>
                <cn cellml:units="dimensionless"> 0.75 </cn>
                <ci> Nao </ci>
              </apply>
            </apply>
            <apply>
              <minus/>
              <apply>
                <exp/>
                <apply>
                  <divide/>
                  <apply>
                    <times/>
                    <ci> V </ci>
                    <ci> F </ci>
                  </apply>
                  <apply>
                    <times/>
                    <ci> R </ci>
                    <ci> T </ci>
                  </apply>
                </apply>
              </apply>
              <cn cellml:units="dimensionless"> 1.0 </cn>
            </apply>
          </apply>
        </apply>
      </apply>     
      
      <apply id="I_ns_K_calculation">
        <eq/>
        <ci> I_ns_K </ci>
        <apply>
          <times/>
          <ci> Pns_K </ci>
          <apply>
            <divide/>
            <apply>
              <times/>
              <ci> V </ci>
              <apply>
                <power/>
                <ci> F </ci>
                <cn cellml:units="dimensionless"> 2.0 </cn>
              </apply>
            </apply>
            <apply>
              <times/>
              <ci> R </ci>
              <ci> T </ci>
            </apply>
          </apply>
          <apply>
            <divide/>
            <apply>
              <minus/>
              <apply>
                <times/>
                <cn cellml:units="dimensionless"> 0.75 </cn>
                <ci> Ki </ci>
                <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>
              <apply>
                <times/>
                <cn cellml:units="dimensionless"> 0.75 </cn>
                <ci> Ko </ci>
              </apply>
            </apply>
            <apply>
              <minus/>
              <apply>
                <exp/>
                <apply>
                  <divide/>
                  <apply>
                    <times/>
                    <ci> V </ci>
                    <ci> F </ci>
                  </apply>
                  <apply>
                    <times/>
                    <ci> R </ci>
                    <ci> T </ci>
                  </apply>
                </apply>
              </apply>
              <cn cellml:units="dimensionless"> 1.0 </cn>
            </apply>
          </apply>
        </apply>
      </apply>     
    </math>
  </component>
  
  <component name="sarcolemmal_calcium_pump">
    <variable units="microA_per_microF" public_interface="out" name="i_p_Ca"/>
    
    <variable units="micromolar" name="K_mpCa" initial_value="0.5"/>
    <variable units="microA_per_microF" name="I_pCa" initial_value="0.575"/>
    
    <variable units="millisecond" public_interface="in" name="time"/>
    <variable units="micromolar" public_interface="in" name="Cai"/>   
 
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="i_p_Ca_calculation">
        <eq/>
        <ci> i_p_Ca </ci>
        <apply>
          <times/>
          <ci> I_pCa </ci>
          <apply>
            <divide/>
            <ci> Cai </ci>
            <apply>
              <plus/>
              <ci> K_mpCa </ci>
              <ci> Cai </ci>
            </apply>
          </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.006032"/>
    <variable units="millivolt" name="E_CaN"/>
    
    <variable units="millisecond" public_interface="in" name="time"/>
    <variable units="millivolt" public_interface="in" name="V"/>
    <variable units="joule_per_mole_kelvin" public_interface="in" name="R"/>
    <variable units="kelvin" public_interface="in" name="T"/>
    <variable units="coulomb_per_millimole" public_interface="in" name="F"/>
    <variable units="micromolar" public_interface="in" name="Cai"/>
    <variable units="millimolar" public_interface="in" name="Cao"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="E_CaN_calculation">
        <eq/>
        <ci> E_CaN </ci>
        <apply>
          <times/>
          <apply>
            <divide/>
            <apply>
              <times/>
              <ci> R </ci>
              <ci> T </ci>
            </apply>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> 2.0 </cn>
              <ci> F </ci>
            </apply>
          </apply>
          <apply>
            <ln/>
            <apply>
              <divide/>
              <ci> Cao </ci>
              <ci> Cai </ci>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply id="i_Ca_b_calculation">
        <eq/>
        <ci> i_Ca_b </ci>
        <apply>
          <times/>
          <ci> g_Cab </ci>
          <apply>
            <minus/>
            <ci> V </ci>
            <ci> E_CaN </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.00141"/>
    <variable units="millivolt" name="E_NaN"/>
    
    <variable units="millisecond" public_interface="in" name="time"/>
    <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 id="E_NaN_calculation">
        <eq/>
        <ci> E_NaN </ci>
        <ci> E_Na </ci>
      </apply>
      
      <apply id="i_Na_b_calculation">
        <eq/>
        <ci> i_Na_b </ci>
        <apply>
          <times/>
          <ci> g_Nab </ci>
          <apply>
            <minus/>
            <ci> V </ci>
            <ci> E_NaN </ci>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="L_type_Ca_channel">
    <variable units="microA_per_microF" public_interface="out" name="i_Ca"/>
    <variable units="microA_per_microF" public_interface="out" name="i_Ca_K"/>
    
    <variable units="cm_per_second" name="P_Ca" initial_value="4.32E-3"/>
    <variable units="cm_per_second" name="P_K" initial_value="1.11E-11"/>
    <variable units="cm_per_second" name="p_k"/>
    <variable units="microA_per_microF" name="i_Ca_half" initial_value="-2.16"/>
    <variable units="microA_per_microF" name="i_Ca_max"/>
    <variable units="dimensionless" name="O" initial_value="4.21370E-19"/>
    <variable units="dimensionless" name="O_Ca" initial_value="0.0"/>
    <variable units="per_millisecond" name="alpha"/>
    <variable units="per_millisecond" name="beta"/>
    <variable units="per_millisecond" name="gamma"/>
    <variable units="per_millisecond" name="alpha_a"/>
    <variable units="per_millisecond" name="beta_b"/>
    <variable units="dimensionless" name="a" initial_value="2.0"/>
    <variable units="dimensionless" name="b" initial_value="2.0"/>
    <variable units="per_millisecond" name="g" initial_value="2.0"/>
    <variable units="per_millisecond" name="f" initial_value="0.3"/>
    <variable units="per_millisecond" name="g_" initial_value="0.0"/>
    <variable units="per_millisecond" name="f_" initial_value="0.0"/>
    <variable units="per_millisecond" name="omega" initial_value="0.01"/>
    <variable units="dimensionless" name="C0" initial_value="0.998666"/>
    <variable units="dimensionless" name="C1" initial_value="1.63594E-4"/>
    <variable units="dimensionless" name="C2" initial_value="1.00495E-8"/>
    <variable units="dimensionless" name="C3" initial_value="2.74371E-13"/>
    <variable units="dimensionless" name="C4" initial_value="2.80908E-18"/>
    <variable units="dimensionless" name="C_Ca0" initial_value="0.00272826"/>
    <variable units="dimensionless" name="C_Ca1" initial_value="7.67674E-7"/>
    <variable units="dimensionless" name="C_Ca2" initial_value="1.88635E-10"/>
    <variable units="dimensionless" name="C_Ca3" initial_value="2.06007E-14"/>
    <variable units="dimensionless" name="C_Ca4" initial_value="8.43672E-19"/>
    
    <variable units="millisecond" 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="Ca_SS"/>
    <variable units="millimolar" public_interface="in" name="Cao"/>
    <variable units="millimolar" public_interface="in" name="Ko"/>
    <variable units="millimolar" public_interface="in" name="Ki"/>
    <variable units="joule_per_mole_kelvin" public_interface="in" name="R"/>
    <variable units="kelvin" public_interface="in" name="T"/>
    <variable units="coulomb_per_millimole" public_interface="in" name="F"/>
    
    <variable units="dimensionless" private_interface="in" name="y"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">    
      <apply id="i_Ca_calculation">
        <eq/>
        <ci> i_Ca </ci>
        <apply>
          <times/>
          <ci> i_Ca_max </ci>
          <ci> y </ci>
          <apply>
            <plus/>
            <ci> O </ci>
            <ci> O_Ca </ci>
          </apply>
        </apply>
      </apply>
      
      <apply id="i_Ca_max_calculation">
        <eq/>
        <ci> i_Ca_max </ci>
        <apply>
          <times/>
          <apply>
            <power/>
            <ci> P_Ca </ci>
            <cn cellml:units="dimensionless"> 4.0 </cn>
          </apply>
          <apply>
            <divide/>
            <apply>
              <times/>
              <ci> V </ci>
              <apply>
                <power/>
                <ci> F </ci>
                <cn cellml:units="dimensionless"> 2.0 </cn>
              </apply>
            </apply>
            <apply>
              <times/>
              <ci> R </ci>
              <ci> T </ci>
            </apply>
          </apply>
          <apply>
            <divide/>
            <apply>
              <minus/>
              <apply>
                <times/>
                <cn cellml:units="dimensionless"> 0.001 </cn>
                <apply>
                  <exp/>
                  <apply>
                    <divide/>
                    <apply>
                      <times/>
                      <cn cellml:units="dimensionless"> 2.0 </cn>
                      <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.341 </cn>
                <ci> Cao </ci>
              </apply>
            </apply> 
            <apply>
              <minus/>
              <apply>
                <exp/>
                <apply>
                  <divide/>
                  <apply>
                    <times/>
                    <cn cellml:units="dimensionless"> 2.0 </cn>
                    <ci> V </ci>
                    <ci> F </ci>
                  </apply>
                  <apply>
                    <times/>
                    <ci> R </ci>
                    <ci> T </ci>
                  </apply>
                </apply>
              </apply>
              <cn cellml:units="dimensionless"> 1.0 </cn>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply id="i_Ca_K_calculation">
        <eq/>
        <ci> i_Ca_K </ci>
        <apply>
          <times/>
          <ci> p_k </ci>
          <ci> y </ci>
          <apply>
            <plus/>
            <ci> O </ci>
            <ci> O_Ca </ci>
          </apply>
          <apply>
            <divide/>
            <apply>
              <times/>
              <ci> V </ci>
              <apply>
                <power/>
                <ci> F </ci>
                <cn cellml:units="dimensionless"> 2.0 </cn>
              </apply>
            </apply>
            <apply>
              <times/>
              <ci> R </ci>
              <ci> T </ci>
            </apply>
          </apply>
          <apply>
            <divide/>
            <apply>
              <minus/>
              <apply>
                <times/>
                <ci> Ki </ci>
                <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>
              <ci> Ko </ci>
            </apply> 
            <apply>
              <minus/>
              <apply>
                <exp/>
                <apply>
                  <divide/>
                  <apply>
                    <times/>
                    <ci> V </ci>
                    <ci> F </ci>
                  </apply>
                  <apply>
                    <times/>
                    <ci> R </ci>
                    <ci> T </ci>
                  </apply>
                </apply>
              </apply>
              <cn cellml:units="dimensionless"> 1.0 </cn>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply id="p_k_calculation">
        <eq/>
        <ci> p_k </ci>
        <apply>
          <divide/>
          <ci> P_K </ci>
          <apply>
            <plus/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
              <divide/>
              <ci> i_Ca_max </ci>
              <ci> i_Ca_half </ci>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply id="alpha_calculation">
        <eq/>
        <ci> alpha </ci>
        <apply>
          <times/>
          <cn cellml:units="per_millisecond"> 0.4 </cn>
          <apply>
            <exp/>
            <apply>
              <divide/>
              <apply>
                <plus/>
                <ci> V </ci>
                <cn cellml:units="millivolt"> 12.0 </cn>
              </apply>
              <cn cellml:units="millivolt"> 10.0 </cn>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply id="beta_calculation">
        <eq/>
        <ci> beta </ci>
        <apply>
          <times/>
          <cn cellml:units="per_millisecond"> 0.05 </cn>
          <apply>
            <exp/>
            <apply>
              <divide/>
              <apply>
                <minus/>
                <apply>
                  <plus/>
                  <ci> V </ci>
                  <cn cellml:units="millivolt"> 12.0 </cn>
                </apply>
              </apply>
              <cn cellml:units="millivolt"> 13.0 </cn>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply id="alpha_a_calculation">
        <eq/>
        <ci> alpha_a </ci>
        <apply>
          <times/>
          <ci> alpha </ci>
          <ci> a </ci>
        </apply>
      </apply>
      
      <apply id="beta_b_calculation">
        <eq/>
        <ci> beta_b </ci>
        <apply>
          <divide/>
          <ci> beta </ci>
          <ci> b </ci>
        </apply>
      </apply>
      
      <apply id="gamma_calculation">
        <eq/>
        <ci> gamma </ci>
        <apply>
          <times/>
          <cn cellml:units="dimensionless"> 0.5625 </cn>
          <ci> Ca_SS </ci>
        </apply>
      </apply>
      
      <apply id="C0_diff_eq">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> C0 </ci>
        </apply>
        <apply>
          <minus/>
          <apply>
            <plus/>
            <apply>
              <times/>
              <ci> beta </ci>
              <ci> C1 </ci>
            </apply>
            <apply>
              <times/>
              <ci> omega </ci>
              <ci> C_Ca0 </ci>
            </apply>
          </apply>
          <apply>
            <times/>
            <apply>
              <plus/>
              <apply>
                <times/>
                <cn cellml:units="dimensionless"> 4.0 </cn>
                <ci> alpha </ci>
              </apply>
              <ci> gamma </ci>
            </apply>
            <ci> C0 </ci>
          </apply>
        </apply>
      </apply> 
      
      <apply id="C1_diff_eq">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> C1 </ci>
        </apply>
        <apply>
          <minus/>
          <apply>
            <plus/>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> 4.0 </cn>
              <ci> alpha </ci>
              <ci> C0 </ci>
            </apply>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> 2.0 </cn>
              <ci> beta </ci>
              <ci> C2 </ci>
            </apply>
            <apply>
              <times/>
              <apply>
                <divide/>
                <ci> omega </ci>
                <ci> b </ci>
              </apply>
              <ci> C_Ca1 </ci>
            </apply>
          </apply>
          <apply>
            <times/>
            <apply>
              <plus/>
              <ci> beta </ci>
              <apply>
                <times/>
                <cn cellml:units="dimensionless"> 3.0 </cn>
                <ci> alpha </ci>
              </apply>
              <apply>
                <times/>
                <ci> gamma </ci>
                <ci> a </ci>
              </apply>
            </apply>
            <ci> C1 </ci>
          </apply>
        </apply>
      </apply>
      
      <apply id="C2_diff_eq">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> C2 </ci>
        </apply>
        <apply>
          <minus/>
          <apply>
            <plus/>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> 3.0 </cn>
              <ci> alpha </ci>
              <ci> C1 </ci>
            </apply>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> 3.0 </cn>
              <ci> beta </ci>
              <ci> C3 </ci>
            </apply>
            <apply>
              <times/>
              <apply>
                <divide/>
                <ci> omega </ci>
                <apply>
                  <power/>
                  <ci> b </ci>
                  <cn cellml:units="dimensionless"> 2.0 </cn>
                </apply>
              </apply>
              <ci> C_Ca2 </ci>
            </apply>
          </apply>
          <apply>
            <times/>
            <apply>
              <plus/>
              <apply>
                <times/>
                <ci> beta </ci>
                <cn cellml:units="dimensionless"> 2.0 </cn>
              </apply>
              <apply>
                <times/>
                <cn cellml:units="dimensionless"> 2.0 </cn>
                <ci> alpha </ci>
              </apply>
              <apply>
                <times/>
                <ci> gamma </ci>
                <apply>
                  <power/>
                  <ci> a </ci>
                  <cn cellml:units="dimensionless"> 2.0 </cn>
                </apply>
              </apply>
            </apply>
            <ci> C2 </ci>
          </apply>
        </apply>
      </apply> 
      
      <apply id="C3_diff_eq">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> C3 </ci>
        </apply>
        <apply>
          <minus/>
          <apply>
            <plus/>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> 2.0 </cn>
              <ci> alpha </ci>
              <ci> C2 </ci>
            </apply>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> 4.0 </cn>
              <ci> beta </ci>
              <ci> C4 </ci>
            </apply>
            <apply>
              <times/>
              <apply>
                <divide/>
                <ci> omega </ci>
                <apply>
                  <power/>
                  <ci> b </ci>
                  <cn cellml:units="dimensionless"> 3.0 </cn>
                </apply>
              </apply>
              <ci> C_Ca3 </ci>
            </apply>
          </apply>
          <apply>
            <times/>
            <apply>
              <plus/>
              <apply>
                <times/>
                <ci> beta </ci>
                <cn cellml:units="dimensionless"> 3.0 </cn>
              </apply>
              <ci> alpha </ci>
              <apply>
                <times/>
                <ci> gamma </ci>
                <apply>
                  <power/>
                  <ci> a </ci>
                  <cn cellml:units="dimensionless"> 3.0 </cn>
                </apply>
              </apply>
            </apply>
            <ci> C3 </ci>
          </apply>
        </apply>
      </apply> 
      
      <apply id="C4_diff_eq">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> C4 </ci>
        </apply>
        <apply>
          <minus/>
          <apply>
            <plus/>
            <apply>
              <times/>
              <ci> alpha </ci>
              <ci> C3 </ci>
            </apply>
            <apply>
              <times/>
              <ci> g </ci>
              <ci> O </ci>
            </apply>
            <apply>
              <times/>
              <apply>
                <divide/>
                <ci> omega </ci>
                <apply>
                  <power/>
                  <ci> b </ci>
                  <cn cellml:units="dimensionless"> 4.0 </cn>
                </apply>
              </apply>
              <ci> C_Ca4 </ci>
            </apply>
          </apply>
          <apply>
            <times/>
            <apply>
              <plus/>
              <apply>
                <times/>
                <ci> beta </ci>
                <cn cellml:units="dimensionless"> 4.0 </cn>
              </apply>
              <ci> f </ci>
              <apply>
                <times/>
                <ci> gamma </ci>
                <apply>
                  <power/>
                  <ci> a </ci>
                  <cn cellml:units="dimensionless"> 4.0 </cn>
                </apply>
              </apply>
            </apply>
            <ci> C4 </ci>
          </apply>
        </apply>
      </apply> 
      
      <apply id="O_diff_eq">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> O </ci>
        </apply>
        <apply>
          <minus/>
          <apply>
            <times/>
            <ci> f </ci>
            <ci> C4 </ci>
          </apply>
          <apply>
            <times/>
            <ci> g </ci>
            <ci> O </ci>
          </apply>
        </apply>
      </apply>
      
      <apply id="C_Ca0_diff_eq">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> C_Ca0 </ci>
        </apply>
        <apply>
          <minus/>
          <apply>
            <plus/>
            <apply>
              <times/>
              <ci> beta_b </ci>
              <ci> C_Ca1 </ci>
            </apply>
            <apply>
              <times/>
              <ci> gamma </ci>
              <ci> C_Ca0 </ci>
            </apply>
          </apply>
          <apply>
            <times/>
            <apply>
              <plus/>
              <apply>
                <times/>
                <cn cellml:units="dimensionless"> 4.0 </cn>
                <ci> alpha_a </ci>
              </apply>
              <ci> omega </ci>
            </apply>
            <ci> C_Ca0 </ci>
          </apply>
        </apply>
      </apply>
      
      <apply id="C_Ca1_diff_eq">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> C_Ca1 </ci>
        </apply>
        <apply>
          <minus/>
          <apply>
            <plus/>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> 4.0 </cn>
              <ci> alpha_a </ci>
              <ci> C_Ca0 </ci>
            </apply>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> 2.0 </cn>
              <ci> beta_b </ci>
              <ci> C_Ca2 </ci>
            </apply>
            <apply>
              <times/>
              <ci> gamma </ci>
              <ci> a </ci>
              <ci> C1 </ci>
            </apply>
          </apply>
          <apply>
            <times/>
            <apply>
              <plus/>
              <ci> beta_b </ci>
              <apply>
                <times/>
                <cn cellml:units="dimensionless"> 3.0 </cn>
                <ci> alpha_a </ci>
              </apply>
              <apply>
                <divide/>
                <ci> omega </ci>
                <ci> b </ci>
              </apply>
            </apply>
            <ci> C_Ca1 </ci>
          </apply>
        </apply>
      </apply>
      
      <apply id="C_Ca2_diff_eq">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> C_Ca2 </ci>
        </apply>
        <apply>
          <minus/>
          <apply>
            <plus/>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> 3.0 </cn>
              <ci> alpha_a </ci>
              <ci> C_Ca1 </ci>
            </apply>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> 3.0 </cn>
              <ci> beta_b </ci>
              <ci> C_Ca3 </ci>
            </apply>
            <apply>
              <times/>
              <ci> gamma </ci>
              <apply>
                <power/>
                <ci> a </ci>
                <cn cellml:units="dimensionless"> 2.0 </cn> 
              </apply>
              <ci> C2 </ci>
            </apply>
          </apply>
          <apply>
            <times/>
            <apply>
              <plus/>
              <apply>
                <times/>
                <ci> beta_b </ci>
                <cn cellml:units="dimensionless"> 2.0 </cn>
              </apply>
              <apply>
                <times/>
                <cn cellml:units="dimensionless"> 2.0 </cn>
                <ci> alpha_a </ci>
              </apply>
              <apply>
                <divide/>
                <ci> omega </ci>
                <apply>
                  <power/>
                  <ci> b </ci>
                  <cn cellml:units="dimensionless"> 2.0 </cn>
                </apply>
              </apply>
            </apply>
            <ci> C_Ca2 </ci>
          </apply>
        </apply>
      </apply>
      
      <apply id="C_Ca3_diff_eq">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> C_Ca3 </ci>
        </apply>
        <apply>
          <minus/>
          <apply>
            <plus/>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> 2.0 </cn>
              <ci> alpha_a </ci>
              <ci> C_Ca2 </ci>
            </apply>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> 4.0 </cn>
              <ci> beta_b </ci>
              <ci> C_Ca4 </ci>
            </apply>
            <apply>
              <times/>
              <ci> gamma </ci>
              <apply>
                <power/>
                <ci> a </ci>
                <cn cellml:units="dimensionless"> 3.0 </cn>
              </apply> 
              <ci> C3 </ci>
            </apply>
          </apply>
          <apply>
            <times/>
            <apply>
              <plus/>
              <apply>
                <times/>
                <ci> beta_b </ci>
                <cn cellml:units="dimensionless"> 3.0 </cn>
              </apply>
              <ci> alpha_a </ci>
              <apply>
                <divide/>
                <ci> omega </ci>
                <apply>
                  <power/>
                  <ci> b </ci>
                  <cn cellml:units="dimensionless"> 3.0 </cn>
                </apply>
              </apply>
            </apply>
            <ci> C_Ca3 </ci>
          </apply>
        </apply>
      </apply>
      
      <apply id="C_Ca4_diff_eq">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> C_Ca4 </ci>
        </apply>
        <apply>
          <minus/>
          <apply>
            <plus/>
            <apply>
              <times/>
              <ci> alpha_a </ci>
              <ci> C_Ca3 </ci>
            </apply>
            <apply>
              <times/>
              <ci> g_ </ci>
              <ci> O_Ca </ci>
            </apply>
            <apply>
              <times/>
              <ci> gamma </ci>
              <apply>
                <power/>
                <ci> a </ci>
                <cn cellml:units="dimensionless"> 4.0 </cn>
              </apply> 
              <ci> C4 </ci>
            </apply>
          </apply>
          <apply>
            <times/>
            <apply>
              <plus/>
              <apply>
                <times/>
                <ci> beta_b </ci>
                <cn cellml:units="dimensionless"> 4.0 </cn>
              </apply>
              <ci> f_ </ci>
              <apply>
                <divide/>
                <ci> omega </ci>
                <apply>
                  <power/>
                  <ci> b </ci>
                  <cn cellml:units="dimensionless"> 4.0 </cn>
                </apply>
              </apply>
            </apply>
            <ci> C_Ca4 </ci>
          </apply>
        </apply>
      </apply>
      
      <apply id="O_Ca_diff_eq">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> O_Ca </ci>
        </apply>
        <apply>
          <minus/>
          <apply>
            <times/>
            <ci> f_ </ci>
            <ci> C_Ca4 </ci>
          </apply>
          <apply>
            <times/>
            <ci> g_ </ci>
            <ci> O_Ca </ci>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="L_type_Ca_channel_y_gate">
    <variable units="dimensionless" public_interface="out" name="y" initial_value="0.862010"/>
    
    <variable units="dimensionless" name="y_infinity"/>
    <variable units="dimensionless" name="tau_y"/>
    
    <variable units="millivolt" public_interface="in" name="V"/>
    <variable units="millisecond" public_interface="in" name="time"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="y_diff_eq">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> y </ci>
        </apply>
        <apply>
          <divide/>
          <apply>
            <minus/>
            <ci> y_infinity </ci>
            <ci> y </ci>
          </apply>
          <ci> tau_y </ci>
        </apply>
      </apply>
      
      <apply id="y_infinity_calculation">
        <eq/>
        <ci> y_infinity </ci>
        <apply>
          <plus/>
          <apply>
            <divide/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
              <plus/>
              <cn cellml:units="dimensionless"> 1.0 </cn>
              <apply>
                <exp/>
                <apply>
                  <divide/>
                  <apply>
                    <plus/>
                    <ci> V </ci>
                    <cn cellml:units="millivolt"> 55.0 </cn>
                  </apply>
                  <cn cellml:units="millivolt"> 7.5 </cn>
                </apply>
              </apply>
            </apply>
          </apply>
          <apply>
            <divide/>
            <cn cellml:units="dimensionless"> 0.1 </cn>
            <apply>
              <plus/>
              <cn cellml:units="dimensionless"> 1.0 </cn>
              <apply>
                <exp/>
                <apply>
                  <divide/>
                  <apply>
                    <plus/>
                    <apply>
                      <minus/>
                      <ci> V </ci>
                    </apply>
                    <cn cellml:units="millivolt"> 21.0 </cn>
                  </apply>
                  <cn cellml:units="millivolt"> 6.0 </cn>
                </apply>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply id="tau_y_calculation">
        <eq/>
        <ci> tau_y </ci>
        <apply>
          <plus/>
          <cn cellml:units="dimensionless"> 20.0 </cn>
          <apply>
            <divide/>
            <cn cellml:units="dimensionless"> 600.0 </cn>
            <apply>
              <plus/>
              <cn cellml:units="dimensionless"> 1.0 </cn>
              <apply>
                <exp/>
                <apply>
                  <divide/>
                  <apply>
                    <plus/>
                    <ci> V </ci>
                    <cn cellml:units="millivolt"> 30.0 </cn>
                  </apply>
                  <cn cellml:units="millivolt"> 9.5 </cn>
                </apply>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="RyR_channel_states">
    <variable units="dimensionless" public_interface="out" name="P_O1" initial_value="0.428917E-3"/>
    <variable units="dimensionless" public_interface="out" name="P_O2" initial_value="0.487455E-9"/>
    
    <variable units="dimensionless" name="P_C1" initial_value="0.637931"/>
    <variable units="dimensionless" name="P_C2" initial_value="0.361643"/>
    <variable units="per_micromolar4_per_millisecond" name="k_a_plus" initial_value="0.01215"/>
    <variable units="per_millisecond" name="k_a_minus" initial_value="0.1425"/>
    <variable units="per_micromolar3_per_millisecond" name="k_b_plus" initial_value="0.00405"/>
    <variable units="per_millisecond" name="k_b_minus" initial_value="1.930"/>
    <variable units="per_millisecond" name="k_c_plus" initial_value="0.010"/>
    <variable units="per_millisecond" name="k_c_minus" initial_value="0.0008"/>
    <variable units="dimensionless" name="n" initial_value="4.0"/>
    <variable units="dimensionless" name="m" initial_value="3.0"/>
    
    <variable units="millisecond" public_interface="in" name="time"/>
    <variable units="millimolar" public_interface="in" name="Ca_SS"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="P_C1_diff_eq">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> P_C1 </ci>
        </apply>
        <apply>
          <plus/>
          <apply>
            <times/>
            <apply>
              <minus/>
              <ci> k_a_plus </ci>
            </apply>
            <apply>
              <power/>
              <ci> Ca_SS </ci>
              <ci> n </ci>
            </apply>
            <ci> P_C1 </ci>
          </apply>
          <apply>
            <times/>
            <ci> k_a_minus </ci>
            <ci> P_O1 </ci>
          </apply>
        </apply>
      </apply>
      
      <apply id="P_O1_diff_eq">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> P_O1 </ci>
        </apply>
        <apply>
          <plus/>
          <apply>
            <minus/>
            <apply>
              <times/>
              <ci> k_a_plus </ci>
              <apply>
                <power/>
                <ci> Ca_SS </ci>
                <ci> n </ci>
              </apply>
              <ci> P_C1 </ci>
            </apply>
            <apply>
              <plus/>
              <apply>
                <times/>
                <ci> k_a_minus </ci>
                <ci> P_O1 </ci>
              </apply>
              <apply>
                <times/>
                <ci> k_b_plus </ci>
                <apply>
                  <power/>
                  <ci> Ca_SS </ci>
                  <ci> m </ci>
                </apply>
                <ci> P_O1 </ci>
              </apply>
              <apply>
                <times/>
                <ci> k_c_plus </ci>
                <ci> P_O1 </ci>
              </apply>
            </apply>
          </apply>
          <apply>
            <times/>
            <ci> k_b_minus </ci>
            <ci> P_O2 </ci>
          </apply>
          <apply>
            <times/>
            <ci> k_c_minus </ci>
            <ci> P_C2 </ci>
          </apply>
        </apply>
      </apply>
      
      <apply id="P_O2_diff_eq">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> P_O2 </ci>
        </apply>
        <apply>
          <minus/>
          <apply>
            <times/>
            <ci> k_b_plus </ci>
            <apply>
              <power/>
              <ci> Ca_SS </ci>
              <ci> m </ci>
            </apply>
            <ci> P_O1 </ci>
          </apply>
          <apply>
            <times/>
            <ci> k_b_minus </ci>
            <ci> P_O2 </ci>
          </apply>
        </apply>
      </apply>
      
      <apply id="P_C2_diff_eq">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> P_C2 </ci>
        </apply>
        <apply>
          <minus/>
          <apply>
            <times/>
            <ci> k_c_plus </ci>
            <ci> P_O1 </ci>
          </apply>
          <apply>
            <times/>
            <ci> k_c_minus </ci>
            <ci> P_C2 </ci>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="SERCA_pump">
    <variable units="micromolar_per_millisecond" public_interface="out" name="J_up"/>
    
    <variable units="micromolar" name="k_fb" initial_value="0.24"/>
    <variable units="millimolar" name="k_rb" initial_value="1.64"/>
    <variable units="micromolar_per_second" name="fb" initial_value="0.292"/>
    <variable units="micromolar_per_second" name="rb" initial_value="0.311"/>
    <variable units="micromolar_per_second" name="Vmax_f"/>
    <variable units="micromolar_per_second" name="Vmax_r"/>
    <variable units="dimensionless" name="n_fb" initial_value="1.4"/>
    <variable units="dimensionless" name="n_rb"/>
    
    <variable units="millisecond" public_interface="in" name="time"/>
    <variable units="micromolar" public_interface="in" name="Ca_NSR"/>
    <variable units="micromolar" public_interface="in" name="Cai"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="J_up_calculation">
        <eq/>
        <ci> J_up </ci>
        <apply>
          <divide/>
          <apply>
            <minus/>
            <apply>
              <times/>
              <ci> Vmax_f </ci>
              <ci> fb </ci>
            </apply>
            <apply>
              <times/>
              <ci> Vmax_r </ci>
              <ci> rb </ci>
            </apply>
          </apply>
          <apply>
            <plus/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <ci> fb </ci>
            <ci> rb </ci>
          </apply>
        </apply>
      </apply>
      
      <apply id="fb_calculation">
        <eq/>
        <ci> fb </ci>
        <apply>
          <power/>
          <apply>
            <divide/>
            <ci> Cai </ci>
            <ci> k_fb </ci>
          </apply> 
          <ci> n_fb </ci>
        </apply> 
      </apply> 
      
      <apply id="rb_calculation">
        <eq/>
        <ci> rb </ci>
        <apply>
          <power/>
          <apply>
            <divide/>
            <ci> Ca_NSR </ci>
            <ci> k_rb </ci>
          </apply> 
          <ci> n_rb </ci>
        </apply> 
      </apply> 
    </math>
  </component>
  
  <component name="intracellular_Ca_fluxes">
    <variable units="micromolar_per_millisecond" public_interface="out" name="J_rel"/>
    <variable units="micromolar_per_millisecond" public_interface="out" name="J_trpn"/>
    <variable units="micromolar_per_millisecond" public_interface="out" name="J_tr"/> 
    <variable units="micromolar_per_millisecond" public_interface="out" name="J_xfer"/>
    
    <variable units="per_millisecond" name="v1" initial_value="3.6"/>
    <variable units="millisecond" name="tau_tr" initial_value="0.5747"/>
    <variable units="millisecond" name="tau_xfer" initial_value="9.09"/>
    
    <variable units="micromolar" public_interface="in" name="HTRPN_tot"/>
    <variable units="micromolar" public_interface="in" name="LTRPN_tot"/>
    <variable units="millimolar" public_interface="in" name="HTRPNCa"/>
    <variable units="millimolar" public_interface="in" name="LTRPNCa"/>
    <variable units="per_micromolar_per_second" public_interface="in" name="k_htrpn_plus"/>
    <variable units="per_second" public_interface="in" name="k_htrpn_minus"/>
    <variable units="per_micromolar_per_second" public_interface="in" name="k_ltrpn_plus"/>
    <variable units="per_micromolar_per_second" public_interface="in" name="k_ltrpn_minus"/>
    <variable units="millisecond" public_interface="in" name="time"/>
    <variable units="dimensionless" public_interface="in" name="P_O1"/>
    <variable units="dimensionless" public_interface="in" name="P_O2"/>
    <variable units="micromolar" public_interface="in" name="Cai"/>
    <variable units="millimolar" public_interface="in" name="Ca_SS"/>
    <variable units="micromolar" public_interface="in" name="Ca_NSR"/>
    <variable units="millimolar" public_interface="in" name="Ca_JSR"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="J_rel_calculation">
        <eq/>
        <ci> J_rel </ci>
        <apply>
          <times/>
          <ci> v1 </ci>
          <apply>
            <plus/>
            <ci> P_O1 </ci>
            <ci> P_O2 </ci>
          </apply>
          <apply>
            <minus/>
            <ci> Ca_JSR </ci>
            <ci> Ca_SS </ci>
          </apply>
        </apply>
      </apply>
      
      <apply id="J_tr_calculation">
        <eq/>
        <ci> J_tr </ci>
        <apply>
          <divide/>
          <apply>
            <minus/>
            <ci> Ca_NSR </ci>
            <ci> Ca_JSR </ci>
          </apply>
          <ci> tau_tr </ci>
        </apply>
      </apply>
      
      <apply id="J_xfer_calculation">
        <eq/>
        <ci> J_xfer </ci>
        <apply>
          <divide/>
          <apply>
            <minus/>
            <ci> Ca_SS </ci>
            <ci> Cai </ci>
          </apply>
          <ci> tau_xfer </ci>
        </apply>
      </apply>
      
      <apply id="J_trpn_calculation">
        <eq/>
        <ci> J_trpn </ci>
        <apply>
          <plus/>
          <apply>
            <times/>
            <ci> k_htrpn_plus </ci>
            <ci> Cai </ci>
            <apply>
              <minus/>
              <apply>
                <minus/>
                <ci> HTRPN_tot </ci>
                <ci> HTRPNCa </ci>
              </apply>
              <apply>
                <times/>
                <ci> k_htrpn_minus </ci>
                <ci> HTRPNCa </ci>
              </apply>
            </apply>
          </apply>
          <apply>
            <times/>
            <ci> k_ltrpn_plus </ci>
            <ci> Cai </ci>
            <apply>
              <minus/>
              <apply>
                <minus/>
                <ci> LTRPN_tot </ci>
                <ci> LTRPNCa </ci>
              </apply>
              <apply>
                <times/>
                <ci> k_ltrpn_minus </ci>
                <ci> LTRPNCa </ci>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="intracellular_ionic_concentrations">
    <variable units="millimolar" public_interface="out" name="Nai"/>
    <variable units="millimolar" public_interface="out" name="Nao" initial_value="140.0"/>
    <variable units="millimolar" public_interface="out" name="Ki"/>
    <variable units="millimolar" public_interface="out" name="Ko" initial_value="5.4"/>
    <variable units="millimolar" public_interface="out" name="Cai"/>
    <variable units="millimolar" public_interface="out" name="Cao" initial_value="1.8"/>
    <variable units="millimolar" public_interface="out" name="Ca_SS"/>
    <variable units="millimolar" public_interface="out" name="Ca_NSR"/>
    <variable units="millimolar" public_interface="out" name="Ca_JSR"/>
              
    <variable units="dimensionless" name="betai"/>
    <variable units="dimensionless" name="beta_SS"/>
    <variable units="dimensionless" name="beta_JSR"/>
    <variable units="micromolar" name="K_mCMDN" initial_value="2.38"/>
    <variable units="millimolar" name="K_mCSQN" initial_value="0.8"/>
    <variable units="millimolar" name="CSQN_tot" initial_value="15.0"/>
    <variable units="micromolar" name="CMDN_tot" initial_value="50.0"/>
              
    <variable units="millisecond" public_interface="in" name="time"/>
    <variable units="micromolar_per_millisecond" public_interface="in" name="J_rel"/>
    <variable units="micromolar_per_millisecond" public_interface="in" name="J_up"/>
    <variable units="micromolar_per_millisecond" public_interface="in" name="J_tr"/>
    <variable units="micromolar_per_millisecond" public_interface="in" name="J_xfer"/>
    <variable units="micromolar_per_millisecond" public_interface="in" name="J_trpn"/>
    <variable units="coulomb_per_millimole" 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_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_Ca"/>
    <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_K"/>
    <variable units="microA_per_microF" public_interface="in" name="i_K"/>
    <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_ns_K"/>
    <variable units="cm2" public_interface="in" name="A_cap"/>
    <variable units="micro_litre" public_interface="in" name="V_myo"/>
    <variable units="micro_litre" public_interface="in" name="V_JSR"/>
    <variable units="micro_litre" public_interface="in" name="V_NSR"/>
    <variable units="micro_litre" public_interface="in" name="V_SS"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="betai_calculation">
        <eq/>
        <ci> betai </ci>
        <apply>
          <divide/>
          <cn cellml:units="dimensionless"> 1.0 </cn>
          <apply>
            <plus/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
              <divide/>
              <apply>
                <times/>
                <ci> CMDN_tot </ci>
                <ci> K_mCMDN </ci>
              </apply>
              <apply>
                <power/>
                <apply>
                  <plus/>
                  <ci> K_mCMDN </ci>
                  <ci> Cai </ci>
                </apply>
                <cn cellml:units="dimensionless"> 2.0 </cn>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply id="beta_SS_calculation">
        <eq/>
        <ci> beta_SS </ci>
        <apply>
          <divide/>
          <cn cellml:units="dimensionless"> 1.0 </cn>
          <apply>
            <plus/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
              <divide/>
              <apply>
                <times/>
                <ci> CMDN_tot </ci>
                <ci> K_mCMDN </ci>
              </apply>
              <apply>
                <power/>
                <apply>
                  <plus/>
                  <ci> K_mCMDN </ci>
                  <ci> Ca_SS </ci>
                </apply>
                <cn cellml:units="dimensionless"> 2.0 </cn>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply id="beta_JSR_calculation">
        <eq/>
        <ci> beta_JSR </ci>
        <apply>
          <divide/>
          <cn cellml:units="dimensionless"> 1.0 </cn>
          <apply>
            <plus/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
              <divide/>
              <apply>
                <times/>
                <ci> CSQN_tot </ci>
                <ci> K_mCSQN </ci>
              </apply>
              <apply>
                <power/>
                <apply>
                  <plus/>
                  <ci> K_mCSQN </ci>
                  <ci> Ca_JSR </ci>
                </apply>
                <cn cellml:units="dimensionless"> 2.0 </cn>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply id="calcium_internal_diff_eq">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> Cai </ci>
        </apply>
        <apply>
          <times/>
          <ci> betai </ci>
          <apply>
            <minus/>
            <apply>
              <plus/>
              <ci> J_xfer </ci>
              <ci> J_up </ci>
              <ci> J_trpn </ci>
              <apply>
                <times/>
                <apply>
                  <plus/>
                  <ci> i_Ca_b </ci>
                  <apply>
                    <minus/>
                    <apply>
                      <times/>
                      <cn cellml:units="dimensionless"> 2.0 </cn>
                      <ci> i_NaCa </ci>
                    </apply>
                  </apply>
                  <ci> i_p_Ca </ci>
                </apply>
                <apply>
                  <divide/>
                  <ci> A_cap </ci>
                  <apply>
                    <times/>
                    <cn cellml:units="dimensionless"> 2.0 </cn>
                    <ci> V_myo </ci>
                    <ci> F </ci>
                  </apply>
                </apply>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply id="calcium_subspace_diff_eq">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> Ca_SS </ci>
        </apply>
        <apply>
          <times/>
          <ci> beta_SS </ci>
          <apply>
            <minus/>
            <apply>
              <minus/>
              <apply>
                <times/>
                <ci> J_rel </ci>
                <apply>
                  <divide/>
                  <ci> V_JSR </ci>
                  <ci> V_SS </ci>
                </apply>
              </apply>
              <apply>
                <times/>
                <ci> J_xfer </ci>
                <apply>
                  <divide/>
                  <ci> V_myo </ci>
                  <ci> V_SS </ci>
                </apply>
              </apply>
            </apply>
            <apply>
              <times/>
              <ci> i_Ca </ci>
              <apply>
                <divide/>
                <ci> A_cap </ci>
                <apply>
                  <times/>
                  <cn cellml:units="dimensionless"> 2.0 </cn>
                  <ci> V_SS </ci>
                  <ci> F </ci>
                </apply>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply id="calcium_JSR_diff_eq">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> Ca_JSR </ci>
        </apply>
        <apply>
          <times/>
          <ci> beta_JSR </ci>
          <apply>
            <minus/>
            <ci> J_tr </ci>
            <ci> J_rel </ci>
          </apply>
        </apply>
      </apply>
      
      <apply id="calcium_NSR_diff_eq">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> Ca_NSR </ci>
        </apply>
        <apply>
          <minus/>
          <apply>
            <times/>
            <ci> J_up </ci>
            <apply>
              <divide/>
              <ci> V_myo </ci>
              <ci> V_NSR </ci>
            </apply>
          </apply>
          <apply>
            <times/>
            <ci> J_tr </ci>
            <apply>
              <divide/>
              <ci> V_JSR </ci>
              <ci> V_NSR </ci>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply id="sodium_diff_eq">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> Nai </ci>
        </apply>
        <apply>
          <times/>
          <apply>
            <minus/>
            <apply>
              <plus/>
              <ci> i_Na </ci>
              <ci> i_Na_b </ci>
              <ci> i_ns_Na </ci>
              <apply>
                <times/>
                <ci> i_NaCa </ci>
                <cn cellml:units="dimensionless"> 3.0 </cn>
              </apply>
              <apply>
                <times/>
                <ci> i_NaK </ci>
                <cn cellml:units="dimensionless"> 3.0 </cn>
              </apply>
            </apply>
          </apply>
          <apply>
            <divide/>
            <ci> A_cap </ci>
            <apply>
              <times/>
              <ci> V_myo </ci>
              <ci> F </ci>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply id="potassium_internal_diff_eq">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> Ki </ci>
        </apply>
        <apply>
          <times/>
          <apply>
            <minus/>
            <apply>
              <plus/>
              <ci> i_Ca_K </ci>
              <ci> i_K </ci>
              <ci> i_K1 </ci>
              <ci> i_Kp </ci>
              <ci> i_ns_K </ci>
              <apply>
                <minus/>
                <apply>
                  <times/>
                  <ci> i_NaK </ci>
                  <cn cellml:units="dimensionless"> 2.0 </cn>
                </apply>
              </apply>
            </apply>
          </apply>
          <apply>
            <divide/>
            <ci> A_cap </ci>
            <apply>
              <times/>
              <ci> V_myo </ci>
              <ci> F </ci>
            </apply>
          </apply>
        </apply>        
      </apply>
    </math>
  </component>
  
  <component name="troponin">
    <variable units="per_micromolar_per_second" public_interface="out" name="k_htrpn_plus" initial_value="1.0E8"/>
    <variable units="per_second" public_interface="out" name="k_htrpn_minus" initial_value="0.33"/>
    <variable units="per_micromolar_per_second" public_interface="out" name="k_ltrpn_plus" initial_value="4.0E7"/>
    <variable units="per_micromolar_per_second" public_interface="out" name="k_ltrpn_minus" initial_value="40.0"/>
    <variable units="micromolar" public_interface="out" name="HTRPN_tot" initial_value="140.0"/>
    <variable units="micromolar" public_interface="out" name="LTRPN_tot" initial_value="70.0"/>
    <variable units="millimolar" public_interface="out" name="HTRPNCa" initial_value="0.132944"/>
    <variable units="millimolar" public_interface="out" name="LTRPNCa" initial_value="0.4068232"/>
    
    <variable units="N_per_mm2" public_interface="in" name="Force_norm"/>
    <variable units="millimolar" public_interface="in" name="Cai"/>
    <variable units="millisecond" public_interface="in" name="time"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
  
      <apply id="HTRPNCa_diff_eq">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> HTRPNCa </ci>
        </apply>
        <apply>
          <minus/>
          <apply>
            <times/>
            <ci> k_htrpn_plus </ci>
            <ci> Cai </ci>
            <apply>
              <minus/>
              <ci> HTRPN_tot </ci>
              <ci> HTRPNCa </ci>
            </apply>
          </apply>
          <apply>
            <times/>
            <ci> k_htrpn_minus </ci>
            <ci> HTRPNCa </ci>
          </apply>
        </apply>
      </apply>
      
      <apply id="LTRPNCa_diff_eq">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> LTRPNCa </ci>
        </apply>
        <apply>
          <minus/>
          <apply>
            <times/>
            <ci> k_ltrpn_plus </ci>
            <ci> Cai </ci>
            <apply>
              <minus/>
              <ci> LTRPN_tot </ci>
              <ci> LTRPNCa </ci>
            </apply>
          </apply>
          <apply>
            <times/>
            <ci> k_ltrpn_minus </ci>
            <ci> LTRPNCa </ci>
            <apply>
              <plus/>
              <cn cellml:units="dimensionless"> 0.333 </cn>
              <apply>
                <times/>
                <cn cellml:units="dimensionless"> 0.667 </cn>
                <apply>
                  <minus/>
                  <cn cellml:units="dimensionless"> 1.0 </cn>
                  <ci> Force_norm </ci>
                </apply>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="tropomyosin_cross_bridges">
    <variable units="per_second" public_interface="out" name="f_01"/>
    <variable units="per_second" public_interface="out" name="f_12"/>
    <variable units="per_second" public_interface="out" name="f_23"/>
    <variable units="per_second" public_interface="out" name="g_XB" initial_value="30.0"/>
    <variable units="micrometre" public_interface="out" name="SL" initial_value="2.15"/>
    <variable units="dimensionless" public_interface="out" name="P1" initial_value="1.47026E-6"/>
    <variable units="dimensionless" public_interface="out" name="P2" initial_value="1.83527E-6"/>
    <variable units="dimensionless" public_interface="out" name="P3" initial_value="1.06786E-6"/>
    <variable units="dimensionless" public_interface="out" name="N1" initial_value="3.88506E-6"/>
    
    <variable units="per_second" name="f_XB" initial_value="10.0"/>
    <variable units="per_second" name="g_XB_SL"/>
    <variable units="per_second" name="g_01_SL"/>
    <variable units="per_second" name="g_12_SL"/>
    <variable units="per_second" name="g_23_SL"/>
    <variable units="dimensionless" name="N0" initial_value="0.99999027"/>
    <variable units="dimensionless" name="P0" initial_value="1.46758E-6"/>
    <variable units="micrometre" name="SL_norm"/>
    <variable units="micrometre" name="N_trop"/> 
    <variable units="micromolar" name="K_trop"/>
    <variable units="micromolar" name="k_trop_pn"/>
    <variable units="micromolar" name="k_trop_np"/>
    <variable units="micromolar" name="K_trop_half"/>
    <variable units="micromolar" name="K_trop_Ca"/>
    
    <variable units="millisecond" public_interface="in" name="time"/>
    <variable units="per_micromolar_per_second" public_interface="in" name="k_ltrpn_plus"/>
    <variable units="per_micromolar_per_second" public_interface="in" name="k_ltrpn_minus"/>
    <variable units="micromolar" public_interface="in" name="LTRPN_tot"/>
    <variable units="millimolar" public_interface="in" name="LTRPNCa"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="f_01_calculation">
        <eq/>
        <ci> f_01 </ci>
        <apply>
          <times/>
          <cn cellml:units="dimensionless"> 3.0 </cn>
          <ci> f_XB </ci>
        </apply>
      </apply>
      
      <apply id="f_12_calculation">
        <eq/>
        <ci> f_12 </ci>
        <apply>
          <times/>
          <cn cellml:units="dimensionless"> 10.0 </cn>
          <ci> f_XB </ci>
        </apply>
      </apply>
      
      <apply id="f_23_calculation">
        <eq/>
        <ci> f_23 </ci>
        <apply>
          <times/>
          <cn cellml:units="dimensionless"> 7.0 </cn>
          <ci> f_XB </ci>
        </apply>
      </apply> 
      
      <apply id="g_01_SL_calculation">
        <eq/>
        <ci> g_01_SL </ci>
        <apply>
          <times/>
          <cn cellml:units="dimensionless"> 1.0 </cn>
          <ci> g_XB_SL </ci>
        </apply>
      </apply>
      
      <apply id="g_12_SL_calculation">
        <eq/>
        <ci> g_12_SL </ci>
        <apply>
          <times/>
          <cn cellml:units="dimensionless"> 2.0 </cn>
          <ci> g_XB_SL </ci>
        </apply>
      </apply>
      
      <apply id="g_23_SL_calculation">
        <eq/>
        <ci> g_23_SL </ci>
        <apply>
          <times/>
          <cn cellml:units="dimensionless"> 3.0 </cn>
          <ci> g_XB_SL </ci>
        </apply>
      </apply>
      
      <apply id="g_XB_SL_calculation">
        <eq/>
        <ci> g_XB_SL </ci>
        <apply>
          <times/>
          <ci> g_XB </ci>
          <apply>
            <plus/>
            <cn cellml:units="micrometre"> 1.0 </cn>
            <apply>
              <power/>
              <apply>
                <minus/>
                <cn cellml:units="micrometre"> 1.0 </cn>
                <ci> SL_norm </ci>
              </apply>
              <cn cellml:units="dimensionless"> 1.6 </cn>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply id="SL_norm_calculation">
        <eq/>
        <ci> SL_norm </ci>
        <apply>
          <divide/>
          <apply>
            <minus/>
            <ci> SL </ci>
            <cn cellml:units="micrometre"> 1.7 </cn>
          </apply>
          <cn cellml:units="micrometre"> 0.7 </cn>
        </apply>
      </apply>
      
      <apply id="k_trop_np_calculation">
        <eq/>
        <ci> k_trop_np </ci>
        <apply>
          <times/>
          <ci> k_trop_pn </ci>
          <apply>
            <power/>
            <apply>
              <divide/>
              <apply>
                <divide/>
                <ci> LTRPNCa </ci>
                <ci> LTRPN_tot </ci>
              </apply>
              <ci> K_trop_half </ci>
            </apply>
            <ci> N_trop </ci>
          </apply>
        </apply>
      </apply>
      
      <apply id="N_trop_calculation">
        <eq/>
        <ci> N_trop </ci>
        <apply>
          <plus/>
          <cn cellml:units="dimensionless"> 3.5 </cn>
          <apply>
            <times/> 
            <cn cellml:units="dimensionless"> 2.5 </cn>
            <ci> SL_norm </ci>
          </apply>
        </apply>
      </apply>
      
      <apply id="K_trop_half_calculation">
        <eq/>
        <ci> K_trop_half </ci>
        <apply>
          <power/>
          <apply>
            <plus/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
              <divide/>
              <ci> K_trop_Ca </ci>
              <apply>
                <minus/>
                <cn cellml:units="micromolar"> 1.7 </cn>
                <apply>
                  <times/>
                  <cn cellml:units="micromolar"> 0.9 </cn>
                  <ci> SL_norm </ci>
                </apply>
              </apply>
            </apply>
          </apply>
          <cn cellml:units="dimensionless"> -1.0 </cn>
        </apply>
      </apply>
      
      <apply id="K_trop_Ca_calculation">
        <eq/>
        <ci> K_trop_Ca </ci>
        <apply>
          <divide/>
          <ci> k_ltrpn_minus </ci>
          <ci> k_ltrpn_plus </ci>
        </apply>
      </apply>
      
      <apply id="N0_diff_eq">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> N0 </ci>
        </apply>
        <apply>
          <plus/>
          <apply>
            <minus/>
            <apply>
              <times/>
              <ci> k_trop_pn </ci>
              <ci> P0 </ci>
            </apply>
            <apply>
              <times/>
              <ci> k_trop_np </ci>
              <ci> N0 </ci>
            </apply>
          </apply>
          <apply>
            <times/>
            <ci> g_01_SL </ci>
            <ci> N1 </ci>
          </apply>
        </apply>
      </apply> 
      
      <apply id="P0_diff_eq">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> P0 </ci>
        </apply>
        <apply>
          <plus/>
          <apply>
            <minus/>
            <apply>
              <times/>
              <apply>
                <plus/>
                <ci> k_trop_pn </ci>
                <ci> f_01 </ci>
              </apply>
              <ci> P0 </ci>
            </apply>
          </apply>
          <apply>
            <times/>
            <ci> k_trop_np </ci>
            <ci> N0 </ci>
          </apply>
          <apply>
            <times/>
            <ci> g_01_SL </ci>
            <ci> P1 </ci>
          </apply>
        </apply>
      </apply>
      
      <apply id="P1_diff_eq">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> P1 </ci>
        </apply>
        <apply>
          <plus/>
          <apply>
            <minus/>
            <apply>
              <times/>
              <apply>
                <plus/>
                <ci> k_trop_pn </ci>
                <ci> f_12 </ci>
                <ci> g_01_SL </ci>
              </apply>
              <ci> P1 </ci>
            </apply>
          </apply>
          <apply>
            <times/>
            <ci> k_trop_np </ci>
            <ci> N1 </ci>
          </apply>
          <apply>
            <times/>
            <ci> k_trop_np </ci>
            <ci> N1 </ci>
          </apply>
          <apply>
            <times/>
            <ci> f_01 </ci>
            <ci> P0 </ci>
          </apply>
          <apply>
            <times/>
            <ci> g_12_SL </ci>
            <ci> P2 </ci>
          </apply>
        </apply>
      </apply>
    
      <apply id="P2_diff_eq">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> P2 </ci>
        </apply>
        <apply>
          <plus/>
          <apply>
            <minus/>
            <apply>
              <times/>
              <apply>
                <plus/>
                <ci> f_23 </ci>
                <ci> g_12_SL </ci>
              </apply>
              <ci> P2 </ci>
            </apply>
          </apply>
          <apply>
            <times/>
            <ci> f_12 </ci>
            <ci> P1 </ci>
          </apply>
          <apply>
            <times/>
            <ci> g_23_SL </ci>
            <ci> P3 </ci>
          </apply>
        </apply>
      </apply>
    
      <apply id="P3_diff_eq">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> P3 </ci>
        </apply>
        <apply>
          <plus/>
          <apply>
            <times/>
            <apply>
              <minus/>
              <ci> g_23_SL </ci>
            </apply>
            <ci> P3 </ci>
          </apply>
          <apply>
            <times/>
            <ci> f_23 </ci>
            <ci> P2 </ci>
          </apply>
        </apply>
      </apply>
    </math>
  </component> 
 
  <component name="force_computation">
    <variable units="N_per_mm2" public_interface="out" name="Force_norm"/>
    
    <variable units="N_per_mm2" name="Force"/>
    <variable units="N_per_mm2" name="Force_max"/>
    <variable units="N_per_mm2" name="zeta" initial_value="0.1"/>
    <variable units="dimensionless" name="phi_SL"/>
    <variable units="per_second" name="P1_max"/>
    <variable units="per_second" name="P2_max"/>
    <variable units="per_second" name="P3_max"/>
    
    <variable units="millisecond" public_interface="in" name="time"/>
    <variable units="micrometre" public_interface="in" name="SL"/>
    <variable units="dimensionless" public_interface="in" name="P1"/>
    <variable units="dimensionless" public_interface="in" name="P2"/>
    <variable units="dimensionless" public_interface="in" name="P3"/>
    <variable units="dimensionless" public_interface="in" name="N1"/>
    <variable units="per_second" public_interface="in" name="f_01"/>
    <variable units="per_second" public_interface="in" name="f_12"/>
    <variable units="per_second" public_interface="in" name="f_23"/>
    <variable units="per_second" public_interface="in" name="g_XB"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="Force_calculation">
        <eq/>
        <ci> Force </ci>
        <apply>
          <times/>
          <ci> zeta </ci>
          <ci> Force_norm </ci>
        </apply>
      </apply>
      
      <apply id="Force_norm_calculation">
        <eq/>
        <ci> Force_norm </ci>
        <apply>
          <times/>
          <ci> phi_SL </ci>
          <apply>
            <divide/>
            <apply>
              <plus/>
              <ci> P1 </ci>
              <ci> N1 </ci>
              <apply>
                <times/>
                <cn cellml:units="dimensionless"> 2.0 </cn>
                <ci> P2 </ci>
              </apply>
              <apply>
                <times/>
                <cn cellml:units="dimensionless"> 3.0 </cn>
                <ci> P3 </ci>
              </apply>
            </apply>
            <ci> Force_max </ci>
          </apply>
        </apply>
      </apply>
      
      <apply id="Force_max_calculation">
        <eq/>
        <ci> Force_max </ci>
        <apply>
          <plus/>
          <ci> P1_max </ci>
          <apply>
            <times/>
            <cn cellml:units="dimensionless"> 2.0 </cn>
            <ci> P2_max </ci>
          </apply>
          <apply>
            <times/>
            <cn cellml:units="dimensionless"> 3.0 </cn>
            <ci> P3_max </ci>
          </apply>
        </apply>
      </apply>
     
      <apply id="phi_SL_calculation">
        <eq/>
        <ci> phi_SL </ci>
        <piecewise>
          <piece>
            <apply>
              <divide/>
              <apply>
                <minus/>
                <ci> SL </ci>
                <cn cellml:units="micrometre"> 0.6 </cn>
              </apply>
              <cn cellml:units="micrometre"> 1.4 </cn>
            </apply>
            <apply>
              <and/>
              <apply>
                <lt/>
                <ci> SL </ci> 
                <cn cellml:units="micrometre"> 2.0 </cn>
              </apply>
              <apply>
                <gt/>
                <ci> SL </ci> 
                <cn cellml:units="micrometre"> 1.7 </cn>
              </apply>
            </apply>
          </piece>
          <piece>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
              <and/>
              <apply>
                <lt/>
                <ci> SL </ci> 
                <cn cellml:units="micrometre"> 2.2 </cn>
              </apply>
              <apply>
                <gt/>
                <ci> SL </ci> 
                <cn cellml:units="micrometre"> 2.0 </cn>
              </apply>
            </apply>
          </piece>
          <piece>
            <apply>
              <divide/>
              <apply>
                <minus/>
                <cn cellml:units="micrometre"> 3.6 </cn>
                <ci> SL </ci>
              </apply>
              <cn cellml:units="micrometre"> 1.4 </cn>
            </apply>
            <apply>
              <and/>
              <apply>
                <lt/>
                <ci> SL </ci> 
                <cn cellml:units="micrometre"> 2.3 </cn>
              </apply>
              <apply>
                <gt/>
                <ci> SL </ci> 
                <cn cellml:units="micrometre"> 2.2 </cn>
              </apply>
            </apply>
          </piece>
        </piecewise>
      </apply>
      
      <apply id="P1_max_calculation">
        <eq/>
        <ci> P1_max </ci>
        <apply>
          <divide/>
          <apply>
            <times/>
            <ci> f_01 </ci>
            <cn cellml:units="dimensionless"> 2.0 </cn>
            <ci> g_XB </ci>
            <cn cellml:units="dimensionless"> 3.0 </cn>
            <ci> g_XB </ci>
          </apply>
          <apply>
            <plus/>
            <apply>
              <times/>
              <ci> g_XB </ci>
              <cn cellml:units="dimensionless"> 2.0 </cn>
              <ci> g_XB </ci>
              <cn cellml:units="dimensionless"> 3.0 </cn>
              <ci> g_XB </ci>
            </apply>
            <apply>
              <times/>
              <ci> f_01 </ci>
              <cn cellml:units="dimensionless"> 2.0 </cn>
              <ci> g_XB </ci>
              <cn cellml:units="dimensionless"> 3.0 </cn>
              <ci> g_XB </ci>
            </apply>
            <apply>
              <times/>
              <ci> f_01 </ci>
              <ci> f_12 </ci>
              <cn cellml:units="dimensionless"> 3.0 </cn>
              <ci> g_XB </ci>
            </apply>
            <apply>
              <times/>
              <ci> f_01 </ci>
              <ci> f_12 </ci>
              <ci> f_23 </ci>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply id="P2_max_calculation">
        <eq/>
        <ci> P2_max </ci>
        <apply>
          <divide/>
          <apply>
            <times/>
            <ci> f_01 </ci>
            <ci> f_12 </ci>
            <cn cellml:units="dimensionless"> 3.0 </cn>
            <ci> g_XB </ci>
          </apply>
          <apply>
            <plus/>
            <apply>
              <times/>
              <ci> g_XB </ci>
              <cn cellml:units="dimensionless"> 2.0 </cn>
              <ci> g_XB </ci>
              <cn cellml:units="dimensionless"> 3.0 </cn>
              <ci> g_XB </ci>
            </apply>
            <apply>
              <times/>
              <ci> f_01 </ci>
              <cn cellml:units="dimensionless"> 2.0 </cn>
              <ci> g_XB </ci>
              <cn cellml:units="dimensionless"> 3.0 </cn>
              <ci> g_XB </ci>
            </apply>
            <apply>
              <times/>
              <ci> f_01 </ci>
              <ci> f_12 </ci>
              <cn cellml:units="dimensionless"> 3.0 </cn>
              <ci> g_XB </ci>
            </apply>
            <apply>
              <times/>
              <ci> f_01 </ci>
              <ci> f_12 </ci>
              <ci> f_23 </ci>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply id="P3_max_calculation">
        <eq/>
        <ci> P3_max </ci>
        <apply>
          <divide/>
          <apply>
            <times/>
            <ci> f_01 </ci>
            <ci> f_12 </ci>
            <ci> f_23 </ci>
          </apply>
          <apply>
            <plus/>
            <apply>
              <times/>
              <ci> g_XB </ci>
              <cn cellml:units="dimensionless"> 2.0 </cn>
              <ci> g_XB </ci>
              <cn cellml:units="dimensionless"> 3.0 </cn>
              <ci> g_XB </ci>
            </apply>
            <apply>
              <times/>
              <ci> f_01 </ci>
              <cn cellml:units="dimensionless"> 2.0 </cn>
              <ci> g_XB </ci>
              <cn cellml:units="dimensionless"> 3.0 </cn>
              <ci> g_XB </ci>
            </apply>
            <apply>
              <times/>
              <ci> f_01 </ci>
              <ci> f_12 </ci>
              <cn cellml:units="dimensionless"> 3.0 </cn>
              <ci> g_XB </ci>
            </apply>
            <apply>
              <times/>
              <ci> f_01 </ci>
              <ci> f_12 </ci>
              <ci> f_23 </ci>
            </apply>
          </apply>
        </apply>
      </apply>
    </math>          
  </component>
  
  <group>
    <relationship_ref relationship="containment"/>
    <component_ref component="membrane">
      <component_ref component="fast_sodium_current">
        <component_ref component="fast_sodium_current_m_gate"/>
        <component_ref component="fast_sodium_current_h_gate"/>
        <component_ref component="fast_sodium_current_j_gate"/>
      </component_ref>
      <component_ref component="L_type_Ca_channel">
        <component_ref component="L_type_Ca_channel_y_gate"/>
      </component_ref>
      <component_ref component="time_dependent_potassium_current">
        <component_ref component="time_dependent_potassium_current_X_gate"/>
        <component_ref component="time_dependent_potassium_current_Xi_gate"/>
      </component_ref>
      <component_ref component="Na_Ca_exchanger"/>
      <component_ref component="time_independent_potassium_current">
        <component_ref component="time_independent_potassium_current_K1_gate"/>
      </component_ref> 
      <component_ref component="plateau_potassium_current"/>
      <component_ref component="sarcolemmal_calcium_pump"/>
      <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="intracellular_Ca_fluxes"/>
      <component_ref component="SERCA_pump"/>
      <component_ref component="RyR_channel_states"/>
      <component_ref component="intracellular_ionic_concentrations"/>
      <component_ref component="troponin"/>
      <component_ref component="tropomyosin_cross_bridges"/>
    </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_y_gate"/>
    </component_ref>
    <component_ref component="time_dependent_potassium_current">
      <component_ref component="time_dependent_potassium_current_X_gate"/>
      <component_ref component="time_dependent_potassium_current_Xi_gate"/>
    </component_ref>
    <component_ref component="time_independent_potassium_current">
      <component_ref component="time_independent_potassium_current_K1_gate"/>
    </component_ref>
  </group>
  
  <connection>
    <map_components component_2="environment" component_1="membrane"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>


  <connection>
    <map_components component_2="environment" component_1="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="time_dependent_potassium_current"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  
  <connection>
    <map_components component_2="environment" component_1="Na_Ca_exchanger"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  
  <connection>
    <map_components component_2="environment" component_1="time_independent_potassium_current"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  
  <connection>
    <map_components component_2="environment" component_1="plateau_potassium_current"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  
  <connection>
    <map_components component_2="environment" component_1="sarcolemmal_calcium_pump"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  
  <connection>
    <map_components component_2="environment" component_1="sodium_background_current"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  
  <connection>
    <map_components component_2="environment" component_1="calcium_background_current"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  
  <connection>
    <map_components component_2="environment" component_1="sodium_potassium_pump"/>
    <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="RyR_channel_states"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="environment" component_1="SERCA_pump"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="environment" component_1="intracellular_Ca_fluxes"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="environment" component_1="intracellular_ionic_concentrations"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="environment" component_1="troponin"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="environment" component_1="tropomyosin_cross_bridges"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="environment" component_1="force_computation"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <!--
    Several variables are passed between the "membrane" and its sub-components.
  -->
  
  <connection>
    <map_components component_2="fast_sodium_current" component_1="membrane"/>
    <map_variables variable_2="V" variable_1="V"/>
     <map_variables variable_2="i_Na" variable_1="i_Na"/>
     <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="membrane"/>
    <map_variables variable_2="V" variable_1="V"/>
    <map_variables variable_2="i_Ca" variable_1="i_Ca"/>
    <map_variables variable_2="i_Ca_K" variable_1="i_Ca_K"/>
    <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_dependent_potassium_current" component_1="membrane"/>
    <map_variables variable_2="V" variable_1="V"/>
    <map_variables variable_2="i_K" variable_1="i_K"/>
    <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="Na_Ca_exchanger" component_1="membrane"/>
    <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="membrane"/>
    <map_variables variable_2="V" variable_1="V"/>
    <map_variables variable_2="i_K1" variable_1="i_K1"/>
    <map_variables variable_2="R" variable_1="R"/>
    <map_variables variable_2="T" variable_1="T"/>
    <map_variables variable_2="F" variable_1="F"/>
  </connection>
  
  
  <connection>
    <map_components component_2="plateau_potassium_current" component_1="membrane"/>
    <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="membrane"/>
    <map_variables variable_2="i_p_Ca" variable_1="i_p_Ca"/>
  </connection>
  
  
  <connection>
    <map_components component_2="sodium_background_current" component_1="membrane"/>
    <map_variables variable_2="V" variable_1="V"/>
     <map_variables variable_2="i_Na_b" variable_1="i_Na_b"/>
  </connection>


  <connection>
    <map_components component_2="calcium_background_current" component_1="membrane"/>
    <map_variables variable_2="V" variable_1="V"/>
    <map_variables variable_2="i_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="membrane"/>
    <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="membrane"/>
    <map_variables variable_2="V" variable_1="V"/>
    <map_variables variable_2="i_ns_Ca" variable_1="i_ns_Ca"/>
    <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="intracellular_ionic_concentrations" component_1="membrane"/>
    <map_variables variable_2="F" variable_1="F"/>
  </connection>
  
  <!--  Several variables are passed between the sibling components. -->
  
  <connection>
    <map_components component_2="intracellular_ionic_concentrations" component_1="fast_sodium_current"/>
    <map_variables variable_2="Nai" variable_1="Nai"/>
    <map_variables variable_2="Nao" variable_1="Nao"/>
    <map_variables variable_2="i_Na" variable_1="i_Na"/>
  </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="intracellular_ionic_concentrations" component_1="L_type_Ca_channel"/>
    <map_variables variable_2="Ca_SS" variable_1="Ca_SS"/>
    <map_variables variable_2="i_Ca" variable_1="i_Ca"/>
    <map_variables variable_2="i_Ca_K" variable_1="i_Ca_K"/>
    <map_variables variable_2="Ki" variable_1="Ki"/>
    <map_variables variable_2="Ko" variable_1="Ko"/>
    <map_variables variable_2="Cao" variable_1="Cao"/>
  </connection>
  
  <connection>
    <map_components component_2="intracellular_ionic_concentrations" component_1="time_dependent_potassium_current"/>
    <map_variables variable_2="Ki" variable_1="Ki"/>
    <map_variables variable_2="Ko" variable_1="Ko"/>
    <map_variables variable_2="Nai" variable_1="Nai"/>
    <map_variables variable_2="Nao" variable_1="Nao"/>
    <map_variables variable_2="i_K" variable_1="i_K"/>
  </connection>
  
  
  <connection>
    <map_components component_2="intracellular_ionic_concentrations" component_1="Na_Ca_exchanger"/>
    <map_variables variable_2="Cao" variable_1="Cao"/>
    <map_variables variable_2="Nai" variable_1="Nai"/>
    <map_variables variable_2="Nao" variable_1="Nao"/>
    <map_variables variable_2="Cai" variable_1="Cai"/>
    <map_variables variable_2="i_NaCa" variable_1="i_NaCa"/>
  </connection>
  
  
  <connection>
    <map_components component_2="intracellular_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="time_independent_potassium_current" component_1="plateau_potassium_current"/>
    <map_variables variable_2="E_K1" variable_1="E_K1"/>
  </connection>
  
  
  <connection>
    <map_components component_2="intracellular_ionic_concentrations" component_1="plateau_potassium_current"/>
    <map_variables variable_2="i_Kp" variable_1="i_Kp"/>
  </connection>
  
  <connection>
    <map_components component_2="intracellular_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="intracellular_ionic_concentrations" 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="intracellular_ionic_concentrations" 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="intracellular_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="intracellular_ionic_concentrations" component_1="non_specific_calcium_activated_current"/>
    <map_variables variable_2="i_ns_Na" variable_1="i_ns_Na"/>
    <map_variables variable_2="i_ns_K" variable_1="i_ns_K"/>
    <map_variables variable_2="Cai" variable_1="Cai"/>
    <map_variables variable_2="Nai" variable_1="Nai"/>
    <map_variables variable_2="Nao" variable_1="Nao"/>
    <map_variables variable_2="Ki" variable_1="Ki"/>
    <map_variables variable_2="Ko" variable_1="Ko"/>
  </connection>
 
  
  <connection>
    <map_components component_2="cell_geometry" component_1="intracellular_ionic_concentrations"/>
    <map_variables variable_2="A_cap" variable_1="A_cap"/>
    <map_variables variable_2="V_myo" variable_1="V_myo"/>
    <map_variables variable_2="V_SS" variable_1="V_SS"/>
    <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="intracellular_Ca_fluxes" component_1="RyR_channel_states"/>
    <map_variables variable_2="P_O1" variable_1="P_O1"/>
    <map_variables variable_2="P_O2" variable_1="P_O2"/>
  </connection>
  
  <connection>
    <map_components component_2="intracellular_ionic_concentrations" component_1="RyR_channel_states"/>
    <map_variables variable_2="Ca_SS" variable_1="Ca_SS"/>
  </connection>
  
  
  <connection>
    <map_components component_2="SERCA_pump" component_1="intracellular_ionic_concentrations"/>
    <map_variables variable_2="J_up" variable_1="J_up"/>
    <map_variables variable_2="Ca_NSR" variable_1="Ca_NSR"/>
    <map_variables variable_2="Cai" variable_1="Cai"/>
  </connection>
  
  
  <connection>
    <map_components component_2="intracellular_ionic_concentrations" component_1="intracellular_Ca_fluxes"/>
    <map_variables variable_2="J_rel" variable_1="J_rel"/>
    <map_variables variable_2="J_trpn" variable_1="J_trpn"/>
    <map_variables variable_2="J_tr" variable_1="J_tr"/>
    <map_variables variable_2="J_xfer" variable_1="J_xfer"/>
    <map_variables variable_2="Ca_SS" variable_1="Ca_SS"/>
    <map_variables variable_2="Ca_JSR" variable_1="Ca_JSR"/>
    <map_variables variable_2="Ca_NSR" variable_1="Ca_NSR"/>
    <map_variables variable_2="Cai" variable_1="Cai"/>
  </connection>
  
  <connection>
    <map_components component_2="troponin" component_1="intracellular_Ca_fluxes"/>
    <map_variables variable_2="HTRPNCa" variable_1="HTRPNCa"/>
    <map_variables variable_2="HTRPN_tot" variable_1="HTRPN_tot"/>
    <map_variables variable_2="LTRPNCa" variable_1="LTRPNCa"/>
    <map_variables variable_2="LTRPN_tot" variable_1="LTRPN_tot"/>
    <map_variables variable_2="k_htrpn_minus" variable_1="k_htrpn_minus"/>
    <map_variables variable_2="k_htrpn_plus" variable_1="k_htrpn_plus"/>
    <map_variables variable_2="k_ltrpn_minus" variable_1="k_ltrpn_minus"/>
    <map_variables variable_2="k_ltrpn_plus" variable_1="k_ltrpn_plus"/>
  </connection>
  
  <connection>
    <map_components component_2="troponin" component_1="force_computation"/>
    <map_variables variable_2="Force_norm" variable_1="Force_norm"/>
  </connection>
  
  <connection>
    <map_components component_2="troponin" component_1="intracellular_ionic_concentrations"/>
    <map_variables variable_2="Cai" variable_1="Cai"/>
  </connection>
  
  <connection>
    <map_components component_2="tropomyosin_cross_bridges" component_1="troponin"/>
    <map_variables variable_2="k_ltrpn_minus" variable_1="k_ltrpn_minus"/>
    <map_variables variable_2="k_ltrpn_plus" variable_1="k_ltrpn_plus"/>
    <map_variables variable_2="LTRPNCa" variable_1="LTRPNCa"/>
    <map_variables variable_2="LTRPN_tot" variable_1="LTRPN_tot"/>      
  </connection>
  
  <connection>
    <map_components component_2="tropomyosin_cross_bridges" component_1="force_computation"/>
    <map_variables variable_2="f_01" variable_1="f_01"/>
    <map_variables variable_2="f_12" variable_1="f_12"/>
    <map_variables variable_2="f_23" variable_1="f_23"/>
    <map_variables variable_2="g_XB" variable_1="g_XB"/>
    <map_variables variable_2="SL" variable_1="SL"/>
    <map_variables variable_2="P1" variable_1="P1"/>
    <map_variables variable_2="P2" variable_1="P2"/>
    <map_variables variable_2="N1" variable_1="N1"/>
    <map_variables variable_2="P3" variable_1="P3"/>
  </connection>
  
  <!--  
    Various variables are passed between parent components and their            
    encapsulated gates.
  -->
  
  <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_y_gate" component_1="L_type_Ca_channel"/>
    <map_variables variable_2="y" variable_1="y"/>
    <map_variables variable_2="time" variable_1="time"/>
    <map_variables variable_2="V" variable_1="V"/>
  </connection>
  
  <connection>
    <map_components component_2="time_dependent_potassium_current_X_gate" component_1="time_dependent_potassium_current"/>
    <map_variables variable_2="X" variable_1="X"/>
    <map_variables variable_2="time" variable_1="time"/>
    <map_variables variable_2="V" variable_1="V"/>
  </connection>
  
  <connection>
    <map_components component_2="time_dependent_potassium_current_Xi_gate" component_1="time_dependent_potassium_current"/>
    <map_variables variable_2="Xi" variable_1="Xi"/>
    <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="time" variable_1="time"/>
    <map_variables variable_2="V" variable_1="V"/>
    <map_variables variable_2="E_K1" variable_1="E_K1"/>
    <map_variables variable_2="K1_infinity" variable_1="K1_infinity"/>
  </connection>
    






<rdf:RDF>
  <rdf:Bag rdf:about="rdf:#3327930a-1af9-4c82-a86d-70774454de3d">
    <rdf:li>calcium dynamics</rdf:li>
    <rdf:li>cardiac myocyte</rdf:li>
    <rdf:li>Cardiac Myocyte</rdf:li>
    <rdf:li>electrophysiology</rdf:li>
    <rdf:li>myofilament mechanics</rdf:li>
    <rdf:li>cardiac</rdf:li>
  </rdf:Bag>
  <rdf:Seq rdf:about="rdf:#citationAuthorsSeq">
    <rdf:li rdf:resource="rdf:#author1Vcard"/>
    <rdf:li rdf:resource="rdf:#author2Vcard"/>
    <rdf:li rdf:resource="rdf:#author3Vcard"/>
  </rdf:Seq>
  <rdf:Description rdf:about="rdf:#cf4e6901-4481-4416-b461-cc264e0a9d59">
    <dcterms:W3CDTF>2002-06-24</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#8dc79e2d-d733-46bf-9823-6cba8b2c2588">
    <dcterms:W3CDTF>2003-04-09</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#cfa41065-f3e9-41d7-b2ae-a47cbc211110">
    <dcterms:modified rdf:resource="rdf:#735565b1-a521-4d7c-b198-09bd6d8f0b7a"/>
    <rdf:value>
          Added more metadata.
        </rdf:value>
    <cmeta:modifier rdf:resource="rdf:#a165ac20-3a7c-4713-8edf-bcb47e57a145"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#7d005e9c-1569-4f07-b235-bc637c2c3012">
    <bqs:subject_type>keyword</bqs:subject_type>
    <rdf:value rdf:resource="rdf:#3327930a-1af9-4c82-a86d-70774454de3d"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#6de2c6e3-a9e5-4924-9298-8704b91b2c4e">
    <dc:title>American Journal of Physiology:  Heart and Circulatory Physiology</dc:title>
  </rdf:Description>
  <rdf:Description rdf:about="#rice_jafri_winslow_2000_version02">
    <dc:title>
        The Rice et al's 2000 detailed model of ventricular cardiac 
        myocytes
      </dc:title>
    <cmeta:bio_entity>Cardiac Myocyte</cmeta:bio_entity>
    <cmeta:comment rdf:resource="rdf:#fdac2647-6378-4d56-96e0-1b7906026c31"/>
    <bqs:reference rdf:resource="rdf:#b3a3df16-9fbf-4039-abc9-e2b61a621801"/>
    <bqs:reference rdf:resource="rdf:#dd6972d1-9747-47f0-9bb1-26d908bbb8e2"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#author3VcardN">
    <vCard:Given>Raimond</vCard:Given>
    <vCard:Family>Winslow</vCard:Family>
    <vCard:Other>L</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="">
    <dc:publisher>
        The University of Auckland, Bioengineering Institute
      </dc:publisher>
    <cmeta:modification rdf:resource="rdf:#24eb0b3f-f381-4ccf-ac5d-70644a657622"/>
    <cmeta:modification rdf:resource="rdf:#cfa41065-f3e9-41d7-b2ae-a47cbc211110"/>
    <dcterms:created rdf:resource="rdf:#cf4e6901-4481-4416-b461-cc264e0a9d59"/>
    <dc:creator rdf:resource="rdf:#2f43be39-8e3c-4605-9fde-02dbe2610013"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#4ca4b792-d074-4e5c-b3f9-041f9bd45d57">
    <vCard:Given>Catherine</vCard:Given>
    <vCard:Family>Lloyd</vCard:Family>
    <vCard:Other>May</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#2f43be39-8e3c-4605-9fde-02dbe2610013">
    <vCard:ORG rdf:resource="rdf:#14d2d51a-cf53-434a-b830-54a20a4b8ac3"/>
    <vCard:EMAIL rdf:resource="rdf:#9260fb1a-dd92-4fb1-b17e-d6545f37316d"/>
    <vCard:N rdf:resource="rdf:#9b89d4f9-1535-466b-8e1e-45d105a5df0e"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#author2Vcard">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#author2VcardN"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#dd6972d1-9747-47f0-9bb1-26d908bbb8e2">
    <dc:subject rdf:resource="rdf:#7d005e9c-1569-4f07-b235-bc637c2c3012"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#author2VcardN">
    <vCard:Given>M</vCard:Given>
    <vCard:Family>Jafri</vCard:Family>
    <vCard:Other>Saleet</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#14d2d51a-cf53-434a-b830-54a20a4b8ac3">
    <vCard:Orgname>The University of Auckland</vCard:Orgname>
    <vCard:Orgunit>The Bioengineering Institute</vCard:Orgunit>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#1be093f0-0715-463d-8a1b-938f856c7139">
    <vCard:FN>Catherine Lloyd</vCard:FN>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#045b84e5-d335-4460-a7e8-db2f57cffc3e">
    <dcterms:W3CDTF>2000-03-01</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#74b49ad5-99d9-486c-a023-ed1073ab0991">
    <vCard:N rdf:resource="rdf:#59c07677-3e2a-4a4d-ada1-48f44959a4f4"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#a165ac20-3a7c-4713-8edf-bcb47e57a145">
    <vCard:N rdf:resource="rdf:#4ca4b792-d074-4e5c-b3f9-041f9bd45d57"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#fb91e637-d3e7-42cc-979f-b922d3b981c0">
    <dc:creator rdf:resource="rdf:#citationAuthorsSeq"/>
    <dc:title>Modeling short-term interval-force relations in cardiac muscle</dc:title>
    <bqs:volume>278</bqs:volume>
    <bqs:first_page>H913</bqs:first_page>
    <bqs:Journal rdf:resource="rdf:#6de2c6e3-a9e5-4924-9298-8704b91b2c4e"/>
    <dcterms:issued rdf:resource="rdf:#045b84e5-d335-4460-a7e8-db2f57cffc3e"/>
    <bqs:last_page>H931</bqs:last_page>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#59c07677-3e2a-4a4d-ada1-48f44959a4f4">
    <vCard:Given>Autumn</vCard:Given>
    <vCard:Family>Cuellar</vCard:Family>
    <vCard:Other>A</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#b3a3df16-9fbf-4039-abc9-e2b61a621801">
    <bqs:Pubmed_id>10710361</bqs:Pubmed_id>
    <bqs:JournalArticle rdf:resource="rdf:#fb91e637-d3e7-42cc-979f-b922d3b981c0"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#fdac2647-6378-4d56-96e0-1b7906026c31">
    <dc:creator rdf:resource="rdf:#1be093f0-0715-463d-8a1b-938f856c7139"/>
    <rdf:value>This is the CellML description of Rice et al's 2000 detailed model of ventricular cardiac myocytes. It simulates action potentials, Ca2+-handling mechanisms and isometric force generation by the myofilaments.</rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#author1Vcard">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#author1VcardN"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#author1VcardN">
    <vCard:Given>John</vCard:Given>
    <vCard:Family>Rice</vCard:Family>
    <vCard:Other>Jeremy</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#9260fb1a-dd92-4fb1-b17e-d6545f37316d">
    <rdf:type rdf:resource="http://imc.org/vCard/3.0#internet"/>
    <rdf:value>c.lloyd@auckland.ac.nz</rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#author3Vcard">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#author3VcardN"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#9b89d4f9-1535-466b-8e1e-45d105a5df0e">
    <vCard:Given>Catherine</vCard:Given>
    <vCard:Family>Lloyd</vCard:Family>
    <vCard:Other>May</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#735565b1-a521-4d7c-b198-09bd6d8f0b7a">
    <dcterms:W3CDTF>2002-07-18</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#24eb0b3f-f381-4ccf-ac5d-70644a657622">
    <dcterms:modified rdf:resource="rdf:#8dc79e2d-d733-46bf-9823-6cba8b2c2588"/>
    <rdf:value>
          Added publication date information.
        </rdf:value>
    <cmeta:modifier rdf:resource="rdf:#74b49ad5-99d9-486c-a023-ed1073ab0991"/>
  </rdf:Description>
</rdf:RDF>
</model>