Location: Noble, Varghese, Kohl, Noble, 1998 @ 0be3bf4e376d / noble_varghese_kohl_noble_1998.cellml

Author:
pmr2.import <nobody@models.cellml.org>
Date:
2006-09-21 04:26:52+12:00
Desc:
committing version05 of noble_varghese_kohl_noble_1998
Permanent Source URI:
https://staging.physiomeproject.org/workspace/noble_varghese_kohl_noble_1998/rawfile/0be3bf4e376ddb8f541244b616c564a27b1028b9/noble_varghese_kohl_noble_1998.cellml

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

CREATED :  18th December 2001

LAST MODIFIED : 22nd July 2002

AUTHOR :  Catherine Lloyd
          Department of Engineering Science
          The University of Auckland
          
MODEL STATUS :  This model conforms to the CellML 1.0 Specification released on
10th August 2001, and the 16/1/02 CellML Metadata 1.0 Specification.

DESCRIPTION :  This file contains a CellML description of Noble's 1998 improved
guinea-pig ventricular cell model.

CHANGES:  
  04/01/2002 - CML - Altered some of the connections.
  21/01/2002 - AAC - Updated metadata to conform to the 16/1/02 CellML Metadata
                     1.0 Specification.
  25/02/2002 - CML - Corrected i_K1 and i_K_ACh calculations.
  26/02/2002 - CML - Corrected several more equations.
  01/03/2002 - CML - Corrected units.
  06/05/2002 - CML - Added some initial values.  Removed the blocked h_gate 
                     component from the fast sodium current as it belongs to a
                     separate model which considers the effect of drugs.
  22/07/2002 - CML - Added more metadata.
  08/02/2003 - DPN - Added encapsulation hierarchy for transient_outward_current.
--><!--
  This version of the Noble 1998 model is extensively modified from the
  version available in the CellML model repository on cellml.org. The default
  variable values are all included in this model description, although this
  should change when we go to CellML 1.1.
--><model xmlns="http://www.cellml.org/cellml/1.0#" xmlns:cmeta="http://www.cellml.org/metadata/1.0#" xmlns:dc="http://purl.org/dc/elements/1.1/" xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#" xmlns:bqs="http://www.cellml.org/bqs/1.0#" xmlns:cellml="http://www.cellml.org/cellml/1.0#" xmlns:dcterms="http://purl.org/dc/terms/" xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" cmeta:id="noble_varghese_kohl_noble_1998_version05" name="noble_varghese_kohl_noble_1998_version05">

<documentation xmlns="http://cellml.org/tmp-documentation">
<article>
  <articleinfo>
  <title>Noble Ventricular Cell Model 1998</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 version of the model has been curated by Penny Noble using
      COR and is known to read in to JSim.
    </para>
  </section>
  <sect1 id="sec_structure">
<title>Model Structure</title>

<para>In their 1998 paper, Denis Noble, Anthony Varghese, Peter Kohl and Penny Noble extended the guinea-pig ventricular cell model (originally developed in 1991 by Noble <emphasis>et al</emphasis>) to include accumulation and depletion of calcium in a diadic space between the sarcolemma and the sarcoplasmic reticulum, where the majority of calcium-induced calcium release is triggered.  The model also includes rapid (i_Kr) and slow (i_Ks) components of the delayed rectifier current, and length- and tension-dependent changes in mechanical and electrophysiological processes have been incorporated (see the figure below).  Drug-receptor interactions have also started to be modeled using the sodium channel as an example.  The drugs used have the effect of blocking the inactivation h-gate of the fast sodium channel.
</para>

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

<para>
<ulink url="http://www.pulsus.com/CARDIOL/14_01/nobl_ed.htm">Improved guinea-pig ventricular cell model incorporating a diadic space, I<subscript>Kr</subscript> and I<subscript>Ks</subscript>, and length- and tension-dependent processes</ulink>, Denis Noble, Anthony Varghese, Peter Kohl and Penelope Noble, 1998, <ulink url="http://www.pulsus.com/CARDIOL/home.htm">
            <emphasis>Can J Cardiol</emphasis>
          </ulink>, 14, 123-134.  <ulink url="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?cmd=Retrieve&amp;db=PubMed&amp;list_uids=9487284&amp;dopt=Abstract">PubMed ID: 9487284</ulink>
</para>

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

</sect1>
</article>
</documentation>
    
  
  
  <!--
    Standard CMISS units
  -->
  <units name="mm">
    <unit units="metre" prefix="milli"/>
  </units>
  <units name="ng">
    <unit units="gram" prefix="nano"/>
  </units>
  <units name="ms">
    <unit units="second" prefix="milli"/>
  </units>
  <units name="uA">
    <unit units="ampere" prefix="micro"/>
  </units>
  <units name="K">
    <unit units="kelvin"/>
  </units>
  <units name="nmol">
    <unit units="mole" prefix="nano"/>
  </units>
  
  <units name="mm2">
    <unit units="mm" exponent="2"/>
  </units>
  <units name="mm3">
    <unit units="mm" exponent="3"/>
  </units>
  <units name="mV">
    <unit units="volt" prefix="milli"/>
  </units>
  <units name="mS">
    <unit units="siemens" prefix="milli"/>
  </units>
  <units name="mS_per_mmsq">
    <unit units="mS"/>
    <unit units="mm" exponent="-2"/>
  </units>
  <units name="mS_per_mm">
    <unit units="mS"/>
    <unit units="mm" exponent="-1"/>
  </units>
  <units name="uA_per_mmsq">
    <unit units="uA"/>
    <unit units="mm" exponent="-2"/>
  </units>
  <units name="uA_per_mmcu">
    <unit units="uA"/>
    <unit units="mm" exponent="-3"/>
  </units>
  <units name="nC">
    <unit units="coulomb" prefix="nano"/>
  </units>
  <units name="uF">
    <unit units="farad" prefix="micro"/>
  </units>
  <units name="uF_per_mmsq">
    <unit units="uF"/>
    <unit units="mm" exponent="-2"/>
  </units>
  <units name="mM">
    <unit units="nmol"/>
    <unit units="mm" exponent="-3"/>
  </units>
  <units name="mM_per_ms">
    <unit units="mM"/>
    <unit units="ms" exponent="-1"/>
  </units>
  <units name="nN">
    <unit units="newton" prefix="nano"/>
  </units>
  <units name="mPa">
    <unit units="pascal" prefix="milli"/>
  </units>
  <units name="pJ">
    <unit units="joule" prefix="pico"/>
  </units>
  <units name="nW">
    <unit units="watt" prefix="nano"/>
  </units>
  
  <units name="per_ms">
    <unit units="ms" exponent="-1"/>
  </units>
  <units name="per_mm">
    <unit units="mm" exponent="-1"/>
  </units>
  <units name="per_mV">
    <unit units="mV" exponent="-1"/>
  </units>
  <units name="per_mV_per_ms">
    <unit units="mV" exponent="-1"/>
    <unit units="per_ms"/>
  </units>
  <units name="per_mM_per_ms">
    <unit units="mM" exponent="-1"/>
    <unit units="per_ms"/>
  </units>
  <units name="uA_per_mmsq_per_mM">
    <unit units="uA"/>
    <unit units="mm" exponent="-2"/>
    <unit units="mM" exponent="-1"/>
  </units>
  
  <units name="gas_constant">
    <unit units="pJ"/>
    <unit units="nmol" exponent="-1"/>
    <unit units="K" exponent="-1"/>
  </units>
  
  <units name="faradays_constant">
    <unit units="nC"/>
    <unit units="nmol" exponent="-1"/>
  </units>
  
  <component name="environment" cmeta:id="environment">
    
    <variable units="ms" public_interface="out" name="time"/>
  </component>

  <component name="membrane" cmeta:id="membrane">
    

    <variable units="mV" public_interface="out" name="V" initial_value="-92.499"/>
    <variable units="gas_constant" public_interface="out" name="R" initial_value="8314.472"/>
    <variable units="K" public_interface="out" name="T" initial_value="310.0"/>
    <variable units="faradays_constant" public_interface="out" name="F" initial_value="96485.0"/>
    <variable units="per_mm" public_interface="out" name="Am" initial_value="193.69"/>
    
    <variable units="uF_per_mmsq" name="Cm" initial_value="0.014651"/>
    
    <variable units="ms" public_interface="in" name="time"/>
    <variable units="uA_per_mmsq" public_interface="in" name="i_K1"/>
    <variable units="uA_per_mmsq" public_interface="in" name="i_to"/>
    <variable units="uA_per_mmsq" public_interface="in" name="i_K"/>
    <variable units="uA_per_mmsq" public_interface="in" name="i_K_ATP"/>
    <variable units="uA_per_mmsq" public_interface="in" name="i_Ca_L_K"/>
    <variable units="uA_per_mmsq" public_interface="in" name="i_b_K"/>
    <variable units="uA_per_mmsq" public_interface="in" name="i_NaK"/>
    <variable units="uA_per_mmsq" public_interface="in" name="i_Na"/>
    <variable units="uA_per_mmsq" public_interface="in" name="i_b_Na"/>
    <variable units="uA_per_mmsq" public_interface="in" name="i_p_Na"/>
    <variable units="uA_per_mmsq" public_interface="in" name="i_Ca_L_Na"/>
    <variable units="uA_per_mmsq" public_interface="in" name="i_NaCa"/>
    <variable units="uA_per_mmsq" public_interface="in" name="i_NaCa_ds"/>
    <variable units="uA_per_mmsq" public_interface="in" name="i_Ca_L_Ca"/>
    <variable units="uA_per_mmsq" public_interface="in" name="i_Ca_L_Ca_ds"/>
    <variable units="uA_per_mmsq" public_interface="in" name="i_Ca_L_K_ds"/>
    <variable units="uA_per_mmsq" public_interface="in" name="i_Ca_L_Na_ds"/>
    <variable units="uA_per_mmsq" public_interface="in" name="i_b_Ca"/>
    <variable units="uA_per_mmsq" public_interface="in" name="i_Na_stretch"/>
    <variable units="uA_per_mmsq" public_interface="in" name="i_K_stretch"/>
    <variable units="uA_per_mmsq" public_interface="in" name="i_Ca_stretch"/>
    <variable units="uA_per_mmsq" public_interface="in" name="i_Ns_stretch"/>
    <variable units="uA_per_mmsq" public_interface="in" name="i_An_stretch"/>
    <variable units="uA_per_mmsq" public_interface="in" name="i_K_ACh"/>
    
    <variable units="uA_per_mmsq" name="IStim"/>
    <variable units="ms" name="stimPeriod" initial_value="750"/>
    <variable units="ms" name="stimDuration" initial_value="1"/>
    <variable units="uA_per_mmcu" name="stimCurrent" initial_value="250"/> 

    
    <math xmlns="http://www.w3.org/1998/Math/MathML" cmeta:id="IStim_for_cmiss_eq">
      <apply id="stimulus_calculation">
        <eq/>
        <ci>IStim</ci>
        <piecewise>
          <piece>
            <apply>
              <divide/>
              <ci>stimCurrent</ci>
              <ci>Am</ci>
            </apply>
            <apply>
              <lt/>
              <apply>
                <rem/>
                <ci>time</ci>
                <ci>stimPeriod</ci>
              </apply>
              <ci>stimDuration</ci>
            </apply>
          </piece>
          <otherwise>
            <cn cellml:units="uA_per_mmsq">0.0</cn>
          </otherwise>
        </piecewise>
      </apply>
    </math>

    
    <math xmlns="http://www.w3.org/1998/Math/MathML" cmeta:id="membrane_voltage_diff_eq">
      <apply id="membrane_voltage_diff_eq">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> V </ci>
        </apply>
        <apply>
          <divide/>
          <apply>
            <minus/>
            <ci>IStim</ci>
            <apply>
              <plus/>
              <ci> i_K1 </ci>
              <ci> i_to </ci>
              <ci> i_K </ci>
              <ci> i_K_ATP </ci>
              <ci> i_b_K </ci>
              <ci> i_NaK </ci>
              <ci> i_Na </ci>
              <ci> i_b_Na </ci>
              <ci> i_p_Na </ci>
              <ci> i_NaCa </ci>
              <ci> i_NaCa_ds </ci>
              <ci> i_Ca_L_Ca </ci>
              <ci> i_Ca_L_Ca_ds </ci>
              <ci> i_Ca_L_K </ci>
              <ci> i_Ca_L_K_ds </ci>
              <ci> i_Ca_L_Na </ci>
              <ci> i_Ca_L_Na_ds </ci>
              <ci> i_b_Ca </ci>
              <ci> i_Na_stretch </ci>
              <ci> i_K_stretch </ci>
              <ci> i_Ca_stretch </ci>
              <ci> i_Ns_stretch </ci>
              <ci> i_An_stretch </ci>
              <ci> i_K_ACh </ci>
            </apply>
            <!--plus-->
          </apply>
          <!--minus-->
          <ci>Cm</ci>
        </apply>
        <!--divide-->
      </apply>
      <!--eq-->
    </math>
  </component>
  
  <component name="reversal_potentials" cmeta:id="reversal_potentials">
    

    <variable units="mV" public_interface="out" name="E_Na"/>
    <variable units="mV" public_interface="out" name="E_K"/>
    <variable units="mV" public_interface="out" name="E_Ks"/>
    <variable units="mV" public_interface="out" name="E_Ca"/>
    <variable units="mV" public_interface="out" name="E_mh"/>
    
    <variable units="dimensionless" name="P_kna" initial_value="0.03"/>
    <variable units="dimensionless" name="P_nak" initial_value="0.12"/>
    
    <variable units="mM" public_interface="in" name="K_o"/>
    <variable units="mM" public_interface="in" name="Na_o"/>
    <variable units="mM" public_interface="in" name="K_i"/>
    <variable units="mM" public_interface="in" name="Na_i"/>
    <variable units="mM" public_interface="in" name="Ca_o"/>
    <variable units="mM" public_interface="in" name="Ca_i"/>
    <variable units="gas_constant" public_interface="in" name="R"/>
    <variable units="faradays_constant" public_interface="in" name="F"/>
    <variable units="K" public_interface="in" name="T"/>
    <variable units="ms" public_interface="in" name="time"/>
    
    
    <math xmlns="http://www.w3.org/1998/Math/MathML" cmeta:id="E_Na_calculation_eq">
      <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> Na_o </ci>
              <ci> Na_i </ci>
            </apply>
          </apply>
        </apply>
      </apply>
    </math>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML" cmeta:id="E_K_calculation_eq">
      <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/>
              <ci> K_o </ci>
              <ci> K_i </ci>
            </apply>
          </apply>
        </apply>
      </apply>
    </math>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML" cmeta:id="E_Ks_calculation_eq">
      <apply id="E_Ks_calculation">
        <eq/>
        <ci> E_Ks </ci>
        <apply>
          <times/>
          <apply>
            <divide/>
            <apply>
              <times/>
              <ci> R </ci>
              <ci> T </ci>
            </apply>
            <ci> F </ci>
          </apply>
          <apply>
            <ln/>
            <apply>
              <divide/>
              <apply>
                <plus/>
                <ci> K_o </ci>
                <apply>
                  <times/>
                  <ci> P_kna </ci>
                  <ci> Na_o </ci>
                </apply>
              </apply>
              <apply>
                <plus/>
                <ci> K_i </ci>
                <apply>
                  <times/>
                  <ci> P_kna </ci>
                  <ci> Na_i </ci>
                </apply>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
    </math>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML" cmeta:id="E_Ca_calculation_eq">
      <apply id="E_Ca_calculation">
        <eq/>
        <ci> E_Ca </ci>
        <apply>
          <times/>
          <cn cellml:units="dimensionless">0.5</cn>
          <apply>
            <divide/>
            <apply>
              <times/>
              <ci> R </ci>
              <ci> T </ci>
            </apply>
            <ci> F </ci>
          </apply>
          <apply>
            <ln/>
            <apply>
              <divide/>
              <ci> Ca_o </ci>
              <ci> Ca_i </ci>
            </apply>
          </apply>
        </apply>
      </apply>
    </math>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML" cmeta:id="E_mh_calculation_eq">
      <apply id="E_mh_calculation">
        <eq/>
        <ci> E_mh </ci>
        <apply>
          <times/>
          <apply>
            <divide/>
            <apply>
              <times/>
              <ci> R </ci>
              <ci> T </ci>
            </apply>
            <ci> F </ci>
          </apply>
          <apply>
            <ln/>
            <apply>
              <divide/>
              <apply>
                <plus/>
                <ci> Na_o </ci>
                <apply>
                  <times/>
                  <ci> P_nak </ci>
                  <ci> K_o </ci>
                </apply>
              </apply>
              <apply>
                <plus/>
                <ci> Na_i </ci>
                <apply>
                  <times/>
                  <ci> P_nak </ci>
                  <ci> K_i </ci>
                </apply>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
      
  <component name="time_independent_potassium_current" cmeta:id="time_independent_potassium_current">
    
    
    <variable units="uA_per_mmsq" public_interface="out" name="i_K1"/>
    <variable units="mM" public_interface="out" name="K_mk1" initial_value="10.0"/> 
    
    <variable units="mS_per_mmsq" name="g_K1" initial_value="77.11e-3"/>
    
    <variable units="mM" public_interface="in" name="K_o"/>
    <variable units="mV" public_interface="in" name="V"/>
    <variable units="mV" public_interface="in" name="E_K"/>
    <variable units="gas_constant" public_interface="in" name="R"/>
    <variable units="faradays_constant" public_interface="in" name="F"/>
    <variable units="K" public_interface="in" name="T"/>
    <variable units="ms" public_interface="in" name="time"/>
    
    
    <math xmlns="http://www.w3.org/1998/Math/MathML" cmeta:id="i_K1_calculation_eq">
      <apply id="i_K1_calculation">
        <eq/>
        <ci> i_K1 </ci>
        <apply>
          <times/>
          <ci> g_K1 </ci>
          <apply>
            <divide/>
            <ci> K_o </ci>
            <apply>
              <plus/>
              <ci> K_o </ci>
              <ci> K_mk1 </ci>
            </apply>
          </apply>
          <apply>
            <divide/>
            <apply>
              <minus/>
              <ci> V </ci>
              <ci> E_K </ci>
            </apply>
            <apply>
              <plus/>
              <cn cellml:units="dimensionless"> 1.0 </cn>
              <apply>
                <exp/>
                <apply>
                  <divide/>
                  <apply>
                    <times/>
                    <cn cellml:units="dimensionless"> 1.25 </cn>
                    <apply>
                      <minus/>
                      <ci> V </ci>
                      <apply>
                        <plus/>
                        <ci> E_K </ci>
                        <cn cellml:units="mV"> 10.0 </cn>
                      </apply>
                    </apply>
                  </apply>
                  <apply>
                    <divide/>
                    <apply>
                      <times/>
                      <ci> R </ci>
                      <ci> T </ci>
                    </apply>
                    <ci> F </ci>
                  </apply>
                </apply>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="total_potassium_current" cmeta:id="total_potassium_current">
    
    <variable units="uA_per_mmsq" public_interface="out" name="i_K"/>
    <variable units="uA_per_mmsq" public_interface="in" name="i_Kr"/>
    <variable units="uA_per_mmsq" public_interface="in" name="i_Ks"/>
    <variable units="uA_per_mmsq" public_interface="in" name="i_KNa"/> 
    <variable units="ms" public_interface="in" name="time"/>
    
    
    <math xmlns="http://www.w3.org/1998/Math/MathML" cmeta:id="i_K_calculation_eq">
      <apply id="i_K_calculation">
        <eq/>
        <ci> i_K </ci>
        <apply>
          <plus/>
          <ci> i_Kr </ci>
          <ci> i_Ks </ci>
          <ci> i_KNa </ci>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="rapid_time_dependent_potassium_current" cmeta:id="rapid_time_dependent_potassium_current">
    

    <variable units="uA_per_mmsq" public_interface="out" name="i_Kr"/>
    
    <variable units="mS_per_mmsq" name="g_Kr1" initial_value="0.32386e-3"/>
    <variable units="mS_per_mmsq" name="g_Kr2" initial_value="0.20049e-3"/>
    
    <variable units="mV" public_interface="in" name="E_K"/>
    <variable units="ms" public_interface="in" private_interface="out" name="time"/>
    <variable units="mV" public_interface="in" private_interface="out" name="V"/>
    
    <variable units="dimensionless" private_interface="in" name="Xr1"/>
    <variable units="dimensionless" private_interface="in" name="Xr2"/>
    
    
    <math xmlns="http://www.w3.org/1998/Math/MathML" cmeta:id="i_Kr_calculation_eq">
      <apply id="i_Kr_calculation">
        <eq/>
        <ci> i_Kr </ci>
        <apply>
          <times/>
          <apply>
            <plus/>
            <apply>
              <times/>
              <ci> g_Kr1 </ci>
              <ci> Xr1 </ci>
            </apply>
            <apply>
              <times/>
              <ci> g_Kr2 </ci>
              <ci> Xr2 </ci>
            </apply>
          </apply>
          <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="mV"> 9.0 </cn>
                  </apply>
                  <cn cellml:units="mV"> 22.4 </cn>
                </apply>
              </apply>
            </apply>
          </apply>
          <apply>
            <minus/>
            <ci> V </ci>
            <ci> E_K </ci>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="rapid_time_dependent_potassium_current_Xr1_gate" cmeta:id="rapid_time_dependent_potassium_current_Xr1_gate">
    

    <variable units="dimensionless" public_interface="out" name="Xr1" initial_value="0.01428"/>
    <variable units="mV" public_interface="in" name="V"/>
    <variable units="ms" public_interface="in" name="time"/>
    
    
    <math xmlns="http://www.w3.org/1998/Math/MathML" cmeta:id="Xr1_diff_eq">
      <apply id="Xr1_diff">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> Xr1 </ci>
        </apply>
        <apply>
          <minus/>
          <apply>
            <times/>
            <apply>
              <divide/>
              <cn cellml:units="per_ms"> 0.05 </cn>
              <apply>
                <plus/>
                <cn cellml:units="dimensionless"> 1.0 </cn>
                <apply>
                  <exp/>
                  <apply>
                    <divide/>
                    <apply>
                      <minus/>
                      <apply>
                        <minus/>
                        <ci> V </ci>
                        <cn cellml:units="mV"> 5.0 </cn>
                      </apply>
                    </apply>
                    <cn cellml:units="mV"> 9.0 </cn>
                  </apply>
                </apply>
              </apply>
            </apply>
            <apply>
              <minus/>
              <cn cellml:units="dimensionless"> 1.0 </cn>
              <ci> Xr1 </ci>
            </apply>
          </apply>
          <apply>
            <times/>
            <cn cellml:units="per_ms"> 0.00005 </cn>
            <apply>
              <exp/>
              <apply>
                <minus/>
                <apply>
                  <divide/>
                  <apply>
                    <minus/>
                    <ci> V </ci>
                    <cn cellml:units="mV"> 20.0 </cn>
                  </apply>
                  <cn cellml:units="mV"> 15.0 </cn>
                </apply>
              </apply>
            </apply>
            <ci> Xr1 </ci>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="rapid_time_dependent_potassium_current_Xr2_gate" cmeta:id="rapid_time_dependent_potassium_current_Xr2_gate">
    

    <variable units="dimensionless" public_interface="out" name="Xr2" initial_value="0.5"/>
    <variable units="mV" public_interface="in" name="V"/>
    <variable units="ms" public_interface="in" name="time"/>
    
    
    <math xmlns="http://www.w3.org/1998/Math/MathML" cmeta:id="Xr2_diff_eq">
      <apply id="Xr2_diff">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> Xr2 </ci>
        </apply>
        <apply>
          <minus/>
          <apply>
            <times/>
            <apply>
              <divide/>
              <cn cellml:units="per_ms"> 0.05 </cn>
              <apply>
                <plus/>
                <cn cellml:units="dimensionless"> 1.0 </cn>
                <apply>
                  <exp/>
                  <apply>
                    <divide/>
                    <apply>
                      <minus/>
                      <apply>
                        <minus/>
                        <ci> V </ci>
                        <cn cellml:units="mV"> 5.0 </cn>
                      </apply>
                    </apply>
                    <cn cellml:units="mV"> 9.0 </cn>
                  </apply>
                </apply>
              </apply>
            </apply>
            <apply>
              <minus/>
              <cn cellml:units="dimensionless"> 1.0 </cn>
              <ci> Xr2 </ci>
            </apply>
          </apply>
          <apply>
            <times/>
            <cn cellml:units="per_ms"> 0.0004 </cn>
            <apply>
              <exp/>
              <apply>
                <minus/>
                <apply>
                  <power/>
                  <apply>
                    <divide/>
                    <apply>
                      <plus/>
                      <ci> V </ci>
                      <cn cellml:units="mV"> 30.0 </cn>
                    </apply>
                    <cn cellml:units="mV"> 30.0 </cn>
                  </apply>
                  <cn cellml:units="dimensionless"> 3.0 </cn>
                </apply>
              </apply>
            </apply>
            <ci> Xr2 </ci>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="slow_time_dependent_potassium_current" cmeta:id="slow_time_dependent_potassium_current">
    

    <variable units="uA_per_mmsq" public_interface="out" name="i_Ks"/>
    
    <variable units="mS_per_mmsq" name="g_Ks" initial_value="0.40097e-3"/>
    
    <variable units="mV" public_interface="in" name="E_Ks"/>
    <variable units="ms" public_interface="in" private_interface="out" name="time"/>
    <variable units="mV" public_interface="in" private_interface="out" name="V"/>
    
    <variable units="dimensionless" private_interface="in" name="Xs"/>
    
    
    <math xmlns="http://www.w3.org/1998/Math/MathML" cmeta:id="i_Ks_calculation_eq">
      <apply id="i_Ks_calculation">
        <eq/>
        <ci> i_Ks </ci>
        <apply>
          <times/>
          <apply>
            <times/>
            <ci> g_Ks </ci>
            <apply>
              <power/>
              <ci> Xs </ci>
              <cn cellml:units="dimensionless"> 2.0 </cn>
            </apply>
          </apply>
          <apply>
            <minus/>
            <ci> V </ci>
            <ci> E_Ks </ci>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="slow_time_dependent_potassium_current_Xs_gate" cmeta:id="slow_time_dependent_potassium_current_Xs_gate">
    

    <variable units="dimensionless" public_interface="out" name="Xs" initial_value="0.5"/>
    <variable units="mV" public_interface="in" name="V"/>
    <variable units="ms" public_interface="in" name="time"/>
    
    
    <math xmlns="http://www.w3.org/1998/Math/MathML" cmeta:id="Xs_diff_eq">
      <apply id="Xs_diff">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> Xs </ci>
        </apply>
        <apply>
          <minus/>
          <apply>
            <times/>
            <apply>
              <divide/>
              <cn cellml:units="per_ms"> 0.014 </cn>
              <apply>
                <plus/>
                <cn cellml:units="dimensionless"> 1.0 </cn>
                <apply>
                  <exp/>
                  <apply>
                    <divide/>
                    <apply>
                      <minus/>
                      <apply>
                        <minus/>
                        <ci> V </ci>
                        <cn cellml:units="mV"> 40.0 </cn>
                      </apply>
                    </apply>
                    <cn cellml:units="mV"> 9.0 </cn>
                  </apply>
                </apply>
              </apply>
            </apply>
            <apply>
              <minus/>
              <cn cellml:units="dimensionless"> 1.0 </cn>
              <ci> Xs </ci>
            </apply>
          </apply>
          <apply>
            <times/>
            <ci> Xs </ci>
            <cn cellml:units="per_ms"> 0.001 </cn>
            <apply>
              <exp/>
              <apply>
                <divide/>
                <apply>
                  <minus/>
                  <ci> V </ci>
                </apply>
                <cn cellml:units="mV"> 45.0 </cn>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="potassium_background_current" cmeta:id="potassium_background_current">
    

    <variable units="uA_per_mmsq" public_interface="out" name="i_b_K"/>
    
    <variable units="mS_per_mmsq" name="g_bk" initial_value="0.0"/>
    
    <variable units="mV" public_interface="in" name="E_K"/>
    <variable units="mV" public_interface="in" name="V"/>
    <variable units="ms" public_interface="in" name="time"/>
    
    
    <math xmlns="http://www.w3.org/1998/Math/MathML" cmeta:id="i_b_K_calculation_eq">
      <apply id="i_b_K_calculation">
        <eq/>
        <ci> i_b_K </ci>
        <apply>
          <times/>
          <ci> g_bk </ci>
          <apply>
            <minus/>
            <ci> V </ci>
            <ci> E_K </ci>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="ATP_dependent_potassium_current" cmeta:id="ATP_dependent_potassium_current">
    

    <variable units="uA_per_mmsq" public_interface="out" name="i_K_ATP"/>
    
    <variable units="mS_per_mmsq" name="g_K_ATP" initial_value="0.0"/>
    <variable units="mM" name="K_ATP" initial_value="0.1"/>
    <variable units="mM" name="ATP" initial_value="5.0"/>
    
    <variable units="mV" public_interface="in" name="V"/>
    <variable units="ms" public_interface="in" name="time"/>
    
    
    <math xmlns="http://www.w3.org/1998/Math/MathML" cmeta:id="i_K_ATP_calculation_eq">
      <apply id="i_K_ATP_calculation">
        <eq/>
        <ci> i_K_ATP </ci>
        <apply>
          <divide/>
          <apply>
            <times/>
            <ci> g_K_ATP </ci>
            <apply>
              <plus/>
              <ci> V </ci>
              <cn cellml:units="mV"> 80.0 </cn>
            </apply>
          </apply>
          <apply>
            <plus/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
              <power/>
              <apply>
                <divide/>
                <ci> ATP </ci>
                <ci> K_ATP </ci>
              </apply>
              <cn cellml:units="dimensionless"> 2.0 </cn>
            </apply>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="sodium_activated_potassium_current" cmeta:id="sodium_activated_potassium_current">
    

    <variable units="uA_per_mmsq" public_interface="out" name="i_KNa"/>
    
    <variable units="mS_per_mmsq" name="g_K_Na" initial_value="0.0"/>
    <variable units="mM" name="K_kna" initial_value="20.0"/>
    
    <variable units="mM" public_interface="in" name="Na_i"/>
    <variable units="mV" public_interface="in" name="E_K"/>
    <variable units="mV" public_interface="in" name="V"/>
    <variable units="ms" public_interface="in" name="time"/>
    
    
    <math xmlns="http://www.w3.org/1998/Math/MathML" cmeta:id="i_KNa_calculation_eq">
      <apply id="i_KNa_calculation">
        <eq/>
        <ci> i_KNa </ci>
        <apply>
          <times/>
          <ci> g_K_Na </ci>
          <apply>
            <divide/>
            <ci> Na_i </ci>
            <apply>
              <plus/>
              <ci> Na_i </ci>
              <ci> K_kna </ci>
            </apply>
          </apply>
          <apply>
            <minus/>
            <ci> V </ci>
            <ci> E_K </ci>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="fast_sodium_current" cmeta:id="fast_sodium_current">
    

    <variable units="uA_per_mmsq" public_interface="out" name="i_Na"/>
    
    <variable units="mS_per_mmsq" name="g_Na" initial_value="385.5e-3"/>
    
    <variable units="mV" public_interface="in" name="E_mh"/>
    <variable units="ms" public_interface="in" private_interface="out" name="time"/>
    <variable units="mV" public_interface="in" private_interface="out" name="V"/>
    
    <variable units="dimensionless" private_interface="in" name="m"/>
    <variable units="dimensionless" private_interface="in" name="h"/>
    
    
    <math xmlns="http://www.w3.org/1998/Math/MathML" cmeta:id="i_Na_calculation_eq">
      <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>
          <apply>
            <minus/>
            <ci> V </ci>
            <ci> E_mh </ci>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="fast_sodium_current_m_gate" cmeta:id="fast_sodium_current_m_gate">
    

    <variable units="dimensionless" public_interface="out" name="m" initial_value="0.0017"/>
    <variable units="mV" public_interface="in" name="V"/>
    <variable units="ms" public_interface="in" name="time"/>
   
    <variable units="per_ms" name="alpha"/>
    <variable units="per_ms" name="beta"/>
    
    
    <math xmlns="http://www.w3.org/1998/Math/MathML" cmeta:id="m_diff_eq">
      <apply id="m_diff">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> m </ci>
        </apply>
        <apply>
          <minus/>
          <apply>
            <times/>
            <ci>alpha</ci>
            <apply>
              <minus/>
              <cn cellml:units="dimensionless"> 1.0 </cn>
              <ci> m </ci>
            </apply>
          </apply>
          <apply>
            <times/>
            <ci>beta</ci>
            <ci>m</ci>
          </apply>
        </apply>
      </apply>
      <!--m_diff_eq-->
    </math>    
    
    <math xmlns="http://www.w3.org/1998/Math/MathML" cmeta:id="alpha_calculation_eq">
      <apply id="alpha_calculation">
        <eq/>
        <ci>alpha</ci>
        <piecewise>
          <piece>
            <cn cellml:units="per_ms">2.0</cn>
            <apply>
              <lt/>
              <apply>
                <abs/>
                <apply>
                  <plus/>
                  <ci>V</ci>
                  <cn cellml:units="mV">41.0</cn>
                </apply>
              </apply>
              <cn cellml:units="mV">0.1e-4</cn>
            </apply>
            <!--lt-->
          </piece>
          <otherwise>
            <apply>
              <divide/>
              <apply>
                <times/>
                <cn cellml:units="per_mV_per_ms"> 0.2 </cn>
                <apply>
                  <plus/>
                  <ci> V </ci>
                  <cn cellml:units="mV"> 41.0 </cn>
                </apply>
              </apply>
              <apply>
                <minus/>
                <cn cellml:units="dimensionless"> 1.0 </cn>
                <apply>
                  <exp/>
                  <apply>
                    <times/>
                    <cn cellml:units="per_mV"> -0.1 </cn>
                    <apply>
                      <plus/>
                      <ci> V </ci>
                      <cn cellml:units="mV"> 41.0 </cn>
                    </apply>
                  </apply>
                </apply>
              </apply>
            </apply>
          </otherwise>
        </piecewise>
      </apply>
      <!--alpha_calculation-->
    </math>    
    
    <math xmlns="http://www.w3.org/1998/Math/MathML" cmeta:id="beta_calculation_eq">
      <apply id="beta_calculation">
        <eq/>
        <ci>beta</ci>
        <apply>
          <times/>
          <cn cellml:units="per_ms"> 8.0 </cn>
          <apply>
            <exp/>
            <apply>
              <times/>
              <cn cellml:units="per_mV"> -0.055556 </cn>
              <apply>
                <plus/>
                <ci> V </ci>
                <cn cellml:units="mV"> 66.0 </cn>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      <!--beta_calculation-->
    </math>
  </component>
  
  <component name="fast_sodium_current_h_gate" cmeta:id="fast_sodium_current_h_gate">
    

    <variable units="dimensionless" public_interface="out" name="h" initial_value="0.9939474"/>
    <variable units="mV" public_interface="in" name="V"/>
    <variable units="ms" public_interface="in" name="time"/>
   
    
    <math xmlns="http://www.w3.org/1998/Math/MathML" cmeta:id="h_diff_eq">
      <apply id="h_diff">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> h </ci>
        </apply>
        <apply>
          <minus/>
          <apply>
            <times/>
            <apply>
              <times/>
              <cn cellml:units="per_ms"> 0.02 </cn>
              <apply>
                <exp/>
                <apply>
                  <times/>
                  <cn cellml:units="per_mV"> -0.125 </cn>
                  <apply>
                    <plus/>
                    <ci> V </ci>
                    <cn cellml:units="mV"> 75.0 </cn>
                  </apply>
                </apply>
              </apply>
            </apply>
            <apply>
              <minus/>
              <cn cellml:units="dimensionless"> 1.0 </cn>
              <ci> h </ci>
            </apply>
          </apply>
          <apply>
            <times/>
            <ci> h </ci>
            <apply>
              <divide/>
              <cn cellml:units="per_ms"> 2.0 </cn>
              <apply>
                <plus/>
                <cn cellml:units="dimensionless"> 1.0 </cn>
                <apply>
                  <times/>
                  <cn cellml:units="dimensionless"> 320.0 </cn>
                  <apply>
                    <exp/>
                    <apply>
                      <times/>
                      <cn cellml:units="per_mV"> -0.1 </cn>
                      <apply>
                        <plus/>
                        <ci> V </ci>
                        <cn cellml:units="mV"> 75.0 </cn>
                      </apply>
                    </apply>
                  </apply>
                </apply>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="persistent_sodium_current" cmeta:id="persistent_sodium_current">
    

    <variable units="uA_per_mmsq" public_interface="out" name="i_p_Na"/>
    
    <variable units="mS_per_mmsq" name="g_pna" initial_value="0.61688e-3"/>
    
    <variable units="mV" public_interface="in" name="E_Na"/>
    <variable units="ms" public_interface="in" name="time"/>
    <variable units="mV" public_interface="in" name="V"/>
    
    
    <math xmlns="http://www.w3.org/1998/Math/MathML" cmeta:id="i_p_Na_calculation_eq">
      <apply id="i_p_Na_calculation">
        <eq/>
        <ci> i_p_Na </ci>
        <apply>
          <times/>
          <ci> g_pna </ci>
          <apply>
            <divide/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
              <plus/>
              <cn cellml:units="dimensionless"> 1.0 </cn>
              <apply>
                <exp/>
                <apply>
                  <minus/>
                  <apply>
                    <divide/>
                    <apply>
                      <plus/>
                      <ci> V </ci> 
                      <cn cellml:units="mV"> 52.0 </cn>
                    </apply>
                    <cn cellml:units="mV"> 8.0 </cn>
                  </apply>
                </apply>
              </apply>
            </apply>
          </apply>
          <apply>
            <minus/>
            <ci> V </ci>
            <ci> E_Na </ci>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="sodium_background_current" cmeta:id="sodium_background_current">
    

    <variable units="uA_per_mmsq" public_interface="out" name="i_b_Na"/>
    
    <variable units="mS_per_mmsq" name="g_bna" initial_value="0.092532e-3"/>
    
    <variable units="mM" public_interface="in" name="Na_i"/>
    <variable units="mV" public_interface="in" name="E_Na"/>
    <variable units="ms" public_interface="in" name="time"/>
    <variable units="mV" public_interface="in" name="V"/>
    
    
    <math xmlns="http://www.w3.org/1998/Math/MathML" cmeta:id="i_b_Na_calculation_eq">
      <apply id="i_b_Na_calculation">
        <eq/>
        <ci> i_b_Na </ci>
        <apply>
          <times/>
          <ci> g_bna </ci>
          <apply>
            <minus/>
            <ci> V </ci>
            <ci> E_Na </ci> 
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  
  <!--  
    The "L_type_Ca_channel" component describes an inward ionic current which is
    the sum of Ca, Na and K ions through the membrane channel.  The channel has 
    one activation gate (d) and two inactivation gates (f and   
    f2 (or f2ds in the diadic space)).
  -->
  <component name="L_type_Ca_channel" cmeta:id="L_type_Ca_channel">
    

    <variable units="uA_per_mmsq" public_interface="out" name="i_Ca_L_Ca"/>
    <variable units="uA_per_mmsq" public_interface="out" name="i_Ca_L_K"/>
    <variable units="uA_per_mmsq" public_interface="out" name="i_Ca_L_Na"/>
    <variable units="uA_per_mmsq" public_interface="out" name="i_Ca_L_Ca_ds"/>
    <variable units="uA_per_mmsq" public_interface="out" name="i_Ca_L_K_ds"/>
    <variable units="uA_per_mmsq" public_interface="out" name="i_Ca_L_Na_ds"/>
    
    <variable units="uA_per_mmsq_per_mM" name="P_ca" initial_value="15.422e-3"/>
    <variable units="dimensionless" name="P_caK" initial_value="0.002"/>
    <variable units="dimensionless" name="P_caNa" initial_value="0.01"/>
    <variable units="dimensionless" name="ICaLfract" initial_value="1.0"/>
    
    <variable units="mM" public_interface="in" name="Ca_o"/>
    <variable units="mM" public_interface="in" private_interface="out" name="Ca_i"/>
    <variable units="mM" public_interface="in" private_interface="out" name="Ca_ds"/>
    <variable units="mM" public_interface="in" name="K_o"/>
    <variable units="mM" public_interface="in" name="K_i"/>
    <variable units="mM" public_interface="in" name="Na_o"/>
    <variable units="mM" public_interface="in" name="Na_i"/>
    <variable units="gas_constant" public_interface="in" name="R"/>
    <variable units="faradays_constant" public_interface="in" name="F"/>
    <variable units="K" public_interface="in" name="T"/>
    <variable units="mM" private_interface="out" name="K_cachoff" initial_value="0.001"/>
    <variable units="ms" public_interface="in" private_interface="out" name="time"/>
    <variable units="mV" public_interface="in" private_interface="out" name="V"/>
        
    <variable units="mM" public_interface="in" name="ACh"/>
    
    <variable units="mM" name="K_AChICaL" initial_value="0.2e-6"/>
    
    <variable units="dimensionless" private_interface="in" name="d"/>
    <variable units="dimensionless" private_interface="in" name="f"/>
    <variable units="dimensionless" private_interface="in" name="f2"/>
    <variable units="dimensionless" private_interface="in" name="f2ds"/>
    
    
    <math xmlns="http://www.w3.org/1998/Math/MathML" cmeta:id="i_Ca_L_Ca_calculation_eq">
      <apply id="i_Ca_L_Ca_calculation">
        <eq/>
        <ci> i_Ca_L_Ca </ci>
        <apply>
          <times/>
          <apply>
            <minus/>
            <cn cellml:units="dimensionless">1.0</cn>
            <ci>ICaLfract</ci>
          </apply>
          <cn cellml:units="dimensionless"> 4.0 </cn>
          <ci> P_ca </ci>
          <ci> d </ci>
          <ci> f </ci>
          <ci> f2 </ci>
          <apply>
            <divide/>
            <apply>
              <divide/>
              <apply>
                <minus/>
                <ci> V </ci>
                <cn cellml:units="mV"> 50.0 </cn>
              </apply>
              <apply>
                <divide/>
                <apply>
                  <times/>
                  <ci> R </ci>
                  <ci> T </ci>
                </apply>
                <ci> F </ci>
              </apply>
            </apply>
            <apply>
              <minus/>
              <cn cellml:units="dimensionless"> 1.0 </cn>
              <apply>
                <exp/>
                <apply>
                  <divide/>
                  <apply>
                    <minus/>
                    <apply>
                      <minus/>
                      <ci> V </ci>
                      <cn cellml:units="mV"> 50.0 </cn>
                    </apply>
                  </apply>
                  <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>
              </apply>
            </apply>
          </apply>
          <apply>
            <minus/>
            <apply>
              <times/>
              <ci> Ca_i </ci>
              <apply>
                <exp/>
                <apply>
                  <divide/>
                  <cn cellml:units="mV"> 50.0 </cn>
                  <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>
              </apply>
            </apply>
            <apply>
              <times/>
              <ci> Ca_o </ci>
              <apply>
                <exp/>
                <apply>
                  <divide/>
                  <apply>
                    <minus/>
                    <apply>
                      <minus/>
                      <ci> V </ci>
                      <cn cellml:units="mV"> 50.0 </cn>
                    </apply>
                  </apply>
                  <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>
              </apply>
            </apply>
          </apply>
          <apply>
            <minus/>
            <cn cellml:units="dimensionless">1.0</cn>
            <apply>
              <divide/>  
              <ci>ACh</ci>
              <apply>
                <plus/>
                <ci>ACh</ci>
                <ci>K_AChICaL</ci>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
    </math>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML" cmeta:id="i_Ca_L_K_calculation_eq">
      <apply id="i_Ca_L_K_calculation">
        <eq/>
        <ci> i_Ca_L_K </ci>
        <apply>
          <times/>
          <apply>
            <minus/>
            <cn cellml:units="dimensionless">1.0</cn>
            <ci>ICaLfract</ci>
          </apply>
          <ci> P_caK </ci>
          <ci> P_ca </ci>
          <ci> d </ci>
          <ci> f </ci>
          <ci> f2 </ci>
          <apply>
            <minus/>
            <apply>
              <times/>
              <ci> K_i </ci>
              <apply>
                <exp/>
                <apply>
                  <divide/>
                  <cn cellml:units="mV"> 50.0 </cn>
                  <apply>
                    <divide/>
                    <apply>
                      <times/>
                      <ci> R </ci>
                      <ci> T </ci>
                    </apply>
                    <ci> F </ci>
                  </apply>
                </apply>
              </apply>
            </apply>
            <apply>
              <times/>
              <ci> K_o </ci>
              <apply>
                <exp/>
                <apply>
                  <divide/>
                  <apply>
                    <minus/>
                    <apply>
                      <minus/>
                      <ci> V </ci>
                      <cn cellml:units="mV"> 50.0 </cn>
                    </apply>
                  </apply>
                  <apply>
                    <divide/>
                    <apply>
                      <times/>
                      <ci> R </ci>
                      <ci> T </ci>
                    </apply>
                    <ci> F </ci>
                  </apply>
                </apply>
              </apply>
            </apply>
          </apply>
          <apply>
            <minus/>
            <cn cellml:units="dimensionless">1.0</cn>
            <apply>
              <divide/>  
              <ci>ACh</ci>
              <apply>
                <plus/>
                <ci>ACh</ci>
                <ci>K_AChICaL</ci>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
    </math>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML" cmeta:id="i_Ca_L_Na_calculation_eq">
      <apply id="i_Ca_L_Na_calculation">
        <eq/>
        <ci> i_Ca_L_Na </ci>
        <apply>
          <times/>
          <apply>
            <minus/>
            <cn cellml:units="dimensionless">1.0</cn>
            <ci>ICaLfract</ci>
          </apply>
          <ci> P_caNa </ci>
          <ci> P_ca </ci>
          <ci> d </ci>
          <ci> f </ci>
          <ci> f2 </ci>
          <apply>
            <minus/>
            <apply>
              <times/>
              <ci> Na_i </ci>
              <apply>
                <exp/>
                <apply>
                  <divide/>
                  <cn cellml:units="mV"> 50.0 </cn>
                  <apply>
                    <divide/>
                    <apply>
                      <times/>
                      <ci> R </ci>
                      <ci> T </ci>
                    </apply>
                    <ci> F </ci>
                  </apply>
                </apply>
              </apply>
            </apply>
            <apply>
              <times/>
              <ci> Na_o </ci>
              <apply>
                <exp/>
                <apply>
                  <divide/>
                  <apply>
                    <minus/>
                    <apply>
                      <minus/>
                      <ci> V </ci>
                      <cn cellml:units="mV"> 50.0 </cn>
                    </apply>
                  </apply>
                  <apply>
                    <divide/>
                    <apply>
                      <times/>
                      <ci> R </ci>
                      <ci> T </ci>
                    </apply>
                    <ci> F </ci>
                  </apply>
                </apply>
              </apply>
            </apply>
          </apply>
          <apply>
            <minus/>
            <cn cellml:units="dimensionless">1.0</cn>
            <apply>
              <divide/>  
              <ci>ACh</ci>
              <apply>
                <plus/>
                <ci>ACh</ci>
                <ci>K_AChICaL</ci>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
    </math>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML" cmeta:id="i_Ca_L_Ca_ds_calculation_eq">
      <apply id="i_Ca_L_Ca_ds_calculation">
        <eq/>
        <ci> i_Ca_L_Ca_ds </ci>
        <apply>
          <times/>
          <ci>ICaLfract</ci>
          <cn cellml:units="dimensionless"> 4.0 </cn>
          <ci> P_ca </ci>
          <ci> d </ci>
          <ci> f </ci>
          <ci> f2ds </ci>
          <apply>
            <divide/>
            <apply>
              <divide/>
              <apply>
                <minus/>
                <ci> V </ci>
                <cn cellml:units="mV"> 50.0 </cn>
              </apply>
              <apply>
                <divide/>
                <apply>
                  <times/>
                  <ci> R </ci>
                  <ci> T </ci>
                </apply>
                <ci> F </ci>
              </apply>
            </apply>
            <apply>
              <minus/>
              <cn cellml:units="dimensionless"> 1.0 </cn>
              <apply>
                <exp/>
                <apply>
                  <divide/>
                  <apply>
                    <minus/>
                    <apply>
                      <minus/>
                      <ci> V </ci>
                      <cn cellml:units="mV"> 50.0 </cn>
                    </apply>
                  </apply>
                  <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>
              </apply>
            </apply>
          </apply>
          <apply>
            <minus/>
            <apply>
              <times/>
              <ci> Ca_ds </ci>
              <apply>
                <exp/>
                <apply>
                  <divide/>
                  <cn cellml:units="mV"> 50.0 </cn>
                  <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>
              </apply>
            </apply>
            <apply>
              <times/>
              <ci> Ca_o </ci>
              <apply>
                <exp/>
                <apply>
                  <divide/>
                  <apply>
                    <minus/>
                    <apply>
                      <minus/>
                      <ci> V </ci>
                      <cn cellml:units="mV"> 50.0 </cn>
                    </apply>
                  </apply>
                  <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>
              </apply>
            </apply>
          </apply>
          <apply>
            <minus/>
            <cn cellml:units="dimensionless">1.0</cn>
            <apply>
              <divide/>  
              <ci>ACh</ci>
              <apply>
                <plus/>
                <ci>ACh</ci>
                <ci>K_AChICaL</ci>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
    </math>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML" cmeta:id="i_Ca_L_K_ds_calculation_eq">
      <apply id="i_Ca_L_K_ds_calculation">
        <eq/>
        <ci> i_Ca_L_K_ds </ci>
        <apply>
          <times/>
          <ci>ICaLfract</ci>
          <ci> P_caK </ci>
          <ci> P_ca </ci>
          <ci> d </ci>
          <ci> f </ci>
          <ci> f2ds </ci>
          <apply>
            <minus/>
            <apply>
              <times/>
              <ci> K_i </ci>
              <apply>
                <exp/>
                <apply>
                  <divide/>
                  <cn cellml:units="mV"> 50.0 </cn>
                  <apply>
                    <divide/>
                    <apply>
                      <times/>
                      <ci> R </ci>
                      <ci> T </ci>
                    </apply>
                    <ci> F </ci>
                  </apply>
                </apply>
              </apply>
            </apply>
            <apply>
              <times/>
              <ci> K_o </ci>
              <apply>
                <exp/>
                <apply>
                  <divide/>
                  <apply>
                    <minus/>
                    <apply>
                      <minus/>
                      <ci> V </ci>
                      <cn cellml:units="mV"> 50.0 </cn>
                    </apply>
                  </apply>
                  <apply>
                    <divide/>
                    <apply>
                      <times/>
                      <ci> R </ci>
                      <ci> T </ci>
                    </apply>
                    <ci> F </ci>
                  </apply>
                </apply>
              </apply>
            </apply>
          </apply>
          <apply>
            <minus/>
            <cn cellml:units="dimensionless">1.0</cn>
            <apply>
              <divide/>  
              <ci>ACh</ci>
              <apply>
                <plus/>
                <ci>ACh</ci>
                <ci>K_AChICaL</ci>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
    </math>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML" cmeta:id="i_Ca_L_Na_ds_calculation_eq">
      <apply id="i_Ca_L_Na_ds_calculation">
        <eq/>
        <ci> i_Ca_L_Na_ds </ci>
        <apply>
          <times/>
          <ci>ICaLfract</ci>
          <ci> P_caNa </ci>
          <ci> P_ca </ci>
          <ci> d </ci>
          <ci> f </ci>
          <ci> f2ds </ci>
          <apply>
            <minus/>
            <apply>
              <times/>
              <ci> Na_i </ci>
              <apply>
                <exp/>
                <apply>
                  <divide/>
                  <cn cellml:units="mV"> 50.0 </cn>
                  <apply>
                    <divide/>
                    <apply>
                      <times/>
                      <ci> R </ci>
                      <ci> T </ci>
                    </apply>
                    <ci> F </ci>
                  </apply>
                </apply>
              </apply>
            </apply>
            <apply>
              <times/>
              <ci> Na_o </ci>
              <apply>
                <exp/>
                <apply>
                  <divide/>
                  <apply>
                    <minus/>
                    <apply>
                      <minus/>
                      <ci> V </ci>
                      <cn cellml:units="mV"> 50.0 </cn>
                    </apply>
                  </apply>
                  <apply>
                    <divide/>
                    <apply>
                      <times/>
                      <ci> R </ci>
                      <ci> T </ci>
                    </apply>
                    <ci> F </ci>
                  </apply>
                </apply>
              </apply>
            </apply>
          </apply>
          <apply>
            <minus/>
            <cn cellml:units="dimensionless">1.0</cn>
            <apply>
              <divide/>  
              <ci>ACh</ci>
              <apply>
                <plus/>
                <ci>ACh</ci>
                <ci>K_AChICaL</ci>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="L_type_Ca_channel_d_gate" cmeta:id="L_type_Ca_channel_d_gate">
    

    <variable units="dimensionless" public_interface="out" name="d" initial_value="0.0"/>
    <variable units="mV" public_interface="in" name="V"/>
    <variable units="ms" public_interface="in" name="time"/>
   
    <variable units="per_ms" name="alpha"/>
    <variable units="per_ms" name="beta"/>

    
    <math xmlns="http://www.w3.org/1998/Math/MathML" cmeta:id="d_diff_eq">
      <apply id="d_diff">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> d </ci>
        </apply>
        <apply>
          <minus/>
          <apply>
            <times/>
            <ci>alpha</ci>
            <apply>
              <minus/>
              <cn cellml:units="dimensionless">1.0</cn>
              <ci>d</ci>
            </apply>
          </apply>
          <apply>
            <times/>
            <ci>beta</ci>
            <ci>d</ci>
          </apply>
        </apply>
      </apply>
      <!--d_diff_eq-->
    </math>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML" cmeta:id="alpha_calculation_eq2">
      <apply id="alpha_calculation">
        <eq/>
        <ci>alpha</ci>
        <piecewise>
          <piece>
            <cn cellml:units="per_ms">0.36</cn>
            <apply>
              <lt/>
              <apply>
                <abs/>
                <apply>
                  <plus/>
                  <ci>V</ci>
                  <cn cellml:units="mV">19.0</cn>
                </apply>
              </apply>
              <cn cellml:units="mV">1.0e-3</cn>
            </apply>
          </piece>
          <otherwise>
            <apply>
              <times/>
              <cn cellml:units="per_ms">1.0</cn>
              <apply>
                <divide/>
                <apply>
                  <times/>
                  <cn cellml:units="per_mV"> 0.09 </cn>
                  <apply>
                    <plus/>
                    <ci> V </ci>
                    <cn cellml:units="mV"> 19.0 </cn>
                  </apply>
                </apply>
                <apply>
                  <minus/>
                  <cn cellml:units="dimensionless"> 1.0 </cn>
                  <apply>
                    <exp/>
                    <apply>
                      <divide/>
                      <apply>
                        <minus/>
                        <apply>
                          <plus/>
                          <ci> V </ci>
                          <cn cellml:units="mV"> 19.0 </cn>
                        </apply>
                      </apply>
                      <cn cellml:units="mV"> 4.0 </cn>
                    </apply>
                  </apply>
                </apply>
              </apply>
            </apply>
          </otherwise>
        </piecewise>
      </apply>
      <!--alpha_calculation-->
    </math>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML" cmeta:id="beta_calculation_eq2">
      <apply id="beta_calculation">
        <eq/>
        <ci>beta</ci>
        <piecewise>
          <piece>
            <cn cellml:units="per_ms">0.36</cn>
            <apply>
              <lt/>
              <apply>
                <abs/>
                <apply>
                  <plus/>
                  <ci>V</ci>
                  <cn cellml:units="mV">19.0</cn>
                </apply>
              </apply>
              <cn cellml:units="mV">1.0e-3</cn>
            </apply>
          </piece>
          <otherwise>
            <apply>
              <times/>
              <cn cellml:units="per_ms">1.0</cn>
              <apply>
                <divide/>
                <apply>
                  <times/>
                  <cn cellml:units="per_mV"> -0.036 </cn>
                  <apply>
                    <plus/>
                    <ci> V </ci>
                    <cn cellml:units="mV"> 19.0 </cn>
                  </apply>
                </apply>
                <apply>
                  <minus/>
                  <cn cellml:units="dimensionless"> 1.0 </cn>
                  <apply>
                    <exp/>
                    <apply>
                      <divide/>
                      <apply>
                        <plus/>
                        <ci> V </ci>
                        <cn cellml:units="mV"> 19.0 </cn>
                      </apply>
                      <cn cellml:units="mV"> 10.0 </cn>
                    </apply>
                  </apply>
                </apply>
              </apply>
            </apply>
          </otherwise>
        </piecewise>
      </apply>
      <!--beta_calculation-->
    </math>
  </component>
  
  <component name="L_type_Ca_channel_f_gate" cmeta:id="L_type_Ca_channel_f_gate">
    

    <variable units="dimensionless" public_interface="out" name="f" initial_value="1.0"/>
    <variable units="mM" public_interface="in" name="K_cachoff"/>
    <variable units="mM" public_interface="in" name="Ca_i"/>
    <variable units="mV" public_interface="in" name="V"/>
    <variable units="ms" public_interface="in" name="time"/>
   
    <variable units="per_ms" name="alpha"/>
    <variable units="per_ms" name="beta"/>
    
    
    <math xmlns="http://www.w3.org/1998/Math/MathML" cmeta:id="f_diff_eq">
      <apply id="f_diff">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> f </ci>
        </apply>
        <apply>
          <minus/>
          <apply>
            <times/>
            <ci>alpha</ci>
            <apply>
              <minus/>
              <cn cellml:units="dimensionless">1.0</cn>
              <ci>f</ci>
            </apply>
          </apply>
          <apply>
            <times/>
            <ci>beta</ci>
            <ci>f</ci>
          </apply>
        </apply>
      </apply>
      <!--f_diff_eq-->
    </math>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML" cmeta:id="alpha_calculation_eq3">
      <apply id="alpha_calculation">
        <eq/>
        <ci>alpha</ci>
        <piecewise>
          <piece>
            <cn cellml:units="per_ms">0.0075</cn>
            <apply>
              <lt/>
              <apply>
                <abs/>
                <apply>
                  <plus/>
                  <ci>V</ci>
                  <cn cellml:units="mV">34.0</cn>
                </apply>
              </apply>
              <cn cellml:units="mV">1.0e-3</cn>
            </apply>
          </piece>
          <otherwise>
            <apply>
              <times/>
              <cn cellml:units="per_ms">1.0</cn>
              <apply>
                <times/>
                <cn cellml:units="per_mV">-0.001875</cn>
                <apply>
                  <divide/>
                  <apply>
                    <plus/>
                    <ci>V</ci>
                    <cn cellml:units="mV">34.0</cn>
                  </apply>
                  <!--plus-->
                  <apply>
                    <minus/>
                    <cn cellml:units="dimensionless">1.0</cn>
                    <apply>
                      <exp/>
                      <apply>
                        <divide/>
                        <apply>
                          <plus/>
                          <ci>V</ci>
                          <cn cellml:units="mV">34.0</cn>
                        </apply>
                        <!--plus-->
                        <cn cellml:units="mV">4.0</cn>
                      </apply>
                      <!--divide-->
                    </apply>
                    <!--exp-->
                  </apply>
                  <!--minus-->
                </apply>
                <!--divide-->
              </apply>
              <!--times-->
            </apply>
          </otherwise>
        </piecewise>
      </apply>
      <!--alpha_calculation-->
    </math>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML" cmeta:id="beta_calculation_eq3">
      <apply id="beta_calculation">
        <eq/>
        <ci>beta</ci>
        <apply>
          <divide/>
          <cn cellml:units="per_ms">0.0036</cn>
          <apply>
            <plus/>
            <cn cellml:units="dimensionless">1.0</cn>
            <apply>
              <exp/>
              <apply>
                <minus/>
                <apply>
                  <divide/>
                  <apply>
                    <plus/>
                    <ci>V</ci>
                    <cn cellml:units="mV">34.0</cn>
                  </apply>
                  <!--plus-->
                  <cn cellml:units="mV">4.0</cn>
                </apply>
                <!--divide-->
              </apply>
              <!--minus-->
            </apply>
            <!--exp-->
          </apply>
          <!--plus-->
        </apply>
        <!--divide-->
      </apply>
      <!--eq-->
    </math>
  </component>
  
  <component name="L_type_Ca_channel_f2_gate" cmeta:id="L_type_Ca_channel_f2_gate">
    

    <variable units="dimensionless" public_interface="out" name="f2" initial_value="1.0"/>
    
    <variable units="per_ms" name="CaInact" initial_value="0.02"/>
    
    <variable units="mM" public_interface="in" name="K_cachoff"/>
    <variable units="mM" public_interface="in" name="Ca_i"/>
    <variable units="mV" public_interface="in" name="V"/>
    <variable units="ms" public_interface="in" name="time"/>
   
    
    <math xmlns="http://www.w3.org/1998/Math/MathML" cmeta:id="f2_diff_eq">
      <apply id="f2_diff">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> f2 </ci>
        </apply>
        <apply>
          <times/>
          <ci>CaInact</ci>
          <apply>
            <minus/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
              <plus/> 
              <apply>
                <divide/>
                <ci> Ca_i </ci>
                <apply>
                  <plus/>
                  <ci> K_cachoff </ci>
                  <ci> Ca_i </ci>
                </apply>
              </apply>
              <ci> f2 </ci>
            </apply>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="L_type_Ca_channel_f2ds_gate" cmeta:id="L_type_Ca_channel_f2ds_gate">
    

    <variable units="dimensionless" public_interface="out" name="f2ds" initial_value="1.0"/>
    <variable units="mM" name="K_dsoff" initial_value="0.001"/>
    <variable units="mM" public_interface="in" name="Ca_ds"/>
    <variable units="mV" public_interface="in" name="V"/>
    <variable units="ms" public_interface="in" name="time"/>
    
    <variable units="per_ms" name="CaInactDS" initial_value="0.02"/>
   
    
    <math xmlns="http://www.w3.org/1998/Math/MathML" cmeta:id="f2ds_diff_eq">
      <apply id="f2ds_diff">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> f2ds </ci>
        </apply>
        <apply>
          <times/>
          <ci> CaInactDS </ci>
          <apply>
            <minus/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
              <plus/> 
              <apply>
                <divide/>
                <ci> Ca_ds </ci>
                <apply>
                  <plus/>
                  <ci> K_dsoff </ci>
                  <ci> Ca_ds </ci>
                </apply>
              </apply>
              <ci> f2ds </ci>
            </apply>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="calcium_background_current" cmeta:id="calcium_background_current">
    

    <variable units="uA_per_mmsq" public_interface="out" name="i_b_Ca"/>
    
    <variable units="mS_per_mmsq" name="g_bca" initial_value="0.038555e-3"/>
    
    <variable units="mV" public_interface="in" name="E_Ca"/>
    <variable units="ms" public_interface="in" name="time"/>
    <variable units="mV" public_interface="in" name="V"/>
    
    
    <math xmlns="http://www.w3.org/1998/Math/MathML" cmeta:id="i_b_Ca_calculation_eq">
      <apply id="i_b_Ca_calculation">
        <eq/>
        <ci> i_b_Ca </ci>
        <apply>
          <times/>
          <ci> g_bca </ci>
          <apply>
            <minus/>
            <ci> V </ci>
            <ci> E_Ca </ci> 
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="transient_outward_current" cmeta:id="transient_outward_current">
    

    <variable units="uA_per_mmsq" public_interface="out" name="i_to"/>
    
    <variable units="mS_per_mmsq" name="g_to" initial_value="0.77100e-3"/>
    <variable units="dimensionless" name="g_tos" initial_value="0.0"/>
    
    <variable units="mV" public_interface="in" name="E_K"/>
    <variable units="mM" public_interface="in" name="Ca_i"/>
    <variable units="mV" public_interface="in" private_interface="out" name="V"/>
    <variable units="ms" public_interface="in" private_interface="out" name="time"/>
    
    <variable units="dimensionless" private_interface="in" name="s"/>
    <variable units="dimensionless" private_interface="in" name="r"/>
    
    
    <math xmlns="http://www.w3.org/1998/Math/MathML" cmeta:id="i_to_calculation_eq">
      <apply id="i_to_calculation">
        <eq/>
        <ci> i_to </ci>
        <apply>
          <times/>
          <ci> g_to </ci>
          <apply>
            <plus/>
            <ci> g_tos </ci>
            <apply>
              <times/>
              <ci> s </ci>
              <apply>
                <minus/>
                <cn cellml:units="dimensionless"> 1.0 </cn>
                <ci> g_tos </ci>
              </apply>
            </apply>
          </apply>
          <ci> r </ci>
          <apply>
            <minus/>
            <ci> V </ci>
            <ci> E_K </ci>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="transient_outward_current_s_gate" cmeta:id="transient_outward_current_s_gate">
    

    <variable units="dimensionless" public_interface="out" name="s" initial_value="0.9950485"/>
    <variable units="mV" public_interface="in" name="V"/>
    <variable units="ms" public_interface="in" name="time"/>
   
    
    <math xmlns="http://www.w3.org/1998/Math/MathML" cmeta:id="s_diff_eq">
      <apply id="s_diff">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> s </ci>
        </apply>
        <apply>
          <minus/>
          <apply>
            <times/>
            <cn cellml:units="per_ms"> 0.000033 </cn>
            <apply>
              <exp/>
              <apply>
                <divide/>
                <apply>
                  <minus/>
                  <ci> V </ci>
                </apply>
                <cn cellml:units="mV"> 17.0 </cn>
              </apply>
            </apply>
            <apply>
              <minus/>
              <cn cellml:units="dimensionless"> 1.0 </cn>
              <ci> s </ci>
            </apply>
          </apply>
          <apply>
            <times/>
            <apply>
              <divide/>
              <cn cellml:units="per_ms"> 0.033 </cn>
              <apply>
                <plus/>
                <cn cellml:units="dimensionless"> 1.0 </cn>
                <apply>
                  <exp/>
                  <apply>
                    <divide/>
                    <apply>
                      <minus/>
                      <apply>
                        <plus/>
                        <ci> V </ci>
                        <cn cellml:units="mV"> 10.0 </cn>
                      </apply>
                    </apply>
                    <cn cellml:units="mV"> 8.0 </cn>
                  </apply>
                </apply>
              </apply>
            </apply>
            <ci> s </ci>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="transient_outward_current_r_gate" cmeta:id="transient_outward_current_r_gate">
    

    <variable units="dimensionless" public_interface="out" name="r" initial_value="0.0"/>
    <variable units="dimensionless" public_interface="out" name="r_ss"/>
    <variable units="mV" public_interface="in" name="V"/>
    <variable units="ms" public_interface="in" name="time"/>
   
    
    <math xmlns="http://www.w3.org/1998/Math/MathML" cmeta:id="r_diff_eq">
      <apply id="r_diff">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> r </ci>
        </apply>
        <apply>
          <times/>
          <cn cellml:units="per_ms"> 0.333 </cn>
          <apply>
            <minus/>
            <ci> r_ss </ci>
            <ci> r </ci>
          </apply>
        </apply>
      </apply>
    </math>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML" cmeta:id="r_ss_calculation_eq">
      <apply id="r_ss_calculation">
        <eq/>
        <ci> r_ss </ci>
        <apply>
          <divide/>
          <cn cellml:units="dimensionless"> 1.0 </cn>
          <apply>
            <plus/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
              <exp/>
              <apply>
                <times/>
                <cn cellml:units="per_mV"> -0.2 </cn>
                <apply>
                  <plus/>
                  <ci> V </ci>
                  <cn cellml:units="mV"> 4.0 </cn>
                </apply>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
    </math>    
  </component>
  
  <component name="ACh_dependent_potassium_current" cmeta:id="ACh_dependent_potassium_current">
    

    <variable units="uA_per_mmsq" public_interface="out" name="i_K_ACh"/>
    
    <variable units="mS_per_mmsq" name="g_KACh" initial_value="0.0"/>
    <variable units="mM" public_interface="out" name="ACh" initial_value="0.0"/>
    <variable units="mM" name="K_ACh" initial_value="0.2795e-6"/>
    <variable units="mM" name="K_KACh" initial_value="10.0"/>
    
    <variable units="mV" public_interface="in" name="E_K"/>
    <variable units="mM" public_interface="in" name="K_o"/>
    <variable units="mM" public_interface="in" name="K_mk1"/>
    <variable units="gas_constant" public_interface="in" name="R"/>
    <variable units="faradays_constant" public_interface="in" name="F"/>
    <variable units="K" public_interface="in" private_interface="out" name="T"/>
    <variable units="mV" public_interface="in" private_interface="out" name="V"/>
    <variable units="ms" public_interface="in" private_interface="out" name="time"/>
    
    <variable units="dimensionless" private_interface="in" name="x_ACh1"/>
    <variable units="dimensionless" private_interface="in" name="x_ACh2"/>
    
    
    <math xmlns="http://www.w3.org/1998/Math/MathML" cmeta:id="i_K_ACh_calculation_eq">
      <apply id="i_K_ACh_calculation">
        <eq/>
        <ci> i_K_ACh </ci>
        <apply>
          <times/>
          <ci> g_KACh </ci>
          <apply>
            <divide/>
            <ci> K_o </ci>
            <apply>
              <plus/>
              <ci> K_o </ci>
              <ci> K_KACh </ci>
            </apply>
          </apply>
          <ci> x_ACh1 </ci>
          <ci> x_ACh2 </ci>
          <apply>
            <divide/>
            <apply>
              <power/>
              <ci>ACh</ci>
              <cn cellml:units="dimensionless">1.4969</cn>
            </apply>
            <apply>
              <plus/>
              <apply>
                <power/>
                <ci>ACh</ci>
                <cn cellml:units="dimensionless">1.4969</cn>
              </apply>
              <apply>
                <power/>
                <ci>K_ACh</ci>
                <cn cellml:units="dimensionless">1.4969</cn>
              </apply>
            </apply>
          </apply>
          <apply>
            <divide/>
            <apply>
              <minus/>
              <ci> V </ci>
              <ci> E_K </ci>
            </apply>
            <apply>
              <plus/>
              <cn cellml:units="dimensionless"> 1.0 </cn>
              <apply>
                <exp/>
                <apply>
                  <divide/>
                  <apply>
                    <times/>
                    <cn cellml:units="dimensionless"> 0.4 </cn>
                    <apply>
                      <minus/>
                      <ci> V </ci>
                      <apply>
                        <plus/>
                        <ci> E_K </ci>
                        <cn cellml:units="mV"> 140.0 </cn>
                      </apply>
                    </apply>
                  </apply>
                  <apply>
                    <divide/>
                    <apply>
                      <times/>
                      <ci> R </ci>
                      <ci> T </ci>
                    </apply>
                    <ci> F </ci>
                  </apply>
                </apply>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="ACh_dependent_potassium_current_xACh1_gate" cmeta:id="ACh_dependent_potassium_current_xACh1_gate">
    

    <variable units="dimensionless" public_interface="out" name="x_ACh1" initial_value="0.5"/>
    <variable units="per_ms" name="alpha_ACh" initial_value="0.003684211"/>
    <variable units="mV" public_interface="in" name="V"/>
    <variable units="ms" public_interface="in" name="time"/>
   
    
    <math xmlns="http://www.w3.org/1998/Math/MathML" cmeta:id="x_ACh1_diff_eq">
      <apply id="x_ACh1_diff">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> x_ACh1 </ci>
        </apply>
        <apply>
          <minus/>
          <apply>
            <times/>
            <ci> alpha_ACh </ci>
            <apply>
              <minus/>
              <cn cellml:units="dimensionless"> 1.0 </cn>
              <ci> x_ACh1 </ci>
            </apply>
          </apply>
          <apply>
            <times/>
            <apply>
              <divide/>
              <cn cellml:units="per_ms">0.00582</cn>
              <apply>
                <plus/>
                <cn cellml:units="dimensionless"> 1.0 </cn>
                <apply>
                  <exp/>
                  <apply>
                    <minus/>
                    <apply>
                      <divide/>
                      <apply>
                        <plus/>
                        <ci>V</ci>
                        <cn cellml:units="mV">50.0</cn>
                      </apply>
                      <cn cellml:units="mV">15.0</cn>
                    </apply>
                  </apply>
                </apply>
              </apply>
            </apply>
            <ci> x_ACh1 </ci>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
          
  <component name="ACh_dependent_potassium_current_xACh2_gate" cmeta:id="ACh_dependent_potassium_current_xACh2_gate">
    

    <variable units="dimensionless" public_interface="out" name="x_ACh2" initial_value="0.5"/>
    <variable units="per_ms" name="alpha_ACh" initial_value="0.07309924"/>
    <variable units="mV" public_interface="in" name="V"/>
    <variable units="ms" public_interface="in" name="time"/>
   
    
    <math xmlns="http://www.w3.org/1998/Math/MathML" cmeta:id="x_ACh2_diff_eq">
      <apply id="x_ACh2_diff">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> x_ACh2 </ci>
        </apply>
        <apply>
          <minus/>
          <apply>
            <times/>
            <ci> alpha_ACh </ci>
            <apply>
              <minus/>
              <cn cellml:units="dimensionless"> 1.0 </cn>
              <ci> x_ACh2 </ci>
            </apply>
          </apply>
          <apply>
            <times/>
            <apply>
              <divide/>
              <cn cellml:units="per_ms">0.12</cn>
              <apply>
                <plus/>
                <cn cellml:units="dimensionless"> 1.0 </cn>
                <apply>
                  <exp/>
                  <apply>
                    <minus/>
                    <apply>
                      <divide/>
                      <apply>
                        <plus/>
                        <ci>V</ci>
                        <cn cellml:units="mV">50.0</cn>
                      </apply>
                      <cn cellml:units="mV">15.0</cn>
                    </apply>
                  </apply>
                </apply>
              </apply>
            </apply>
            <ci> x_ACh2 </ci>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
          
  <component name="sodium_potassium_pump" cmeta:id="sodium_potassium_pump">
    

    <variable units="uA_per_mmsq" public_interface="out" name="i_NaK"/>
    
    <variable units="uA_per_mmsq" name="i_NaK_max" initial_value="107.95e-3"/>
    <variable units="mM" name="K_mk" initial_value="1.0"/>
    <variable units="mM" name="K_mNa" initial_value="40.0"/>
    
    <variable units="mM" public_interface="in" name="K_o"/>
    <variable units="mM" public_interface="in" name="Na_i"/>
    <variable units="ms" public_interface="in" name="time"/>
    
    
    <math xmlns="http://www.w3.org/1998/Math/MathML" cmeta:id="i_NaK_calculation_eq">
      <apply id="i_NaK_calculation">
        <eq/>
        <ci> i_NaK </ci>
        <apply>
          <times/>
          <ci> i_NaK_max </ci>
          <apply>
            <divide/>
            <ci> K_o </ci>
            <apply>
              <plus/>
              <ci> K_o </ci>
              <ci> K_mk </ci>
            </apply>
          </apply>
          <apply>
            <divide/>
            <ci> Na_i </ci>
            <apply>
              <plus/>
              <ci> Na_i </ci>
              <ci> K_mNa </ci>
            </apply>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="sodium_calcium_pump" cmeta:id="sodium_calcium_pump">
    
  
    <units name="uA_per_mmsq_per_mM4">
      <unit units="uA"/>
      <unit units="mm" exponent="-2"/>
      <unit units="mM" exponent="-4"/>
    </units>
    <units name="per_mM4">
      <unit units="mM" exponent="-4"/>
    </units>
    
    <variable units="uA_per_mmsq" public_interface="out" name="i_NaCa"/>
    <variable units="uA_per_mmsq" public_interface="out" name="i_NaCa_ds"/>
    
    <variable units="uA_per_mmsq_per_mM4" name="K_NaCa" initial_value="0.07711e-3"/>
    <variable units="per_mM4" name="d_NaCa" initial_value="0.0"/>
    <variable units="dimensionless" name="gamma" initial_value="0.5"/>
    
    <variable units="dimensionless" name="INaCaFract" initial_value="0.001"/>
    
    <variable units="mM" public_interface="in" name="Na_i"/>
    <variable units="mM" public_interface="in" name="Na_o"/>
    <variable units="mM" public_interface="in" name="Ca_i"/>  
    <variable units="mM" public_interface="in" name="Ca_o"/>
    <variable units="mM" public_interface="in" name="Ca_ds"/>
    <variable units="gas_constant" public_interface="in" name="R"/>
    <variable units="faradays_constant" public_interface="in" name="F"/>
    <variable units="K" public_interface="in" name="T"/>
    <variable units="mV" public_interface="in" name="V"/>
    <variable units="ms" public_interface="in" name="time"/>
    
    
    <math xmlns="http://www.w3.org/1998/Math/MathML" cmeta:id="i_NaCa_calculation_eq">
      <apply id="i_NaCa_calculation">
        <eq/>
        <ci> i_NaCa </ci>
        <apply>
          <times/>
          <apply>
            <minus/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <ci>INaCaFract</ci>
          </apply>
          <ci> K_NaCa </ci>
          <apply>
            <divide/>
            <apply>
              <minus/>
              <apply>
                <times/>
                <apply>
                  <exp/>
                  <apply>
                    <times/>
                    <ci> gamma </ci>
                    <apply>
                      <divide/>
                      <ci> V </ci>
                      <apply>
                        <divide/>
                        <apply>
                          <times/>
                          <ci> R </ci>
                          <ci> T </ci>
                        </apply>
                        <ci> F </ci>
                      </apply>
                    </apply>
                  </apply>
                </apply>
                <apply>
                  <power/>
                  <ci> Na_i </ci>
                  <cn cellml:units="dimensionless"> 3.0 </cn>
                </apply>
                <ci> Ca_o </ci>
              </apply>
              <apply>
                <times/>
                <apply>
                  <exp/>
                  <apply>
                    <times/>
                    <apply>
                      <minus/>
                      <ci> gamma </ci>
                      <cn cellml:units="dimensionless"> 1.0 </cn>
                    </apply>
                    <apply>
                      <divide/>
                      <ci> V </ci>
                      <apply>
                        <divide/>
                        <apply>
                          <times/>
                          <ci> R </ci>
                          <ci> T </ci>
                        </apply>
                        <ci> F </ci>
                      </apply>
                    </apply>
                  </apply>
                </apply>
                <apply>
                  <power/>
                  <ci> Na_o </ci>
                  <cn cellml:units="dimensionless"> 3.0 </cn>
                </apply>
                <ci> Ca_i </ci>
              </apply>
            </apply>
            <apply>
              <plus/>
              <cn cellml:units="dimensionless"> 1.0 </cn>
              <apply>
                <times/>
                <ci> d_NaCa </ci>
                <apply>
                  <plus/>
                  <apply>
                    <times/>
                    <ci> Ca_i </ci>
                    <apply>
                      <power/>
                      <ci> Na_o </ci>
                      <cn cellml:units="dimensionless"> 3.0 </cn>
                    </apply>
                  </apply>
                  <apply>
                    <times/>
                    <ci> Ca_o </ci>
                    <apply>
                      <power/>
                      <ci> Na_i </ci>
                      <cn cellml:units="dimensionless"> 3.0 </cn>
                    </apply>
                  </apply>
                </apply>
              </apply>
            </apply>
          </apply>
          <apply>
            <divide/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
              <plus/>
              <cn cellml:units="dimensionless"> 1.0 </cn>
              <apply>
                <divide/>
                <ci>Ca_i</ci>
                <cn cellml:units="mM"> 0.0069 </cn>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
    </math>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML" cmeta:id="i_NaCa_ds_calculation_eq">
      <apply id="i_NaCa_ds_calculation">
        <eq/>
        <ci> i_NaCa_ds </ci>
        <apply>
          <times/>
          <ci>INaCaFract</ci>
          <ci> K_NaCa </ci>
          <apply>
            <divide/>
            <apply>
              <minus/>
              <apply>
                <times/>
                <apply>
                  <exp/>
                  <apply>
                    <times/>
                    <ci> gamma </ci>
                    <apply>
                      <divide/>
                      <ci> V </ci>
                      <apply>
                        <divide/>
                        <apply>
                          <times/>
                          <ci> R </ci>
                          <ci> T </ci>
                        </apply>
                        <ci> F </ci>
                      </apply>
                    </apply>
                  </apply>
                </apply>
                <apply>
                  <power/>
                  <ci> Na_i </ci>
                  <cn cellml:units="dimensionless"> 3.0 </cn>
                </apply>
                <ci> Ca_o </ci>
              </apply>
              <apply>
                <times/>
                <apply>
                  <exp/>
                  <apply>
                    <times/>
                    <apply>
                      <minus/>
                      <ci> gamma </ci>
                      <cn cellml:units="dimensionless"> 1.0 </cn>
                    </apply>
                    <apply>
                      <divide/>
                      <ci> V </ci>
                      <apply>
                        <divide/>
                        <apply>
                          <times/>
                          <ci> R </ci>
                          <ci> T </ci>
                        </apply>
                        <ci> F </ci>
                      </apply>
                    </apply>
                  </apply>
                </apply>
                <apply>
                  <power/>
                  <ci> Na_o </ci>
                  <cn cellml:units="dimensionless"> 3.0 </cn>
                </apply>
                <ci> Ca_ds </ci>
              </apply>
            </apply>
            <apply>
              <plus/>
              <cn cellml:units="dimensionless"> 1.0 </cn>
              <apply>
                <times/>
                <ci> d_NaCa </ci>
                <apply>
                  <plus/>
                  <apply>
                    <times/>
                    <ci> Ca_ds </ci>
                    <apply>
                      <power/>
                      <ci> Na_o </ci>
                      <cn cellml:units="dimensionless"> 3.0 </cn>
                    </apply>
                  </apply>
                  <apply>
                    <times/>
                    <ci> Ca_o </ci>
                    <apply>
                      <power/>
                      <ci> Na_i </ci>
                      <cn cellml:units="dimensionless"> 3.0 </cn>
                    </apply>
                  </apply>
                </apply>
              </apply>
            </apply>
          </apply>
          <apply>
            <divide/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
              <plus/>
              <cn cellml:units="dimensionless"> 1.0 </cn>
              <apply>
                <divide/>
                <ci>Ca_ds</ci>
                <cn cellml:units="mM"> 0.0069 </cn>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="calcium_uptake_from_myoplasm_to_NSR" cmeta:id="calcium_uptake_from_myoplasm_to_NSR">
    

    <variable units="mM_per_ms" public_interface="out" name="i_up"/>
    
    <variable units="mM" name="K_cyca" initial_value="0.0003"/>
    <variable units="dimensionless" name="K_xcs" initial_value="0.4"/>
    <variable units="mM" name="K_srca" initial_value="0.5"/>
    
    <variable units="mM" public_interface="in" name="Ca_i"/>
    <variable units="mM" public_interface="in" name="Ca_up"/>
    <variable units="ms" public_interface="in" name="time"/>
    
    
    <math xmlns="http://www.w3.org/1998/Math/MathML" cmeta:id="i_up_calculation_eq">      
      <apply id="i_up_calculation">
        <eq/>
        <ci> i_up </ci>
        <apply>
          <divide/>
          <apply>
            <minus/>
            <apply>
              <times/>
              <cn cellml:units="mM_per_ms"> 0.0004 </cn>
              <ci> Ca_i </ci>
            </apply>
            <apply>
              <times/>
              <cn cellml:units="mM_per_ms"> 0.00003 </cn>
              <ci> Ca_up </ci>
              <apply>
                <divide/>
                <apply>
                  <times/>
                  <ci> K_cyca </ci>
                  <ci> K_xcs </ci>
                </apply>
                <ci> K_srca </ci>
              </apply>
            </apply>
          </apply>
          <apply>
            <plus/>
            <ci> Ca_i </ci>
            <apply>
              <times/>
              <ci> Ca_up </ci>
              <apply>
                <divide/>
                <apply>
                  <times/>
                  <ci> K_cyca </ci>
                  <ci> K_xcs </ci>
                </apply>
                <ci> K_srca </ci>
              </apply>
            </apply>
            <apply>
              <times/>
              <ci> K_cyca </ci>
              <ci> K_xcs </ci>
            </apply>
            <ci> K_cyca </ci>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="calcium_transfer_from_NSR_to_JSR" cmeta:id="calcium_transfer_from_NSR_to_JSR">
    

    <variable units="mM_per_ms" public_interface="out" name="i_tr"/>
    
    <variable units="mM" public_interface="in" name="Ca_rel"/>
    <variable units="mM" public_interface="in" name="Ca_up"/>
    <variable units="ms" public_interface="in" name="time"/>
    
    
    <math xmlns="http://www.w3.org/1998/Math/MathML" cmeta:id="i_tr_calculation_eq">
      <apply id="i_tr_calculation">
        <eq/>
        <ci> i_tr </ci>
        <apply>
          <times/>
          <cn cellml:units="per_ms"> 0.05 </cn>
          <apply>
            <minus/>
            <ci> Ca_up </ci>
            <ci> Ca_rel </ci>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="length_dependence" cmeta:id="length_dependence">
    

    <variable units="dimensionless" public_interface="out" name="f_stretch"/>
    
    <variable units="mm" public_interface="out" name="SL" initial_value="0.0"/>
    <variable units="mm" public_interface="out" name="SLRef" initial_value="2.0e-3"/>
    <variable units="per_mm" public_interface="out" name="gamma_SR_SL" initial_value="1.5e+3"/>
    <variable units="per_mm" name="gamma_SAC_SL" initial_value="2.5e+3"/>
    
    
    <math xmlns="http://www.w3.org/1998/Math/MathML" cmeta:id="f_stretch_calculation_eq">
      <apply id="f_stretch_calculation">
        <eq/>
        <ci> f_stretch </ci>
        <apply>
          <divide/>
          <cn cellml:units="dimensionless"> 1.0 </cn>
          <apply>
            <plus/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
              <exp/>
              <apply>
                <times/>
                <cn cellml:units="dimensionless"> -2.0 </cn> 
                <ci> gamma_SAC_SL </ci>
                <apply>
                  <minus/>
                  <ci> SL </ci>
                  <ci> SLRef </ci>
                </apply>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="calcium_leak_from_SR" cmeta:id="calcium_leak_from_SR">
    

    <variable units="mM_per_ms" public_interface="out" name="i_leak"/>
    
    <variable units="mm" public_interface="in" name="SL"/>
    <variable units="per_mm" public_interface="in" name="gamma_SR_SL"/>
    <variable units="mM" public_interface="in" name="Ca_rel"/>
    
    <variable units="per_ms" name="alpha_SRleak" initial_value="0.00005"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">      
      <apply id="i_leak_calculation">
        <eq/>
        <ci> i_leak </ci>
        <apply>
          <times/>
          <apply>
            <exp/>
            <apply>
              <times/>
              <ci>gamma_SR_SL</ci>
              <ci>SL</ci>
            </apply>
          </apply>
          <ci>alpha_SRleak</ci>
          <ci>Ca_rel</ci>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="ds_decay_flux" cmeta:id="ds_decay_flux">
    

    <variable units="mM_per_ms" public_interface="out" name="i_decay"/>
    
    <variable units="mM" public_interface="in" name="Ca_ds"/>
    
    <variable units="per_ms" name="alpha_ca_ds_decay" initial_value="0.01"/>
    
    
    <math xmlns="http://www.w3.org/1998/Math/MathML" cmeta:id="i_leak_calculation_eq">
      <apply id="i_leak_calculation">
        <eq/>
        <ci> i_decay</ci>
        <apply>
          <times/>
          <ci>alpha_ca_ds_decay</ci>
          <ci>Ca_ds</ci>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="calcium_release_from_SR" cmeta:id="calcium_release_from_SR">
    
    
    <variable units="mM" public_interface="in" name="Ca_rel"/>
    <variable units="mM" public_interface="in" name="Ca_i"/>
    <variable units="mM" public_interface="in" name="Ca_ds"/>
    <variable units="ms" public_interface="in" name="time"/>
    <variable units="mV" public_interface="in" name="V"/>
    
    <variable units="mM_per_ms" public_interface="out" name="i_rel"/>

    <variable units="dimensionless" name="ActivatorFract" initial_value="0.00191"/>
    <variable units="dimensionless" name="ProductFract" initial_value="0.28546"/>
    <variable units="dimensionless" name="OpenReleaseChannelFract"/>
    <variable units="dimensionless" name="RegulatoryBindingSite"/>
    <variable units="dimensionless" name="PrecursorFract"/>
    <variable units="per_ms" name="alpha_act"/>
    <variable units="per_ms" name="alpha_inact"/>
    <variable units="dimensionless" name="SR_speed_factor"/>
    
    <variable units="mM" name="K_mca" initial_value="0.0005"/>
    <variable units="mM" name="K_mca_ds" initial_value="0.01"/>
    <variable units="per_ms" name="alpha_rel" initial_value="0.25"/>
       
    
    <math xmlns="http://www.w3.org/1998/Math/MathML" cmeta:id="OpenReleaseChannelFract_calculation_eq">
      <apply id="OpenReleaseChannelFract_calculation">
        <eq/>
        <ci>OpenReleaseChannelFract</ci>
        <apply>
          <power/>
          <apply>
            <divide/>
            <ci>ActivatorFract</ci>
            <apply>
              <plus/>
              <ci>ActivatorFract</ci>
              <cn cellml:units="dimensionless">0.25</cn>
            </apply>
          </apply>
          <cn cellml:units="dimensionless">2.0</cn>
        </apply>
      </apply>
      <!--OpenReleaseChannelFract_calculation-->
    </math>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML" cmeta:id="i_rel_calculation_eq">
      <apply id="i_rel_calculation">
        <eq/>
        <ci>i_rel</ci>
        <apply>
          <times/>
          <ci>alpha_rel</ci>
          <ci>OpenReleaseChannelFract</ci>
          <ci>Ca_rel</ci>
        </apply>
      </apply>
      <!--i_rel_calculation-->
    </math>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML" cmeta:id="RegulatoryBindingSite_calculation_eq">
      <apply id="RegulatoryBindingSite_calculation">
        <eq/>
        <ci>RegulatoryBindingSite</ci>
        <apply>
          <power/>
          <apply>
            <plus/>
            <apply>
              <divide/>
              <ci>Ca_i</ci>
              <apply>
                <plus/>
                <ci>Ca_i</ci>
                <ci>K_mca</ci>
              </apply>
            </apply>
            <apply>
              <times/>
              <apply>
                <minus/>
                <cn cellml:units="dimensionless">1.0</cn>
                <apply>
                  <divide/>
                  <ci>Ca_i</ci>
                  <apply>
                    <plus/>
                    <ci>Ca_i</ci>
                    <ci>K_mca</ci>
                  </apply>
                </apply>
              </apply>
              <apply>
                <divide/>
                <ci>Ca_ds</ci>
                <apply>
                  <plus/>
                  <ci>Ca_ds</ci>
                  <ci>K_mca_ds</ci>
                </apply>
              </apply>
            </apply>
            <!--times-->
          </apply>
          <!--plus-->
          <cn cellml:units="dimensionless">2.0</cn>
        </apply>
        <!--power-->
      </apply>
      <!--RegulatoryBindingSite_calculation-->
    </math>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML" cmeta:id="PrecursorFract_calculation_eq">
      <apply id="PrecursorFract_calculation">
        <eq/>
        <ci>PrecursorFract</ci>
        <apply>
          <minus/>
          <cn cellml:units="dimensionless">1.0</cn>
          <apply>
            <plus/>
            <ci>ActivatorFract</ci>
            <ci>ProductFract</ci>
          </apply>
        </apply>
      </apply>
      <!--PrecursorFract_calculation-->
    </math>
    <math xmlns="http://www.w3.org/1998/Math/MathML" cmeta:id="alpha_act_calculation_eq">
      <apply id="alpha_act_calculation">
        <eq/>
        <ci>alpha_act</ci>
        <apply>
          <times/>
          <cn cellml:units="per_ms">0.5</cn>
          <ci>RegulatoryBindingSite</ci>
        </apply>
      </apply>
      <!--alpha_act_calculation-->
    </math>
    <math xmlns="http://www.w3.org/1998/Math/MathML" cmeta:id="alpha_inact_calculation_eq">
      <apply id="alpha_inact_calculation">
        <eq/>
        <ci>alpha_inact</ci>
        <apply>
          <plus/>
          <cn cellml:units="per_ms">0.06</cn>
          <apply>
            <times/>
            <cn cellml:units="per_ms">0.5</cn>
            <ci>RegulatoryBindingSite</ci>
          </apply>
        </apply>
      </apply>
      <!--alpha_inact_calculation-->
    </math>
    <math xmlns="http://www.w3.org/1998/Math/MathML" cmeta:id="SR_speed_factor_calculation_eq">
      <apply id="SR_speed_factor_calculation">
        <eq/>
        <ci>SR_speed_factor</ci>
        <piecewise>
          <piece>
            <cn cellml:units="dimensionless">5.0</cn>
            <apply>
              <lt/>
              <ci>V</ci>
              <cn cellml:units="mV">-50.0</cn>
            </apply>
          </piece>
          <otherwise>
            <cn cellml:units="dimensionless">1.0</cn>
          </otherwise>
        </piecewise>
      </apply>
      <!--SR_speed_factor_calculation-->
    </math>
    <math xmlns="http://www.w3.org/1998/Math/MathML" cmeta:id="ActivatorFract_diff_eq">
      <apply id="ActivatorFract_diff">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci>time</ci>
          </bvar>
          <ci>ActivatorFract</ci>
        </apply>
        <apply>
          <times/>
          <ci>SR_speed_factor</ci>
          <apply>
            <minus/>
            <apply>
              <times/>
              <ci>alpha_act</ci>
              <ci>PrecursorFract</ci>
            </apply>
            <apply>
              <times/>
              <ci>alpha_inact</ci>
              <ci>ActivatorFract</ci>
            </apply>
          </apply>
        </apply>
      </apply>
      <!--ActivatorFract_diff_eq-->
    </math>
    <math xmlns="http://www.w3.org/1998/Math/MathML" cmeta:id="ProductFract_diff_eq">
      <apply id="ProductFract_diff">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci>time</ci>
          </bvar>
          <ci>ProductFract</ci>
        </apply>
        <apply>
          <times/>
          <ci>SR_speed_factor</ci>
          <apply>
            <minus/>
            <apply>
              <times/>
              <ci>alpha_inact</ci>
              <ci>ActivatorFract</ci>
            </apply>
            <apply>
              <times/>
              <cn cellml:units="per_ms">0.001</cn>
              <ci>ProductFract</ci>
            </apply>
          </apply>
        </apply>
      </apply>
      <!--ProductFract_diff_eq-->
    </math>
  </component>
  
  <component name="stretch_dependent_Ca_current" cmeta:id="stretch_dependent_Ca_current">
    

    <variable units="uA_per_mmsq" public_interface="out" name="i_Ca_stretch"/>
    
    <variable units="mS_per_mmsq" name="g_Ca_stretch" initial_value="154.22e-3"/>
    
    <variable units="mV" public_interface="in" name="E_Ca"/>
    <variable units="ms" public_interface="in" name="time"/>
    <variable units="dimensionless" public_interface="in" name="f_stretch"/>
    <variable units="mV" public_interface="in" name="V"/>
    
    
    <math xmlns="http://www.w3.org/1998/Math/MathML" cmeta:id="i_Ca_stretch_calculation_eq">
      <apply id="i_Ca_stretch_calculation">
        <eq/>
        <ci> i_Ca_stretch </ci>
        <apply>
          <times/>
          <ci> g_Ca_stretch </ci>
          <ci> f_stretch </ci>
          <apply>
            <minus/>
            <ci> V </ci>
            <ci> E_Ca </ci>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="stretch_dependent_K_current" cmeta:id="stretch_dependent_K_current">
    

    <variable units="uA_per_mmsq" public_interface="out" name="i_K_stretch"/>
    
    <variable units="mS_per_mmsq" name="g_K_stretch" initial_value="154.22e-3"/>
    
    <variable units="mV" public_interface="in" name="E_K"/>
    <variable units="ms" public_interface="in" name="time"/>
    <variable units="dimensionless" public_interface="in" name="f_stretch"/>
    <variable units="mV" public_interface="in" name="V"/>
    
    
    <math xmlns="http://www.w3.org/1998/Math/MathML" cmeta:id="i_K_stretch_calculation_eq">
      <apply id="i_K_stretch_calculation">
        <eq/>
        <ci> i_K_stretch </ci>
        <apply>
          <times/>
          <ci> g_K_stretch </ci>
          <ci> f_stretch </ci>
          <apply>
            <minus/>
            <ci> V </ci>
            <ci> E_K </ci>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="stretch_dependent_Na_current" cmeta:id="stretch_dependent_Na_current">
    

    <variable units="uA_per_mmsq" public_interface="out" name="i_Na_stretch"/>
    
    <variable units="mS_per_mmsq" name="g_Na_stretch" initial_value="154.22e-3"/>
    
    <variable units="mV" public_interface="in" name="E_Na"/>
    <variable units="ms" public_interface="in" name="time"/>
    <variable units="dimensionless" public_interface="in" name="f_stretch"/>
    <variable units="mV" public_interface="in" name="V"/>
    
    
    <math xmlns="http://www.w3.org/1998/Math/MathML" cmeta:id="i_Na_stretch_calculation_eq">
      <apply id="i_Na_stretch_calculation">
        <eq/>
        <ci> i_Na_stretch </ci>
        <apply>
          <times/>
          <ci> g_Na_stretch </ci>
          <ci> f_stretch </ci>
          <apply>
            <minus/>
            <ci> V </ci>
            <ci> E_Na </ci>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="stretch_dependent_Ns_current" cmeta:id="stretch_dependent_Ns_current">
    

    <variable units="uA_per_mmsq" public_interface="out" name="i_Ns_stretch"/>
   
    <variable units="mS_per_mmsq" name="g_Ns_stretch" initial_value="154.22e-3"/>
    <variable units="mV" name="E_Ns_stretch" initial_value="-30.0"/>
   
    <variable units="ms" public_interface="in" name="time"/>
    <variable units="dimensionless" public_interface="in" name="f_stretch"/>
    <variable units="mV" public_interface="in" name="V"/>
    
    
    <math xmlns="http://www.w3.org/1998/Math/MathML" cmeta:id="i_Ns_stretch_calculation_eq">
      <apply id="i_Ns_stretch_calculation">
        <eq/>
        <ci> i_Ns_stretch </ci>
        <apply>
          <times/>
          <ci> g_Ns_stretch </ci>
          <ci> f_stretch </ci>
          <apply>
            <minus/>
            <ci> V </ci>
            <ci> E_Ns_stretch </ci>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="stretch_dependent_An_current" cmeta:id="stretch_dependent_An_current">
    

    <variable units="uA_per_mmsq" public_interface="out" name="i_An_stretch"/>
   
    <variable units="mS_per_mmsq" name="g_An_stretch" initial_value="154.22e-3"/>
    <variable units="mV" name="E_An_stretch" initial_value="-20.0"/>
   
    <variable units="ms" public_interface="in" name="time"/>
    <variable units="dimensionless" public_interface="in" name="f_stretch"/>
    <variable units="mV" public_interface="in" name="V"/>
    
    
    <math xmlns="http://www.w3.org/1998/Math/MathML" cmeta:id="i_An_stretch_calculation_eq">
      <apply id="i_An_stretch_calculation">
        <eq/>
        <ci> i_An_stretch </ci>
        <apply>
          <times/>
          <ci> g_An_stretch </ci>
          <ci> f_stretch </ci>
          <apply>
            <minus/>
            <ci> V </ci>
            <ci> E_An_stretch </ci>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="ionic_concentrations" cmeta:id="ionic_concentrations">
    

    <variable units="mM" public_interface="out" name="Na_i" initial_value="5.5997"/>
    <variable units="mM" public_interface="out" name="Na_o" initial_value="140.0"/>
    <variable units="mM" public_interface="out" name="K_i" initial_value="139.3050"/>
    <variable units="mM" public_interface="out" name="K_o" initial_value="4.0"/>
    <variable units="mM" public_interface="out" name="Ca_i" initial_value="0.0000082"/>
    <variable units="mM" public_interface="out" name="Ca_o" initial_value="2.0"/>
    <variable units="mM" public_interface="out" name="Ca_ds" initial_value="0.0000171"/>
    <variable units="mM" public_interface="out" name="Ca_up" initial_value="0.2872393"/>
    <variable units="mM" public_interface="out" name="Ca_rel" initial_value="0.2846761"/>
    
    <variable units="mM" name="Ca_calmod" initial_value="0.0003258"/>
    <variable units="mM_per_ms" name="dy_Ca_calmod"/>
    <variable units="mM" name="Ca_troponin" initial_value="0.0002060"/>
    <variable units="mM_per_ms" name="dy_Ca_troponin"/>
    <variable units="mM" name="CALM" initial_value="0.02"/>
    <variable units="mM" name="TROP" initial_value="0.05"/>
    <variable units="per_mM_per_ms" name="alpha_trop" initial_value="100.0"/>
    <variable units="per_ms" name="beta_trop" initial_value="0.2"/>
    <variable units="per_mm" name="gamma_trop_SL" initial_value="2.5e+3"/>
    <variable units="dimensionless" name="V_i" initial_value="0.49"/>
    <variable units="dimensionless" name="V_rel" initial_value="0.1"/>
    <variable units="dimensionless" name="V_ds" initial_value="0.1"/>
    <variable units="dimensionless" name="V_up" initial_value="0.01"/>
    <variable units="dimensionless" name="V_ecs" initial_value="0.3"/>
    
    <variable units="per_mm" public_interface="in" name="Am"/>
    <variable units="mm" public_interface="in" name="SL"/>
    <variable units="mM_per_ms" public_interface="in" name="i_leak"/>
    <variable units="mM_per_ms" public_interface="in" name="i_decay"/>
    <variable units="mM_per_ms" public_interface="in" name="i_up"/>
    <variable units="mM_per_ms" public_interface="in" name="i_tr"/>
    <variable units="mM_per_ms" public_interface="in" name="i_rel"/>
    <variable units="uA_per_mmsq" public_interface="in" name="i_Na"/>
    <variable units="uA_per_mmsq" public_interface="in" name="i_p_Na"/>
    <variable units="uA_per_mmsq" public_interface="in" name="i_b_Na"/>
    <variable units="uA_per_mmsq" public_interface="in" name="i_NaK"/>
    <variable units="uA_per_mmsq" public_interface="in" name="i_NaCa"/>
    <variable units="uA_per_mmsq" public_interface="in" name="i_NaCa_ds"/>
    <variable units="uA_per_mmsq" public_interface="in" name="i_Ca_L_Na"/>
    <variable units="uA_per_mmsq" public_interface="in" name="i_Ca_L_Na_ds"/>
    <variable units="uA_per_mmsq" public_interface="in" name="i_Na_stretch"/>
    <variable units="uA_per_mmsq" public_interface="in" name="i_to"/>
    <variable units="uA_per_mmsq" public_interface="in" name="i_K_ACh"/>
    <variable units="uA_per_mmsq" public_interface="in" name="i_K_ATP"/>
    <variable units="uA_per_mmsq" public_interface="in" name="i_K1"/>
    <variable units="uA_per_mmsq" public_interface="in" name="i_Ca_L_K"/>
    <variable units="uA_per_mmsq" public_interface="in" name="i_Ca_L_K_ds"/>
    <variable units="uA_per_mmsq" public_interface="in" name="i_b_K"/>
    <variable units="uA_per_mmsq" public_interface="in" name="i_K"/>
    <variable units="uA_per_mmsq" public_interface="in" name="i_K_stretch"/>
    <variable units="uA_per_mmsq" public_interface="in" name="i_Ca_L_Ca"/>
    <variable units="uA_per_mmsq" public_interface="in" name="i_Ca_L_Ca_ds"/>
    <variable units="uA_per_mmsq" public_interface="in" name="i_b_Ca"/>
    <variable units="uA_per_mmsq" public_interface="in" name="i_Ca_stretch"/>
    <variable units="faradays_constant" public_interface="in" name="F"/>
    <variable units="ms" public_interface="in" name="time"/>
    
    
    <math xmlns="http://www.w3.org/1998/Math/MathML" cmeta:id="Na_i_diff_eq">
      <apply id="Na_i_diff">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> Na_i </ci>
        </apply>
        <apply>
          <times/>
          <apply>
            <minus/>
            <apply>
              <divide/>
              <ci> Am </ci>
              <apply>
                <times/>
                <ci> V_i </ci>
                <ci> F </ci>
              </apply>
            </apply>
          </apply>
          <apply>
            <plus/>
            <ci> i_Na </ci>
            <ci> i_p_Na </ci>
            <apply>
              <times/>
              <ci> i_b_Na </ci>
              <apply>
                <divide/>
                <ci>Na_o</ci>
                <cn cellml:units="mM">140.0</cn>
              </apply>
            </apply>
            <apply>
              <times/>
              <ci> i_NaK </ci>
              <cn cellml:units="dimensionless"> 3.0 </cn>
            </apply>
            <apply>
              <times/>
              <ci> i_NaCa </ci>
              <cn cellml:units="dimensionless"> 3.0 </cn>
            </apply>
            <apply>
              <times/>
              <ci> i_NaCa_ds </ci>
              <cn cellml:units="dimensionless"> 3.0 </cn>
            </apply>
            <ci> i_Ca_L_Na </ci>
            <ci> i_Ca_L_Na_ds </ci>
            <ci> i_Na_stretch </ci>
          </apply>
        </apply>
      </apply>
    </math>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML" cmeta:id="K_i_diff_eq">
      <apply id="K_i_diff">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> K_i </ci>
        </apply>
        <apply>
          <times/>
          <apply>
            <minus/>
            <apply>
              <divide/>
              <ci> Am </ci>
              <apply>
                <times/>
                <ci> V_i </ci>
                <ci> F </ci>
              </apply>
            </apply>
          </apply>
          <apply>
            <minus/>
            <apply>
              <plus/>
              <ci> i_K1 </ci>
              <ci> i_to </ci>
              <ci> i_b_K </ci>
              <ci> i_K </ci>
              <ci> i_K_stretch </ci>
              <ci> i_Ca_L_K </ci>
              <ci> i_Ca_L_K_ds </ci>
              <ci> i_K_ACh </ci>
              <ci> i_K_ATP </ci>
            </apply>
            <apply>
              <times/>
              <ci> i_NaK </ci>
              <cn cellml:units="dimensionless"> 2.0 </cn>
            </apply>
          </apply>
        </apply>
      </apply>
    </math>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML" cmeta:id="Ca_i_diff_eq">
      <apply id="Ca_i_diff">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> Ca_i </ci>
        </apply>
        <apply>
          <minus/>
          <apply>
            <plus/>
            <apply>
              <times/>
              <apply>
                <minus/>
                <apply>
                  <divide/>
                  <ci> Am </ci>
                  <apply>
                    <times/>
                    <cn cellml:units="dimensionless"> 2.0 </cn>
                    <ci> V_i </ci>
                    <ci> F </ci>
                  </apply>
                </apply>
              </apply>
              <apply>
                <minus/>
                <apply>
                  <plus/>
                  <ci> i_Ca_L_Ca </ci>
                  <ci> i_b_Ca </ci>
                  <ci> i_Ca_stretch </ci>
                </apply>
                <apply>
                  <times/>
                  <ci> i_NaCa </ci>
                  <cn cellml:units="dimensionless"> 2.0 </cn>
                </apply>
              </apply>
            </apply>
            <!--times-->
            <apply>
              <times/>
              <ci> i_rel </ci>
              <apply>
                <divide/>
                <ci> V_rel </ci>
                <ci> V_i </ci>
              </apply>
            </apply>
            <apply>
              <times/>
              <ci> i_leak </ci>
              <apply>
                <divide/>
                <ci> V_rel </ci>
                <ci> V_i </ci>
              </apply>
            </apply>
            <apply>
              <times/>
              <ci> i_decay </ci>
              <apply>
                <divide/>
                <ci> V_ds </ci>
                <ci> V_i </ci>
              </apply>
            </apply>
          </apply>
          <!--plus-->
          <apply>
            <plus/>
            <ci>i_up</ci>
            <ci>dy_Ca_calmod</ci>
            <ci>dy_Ca_troponin</ci>
          </apply>
          <!--plus-->
        </apply>
        <!--minus-->
      </apply>
    </math>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML" cmeta:id="Ca_ds_diff_eq">
      <apply id="Ca_ds_diff">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> Ca_ds </ci>
        </apply>
        <apply>
          <minus/>
          <apply>
            <times/>
            <apply>
              <minus/>
              <apply>
                <divide/>
                <ci> Am </ci>
                <apply>
                  <times/>
                  <cn cellml:units="dimensionless"> 2.0 </cn>
                  <ci> V_ds </ci>
                  <ci> F </ci>
                </apply>
              </apply>
            </apply>
            <apply>
              <minus/>
              <ci> i_Ca_L_Ca_ds </ci>
              <apply>
                <times/>
                <ci> i_NaCa_ds </ci>
                <cn cellml:units="dimensionless"> 2.0 </cn>
              </apply>
            </apply>
          </apply>
          <!--times-->
          <ci>i_decay</ci>
        </apply>
        <!--minus-->
      </apply>
    </math>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML" cmeta:id="Ca_up_diff_eq">
      <apply id="Ca_up_diff">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> Ca_up </ci>
        </apply>
        <apply>
          <minus/>
          <apply>
            <times/>
            <apply>
              <divide/>
              <ci> V_i </ci>
              <ci> V_up </ci>
            </apply>
            <ci> i_up </ci>
          </apply>
          <ci> i_tr </ci>
        </apply>
      </apply>
    </math>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML" cmeta:id="Ca_rel_diff_eq">
      <apply id="Ca_rel_diff">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> Ca_rel </ci>
        </apply>
        <apply>
          <minus/>
          <apply>
            <times/>
            <apply>
              <divide/>
              <ci> V_up </ci>
              <ci> V_rel </ci>
            </apply>
            <ci> i_tr </ci>
          </apply>
          <apply>
            <plus/>
            <ci> i_rel </ci>
            <ci>i_leak</ci>
          </apply>
        </apply>
      </apply>
    </math>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML" cmeta:id="dy_Ca_calmod_calculation_eq">
      <apply id="dy_Ca_calmod_calculation">
        <eq/>
        <ci>dy_Ca_calmod</ci>
        <apply>
          <minus/>
          <apply>
            <times/>
            <cn cellml:units="per_mM_per_ms"> 100.0 </cn>
            <ci> Ca_i </ci>
            <apply>
              <minus/>
              <ci> CALM </ci>
              <ci> Ca_calmod </ci>
            </apply>
          </apply>
          <apply>
            <times/>
            <cn cellml:units="per_ms"> 0.05 </cn>
            <ci> Ca_calmod </ci>
          </apply>
        </apply>
      </apply>
    </math>
    <math xmlns="http://www.w3.org/1998/Math/MathML" cmeta:id="Ca_calmod_diff_eq">
      <apply id="Ca_calmod_diff">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> Ca_calmod </ci>
        </apply>
        <ci>dy_Ca_calmod</ci>
      </apply>
    </math>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML" cmeta:id="dy_Ca_troponoin_calculation_eq">
      <apply id="dy_Ca_troponoin_calculation">
        <eq/>
        <ci>dy_Ca_troponin</ci>
        <apply>
          <minus/>
          <apply>
            <times/>
            <ci> alpha_trop </ci>
            <ci> Ca_i </ci>
            <apply>
              <minus/>
              <ci> TROP </ci>
              <ci> Ca_troponin </ci>
            </apply>
          </apply>
          <apply>
            <times/>
            <ci> beta_trop </ci>
            <ci> Ca_troponin </ci>
          </apply>
        </apply>
      </apply>
    </math>
    <math xmlns="http://www.w3.org/1998/Math/MathML" cmeta:id="Ca_troponin_diff_eq">
      <apply id="Ca_troponin_diff">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> Ca_troponin </ci>
        </apply>
        <ci>dy_Ca_troponin</ci>
      </apply>
    </math>
  </component>
    
  <group>
    <relationship_ref relationship="containment"/>
    <component_ref component="membrane">
      <component_ref component="fast_sodium_current">
        <component_ref component="fast_sodium_current_m_gate"/>
        <component_ref component="fast_sodium_current_h_gate"/>
      </component_ref>
      <component_ref component="L_type_Ca_channel">
        <component_ref component="L_type_Ca_channel_d_gate"/>
        <component_ref component="L_type_Ca_channel_f_gate"/>
        <component_ref component="L_type_Ca_channel_f2_gate"/>
        <component_ref component="L_type_Ca_channel_f2ds_gate"/>
      </component_ref>
      <component_ref component="rapid_time_dependent_potassium_current">
        <component_ref component="rapid_time_dependent_potassium_current_Xr1_gate"/>
        <component_ref component="rapid_time_dependent_potassium_current_Xr2_gate"/>
      </component_ref>
      <component_ref component="slow_time_dependent_potassium_current">
        <component_ref component="slow_time_dependent_potassium_current_Xs_gate"/>
      </component_ref>
      <component_ref component="total_potassium_current"/>
      <component_ref component="potassium_background_current"/>
      <component_ref component="time_independent_potassium_current"/>
      <component_ref component="sodium_activated_potassium_current"/>
      <component_ref component="ATP_dependent_potassium_current"/>
      <component_ref component="persistent_sodium_current"/>
      <component_ref component="sodium_background_current"/>
      <component_ref component="calcium_background_current"/>
      <component_ref component="transient_outward_current">
        <component_ref component="transient_outward_current_s_gate"/>
        <component_ref component="transient_outward_current_r_gate"/>
      </component_ref>
      <component_ref component="ACh_dependent_potassium_current">
        <component_ref component="ACh_dependent_potassium_current_xACh1_gate"/>
        <component_ref component="ACh_dependent_potassium_current_xACh2_gate"/>
      </component_ref>
      <component_ref component="sodium_potassium_pump"/>
      <component_ref component="sodium_calcium_pump"/>
      <component_ref component="calcium_uptake_from_myoplasm_to_NSR"/>
      <component_ref component="calcium_transfer_from_NSR_to_JSR"/>
      <component_ref component="calcium_leak_from_SR"/>
      <component_ref component="calcium_release_from_SR"/>
      <component_ref component="ds_decay_flux"/>
      <component_ref component="length_dependence"/>
      <component_ref component="stretch_dependent_Ca_current"/>
      <component_ref component="stretch_dependent_K_current"/>
      <component_ref component="stretch_dependent_Na_current"/>
      <component_ref component="stretch_dependent_Ns_current"/>
      <component_ref component="stretch_dependent_An_current"/>
      <component_ref component="ionic_concentrations"/>
      <component_ref component="reversal_potentials"/>
    </component_ref>
  </group>
  
  <group>
    <relationship_ref relationship="encapsulation"/>
    <component_ref component="fast_sodium_current">
      <component_ref component="fast_sodium_current_m_gate"/>
      <component_ref component="fast_sodium_current_h_gate"/>
    </component_ref>
    <component_ref component="L_type_Ca_channel">
      <component_ref component="L_type_Ca_channel_d_gate"/>
      <component_ref component="L_type_Ca_channel_f_gate"/>
      <component_ref component="L_type_Ca_channel_f2_gate"/>
      <component_ref component="L_type_Ca_channel_f2ds_gate"/>
    </component_ref>
    <component_ref component="rapid_time_dependent_potassium_current">
      <component_ref component="rapid_time_dependent_potassium_current_Xr1_gate"/>
      <component_ref component="rapid_time_dependent_potassium_current_Xr2_gate"/>
    </component_ref>
    <component_ref component="slow_time_dependent_potassium_current">
      <component_ref component="slow_time_dependent_potassium_current_Xs_gate"/>
    </component_ref>
    <component_ref component="ACh_dependent_potassium_current">
      <component_ref component="ACh_dependent_potassium_current_xACh1_gate"/>
      <component_ref component="ACh_dependent_potassium_current_xACh2_gate"/>
    </component_ref>
    <component_ref component="transient_outward_current">
      <component_ref component="transient_outward_current_s_gate"/>
      <component_ref component="transient_outward_current_r_gate"/>
    </component_ref>
  </group>
  
  <!--
    "Time" is passed from the "environment" component into the
    "membrane" and current components.
  -->
  
  <connection>
    <map_components component_2="environment" component_1="membrane"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="environment" component_1="reversal_potentials"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="environment" component_1="fast_sodium_current"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>

  <connection>
    <map_components component_2="environment" component_1="L_type_Ca_channel"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="environment" component_1="rapid_time_dependent_potassium_current"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="environment" component_1="slow_time_dependent_potassium_current"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="environment" component_1="time_independent_potassium_current"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="environment" component_1="sodium_activated_potassium_current"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="environment" component_1="ATP_dependent_potassium_current"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="environment" component_1="ACh_dependent_potassium_current"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="environment" component_1="persistent_sodium_current"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="environment" component_1="sodium_background_current"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="environment" component_1="potassium_background_current"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="environment" component_1="calcium_background_current"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="environment" component_1="transient_outward_current"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="environment" component_1="stretch_dependent_Ca_current"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="environment" component_1="stretch_dependent_K_current"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="environment" component_1="stretch_dependent_Na_current"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="environment" component_1="stretch_dependent_Ns_current"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="environment" component_1="stretch_dependent_An_current"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="environment" component_1="sodium_potassium_pump"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="environment" component_1="sodium_calcium_pump"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="environment" component_1="ionic_concentrations"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="environment" component_1="total_potassium_current"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="environment" component_1="calcium_uptake_from_myoplasm_to_NSR"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="environment" component_1="calcium_transfer_from_NSR_to_JSR"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="environment" component_1="calcium_release_from_SR"/>
    <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="reversal_potentials" component_1="membrane"/>
    <map_variables variable_2="R" variable_1="R"/>
    <map_variables variable_2="T" variable_1="T"/>
    <map_variables variable_2="F" variable_1="F"/>
  </connection>
  
  <connection>
    <map_components component_2="fast_sodium_current" component_1="membrane"/>
    <map_variables variable_2="V" variable_1="V"/>
    <map_variables variable_2="i_Na" variable_1="i_Na"/>
  </connection>
  
  <connection>
    <map_components component_2="L_type_Ca_channel" component_1="membrane"/>
    <map_variables variable_2="V" variable_1="V"/>
    <map_variables variable_2="i_Ca_L_Ca" variable_1="i_Ca_L_Ca"/>
    <map_variables variable_2="i_Ca_L_Na" variable_1="i_Ca_L_Na"/>
    <map_variables variable_2="i_Ca_L_K" variable_1="i_Ca_L_K"/>
    <map_variables variable_2="i_Ca_L_Ca_ds" variable_1="i_Ca_L_Ca_ds"/>
    <map_variables variable_2="i_Ca_L_Na_ds" variable_1="i_Ca_L_Na_ds"/>
    <map_variables variable_2="i_Ca_L_K_ds" variable_1="i_Ca_L_K_ds"/>
    <map_variables variable_2="R" variable_1="R"/>
    <map_variables variable_2="T" variable_1="T"/>
    <map_variables variable_2="F" variable_1="F"/>
  </connection>
  
  <connection>
    <map_components component_2="total_potassium_current" component_1="membrane"/>
    <map_variables variable_2="i_K" variable_1="i_K"/>
  </connection>
  
  <connection>
    <map_components component_2="rapid_time_dependent_potassium_current" component_1="membrane"/>
    <map_variables variable_2="V" variable_1="V"/>
  </connection>
  
  <connection>
    <map_components component_2="slow_time_dependent_potassium_current" component_1="membrane"/>
    <map_variables variable_2="V" variable_1="V"/>
  </connection>
  
  <connection>
    <map_components component_2="transient_outward_current" component_1="membrane"/>
    <map_variables variable_2="V" variable_1="V"/>
    <map_variables variable_2="i_to" variable_1="i_to"/>
  </connection>
  
  <connection>
    <map_components component_2="sodium_activated_potassium_current" component_1="membrane"/>
    <map_variables variable_2="V" variable_1="V"/>
  </connection>
  
  <connection>
    <map_components component_2="ATP_dependent_potassium_current" component_1="membrane"/>
    <map_variables variable_2="V" variable_1="V"/>
    <map_variables variable_2="i_K_ATP" variable_1="i_K_ATP"/>
  </connection>
  
  <connection>
    <map_components component_2="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="persistent_sodium_current" component_1="membrane"/>
    <map_variables variable_2="V" variable_1="V"/>
    <map_variables variable_2="i_p_Na" variable_1="i_p_Na"/>
  </connection>
  
  <connection>
    <map_components component_2="sodium_background_current" component_1="membrane"/>
    <map_variables variable_2="V" variable_1="V"/>
     <map_variables variable_2="i_b_Na" variable_1="i_b_Na"/>
  </connection>
  
  <connection>
    <map_components component_2="calcium_background_current" component_1="membrane"/>
    <map_variables variable_2="V" variable_1="V"/>
     <map_variables variable_2="i_b_Ca" variable_1="i_b_Ca"/>
  </connection>
  
  <connection>
    <map_components component_2="potassium_background_current" component_1="membrane"/>
    <map_variables variable_2="V" variable_1="V"/>
     <map_variables variable_2="i_b_K" variable_1="i_b_K"/>
  </connection>
  
  <connection>
    <map_components component_2="sodium_potassium_pump" component_1="membrane"/>
    <map_variables variable_2="i_NaK" variable_1="i_NaK"/>
  </connection>
  
  <connection>
    <map_components component_2="sodium_calcium_pump" 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="i_NaCa_ds" variable_1="i_NaCa_ds"/>
    <map_variables variable_2="R" variable_1="R"/>
    <map_variables variable_2="T" variable_1="T"/>
    <map_variables variable_2="F" variable_1="F"/>
  </connection>
  
  <connection>
    <map_components component_2="ACh_dependent_potassium_current" component_1="membrane"/>
    <map_variables variable_2="V" variable_1="V"/>
    <map_variables variable_2="i_K_ACh" variable_1="i_K_ACh"/>
    <map_variables variable_2="R" variable_1="R"/>
    <map_variables variable_2="T" variable_1="T"/>
    <map_variables variable_2="F" variable_1="F"/>
  </connection>
  
  <connection>
    <map_components component_2="stretch_dependent_Ca_current" component_1="membrane"/>
    <map_variables variable_2="i_Ca_stretch" variable_1="i_Ca_stretch"/>
    <map_variables variable_2="V" variable_1="V"/>
  </connection>
  
  <connection>
    <map_components component_2="stretch_dependent_K_current" component_1="membrane"/>
    <map_variables variable_2="i_K_stretch" variable_1="i_K_stretch"/>
    <map_variables variable_2="V" variable_1="V"/>
  </connection>
  
  <connection>
    <map_components component_2="stretch_dependent_Na_current" component_1="membrane"/>
    <map_variables variable_2="i_Na_stretch" variable_1="i_Na_stretch"/>
    <map_variables variable_2="V" variable_1="V"/>
  </connection>
  
  <connection>
    <map_components component_2="stretch_dependent_Ns_current" component_1="membrane"/>
    <map_variables variable_2="i_Ns_stretch" variable_1="i_Ns_stretch"/>
    <map_variables variable_2="V" variable_1="V"/>
  </connection>
  
  <connection>
    <map_components component_2="stretch_dependent_An_current" component_1="membrane"/>
    <map_variables variable_2="i_An_stretch" variable_1="i_An_stretch"/>
    <map_variables variable_2="V" variable_1="V"/>
  </connection>
  
  <connection>
    <map_components component_2="ionic_concentrations" component_1="membrane"/>
    <map_variables variable_2="F" variable_1="F"/>
    <map_variables variable_2="Am" variable_1="Am"/>
  </connection>
          
  <!--  Several variables are passed between the sibling components. -->        
  
  <connection>
    <map_components component_2="sodium_background_current" component_1="reversal_potentials"/>
    <map_variables variable_2="E_Na" variable_1="E_Na"/>
  </connection>
  
  <connection>
    <map_components component_2="persistent_sodium_current" component_1="reversal_potentials"/>
    <map_variables variable_2="E_Na" variable_1="E_Na"/>
  </connection>
  
  <connection>
    <map_components component_2="stretch_dependent_Na_current" component_1="reversal_potentials"/>
    <map_variables variable_2="E_Na" variable_1="E_Na"/>
  </connection>
  
  <connection>
    <map_components component_2="stretch_dependent_K_current" component_1="reversal_potentials"/>
    <map_variables variable_2="E_K" variable_1="E_K"/>
  </connection>
  
  <connection>
    <map_components component_2="potassium_background_current" component_1="reversal_potentials"/>
    <map_variables variable_2="E_K" variable_1="E_K"/>
  </connection>
  
  <connection>
    <map_components component_2="sodium_activated_potassium_current" component_1="reversal_potentials"/>
    <map_variables variable_2="E_K" variable_1="E_K"/>
  </connection>
  
  <connection>
    <map_components component_2="time_independent_potassium_current" component_1="reversal_potentials"/>
    <map_variables variable_2="E_K" variable_1="E_K"/>
  </connection>
  
  <connection>
    <map_components component_2="rapid_time_dependent_potassium_current" component_1="reversal_potentials"/>
    <map_variables variable_2="E_K" variable_1="E_K"/>
  </connection>
  
  <connection>
    <map_components component_2="transient_outward_current" component_1="reversal_potentials"/>
    <map_variables variable_2="E_K" variable_1="E_K"/>
  </connection>
  
  <connection>
    <map_components component_2="ACh_dependent_potassium_current" component_1="reversal_potentials"/>
    <map_variables variable_2="E_K" variable_1="E_K"/>
  </connection>
  
  <connection>
    <map_components component_2="slow_time_dependent_potassium_current" component_1="reversal_potentials"/>
    <map_variables variable_2="E_Ks" variable_1="E_Ks"/>
  </connection>
  
  <connection>
    <map_components component_2="calcium_background_current" component_1="reversal_potentials"/>
    <map_variables variable_2="E_Ca" variable_1="E_Ca"/>
  </connection>
  
  <connection>
    <map_components component_2="stretch_dependent_Ca_current" component_1="reversal_potentials"/>
    <map_variables variable_2="E_Ca" variable_1="E_Ca"/>
  </connection>
  
  <connection>
    <map_components component_2="fast_sodium_current" component_1="reversal_potentials"/>
    <map_variables variable_2="E_mh" variable_1="E_mh"/>
  </connection>
  
  <connection>
    <map_components component_2="ionic_concentrations" component_1="reversal_potentials"/>
    <map_variables variable_2="K_o" variable_1="K_o"/>
    <map_variables variable_2="Na_o" variable_1="Na_o"/>
    <map_variables variable_2="K_i" variable_1="K_i"/>
    <map_variables variable_2="Na_i" variable_1="Na_i"/>
    <map_variables variable_2="Ca_o" variable_1="Ca_o"/>
    <map_variables variable_2="Ca_i" variable_1="Ca_i"/>
  </connection>
  
  <connection>
    <map_components component_2="ionic_concentrations" component_1="time_independent_potassium_current"/>
    <map_variables variable_2="K_o" variable_1="K_o"/>
    <map_variables variable_2="i_K1" variable_1="i_K1"/>
  </connection>
  
  <connection>
    <map_components component_2="ACh_dependent_potassium_current" component_1="time_independent_potassium_current"/>
    <map_variables variable_2="K_mk1" variable_1="K_mk1"/>
  </connection>
  
  <connection>
    <map_components component_2="ionic_concentrations" component_1="total_potassium_current"/>
    <map_variables variable_2="i_K" variable_1="i_K"/>
  </connection>
  
  <connection>
    <map_components component_2="rapid_time_dependent_potassium_current" component_1="total_potassium_current"/>
    <map_variables variable_2="i_Kr" variable_1="i_Kr"/>
  </connection>
  
  <connection>
    <map_components component_2="slow_time_dependent_potassium_current" component_1="total_potassium_current"/>
    <map_variables variable_2="i_Ks" variable_1="i_Ks"/>
  </connection>
  
  <connection>
    <map_components component_2="sodium_activated_potassium_current" component_1="total_potassium_current"/>
    <map_variables variable_2="i_KNa" variable_1="i_KNa"/>
  </connection>
  
  <connection>
    <map_components component_2="ionic_concentrations" component_1="potassium_background_current"/>
    <map_variables variable_2="i_b_K" variable_1="i_b_K"/>
  </connection>
  
  <connection>
    <map_components component_2="ionic_concentrations" component_1="sodium_activated_potassium_current"/>
    <map_variables variable_2="Na_i" variable_1="Na_i"/>
  </connection>
  
  <connection>
    <map_components component_2="ionic_concentrations" component_1="fast_sodium_current"/>
    <map_variables variable_2="i_Na" variable_1="i_Na"/>
  </connection>
  
  <connection>
    <map_components component_2="ionic_concentrations" component_1="persistent_sodium_current"/>
    <map_variables variable_2="i_p_Na" variable_1="i_p_Na"/>
  </connection>
  
  <connection>
    <map_components component_2="ionic_concentrations" component_1="sodium_background_current"/>
    <map_variables variable_2="i_b_Na" variable_1="i_b_Na"/>
    <map_variables variable_2="Na_i" variable_1="Na_i"/>
  </connection>
  
  <connection>
    <map_components component_2="ionic_concentrations" component_1="L_type_Ca_channel"/>
    <map_variables variable_2="i_Ca_L_Ca" variable_1="i_Ca_L_Ca"/>
    <map_variables variable_2="i_Ca_L_K" variable_1="i_Ca_L_K"/>
    <map_variables variable_2="i_Ca_L_Na" variable_1="i_Ca_L_Na"/>
    <map_variables variable_2="i_Ca_L_Ca_ds" variable_1="i_Ca_L_Ca_ds"/>
    <map_variables variable_2="i_Ca_L_Na_ds" variable_1="i_Ca_L_Na_ds"/>
    <map_variables variable_2="i_Ca_L_K_ds" variable_1="i_Ca_L_K_ds"/>
    <map_variables variable_2="Ca_o" variable_1="Ca_o"/>
    <map_variables variable_2="Na_o" variable_1="Na_o"/>
    <map_variables variable_2="K_o" variable_1="K_o"/>
    <map_variables variable_2="Ca_i" variable_1="Ca_i"/>
    <map_variables variable_2="Na_i" variable_1="Na_i"/>
    <map_variables variable_2="K_i" variable_1="K_i"/>
    <map_variables variable_2="Ca_ds" variable_1="Ca_ds"/>
  </connection>
  
  <connection>
    <map_components component_2="ionic_concentrations" component_1="calcium_background_current"/>
    <map_variables variable_2="i_b_Ca" variable_1="i_b_Ca"/>
  </connection>
  
  <connection>
    <map_components component_2="ionic_concentrations" component_1="transient_outward_current"/>
    <map_variables variable_2="Ca_i" variable_1="Ca_i"/>
    <map_variables variable_2="i_to" variable_1="i_to"/>
  </connection>
  
  <connection>
    <map_components component_2="ionic_concentrations" component_1="ACh_dependent_potassium_current"/>
    <map_variables variable_2="K_o" variable_1="K_o"/>
    <map_variables variable_2="i_K_ACh" variable_1="i_K_ACh"/>
  </connection>
  
  <connection>
    <map_components component_2="ionic_concentrations" component_1="ATP_dependent_potassium_current"/>
    <map_variables variable_2="i_K_ATP" variable_1="i_K_ATP"/>
  </connection>
  
  <connection>
    <map_components component_2="ionic_concentrations" component_1="sodium_potassium_pump"/>
    <map_variables variable_2="i_NaK" variable_1="i_NaK"/>
    <map_variables variable_2="K_o" variable_1="K_o"/>
    <map_variables variable_2="Na_i" variable_1="Na_i"/>
  </connection> 
  
  <connection>
    <map_components component_2="ionic_concentrations" component_1="sodium_calcium_pump"/>
    <map_variables variable_2="i_NaCa" variable_1="i_NaCa"/>
    <map_variables variable_2="i_NaCa_ds" variable_1="i_NaCa_ds"/>
    <map_variables variable_2="Ca_i" variable_1="Ca_i"/>
    <map_variables variable_2="Na_i" variable_1="Na_i"/>
    <map_variables variable_2="Ca_o" variable_1="Ca_o"/>
    <map_variables variable_2="Na_o" variable_1="Na_o"/>
    <map_variables variable_2="Ca_ds" variable_1="Ca_ds"/>
  </connection>
  
  <connection>
    <map_components component_2="ionic_concentrations" component_1="calcium_uptake_from_myoplasm_to_NSR"/>
    <map_variables variable_2="i_up" variable_1="i_up"/>
    <map_variables variable_2="Ca_i" variable_1="Ca_i"/>
    <map_variables variable_2="Ca_up" variable_1="Ca_up"/>
  </connection>
  
  <connection>
    <map_components component_2="ionic_concentrations" component_1="calcium_transfer_from_NSR_to_JSR"/>
    <map_variables variable_2="i_tr" variable_1="i_tr"/>
    <map_variables variable_2="Ca_rel" variable_1="Ca_rel"/>
    <map_variables variable_2="Ca_up" variable_1="Ca_up"/>
  </connection>
  
  <connection>
    <map_components component_2="ionic_concentrations" component_1="calcium_release_from_SR"/>
    <map_variables variable_2="i_rel" variable_1="i_rel"/>
    <map_variables variable_2="Ca_rel" variable_1="Ca_rel"/>
  </connection>
  
  <connection>
    <map_components component_2="ionic_concentrations" component_1="calcium_leak_from_SR"/>
    <map_variables variable_2="i_leak" variable_1="i_leak"/>
    <map_variables variable_2="Ca_rel" variable_1="Ca_rel"/>
  </connection>
  
  <connection>
    <map_components component_2="length_dependence" component_1="calcium_leak_from_SR"/>
    <map_variables variable_2="gamma_SR_SL" variable_1="gamma_SR_SL"/>
    <map_variables variable_2="SL" variable_1="SL"/>
  </connection>
  
  <connection>
    <map_components component_2="ionic_concentrations" component_1="ds_decay_flux"/>
    <map_variables variable_2="i_decay" variable_1="i_decay"/>
    <map_variables variable_2="Ca_ds" variable_1="Ca_ds"/>
  </connection>
  
  <connection>
    <map_components component_2="ionic_concentrations" component_1="calcium_release_from_SR"/>
    <map_variables variable_2="Ca_i" variable_1="Ca_i"/>
    <map_variables variable_2="Ca_ds" variable_1="Ca_ds"/>
  </connection> 
  
  <connection>
    <map_components component_2="membrane" component_1="calcium_release_from_SR"/>
    <map_variables variable_2="V" variable_1="V"/>
  </connection> 
  
  <connection>
    <map_components component_2="ionic_concentrations" component_1="length_dependence"/>
    <map_variables variable_2="SL" variable_1="SL"/>
  </connection>
  
  <connection>
    <map_components component_2="stretch_dependent_Ca_current" component_1="length_dependence"/>
    <map_variables variable_2="f_stretch" variable_1="f_stretch"/>
  </connection>
  
  <connection>
    <map_components component_2="stretch_dependent_Na_current" component_1="length_dependence"/>
    <map_variables variable_2="f_stretch" variable_1="f_stretch"/>
  </connection>
  
  <connection>
    <map_components component_2="stretch_dependent_K_current" component_1="length_dependence"/>
    <map_variables variable_2="f_stretch" variable_1="f_stretch"/>
  </connection> 
  
  <connection>
    <map_components component_2="stretch_dependent_Ns_current" component_1="length_dependence"/>
    <map_variables variable_2="f_stretch" variable_1="f_stretch"/>
  </connection>
  
  <connection>
    <map_components component_2="stretch_dependent_An_current" component_1="length_dependence"/>
    <map_variables variable_2="f_stretch" variable_1="f_stretch"/>
  </connection>
  
  <connection>
    <map_components component_2="ionic_concentrations" component_1="stretch_dependent_Ca_current"/>
    <map_variables variable_2="i_Ca_stretch" variable_1="i_Ca_stretch"/>
  </connection>
  
  <connection>
    <map_components component_2="ionic_concentrations" component_1="stretch_dependent_K_current"/>
    <map_variables variable_2="i_K_stretch" variable_1="i_K_stretch"/>
  </connection>
  
  <connection>
    <map_components component_2="ionic_concentrations" component_1="stretch_dependent_Na_current"/>
    <map_variables variable_2="i_Na_stretch" variable_1="i_Na_stretch"/>
  </connection>
  
  <!--  
    Various variables are passed between parent components and their            
    encapsulated gates.
  -->
  
  <connection>
    <map_components component_2="rapid_time_dependent_potassium_current_Xr1_gate" component_1="rapid_time_dependent_potassium_current"/>
    <map_variables variable_2="Xr1" variable_1="Xr1"/>
    <map_variables variable_2="time" variable_1="time"/>
    <map_variables variable_2="V" variable_1="V"/>
  </connection>
  
  <connection>
    <map_components component_2="rapid_time_dependent_potassium_current_Xr2_gate" component_1="rapid_time_dependent_potassium_current"/>
    <map_variables variable_2="Xr2" variable_1="Xr2"/>
    <map_variables variable_2="time" variable_1="time"/>
    <map_variables variable_2="V" variable_1="V"/>
  </connection>
  
  <connection>
    <map_components component_2="slow_time_dependent_potassium_current_Xs_gate" component_1="slow_time_dependent_potassium_current"/>
    <map_variables variable_2="Xs" variable_1="Xs"/>
    <map_variables variable_2="time" variable_1="time"/>
    <map_variables variable_2="V" variable_1="V"/>
  </connection>
  
  <connection>
    <map_components component_2="fast_sodium_current_m_gate" component_1="fast_sodium_current"/>
    <map_variables variable_2="m" variable_1="m"/>
    <map_variables variable_2="time" variable_1="time"/>
    <map_variables variable_2="V" variable_1="V"/>
  </connection>

  <connection>
    <map_components component_2="fast_sodium_current_h_gate" component_1="fast_sodium_current"/>
    <map_variables variable_2="h" variable_1="h"/>
    <map_variables variable_2="time" variable_1="time"/>
    <map_variables variable_2="V" variable_1="V"/>
  </connection>
  
  <connection>
    <map_components component_2="L_type_Ca_channel_d_gate" component_1="L_type_Ca_channel"/>
    <map_variables variable_2="d" variable_1="d"/>
    <map_variables variable_2="time" variable_1="time"/>
    <map_variables variable_2="V" variable_1="V"/>
  </connection>
  
  <connection>
    <map_components component_2="L_type_Ca_channel_f_gate" component_1="L_type_Ca_channel"/>
    <map_variables variable_2="f" variable_1="f"/>
    <map_variables variable_2="time" variable_1="time"/>
    <map_variables variable_2="Ca_i" variable_1="Ca_i"/>
    <map_variables variable_2="V" variable_1="V"/>
    <map_variables variable_2="K_cachoff" variable_1="K_cachoff"/>
  </connection>
  
  <connection>
    <map_components component_2="L_type_Ca_channel_f2_gate" component_1="L_type_Ca_channel"/>
    <map_variables variable_2="f2" variable_1="f2"/>
    <map_variables variable_2="time" variable_1="time"/>
    <map_variables variable_2="Ca_i" variable_1="Ca_i"/>
    <map_variables variable_2="V" variable_1="V"/>
    <map_variables variable_2="K_cachoff" variable_1="K_cachoff"/>
  </connection>
  
  <connection>
    <map_components component_2="L_type_Ca_channel_f2ds_gate" component_1="L_type_Ca_channel"/>
    <map_variables variable_2="f2ds" variable_1="f2ds"/>
    <map_variables variable_2="time" variable_1="time"/>
    <map_variables variable_2="Ca_ds" variable_1="Ca_ds"/>
    <map_variables variable_2="V" variable_1="V"/>
  </connection>
  
  <connection>
    <map_components component_2="ACh_dependent_potassium_current" component_1="L_type_Ca_channel"/>
    <map_variables variable_2="ACh" variable_1="ACh"/>
  </connection>
  
  <connection>
    <map_components component_2="transient_outward_current_s_gate" component_1="transient_outward_current"/>
    <map_variables variable_2="s" variable_1="s"/>
    <map_variables variable_2="time" variable_1="time"/>
    <map_variables variable_2="V" variable_1="V"/>
  </connection>
  
  <connection>
    <map_components component_2="transient_outward_current_r_gate" component_1="transient_outward_current"/>
    <map_variables variable_2="r" variable_1="r"/>
    <map_variables variable_2="time" variable_1="time"/>
    <map_variables variable_2="V" variable_1="V"/>
  </connection>
  
  <connection>
    <map_components component_2="ACh_dependent_potassium_current_xACh1_gate" component_1="ACh_dependent_potassium_current"/>
    <map_variables variable_2="x_ACh1" variable_1="x_ACh1"/>
    <map_variables variable_2="time" variable_1="time"/>
    <map_variables variable_2="V" variable_1="V"/>
  </connection>
  <connection>
    <map_components component_2="ACh_dependent_potassium_current_xACh2_gate" component_1="ACh_dependent_potassium_current"/>
    <map_variables variable_2="x_ACh2" variable_1="x_ACh2"/>
    <map_variables variable_2="time" variable_1="time"/>
    <map_variables variable_2="V" variable_1="V"/>
  </connection>
      





<rdf:RDF>
  <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:li rdf:resource="rdf:#author4Vcard"/>
  </rdf:Seq>
  <rdf:Description rdf:about="#calcium_release_from_SR">
    <cmeta:comment rdf:resource="rdf:#2e14661a-49c2-4767-86d4-7fcc7027eba9"/>
  </rdf:Description>
  <rdf:Description rdf:about="#i_to_calculation_eq">
    <cmeta:comment rdf:resource="rdf:#af9c3c31-c6bb-4c9d-9c27-2e904adcec46"/>
  </rdf:Description>
  <rdf:Description rdf:about="#membrane">
    <cmeta:comment rdf:resource="rdf:#f75a5fb6-2029-43ca-b9cb-5ffa72ca92f1"/>
  </rdf:Description>
  <rdf:Description rdf:about="#ACh_dependent_potassium_current_xACh1_gate">
    <cmeta:comment rdf:resource="rdf:#ebd58747-9f34-4a1e-acca-6705214db1ed"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#88eca42f-2bfb-4131-8801-7b68b13ba881">
    <dcterms:W3CDTF>2002-03-01</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#67a795b5-d23c-4c1d-8107-5efc975ec984">
    <rdf:value>
            The sodium ion component of the total current through the channel into
            the diadic subspace.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#66b34f79-38b2-4bd8-8152-811b026043ff">
    <rdf:value>
            Regulatory binding sites.
          </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:#b06662dd-5000-4cf2-91b4-e8e57ca5c7a9">
    <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:#f06d31ad-ec7c-49a0-8b2e-517b0ec3446f">
    <dcterms:W3CDTF>2002-01-04</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#2b9d3310-dbad-4680-82f8-5286f059a431">
    <vCard:N rdf:resource="rdf:#735dee69-ddee-475c-a8e0-2a57e36428c3"/>
  </rdf:Description>
  <rdf:Description rdf:about="#i_leak_calculation_eq">
    <cmeta:comment rdf:resource="rdf:#1a60361b-64f9-478a-a87e-657e818183ad"/>
  </rdf:Description>
  <rdf:Description rdf:about="#Ca_i_diff_eq">
    <cmeta:comment rdf:resource="rdf:#a0ec4572-a007-4ba3-9af2-bf9984fdb07f"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#author3VcardN">
    <vCard:Given>Peter</vCard:Given>
    <vCard:Family>Kohl</vCard:Family>
  </rdf:Description>
  <rdf:Description rdf:about="#i_KNa_calculation_eq">
    <cmeta:comment rdf:resource="rdf:#7eb621bb-57af-4b0e-a58c-0cf72984b47e"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#04ee1311-52db-48a6-be5a-39b8b796aad9">
    <vCard:N rdf:resource="rdf:#9f42bedb-4cb8-40ae-9901-b680e66a8315"/>
  </rdf:Description>
  <rdf:Description rdf:about="#transient_outward_current_s_gate">
    <cmeta:comment rdf:resource="rdf:#35fce50d-7b61-4963-be42-3bfee47a2a4b"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#2e14661a-49c2-4767-86d4-7fcc7027eba9">
    <rdf:value>
            The calcium induced calcium release from the junctional SR into 
            the bulk myoplasm. Assumed that activation of calcium release sites
            by the diadic subspace calcium triggers calcium release in the whole
            cytosol.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#8bcca774-860a-452a-b77b-2a904ec12f4d">
    <dcterms:modified rdf:resource="rdf:#c69bffaa-0405-495e-8b96-04d0ffa3065e"/>
    <rdf:value>
          Fixed up the CellML/MathML to enable the model to be solved in CMISS. The
          model is no longer a direct translation of the Noble 1998 article, but
          I believe it accurately represents the intent of the model.
        </rdf:value>
    <cmeta:modifier rdf:resource="rdf:#c1ca9719-5a5b-448c-ab7c-1da5efa5ce41"/>
  </rdf:Description>
  <rdf:Description rdf:about="#rapid_time_dependent_potassium_current_Xr1_gate">
    <cmeta:comment rdf:resource="rdf:#9cd10f81-fce1-49d5-9532-f6d4663e72ac"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#97d9bd80-107a-46a3-beee-86ce3a643138">
    <rdf:value>
            Evaluation of the ACh dependent potassium current.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="#dy_Ca_troponoin_calculation_eq">
    <cmeta:comment rdf:resource="rdf:#114cc6eb-6b35-4913-aadd-dc9116e42723"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#c2b80383-611d-48cf-96a2-5a49f681f64d">
    <rdf:value>
            This component describes the addition of
            length dependence to the Noble models. However, in this implementation
            of the model, it is unused - sarcomere length (SL) is a constant 0, 
            resulting in the removal of length dependence from the model.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="#i_K1_calculation_eq">
    <cmeta:comment rdf:resource="rdf:#95d0a48a-e028-4234-8cd5-019fdfaf4520"/>
  </rdf:Description>
  <rdf:Description rdf:about="#i_An_stretch_calculation_eq">
    <cmeta:comment rdf:resource="rdf:#66126d5a-6164-4dd8-83ea-90e9aeacd05b"/>
  </rdf:Description>
  <rdf:Description rdf:about="#beta_calculation_eq">
    <cmeta:comment rdf:resource="rdf:#eec8206d-c4e4-416b-80d3-fad37ed38c45"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#5f9ea8c7-0b5a-4647-be0b-65c545740229">
    <rdf:value>
            Stretch activated potassium current - unused and untested, but 
            included for completeness.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="">
    <dc:publisher>
        Bioengineering Institute, The University of Auckland.
      </dc:publisher>
    <cmeta:modification rdf:resource="rdf:#1fd368a3-3129-4d3b-a3b5-2e48aa072ea3"/>
    <cmeta:modification rdf:resource="rdf:#766da4ed-780a-4312-8582-6388bff32a7e"/>
    <cmeta:modification rdf:resource="rdf:#800d17b7-05b5-4b9c-8a2c-c925d29b0473"/>
    <cmeta:modification rdf:resource="rdf:#8bcca774-860a-452a-b77b-2a904ec12f4d"/>
    <cmeta:modification rdf:resource="rdf:#aa513c36-db35-4d5b-86ef-702a30f7592c"/>
    <cmeta:modification rdf:resource="rdf:#b556fcf5-318c-480e-b663-26528d76a3bf"/>
    <cmeta:modification rdf:resource="rdf:#cbda28e7-64b6-4b77-9f2c-21ee2283aa25"/>
    <cmeta:modification rdf:resource="rdf:#d32848b6-2eca-45cc-bc50-caf4195c94e2"/>
    <cmeta:modification rdf:resource="rdf:#db38540d-906b-403c-ad6a-ed863eeaea8b"/>
    <dcterms:created rdf:resource="rdf:#e09722d4-f99e-4416-acd3-14f779c42639"/>
    <dc:creator rdf:resource="rdf:#d77f81a0-2293-4d93-9390-c0c2ccf51d5b"/>
  </rdf:Description>
  <rdf:Description rdf:about="#persistent_sodium_current">
    <cmeta:comment rdf:resource="rdf:#19230fce-5835-4e36-8f52-534df5545598"/>
  </rdf:Description>
  <rdf:Description rdf:about="#f2_diff_eq">
    <cmeta:comment rdf:resource="rdf:#18f889e5-67e6-45e3-b0aa-a298e71ee95c"/>
  </rdf:Description>
  <rdf:Description rdf:about="#stretch_dependent_K_current">
    <cmeta:comment rdf:resource="rdf:#5f9ea8c7-0b5a-4647-be0b-65c545740229"/>
  </rdf:Description>
  <rdf:Description rdf:about="#i_Ca_L_Ca_calculation_eq">
    <cmeta:comment rdf:resource="rdf:#5f298cf4-b690-4df0-9c30-4044222fef23"/>
  </rdf:Description>
  <rdf:Description rdf:about="#sodium_potassium_pump">
    <cmeta:comment rdf:resource="rdf:#d6b8714a-cc64-438b-9065-0bed6400475e"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#55610f93-0041-4d52-8bab-e3de308abcf3">
    <rdf:value>
            Calculation of the current.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="#transient_outward_current">
    <cmeta:comment rdf:resource="rdf:#4e3d90d9-1ed8-4f47-90db-d88fe1e614aa"/>
  </rdf:Description>
  <rdf:Description rdf:about="#stretch_dependent_Na_current">
    <cmeta:comment rdf:resource="rdf:#0b58ad15-9be2-4e72-9f75-b9d5b1d1b1d0"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#444539dc-412c-4815-86d3-5c999e8f2486">
    <rdf:value>
            Calculation of the background calcium current.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="#i_Ca_L_K_calculation_eq">
    <cmeta:comment rdf:resource="rdf:#206027af-a204-4097-a83a-018a9f63a599"/>
  </rdf:Description>
  <rdf:Description rdf:about="#E_Ca_calculation_eq">
    <cmeta:comment rdf:resource="rdf:#5e78a3b9-abc8-41fa-9c71-bb3810a87bd1"/>
  </rdf:Description>
  <rdf:Description rdf:about="#f_stretch_calculation_eq">
    <cmeta:comment rdf:resource="rdf:#6b371cd2-85de-461a-94f6-f5b26e143f1b"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#49d37358-c2fb-45d5-b6a7-3c7d3cde6993">
    <rdf:value>
            The ODE governing the membrane potential -- given by the summation of the
            ionic currents and an applied stimulus.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#6eb2c13b-67b7-4401-8b5e-34b115246aa8">
    <rdf:value>
            The background sodium current.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#d32848b6-2eca-45cc-bc50-caf4195c94e2">
    <dcterms:modified rdf:resource="rdf:#eeaadc01-1293-4b5d-9f86-f9da93d8efb5"/>
    <rdf:value>
          Corrected several equations.
        </rdf:value>
    <cmeta:modifier rdf:resource="rdf:#70e8482b-454a-4405-9adc-1865f5c524ec"/>
  </rdf:Description>
  <rdf:Description rdf:about="#Na_i_diff_eq">
    <cmeta:comment rdf:resource="rdf:#b258a496-ef8c-4ec9-970e-e0ecbec833f2"/>
  </rdf:Description>
  <rdf:Description rdf:about="#L_type_Ca_channel_d_gate">
    <cmeta:comment rdf:resource="rdf:#f3ff6ca9-1cce-406f-8134-08d2f6ea1985"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#20ea7600-7a00-4f13-a11c-666bc807c6be">
    <rdf:value>
            Background potassium current.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#4f2ce802-cd18-45c2-aafe-ce7606667950">
    <rdf:value>
            The slow gate of the rapid component of the delayed rectifier.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="#transient_outward_current_r_gate">
    <cmeta:comment rdf:resource="rdf:#8537a9b9-1340-4295-b7dd-c2b807748695"/>
  </rdf:Description>
  <rdf:Description rdf:about="#stretch_dependent_An_current">
    <cmeta:comment rdf:resource="rdf:#405f02ef-4208-4a2f-87f4-3f682ef64dd1"/>
  </rdf:Description>
  <rdf:Description rdf:about="#fast_sodium_current">
    <cmeta:comment rdf:resource="rdf:#300f7303-a906-4ac5-a80a-5ae2ec7f2544"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#4df4358d-ca8f-415d-8abe-f38f5ca77a67">
    <bqs:Pubmed_id>9487284</bqs:Pubmed_id>
    <bqs:JournalArticle rdf:resource="rdf:#a4fc8d7c-7913-4052-8c0f-88a9bb4a0b8f"/>
  </rdf:Description>
  <rdf:Description rdf:about="#total_potassium_current">
    <cmeta:comment rdf:resource="rdf:#49c210ea-fa2f-4167-88e9-1954148e3a2d"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#db38540d-906b-403c-ad6a-ed863eeaea8b">
    <dcterms:modified rdf:resource="rdf:#4ac98732-130d-4d22-9430-46d4c7cb60e0"/>
    <rdf:value>
          Updated metadata to conform to the 16/1/02 CellML Metadata 1.0 
          Specification.
        </rdf:value>
    <cmeta:modifier rdf:resource="rdf:#03f80f40-5346-48a0-a00a-a61a751d341d"/>
  </rdf:Description>
  <rdf:Description rdf:about="#i_Ca_L_Ca_ds_calculation_eq">
    <cmeta:comment rdf:resource="rdf:#39d41c49-52a7-4cd9-909d-78b6fe6fcbdd"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#cbda28e7-64b6-4b77-9f2c-21ee2283aa25">
    <dcterms:modified rdf:resource="rdf:#f06d31ad-ec7c-49a0-8b2e-517b0ec3446f"/>
    <rdf:value>
          Corrected the i_K1 and i_K_ACh calculations.
        </rdf:value>
    <cmeta:modifier rdf:resource="rdf:#441be557-8978-4445-bae1-33899518cb25"/>
  </rdf:Description>
  <rdf:Description rdf:about="#i_NaCa_calculation_eq">
    <cmeta:comment rdf:resource="rdf:#1aaa3014-e772-43c7-a96c-d58e94931743"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#e1ffd7e5-659f-4c7c-912f-f2988cdcede1">
    <rdf:value>
            Calculation of the current.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="#i_b_Na_calculation_eq">
    <cmeta:comment rdf:resource="rdf:#19f86bf2-1bb6-42fa-81ae-312ad19b9eea"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#8103cce0-3535-4af4-8d30-034f88c6eaef">
    <rdf:value>
            Calculation of the calcium release flux from the JSR into the cytosol.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#b52a9248-2961-4de0-8968-2486025ee2ef">
    <rdf:value>
            The gating kinetics for the inactivation gate.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="#i_Kr_calculation_eq">
    <cmeta:comment rdf:resource="rdf:#0939ce95-5b99-40c5-9e3e-3310386a2aee"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#b97d20cb-6a07-4ee9-a8ea-d39e48d41db0">
    <rdf:value>
            The second inactivation gate for L-type calcium channels in the
            membrane when open into the diadic subspace.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="#OpenReleaseChannelFract_calculation_eq">
    <cmeta:comment rdf:resource="rdf:#038e9bf2-a483-4d58-9dab-22114a2d7170"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#7bbe65da-671b-4fba-a52e-e8a49e7177bf">
    <dcterms:W3CDTF>2002-01-04</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#d8520fa2-0577-4759-833f-e53eb377ba4e">
    <rdf:value>
            This component describes an inward ionic current which is the sum 
            of calcium, sodium and potassium ions through the membrane channel. 
            The channel has one activation gate (d) and two inactivation gates 
            (f and f2 (or f2ds in the diadic space)). A fraction of these channels
            (ICaLfract) open into the diadic subspace, with the remainder being
            evenly distributed throughout the membrane area.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#aa513c36-db35-4d5b-86ef-702a30f7592c">
    <dcterms:modified rdf:resource="rdf:#88eca42f-2bfb-4131-8801-7b68b13ba881"/>
    <rdf:value>
          Corrected units.
        </rdf:value>
    <cmeta:modifier rdf:resource="rdf:#401ef81e-e604-4e1d-80c5-7e39a8587126"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#c5722b0c-8e2e-4b72-8cda-47211a472c4a">
    <rdf:value>
            The rate of change of calcium concentration in the diadic subspace.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#18f889e5-67e6-45e3-b0aa-a298e71ee95c">
    <rdf:value>
            The gating kinetics for the inactivation gate.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#043b5a5c-6c69-4fd3-8731-a6dc58ceebd5">
    <rdf:value>
            Gating kinetics for the first gate.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="#reversal_potentials">
    <cmeta:comment rdf:resource="rdf:#79dee674-8cab-48e6-88cb-19b716608b61"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#19f86bf2-1bb6-42fa-81ae-312ad19b9eea">
    <rdf:value>
            Calculation of the sodium background current.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#af9c3c31-c6bb-4c9d-9c27-2e904adcec46">
    <rdf:value>
            Evaluation of the transient outward current.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="#K_i_diff_eq">
    <cmeta:comment rdf:resource="rdf:#8e5ba07b-8d9e-4238-a928-635c10690f07"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#ac3311fa-7291-47e3-a1e5-34e6a362101f">
    <rdf:value>
            Calculation of the transfer flux between the uptake and release
            calcium stores.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#64c926e1-5129-407b-b860-72039b6185c0">
    <rdf:value>
            Calculation of the reversal potential for the fast sodium channel.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#5c4f64e1-0d79-452a-a78c-d7ecdeb65247">
    <vCard:Given>Autumn</vCard:Given>
    <vCard:Family>Cuellar</vCard:Family>
    <vCard:Other>A.</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#b8f6beea-9218-44ca-95ba-c49b8d1466ba">
    <rdf:value>
            The gating kinetics for the r gate.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="#calcium_background_current">
    <cmeta:comment rdf:resource="rdf:#7268b213-457a-48b4-b94e-1bf9a18abbc0"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#206027af-a204-4097-a83a-018a9f63a599">
    <rdf:value>
            The potassium ion component of the total current through the channel into
            the bulk myoplasm.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="#potassium_background_current">
    <cmeta:comment rdf:resource="rdf:#20ea7600-7a00-4f13-a11c-666bc807c6be"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#17e0ff94-c4e4-439b-a04e-ef6a7221cb00">
    <rdf:value>
            The gating kinetics.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#29f2f1d5-1bbe-4306-a3e7-29ae02c4fdc3">
    <rdf:value>
            Calculation of the sodium calcium exchanger current for the channels
            opening into the diadic subspace.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#a55b5fe6-c072-4f9a-8448-af03be604128">
    <rdf:value>
            Calculation of the background potassium current.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#a039f8b0-0ee6-49b6-bec3-f3e1fb6ad519">
    <rdf:value>
            The gating kinetics for the fast gate of the rapid component
            of the delayed rectifier.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#5bb55bd1-b6f5-439e-952e-aafee108fc34">
    <rdf:value>
            Gating kinetics for the slow gate of the rapid component of
            the delayed rectifier.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#5716b5f6-bcb3-4282-a548-5dcfec1c30ef">
    <rdf:value>
            A convenient grouping of all the ionic concentration differential
            equations.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="#m_diff_eq">
    <cmeta:comment rdf:resource="rdf:#ff8f75b2-4f72-4b0a-89a0-9a6632cd9044"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#9cd10f81-fce1-49d5-9532-f6d4663e72ac">
    <rdf:value>
            The fast activation gate for the rapid component of the delayed 
            rectifier.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#6b371cd2-85de-461a-94f6-f5b26e143f1b">
    <rdf:value>
            Calculation of the scaling factor for the stretch activated currents
            based on the change in length from the reference state.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#3b0f2032-91d9-4c49-b685-9bbf853d829c">
    <vCard:Given>Catherine</vCard:Given>
    <vCard:Family>Lloyd</vCard:Family>
    <vCard:Other>May</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="#i_up_calculation_eq">
    <cmeta:comment rdf:resource="rdf:#a22f31e1-9f2e-4112-bd56-59f4b3c3d380"/>
  </rdf:Description>
  <rdf:Description rdf:about="#L_type_Ca_channel_f2ds_gate">
    <cmeta:comment rdf:resource="rdf:#b97d20cb-6a07-4ee9-a8ea-d39e48d41db0"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#4a4346ca-68dd-442e-871c-6e2e20041d6c">
    <rdf:value>
            Simple summation of the individual potassium currents into a 
            single current.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#038e9bf2-a483-4d58-9dab-22114a2d7170">
    <rdf:value>
            The fraction of release sites in the open state.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="#E_Na_calculation_eq">
    <cmeta:comment rdf:resource="rdf:#94795eca-34f1-40bb-b870-f6b2b5a12a71"/>
  </rdf:Description>
  <rdf:Description rdf:about="#i_tr_calculation_eq">
    <cmeta:comment rdf:resource="rdf:#ac3311fa-7291-47e3-a1e5-34e6a362101f"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#0939ce95-5b99-40c5-9e3e-3310386a2aee">
    <rdf:value>
            Calculation of the rapid component of the delayed rectifier current.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="#i_Ks_calculation_eq">
    <cmeta:comment rdf:resource="rdf:#d67b415f-6f89-4174-9e02-9ec766d08233"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#d75e9429-1273-45af-b812-5136d6476adb">
    <vCard:Given>Catherine</vCard:Given>
    <vCard:Family>Lloyd</vCard:Family>
    <vCard:Other>May</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#eeaadc01-1293-4b5d-9f86-f9da93d8efb5">
    <dcterms:W3CDTF>2002-02-26</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="#ACh_dependent_potassium_current">
    <cmeta:comment rdf:resource="rdf:#112b0393-d9a3-4e7f-8f11-05c57d90ed1f"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#25105cf3-7be9-4a63-8fcb-8e094e573b71">
    <vCard:N rdf:resource="rdf:#d75e9429-1273-45af-b812-5136d6476adb"/>
  </rdf:Description>
  <rdf:Description rdf:about="#i_K_ACh_calculation_eq">
    <cmeta:comment rdf:resource="rdf:#97d9bd80-107a-46a3-beee-86ce3a643138"/>
  </rdf:Description>
  <rdf:Description rdf:about="#h_diff_eq">
    <cmeta:comment rdf:resource="rdf:#68d7f813-7232-47a0-9082-d22a0569b922"/>
  </rdf:Description>
  <rdf:Description rdf:about="#length_dependence">
    <cmeta:comment rdf:resource="rdf:#c2b80383-611d-48cf-96a2-5a49f681f64d"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#a44d364d-b7a5-4820-b73c-53917ef4f283">
    <rdf:value>
            The gating kinetics for the slow component of the delayed 
            rectifier current.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#35fce50d-7b61-4963-be42-3bfee47a2a4b">
    <rdf:value>
            The first gate of the transient outward current.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="#x_ACh1_diff_eq">
    <cmeta:comment rdf:resource="rdf:#043b5a5c-6c69-4fd3-8731-a6dc58ceebd5"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#5e78a3b9-abc8-41fa-9c71-bb3810a87bd1">
    <rdf:value>
            Calculation of the calcium reversal potential.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#5f75dfbd-4ef6-4873-bb98-1e8bee853f77">
    <rdf:value>
            Calculation of the ATP dependent potassium current.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#70b50f28-9684-43ef-99b0-b116e262e1f8">
    <dcterms:W3CDTF>2002-07-22</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#c69bffaa-0405-495e-8b96-04d0ffa3065e">
    <dcterms:W3CDTF>2003-04-05</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="#noble_varghese_kohl_noble_1998_version05">
    <dc:title>
        The Noble 1998 Improved Guinea-Pig Ventricular Cell Model
      </dc:title>
    <cmeta:bio_entity>Ventricular Myocyte</cmeta:bio_entity>
    <cmeta:comment rdf:resource="rdf:#b5cb7c53-1f73-4b8e-9b8e-e1f6c31202db"/>
    <bqs:reference rdf:resource="rdf:#4df4358d-ca8f-415d-8abe-f38f5ca77a67"/>
    <cmeta:species>Guinea-Pig</cmeta:species>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#d505b2e7-3edf-4934-aa21-73df5a31d19a">
    <rdf:value>
            The first inactivation gate for the L-type calcium channel.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#author4Vcard">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#author4VcardN"/>
  </rdf:Description>
  <rdf:Description rdf:about="#ACh_dependent_potassium_current_xACh2_gate">
    <cmeta:comment rdf:resource="rdf:#56b39c88-f35f-4e73-ba8d-d3b7f37550fe"/>
  </rdf:Description>
  <rdf:Description rdf:about="#slow_time_dependent_potassium_current_Xs_gate">
    <cmeta:comment rdf:resource="rdf:#8ab72fb0-2050-4658-bfa1-6ce7e591ef02"/>
  </rdf:Description>
  <rdf:Description rdf:about="#rapid_time_dependent_potassium_current_Xr2_gate">
    <cmeta:comment rdf:resource="rdf:#4f2ce802-cd18-45c2-aafe-ce7606667950"/>
  </rdf:Description>
  <rdf:Description rdf:about="#calcium_uptake_from_myoplasm_to_NSR">
    <cmeta:comment rdf:resource="rdf:#f930f2d5-540c-476d-9146-c279daee3063"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#860cd03d-5ddc-42b1-b1cf-7afe53f9a639">
    <rdf:value>
            A sodium activated potassium current. Again, included for completeness
            but generally not used and untested.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#a4fc8d7c-7913-4052-8c0f-88a9bb4a0b8f">
    <dc:creator rdf:resource="rdf:#citationAuthorsSeq"/>
    <dc:title>Improved guinea-pig ventricular cell model incorporating a diadic space, IKr and IKs, and length- and tension-dependent processes</dc:title>
    <bqs:volume>14</bqs:volume>
    <bqs:first_page>123</bqs:first_page>
    <bqs:Journal rdf:resource="rdf:#547bc783-9328-482b-afe2-c2c088232162"/>
    <dcterms:issued rdf:resource="rdf:#dfcf0375-2107-4d59-a1e2-f385fd660463"/>
    <bqs:last_page>134</bqs:last_page>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#db588734-6dbc-4a22-99d7-8175bb9e26bd">
    <rdf:value>
            Control the stimulus current.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#8537a9b9-1340-4295-b7dd-c2b807748695">
    <rdf:value>
            The second of the two gates for the transient outward current.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#490e2392-5a76-46fc-9207-151d4f6c4f7c">
    <rdf:value>
            The sodium calcium exchanger. A fraction of the Na-Ca exchangers
            (INaCaFract) are assumed to emmpty into the diadic subspace, while
            the remainder of the channels open into the bulk cytosol.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="#rapid_time_dependent_potassium_current">
    <cmeta:comment rdf:resource="rdf:#73dc6f46-b807-4b30-b31f-fa158dd540a7"/>
  </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:#face12b6-70ca-4c84-a329-0bcd0588580e">
    <rdf:value>
            The gating kinetics for the activation gate.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#4e3d90d9-1ed8-4f47-90db-d88fe1e614aa">
    <rdf:value>
            The doubly gated transient outward current.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#401ef81e-e604-4e1d-80c5-7e39a8587126">
    <vCard:N rdf:resource="rdf:#3b0f2032-91d9-4c49-b685-9bbf853d829c"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#94795eca-34f1-40bb-b870-f6b2b5a12a71">
    <rdf:value>
            Calculation of the sodium reversal potential from the current intra- and
            extracellular concentration of sodium ions.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#b5cb7c53-1f73-4b8e-9b8e-e1f6c31202db">
    <dc:creator rdf:resource="rdf:#21075bcd-777d-4bb3-8029-c20a45ff1333"/>
    <rdf:value>This is the CellML description of Noble's 1998 improved guinea-pig ventricular cell model. It incorporates a diadic space, rapid and slow potassium currents and length- and tension-dependent processes.</rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#7eb621bb-57af-4b0e-a58c-0cf72984b47e">
    <rdf:value>
            Calculation of the sodium activated potassium current.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="#r_ss_calculation_eq">
    <cmeta:comment rdf:resource="rdf:#d540d5d1-8a44-4bcd-9732-4c8a2783987f"/>
  </rdf:Description>
  <rdf:Description rdf:about="#f2ds_diff_eq">
    <cmeta:comment rdf:resource="rdf:#b52a9248-2961-4de0-8968-2486025ee2ef"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#d67b415f-6f89-4174-9e02-9ec766d08233">
    <rdf:value>
            Calculation of the slow component current of the delayed rectifier.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="#i_b_Ca_calculation_eq">
    <cmeta:comment rdf:resource="rdf:#444539dc-412c-4815-86d3-5c999e8f2486"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#c8ac0dc5-e94f-40ce-8e3f-32d330938520">
    <rdf:value>
            The inactivation gate for the fast sodium current.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#7268b213-457a-48b4-b94e-1bf9a18abbc0">
    <rdf:value>
            The background calcium current.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#66126d5a-6164-4dd8-83ea-90e9aeacd05b">
    <rdf:value>
            Calculation of the current.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="#ionic_concentrations">
    <cmeta:comment rdf:resource="rdf:#5716b5f6-bcb3-4282-a548-5dcfec1c30ef"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#95d0a48a-e028-4234-8cd5-019fdfaf4520">
    <rdf:value>
            Calculation of the time independent potassium current.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#8e5ba07b-8d9e-4238-a928-635c10690f07">
    <rdf:value>
            The rate of change of intracellular potassium concentration.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="#Xr1_diff_eq">
    <cmeta:comment rdf:resource="rdf:#a039f8b0-0ee6-49b6-bec3-f3e1fb6ad519"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#3aa64ead-ec6b-47bd-abab-3eb860e6537a">
    <rdf:value>
            The sodium ion component of the total current through the channel into
            the bulk myoplasm.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#3277481d-39d4-4513-bd18-a000d7cb3b93">
    <rdf:value>
            Calculation of the potassium reversal potential for the slow component of
            the delayed rectifier.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="#i_rel_calculation_eq">
    <cmeta:comment rdf:resource="rdf:#8103cce0-3535-4af4-8d30-034f88c6eaef"/>
  </rdf:Description>
  <rdf:Description rdf:about="#environment">
    <cmeta:comment rdf:resource="rdf:#2a3927b4-9421-4fdc-946e-02bd99cc40b1"/>
  </rdf:Description>
  <rdf:Description rdf:about="#fast_sodium_current_h_gate">
    <cmeta:comment rdf:resource="rdf:#c8ac0dc5-e94f-40ce-8e3f-32d330938520"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#b4461d55-bc4d-4f44-84f9-e66b11db72fa">
    <dcterms:W3CDTF>2004-12-23</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="#sodium_activated_potassium_current">
    <cmeta:comment rdf:resource="rdf:#860cd03d-5ddc-42b1-b1cf-7afe53f9a639"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#39d41c49-52a7-4cd9-909d-78b6fe6fcbdd">
    <rdf:value>
            The calcium ion component of the total current through the channel into
            the diadic subspace.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="#i_NaCa_ds_calculation_eq">
    <cmeta:comment rdf:resource="rdf:#29f2f1d5-1bbe-4306-a3e7-29ae02c4fdc3"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#205fcb00-ae84-4785-b66d-5c413be536d5">
    <rdf:value>
            The length dependent leakage of calcium from the junctional
            SR into the cytolsol.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#ed46d659-b810-4035-9e3c-56f12d1efaf2">
    <rdf:value>
            The activation gate for the fast sodium current.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#0ed4ead2-b6bf-455d-8923-4c2d011b2c5d">
    <vCard:Given>David</vCard:Given>
    <vCard:Family>Nickerson</vCard:Family>
    <vCard:Other>P</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#5aea6f31-73e0-4fdb-93ce-f7367b9de317">
    <rdf:value>
            The rate of change of calcium concentration in the junctional SR
            release store.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#f3a3f9ed-3ea5-4651-8607-8130af0fec4d">
    <rdf:value>
            The opening rate for the activation gate.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#6ed05548-2c09-4964-bc90-8cc357331048">
    <rdf:value>
            The rate of change of calcium bound to calmodulin.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="#alpha_calculation_eq">
    <cmeta:comment rdf:resource="rdf:#c93d31a6-2e13-489b-a824-df76bc67a294"/>
  </rdf:Description>
  <rdf:Description rdf:about="#IStim_for_cmiss_eq">
    <cmeta:comment rdf:resource="rdf:#db588734-6dbc-4a22-99d7-8175bb9e26bd"/>
  </rdf:Description>
  <rdf:Description rdf:about="#r_diff_eq">
    <cmeta:comment rdf:resource="rdf:#b8f6beea-9218-44ca-95ba-c49b8d1466ba"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#author2VcardN">
    <vCard:Given>Anthony</vCard:Given>
    <vCard:Family>Varghese</vCard:Family>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#f75a5fb6-2029-43ca-b9cb-5ffa72ca92f1">
    <rdf:value>
            This component is the `root' node of our model.
            It defines the action potential variable `V'.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="#ds_decay_flux">
    <cmeta:comment rdf:resource="rdf:#fc766bdd-33c1-4b71-a94c-0e3d0644fc1a"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#1fd368a3-3129-4d3b-a3b5-2e48aa072ea3">
    <dcterms:modified rdf:resource="rdf:#7bbe65da-671b-4fba-a52e-e8a49e7177bf"/>
    <rdf:value>
          Altered some of the connections.
        </rdf:value>
    <cmeta:modifier rdf:resource="rdf:#25105cf3-7be9-4a63-8fcb-8e094e573b71"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#3a2f42cd-54ef-4b5b-aa2f-e40eeebb11b1">
    <rdf:value>
            The second inactivation gate for L-type calcium channels in the
            membrane outside the diadic subspace.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="#L_type_Ca_channel">
    <cmeta:comment rdf:resource="rdf:#d8520fa2-0577-4759-833f-e53eb377ba4e"/>
  </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:#48185889-f99e-48aa-8ed3-35c4534da3cd">
    <vCard:Given>Catherine</vCard:Given>
    <vCard:Family>Lloyd</vCard:Family>
    <vCard:Other>May</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#114cc6eb-6b35-4913-aadd-dc9116e42723">
    <rdf:value>
            The rate of change of calcium bound to troponin.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#1c949c9e-04aa-43b1-86d4-cf249ad72903">
    <vCard:N rdf:resource="rdf:#48185889-f99e-48aa-8ed3-35c4534da3cd"/>
  </rdf:Description>
  <rdf:Description rdf:about="#i_K_ATP_calculation_eq">
    <cmeta:comment rdf:resource="rdf:#5f75dfbd-4ef6-4873-bb98-1e8bee853f77"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#0e3bf76f-c0ac-4ab5-af3c-2723f6f03d5c">
    <rdf:value>
            Diffusion down the concentration gradient to transfer calcium from
            the uptake stores in the network SR to the release stores in the
            junctional SR.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#112b0393-d9a3-4e7f-8f11-05c57d90ed1f">
    <rdf:value>
            Acetylcholine (ACh) dependent potassium current. Once more, 
            this current has been included for completeness with respect 
            to the original publication but is usually left out of simulations.
            It is useful for simulating specific cellular conditions, but
            this formulation has not been tested.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#b258a496-ef8c-4ec9-970e-e0ecbec833f2">
    <rdf:value>
            The rate of change of intracellular sodium concentration.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="#E_mh_calculation_eq">
    <cmeta:comment rdf:resource="rdf:#64c926e1-5129-407b-b860-72039b6185c0"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#156d8060-d5b8-4c48-9cf0-a91b2aaaa209">
    <rdf:value>
            The gating kinetics for the first inactivation gate.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#faa723c7-893a-4ea8-87b9-94f7e9ca80e3">
    <rdf:value>
            The rate of change of calcium concentration in the network SR
            uptake store.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#f930f2d5-540c-476d-9146-c279daee3063">
    <rdf:value>
            The sarcoplasmic reticulum calcium pump which transports calcium
            from the bulk myoplasm into the network SR.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="#f_diff_eq">
    <cmeta:comment rdf:resource="rdf:#156d8060-d5b8-4c48-9cf0-a91b2aaaa209"/>
  </rdf:Description>
  <rdf:Description rdf:about="#i_Ca_stretch_calculation_eq">
    <cmeta:comment rdf:resource="rdf:#e8400ef0-576b-4889-8a84-e2309cb0b0de"/>
  </rdf:Description>
  <rdf:Description rdf:about="#Xs_diff_eq">
    <cmeta:comment rdf:resource="rdf:#a44d364d-b7a5-4820-b73c-53917ef4f283"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#17cadbf6-886e-4f7f-ac5a-ddfccd302fc1">
    <rdf:value>
            Calculation of the current.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#eec8206d-c4e4-416b-80d3-fad37ed38c45">
    <rdf:value>
            Calculation of the closing rate of the activation gate for the fast
            sodium current.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#ff8ce98e-3b9c-4a7a-80b1-513350b95ce7">
    <vCard:Given>Catherine</vCard:Given>
    <vCard:Family>Lloyd</vCard:Family>
    <vCard:Other>May</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#c1ca9719-5a5b-448c-ab7c-1da5efa5ce41">
    <vCard:N rdf:resource="rdf:#0ed4ead2-b6bf-455d-8923-4c2d011b2c5d"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#de79743f-0a6f-4020-a40c-fce8d4917d7e">
    <rdf:value>
            A non-specific stretch activated current. Not used and untested, but
            included for completeness.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="#i_K_stretch_calculation_eq">
    <cmeta:comment rdf:resource="rdf:#17cadbf6-886e-4f7f-ac5a-ddfccd302fc1"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#ce8d92dc-1085-4076-b1ff-7de935b7be77">
    <rdf:value>
            Calculation of the potassium reversal potential from the current intra- and
            extracellular concentration of potassium ions.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#45cc9575-7a50-4072-bffc-7e66fac2c375">
    <rdf:value>
            Calculation of the persistent sodium current.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#1aaa3014-e772-43c7-a96c-d58e94931743">
    <rdf:value>
            Calculation of the sodium calcium exchanger current for the channels
            exposed to the bulk myoplasm.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="#membrane_voltage_diff_eq">
    <cmeta:comment rdf:resource="rdf:#49d37358-c2fb-45d5-b6a7-3c7d3cde6993"/>
  </rdf:Description>
  <rdf:Description rdf:about="#i_Ca_L_K_ds_calculation_eq">
    <cmeta:comment rdf:resource="rdf:#b972c4d4-f35a-468f-8a33-6da921be1fc7"/>
  </rdf:Description>
  <rdf:Description rdf:about="#stretch_dependent_Ca_current">
    <cmeta:comment rdf:resource="rdf:#3ef03e92-b7ba-47c5-862b-c5ba167002aa"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#a34a521e-1b23-47dc-9692-238f091bcbfd">
    <rdf:value>
            Description of the time independent inward potassium current.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#9f42bedb-4cb8-40ae-9901-b680e66a8315">
    <vCard:Given>David</vCard:Given>
    <vCard:Family>Nickerson</vCard:Family>
    <vCard:Other>P</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="#i_Na_stretch_calculation_eq">
    <cmeta:comment rdf:resource="rdf:#55610f93-0041-4d52-8bab-e3de308abcf3"/>
  </rdf:Description>
  <rdf:Description rdf:about="#PrecursorFract_calculation_eq">
    <cmeta:comment rdf:resource="rdf:#72764e6a-d7c0-4305-8750-1c212a670cf8"/>
  </rdf:Description>
  <rdf:Description rdf:about="#beta_calculation_eq2">
    <cmeta:comment rdf:resource="rdf:#1d88ae98-a642-4872-8f26-7023028082e6"/>
  </rdf:Description>
  <rdf:Description rdf:about="#beta_calculation_eq3">
    <cmeta:comment rdf:resource="rdf:#71b11fbf-9aa0-47d4-a145-5a9de06eec9c"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#ebd58747-9f34-4a1e-acca-6705214db1ed">
    <rdf:value>
            The first of two gates for the ACh dependent potassium current.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#405f02ef-4208-4a2f-87f4-3f682ef64dd1">
    <rdf:value>
            A anion specific stretch ativated current - unused and untested, but 
            included for completeness.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#19230fce-5835-4e36-8f52-534df5545598">
    <rdf:value>
            The persistent sodium current.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#8ab72fb0-2050-4658-bfa1-6ce7e591ef02">
    <rdf:value>
            The gate for the slow component of the delayed rectifier current.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#5448d293-03cf-4519-804f-43965bce91d7">
    <rdf:value>
            Calculation of the sodium potassium pump current.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="#alpha_calculation_eq3">
    <cmeta:comment rdf:resource="rdf:#5d46daf9-51e6-4b05-97c5-686d2db03e8e"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#70e8482b-454a-4405-9adc-1865f5c524ec">
    <vCard:N rdf:resource="rdf:#ff8ce98e-3b9c-4a7a-80b1-513350b95ce7"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#800d17b7-05b5-4b9c-8a2c-c925d29b0473">
    <dcterms:modified rdf:resource="rdf:#b07f9f6d-b894-46c8-8217-6bdccc865d0a"/>
    <rdf:value>
          Added some initial values from Penny Noble's documentation.  Removed 
          the blocked h_gate component from the fast sodium current as it 
          belongs to a separate model which considers the effect of drugs.
        </rdf:value>
    <cmeta:modifier rdf:resource="rdf:#1c949c9e-04aa-43b1-86d4-cf249ad72903"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#21075bcd-777d-4bb3-8029-c20a45ff1333">
    <vCard:FN>Catherine Lloyd</vCard:FN>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#5f298cf4-b690-4df0-9c30-4044222fef23">
    <rdf:value>
            The calcium ion component of the total current through the channel into
            the bulk myoplasm.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#a0ec4572-a007-4ba3-9af2-bf9984fdb07f">
    <rdf:value>
            The rate of change of intracellular calcium concentration.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="#slow_time_dependent_potassium_current">
    <cmeta:comment rdf:resource="rdf:#65a1807b-0f20-448d-9663-4da7ffefcca8"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#a22f31e1-9f2e-4112-bd56-59f4b3c3d380">
    <rdf:value>
            Calculation of the calcium uptake flux.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#79dee674-8cab-48e6-88cb-19b716608b61">
    <rdf:value>
            A component that conveniently keeps the reversal potential calculations
            together.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="#i_K_calculation_eq">
    <cmeta:comment rdf:resource="rdf:#4a4346ca-68dd-442e-871c-6e2e20041d6c"/>
  </rdf:Description>
  <rdf:Description rdf:about="#i_NaK_calculation_eq">
    <cmeta:comment rdf:resource="rdf:#5448d293-03cf-4519-804f-43965bce91d7"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#5d46daf9-51e6-4b05-97c5-686d2db03e8e">
    <rdf:value>
            The opening rate for the inactivation gate.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#766da4ed-780a-4312-8582-6388bff32a7e">
    <dcterms:modified rdf:resource="rdf:#b4461d55-bc4d-4f44-84f9-e66b11db72fa"/>
    <rdf:value>
          Encapsulated the model under an interface component for use when
          importing the model into a specific experiment or simulation.
        </rdf:value>
    <cmeta:modifier rdf:resource="rdf:#04ee1311-52db-48a6-be5a-39b8b796aad9"/>
  </rdf:Description>
  <rdf:Description rdf:about="#L_type_Ca_channel_f_gate">
    <cmeta:comment rdf:resource="rdf:#d505b2e7-3edf-4934-aa21-73df5a31d19a"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#441be557-8978-4445-bae1-33899518cb25">
    <vCard:N rdf:resource="rdf:#8e23cfe8-7f55-4499-b523-d56719c74daa"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#e8400ef0-576b-4889-8a84-e2309cb0b0de">
    <rdf:value>
            Calculation of the current.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#04dad8f9-a202-4320-9011-638bb5a3de91">
    <rdf:value>
            An ATP dependent potassium current. Included here for completeness, but 
            this current is only used when modelling ischaemia - and this version
            of the model has not been tested for this.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="#i_Ca_L_Na_calculation_eq">
    <cmeta:comment rdf:resource="rdf:#3aa64ead-ec6b-47bd-abab-3eb860e6537a"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#0b58ad15-9be2-4e72-9f75-b9d5b1d1b1d0">
    <rdf:value>
            A stretch activated sodium current - unused and untested, but
            included for completeness.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="#stretch_dependent_Ns_current">
    <cmeta:comment rdf:resource="rdf:#de79743f-0a6f-4020-a40c-fce8d4917d7e"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#300f7303-a906-4ac5-a80a-5ae2ec7f2544">
    <rdf:value>
            The fast sodium current - the driving force of the upstroke of the
            cardiac action potential.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="#s_diff_eq">
    <cmeta:comment rdf:resource="rdf:#17e0ff94-c4e4-439b-a04e-ef6a7221cb00"/>
  </rdf:Description>
  <rdf:Description rdf:about="#L_type_Ca_channel_f2_gate">
    <cmeta:comment rdf:resource="rdf:#3a2f42cd-54ef-4b5b-aa2f-e40eeebb11b1"/>
  </rdf:Description>
  <rdf:Description rdf:about="#calcium_transfer_from_NSR_to_JSR">
    <cmeta:comment rdf:resource="rdf:#0e3bf76f-c0ac-4ab5-af3c-2723f6f03d5c"/>
  </rdf:Description>
  <rdf:Description rdf:about="#E_Ks_calculation_eq">
    <cmeta:comment rdf:resource="rdf:#3277481d-39d4-4513-bd18-a000d7cb3b93"/>
  </rdf:Description>
  <rdf:Description rdf:about="#d_diff_eq">
    <cmeta:comment rdf:resource="rdf:#face12b6-70ca-4c84-a329-0bcd0588580e"/>
  </rdf:Description>
  <rdf:Description rdf:about="#dy_Ca_calmod_calculation_eq">
    <cmeta:comment rdf:resource="rdf:#6ed05548-2c09-4964-bc90-8cc357331048"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#49c210ea-fa2f-4167-88e9-1954148e3a2d">
    <rdf:value>
            A component to group the potassium currents.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="#Ca_rel_diff_eq">
    <cmeta:comment rdf:resource="rdf:#5aea6f31-73e0-4fdb-93ce-f7367b9de317"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#c93d31a6-2e13-489b-a824-df76bc67a294">
    <rdf:value>
            Calculation of the opening rate of the activation gate for the fast
            sodium current.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#1d88ae98-a642-4872-8f26-7023028082e6">
    <rdf:value>
            The closing rate for the activation gate.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#03f80f40-5346-48a0-a00a-a61a751d341d">
    <vCard:N rdf:resource="rdf:#5c4f64e1-0d79-452a-a78c-d7ecdeb65247"/>
  </rdf:Description>
  <rdf:Description rdf:about="#x_ACh2_diff_eq">
    <cmeta:comment rdf:resource="rdf:#4cf08dd3-41fb-46e2-92cb-2d1246bfa448"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#fc766bdd-33c1-4b71-a94c-0e3d0644fc1a">
    <rdf:value>
            Description of the diffusion of calcium out of the diadic
            subspace into the bulk cytosol.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="#fast_sodium_current_m_gate">
    <cmeta:comment rdf:resource="rdf:#ed46d659-b810-4035-9e3c-56f12d1efaf2"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#4cf08dd3-41fb-46e2-92cb-2d1246bfa448">
    <rdf:value>
            Gating kinetics for the second gate.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#547bc783-9328-482b-afe2-c2c088232162">
    <dc:title>Canadian Journal of Cardiology</dc:title>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#b07f9f6d-b894-46c8-8217-6bdccc865d0a">
    <dcterms:W3CDTF>2002-05-06</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#71b11fbf-9aa0-47d4-a145-5a9de06eec9c">
    <rdf:value>
            The closing rate for the inactivation gate.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="#Ca_ds_diff_eq">
    <cmeta:comment rdf:resource="rdf:#c5722b0c-8e2e-4b72-8cda-47211a472c4a"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#ff8f75b2-4f72-4b0a-89a0-9a6632cd9044">
    <rdf:value>
            The gating kinetics for the activation gate of the fast sodium current.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#68d7f813-7232-47a0-9082-d22a0569b922">
    <rdf:value>
            The gating kinetics for the inactivation gate of the fast sodium 
            current.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#4ac98732-130d-4d22-9430-46d4c7cb60e0">
    <dcterms:W3CDTF>2002-01-21</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="#i_Na_calculation_eq">
    <cmeta:comment rdf:resource="rdf:#84e294b7-8c90-4dc6-8fd1-a119778973f1"/>
  </rdf:Description>
  <rdf:Description rdf:about="#time_independent_potassium_current">
    <cmeta:comment rdf:resource="rdf:#a34a521e-1b23-47dc-9692-238f091bcbfd"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#f3ff6ca9-1cce-406f-8134-08d2f6ea1985">
    <rdf:value>
            The activation gate of the L-type calcium channel.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="#Xr2_diff_eq">
    <cmeta:comment rdf:resource="rdf:#5bb55bd1-b6f5-439e-952e-aafee108fc34"/>
  </rdf:Description>
  <rdf:Description rdf:about="#i_p_Na_calculation_eq">
    <cmeta:comment rdf:resource="rdf:#45cc9575-7a50-4072-bffc-7e66fac2c375"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#8e23cfe8-7f55-4499-b523-d56719c74daa">
    <vCard:Given>Catherine</vCard:Given>
    <vCard:Family>Lloyd</vCard:Family>
    <vCard:Other>May</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="#i_Ns_stretch_calculation_eq">
    <cmeta:comment rdf:resource="rdf:#e1ffd7e5-659f-4c7c-912f-f2988cdcede1"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#735dee69-ddee-475c-a8e0-2a57e36428c3">
    <vCard:Given>Catherine</vCard:Given>
    <vCard:Family>Lloyd</vCard:Family>
    <vCard:Other>May</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#84e294b7-8c90-4dc6-8fd1-a119778973f1">
    <rdf:value>
            Calculation of the fast sodium current.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#author4VcardN">
    <vCard:Given>Penny</vCard:Given>
    <vCard:Family>Noble</vCard:Family>
  </rdf:Description>
  <rdf:Description rdf:about="#i_b_K_calculation_eq">
    <cmeta:comment rdf:resource="rdf:#a55b5fe6-c072-4f9a-8448-af03be604128"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#d540d5d1-8a44-4bcd-9732-4c8a2783987f">
    <rdf:value>
            Calculation of the steady-state component of the gating kinetics
            for the r gate.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#author1VcardN">
    <vCard:Given>Denis</vCard:Given>
    <vCard:Family>Noble</vCard:Family>
  </rdf:Description>
  <rdf:Description rdf:about="#sodium_background_current">
    <cmeta:comment rdf:resource="rdf:#6eb2c13b-67b7-4401-8b5e-34b115246aa8"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#2a3927b4-9421-4fdc-946e-02bd99cc40b1">
    <rdf:value>
            This component is used to declare variables that are used by
            all or most of the other components, in this case just `time'.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="#Ca_up_diff_eq">
    <cmeta:comment rdf:resource="rdf:#faa723c7-893a-4ea8-87b9-94f7e9ca80e3"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#6dce197b-096b-4ad9-9eb9-0e6f7d5e2600">
    <vCard:Given>Catherine</vCard:Given>
    <vCard:Family>Lloyd</vCard:Family>
    <vCard:Other>May</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#d77f81a0-2293-4d93-9390-c0c2ccf51d5b">
    <vCard:ORG rdf:resource="rdf:#d55021bd-737c-46f2-ad6f-709faf7bde38"/>
    <vCard:EMAIL rdf:resource="rdf:#b06662dd-5000-4cf2-91b4-e8e57ca5c7a9"/>
    <vCard:N rdf:resource="rdf:#6dce197b-096b-4ad9-9eb9-0e6f7d5e2600"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#3ef03e92-b7ba-47c5-862b-c5ba167002aa">
    <rdf:value>
            A stretch activated calcium current. Unused and untested in this 
            implementation, but included for completeness.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#e09722d4-f99e-4416-acd3-14f779c42639">
    <dcterms:W3CDTF>2001-12-18</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#b972c4d4-f35a-468f-8a33-6da921be1fc7">
    <rdf:value>
            The potassium ion component of the total current through the channel into
            the diadic subspace.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="#alpha_calculation_eq2">
    <cmeta:comment rdf:resource="rdf:#f3a3f9ed-3ea5-4651-8607-8130af0fec4d"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#d55021bd-737c-46f2-ad6f-709faf7bde38">
    <vCard:Orgname>The University of Auckland</vCard:Orgname>
    <vCard:Orgunit>The Bioengineering Research Group</vCard:Orgunit>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#73dc6f46-b807-4b30-b31f-fa158dd540a7">
    <rdf:value>
            The rapid component of the delayed rectifier current.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#dfcf0375-2107-4d59-a1e2-f385fd660463">
    <dcterms:W3CDTF>1998-01-01</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="#E_K_calculation_eq">
    <cmeta:comment rdf:resource="rdf:#ce8d92dc-1085-4076-b1ff-7de935b7be77"/>
  </rdf:Description>
  <rdf:Description rdf:about="#i_Ca_L_Na_ds_calculation_eq">
    <cmeta:comment rdf:resource="rdf:#67a795b5-d23c-4c1d-8107-5efc975ec984"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#b556fcf5-318c-480e-b663-26528d76a3bf">
    <dcterms:modified rdf:resource="rdf:#70b50f28-9684-43ef-99b0-b116e262e1f8"/>
    <rdf:value>
          Added more metadata.
        </rdf:value>
    <cmeta:modifier rdf:resource="rdf:#2b9d3310-dbad-4680-82f8-5286f059a431"/>
  </rdf:Description>
  <rdf:Description rdf:about="#sodium_calcium_pump">
    <cmeta:comment rdf:resource="rdf:#490e2392-5a76-46fc-9207-151d4f6c4f7c"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#56b39c88-f35f-4e73-ba8d-d3b7f37550fe">
    <rdf:value>
            The second of the two gates for the ACh dependent potassium current.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#65a1807b-0f20-448d-9663-4da7ffefcca8">
    <rdf:value>
            The slow component of the delayed rectifier current.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#d6b8714a-cc64-438b-9065-0bed6400475e">
    <rdf:value>
            The sodium potassium pump.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="#calcium_leak_from_SR">
    <cmeta:comment rdf:resource="rdf:#205fcb00-ae84-4785-b66d-5c413be536d5"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#72764e6a-d7c0-4305-8750-1c212a670cf8">
    <rdf:value>
            Fraction of release channels in the precursor state.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="#ATP_dependent_potassium_current">
    <cmeta:comment rdf:resource="rdf:#04dad8f9-a202-4320-9011-638bb5a3de91"/>
  </rdf:Description>
  <rdf:Description rdf:about="#RegulatoryBindingSite_calculation_eq">
    <cmeta:comment rdf:resource="rdf:#66b34f79-38b2-4bd8-8152-811b026043ff"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#1a60361b-64f9-478a-a87e-657e818183ad">
    <rdf:value>
            Calculation of the flux of calcium out of the diadic subspace.
          </rdf:value>
  </rdf:Description>
</rdf:RDF>
</model>