Location: Matsuoka, Sarai, Kuratomi, Ono, Noma, 2003 @ ed31cab29fdd / matsuoka_sarai_kuratomi_ono_noma_2003.cellml

Author:
pmr2.import <nobody@models.cellml.org>
Date:
2006-09-04 03:11:58+12:00
Desc:
committing version01 of matsuoka_sarai_kuratomi_ono_noma_2003
Permanent Source URI:
https://models.cellml.org/workspace/matsuoka_sarai_kuratomi_ono_noma_2003/rawfile/ed31cab29fdd9bfc12b2ebae4e4dc76a7f473072/matsuoka_sarai_kuratomi_ono_noma_2003.cellml

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

CREATED :  14th March 2003

LAST MODIFIED : 20th April 2005

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

DESCRIPTION :  This file contains a CellML description of Matsuoka et al's 2003 mathematical model of the roles of individual ion currents in sinoatrial node cells. 

CHANGES:  
21/07/2003 - CML - Made changes to several equations on author recommendation.
13/08/2003 - CML - Made more changes to several equations on author recommendation. 
20/04/2005 - PJV - Made MathML id's unique 
--><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="matsuoka_sarai_kuratomi_ono_noma_2003_version02" name="matsuoka_sarai_kuratomi_ono_noma_2003_version02">
<documentation xmlns="http://cellml.org/tmp-documentation">
<article>
  <articleinfo>
  <title>Modelling the Roles of Individual Current Systems in Ventricular Cells and in the SA Node</title>
  <author>
    <firstname>Catherine</firstname>
          <surname>Lloyd</surname>
    <affiliation>
      <shortaffil>Bioengineering Institute, University of Auckland</shortaffil>
    </affiliation>
  </author>
</articleinfo>
  <section id="sec_status">
    <title>Model Status</title>
    <para>
            This is the original unchecked version of the model imported from the previous
            CellML model repository, 24-Jan-2006.
          </para>
  </section>
  <sect1 id="sec_structure">
<title>Model Structure</title>

<para>
Different types of ion channel, with their distinct gating and conductance properties, play a unique role in generating membrane excitation.  Although the behaviours of single ion channels have been extensively studied experimentally, ionic current interactions make it difficult to deduce their individual roles in membrane excitability.  In order to better understand the roles of the individual ion currents, Satoshi Matsuoka <emphasis>et al.</emphasis> have developed two mathematical models of membrane excitation for the sinoatrial node pacemaker cells and for ventricular cells.  The two models are based on essentially the same set of equations, with cell specificity expressed as different current densities, as well as variable cellular components (see <xref linkend="fig_ventricular_cell_diagram"/> and <xref linkend="fig_SA_node_cell_diagram"/> below for a visual comparison of the two different models). 
</para>

<para>
The Matsuoka <emphasis>et al.</emphasis> mathematical models are a development of previously published cardiac cellular models.  The group's general approach to modelling was based on the methods used in the <ulink url="${HTML_EXMPL_DFN_MODEL}">Difrancesco-Noble, 1985 Model</ulink> and the <ulink url="${HTML_EXMPL_LR_II_MODEL}">Luo-Rudy Ventricular Model II (dynamic), 1994</ulink>.  In addition, the compound <emphasis>Kyoto Model</emphasis> of Matsuoka <emphasis>et al.</emphasis> incorporates the Negroni and Lascano (1996) model of cardiac muscle contraction, allowing simulation of sarcomere shortening, and also, the gating model of Shirokov <emphasis>et al.</emphasis> is used for both the fast Na<superscript>+</superscript> and L-type Ca<superscript>2+</superscript> channels.   
</para>

<para>
New experimental data are included, such as new kinetics of the inward rectifier K<superscript>+</superscript> channel, the delayed rectifier K<superscript>+</superscript> channel, and the sustained inward current (for more details on the reaction kinetics, see <xref linkend="fig_reaction_kinetics"/> below).  Model simulations showed that ion conductance by the L-type Ca<superscript>2+</superscript> channel (I<subscript>CaL</subscript>) dominates the membrane conductance during the action potential.  Repolarisation occurs as the following currents are sequentially activated: I<subscript>Ks</subscript>; I<subscript>Kr</subscript>; and I<subscript>K1</subscript>.  Agreement between model simulation results and experimental data was taken as supporting evidence for the validity of the model. 
</para>

<para>
The complete original paper references are cited below:
</para>

<para>
<ulink url="http://www.jstage.jst.go.jp/article/jjphysiol/53/2/53_105/_article">Role of Individual Ionic Current Systems in Ventricular Cells Hypothesized by a Model Study</ulink>, Matsuoka S, Sarai N, Kuratomi S, Ono K, and Noma A, 2003, <ulink url="http://wwwsoc.nii.ac.jp/psj/jjp/jjphome.html">
            <emphasis>The Japanese Journal of Physiology</emphasis>
          </ulink>, 53, 105-123.  (A  PDF version of the article is available to subscribers on the <emphasis>Japanese Journal of Physiology</emphasis> website.)  <ulink url="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?cmd=Retrieve&amp;db=PubMed&amp;list_uids=12877767&amp;dopt=Abstract">PubMed ID: 12877767</ulink>

<ulink url="http://www.jstage.jst.go.jp/article/jjphysiol/53/2/53_125/_article">Role of Individual Ionic Current Systems in the SA Node Hypothesized by a Model Study</ulink>, Sarai N, Matsuoka S, Kuratomi S, Ono K, and Noma A, 2003, <ulink url="http://wwwsoc.nii.ac.jp/psj/jjp/jjphome.html">
            <emphasis>The Japanese Journal of Physiology</emphasis>
          </ulink>, 53, 125-134.  (A PDF version of the article is available to subscribers on the <emphasis>Japanese Journal of Physiology</emphasis> website.)  <ulink url="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?cmd=Retrieve&amp;db=PubMed&amp;list_uids=12877768&amp;dopt=Abstract">PubMed ID: 12877768</ulink>
</para>

<informalfigure float="0" id="fig_ventricular_cell_diagram">
<mediaobject>
  <imageobject>
    <objectinfo>
      <title>ventricular cell diagram</title>
    </objectinfo>
    <imagedata fileref="matsuoka_2003.png"/>
  </imageobject>
</mediaobject>
<caption>A schematic diagram describing the ionic components of the Matsuoka <emphasis>et al.</emphasis> 2003 ventricular cell mathematical model.  Ions are exchanged between the intracellular and extracellular environments, and between the cytosol and the sarcoplasmic reticulum (SR) via channels and pumps.  Unique to the ventricular cells are the transient outward current, I<subscript>to</subscript>, and the slow component of the delayed rectifier K<superscript>+</superscript> current, I<subscript>Ks</subscript>.</caption>
</informalfigure>

<informalfigure float="0" id="fig_SA_node_cell_diagram">
<mediaobject>
  <imageobject>
    <objectinfo>
      <title>SA node cell diagram</title>
    </objectinfo>
    <imagedata fileref="matsuoka_san_2003.png"/>
  </imageobject>
</mediaobject>
<caption>A schematic diagram describing the ionic components of the Matsuoka <emphasis>et al.</emphasis> 2003 SA node cell mathematical model.  Unique to the SA node cells are the acetylecholine-activated K<superscript>+</superscript> current, I<subscript>KACh</subscript>, the hyperpolarisation-activated cation current, I<subscript>ha</subscript>, and the sustained inward current I<subscript>st</subscript>.</caption>
</informalfigure>

<informalfigure float="0" id="fig_reaction_kinetics">
<mediaobject>
  <imageobject>
    <objectinfo>
      <title>reaction kinetics diagram</title>
    </objectinfo>
    <imagedata fileref="matsuoka_3_2003.png"/>
  </imageobject>
</mediaobject>
<caption>A schematic diagram describing the reaction kinetics of the ion channel gates in the Matsuoka <emphasis>et al.</emphasis> 2003 mathematical models.  A) The kinetics of the sodium channel voltage-dependent gate.  B) The kinetics of the Na-Ca exchange pump.  C) The kinetics of the Na-K pump.  D) The kinetics of the inward rectifier K<superscript>+</superscript> current.  E)  The kinetics of the Ca<superscript>2+</superscript>-dependent gate of the L-type Ca<superscript>2+</superscript> channel.  F) The kinetics of the ultra-slow gates of several ion channels.  G) The kinetics of the SR Ca<superscript>2+</superscript> pump.  H) The kinetics of the hyperpolarisation-activated cation current.  I) The kinetics of the RyR channel in the SR.</caption>
</informalfigure>

</sect1>
</article>
</documentation>
  
  
  <!--
    Below, we define some additional units for association with variables and
    constants within the model.
  -->
  
  <units name="millivolt">
    <unit units="volt" prefix="milli"/>
  </units>
  
  <units name="millimole">
    <unit units="mole" prefix="milli"/>
  </units>

  <units name="millimolar">
    <unit units="mole" prefix="milli"/>
    <unit units="litre" exponent="-1"/>
  </units>
  
  <units name="micromolar">
    <unit units="mole" prefix="micro"/>
    <unit units="litre" exponent="-1"/>
  </units>
  
  <units name="millimolar_per_picoA">
    <unit units="millimolar"/>
    <unit units="picoA" exponent="-1"/>
  </units>
  
  <units name="per_micromolar2">
    <unit units="micromolar" exponent="-2"/>
  </units>
  
  <units name="second_per_micromolar2">
    <unit units="second"/>
    <unit units="micromolar" exponent="-2"/>
  </units>
  
  <units name="second_per_micrometre2">
    <unit units="second"/>
    <unit units="micrometre" exponent="-2"/>
  </units>
  
  <units name="mN_per_mm2">
    <unit units="newton" prefix="milli"/>
    <unit units="metre" prefix="milli" exponent="-2"/>
  </units>
  
  <units name="mN_per_mm2_micrometre_micromolar">
    <unit units="newton" prefix="milli"/>
    <unit units="metre" prefix="milli" exponent="-2"/>
    <unit units="micrometre" exponent="-1"/>
    <unit units="micromolar" exponent="-1"/>
  </units>
  
  <units name="per_micrometre2">
    <unit units="micrometre" exponent="-2"/>
  </units>
  
  <units name="millisecond">
    <unit units="second" prefix="milli"/>
  </units>
  
  <units name="flux">
    <unit units="micromolar"/>
    <unit units="second" exponent="-1"/>
  </units>
  
  <units name="first_order_rate_constant">
    <unit units="millisecond" exponent="-1"/>
  </units>
  
  <units name="second_order_rate_constant">
    <unit units="millimolar" exponent="-1"/>
    <unit units="millisecond" exponent="-1"/>
  </units>
  
  <units name="picoF">
    <unit units="farad" prefix="pico"/>
  </units>
  
  <units name="picoF_first_order_rate_constant">
    <unit units="picoF"/>
    <unit units="millisecond" exponent="-1"/>
  </units>
  
  <units name="picoA">
    <unit units="ampere" prefix="pico"/>
  </units>
  
  <units name="micrometre">
    <unit units="metre" prefix="micro"/>
  </units>
  
  <units name="micrometre3">
    <unit units="metre" prefix="micro" exponent="3"/>
  </units>
  
  <units name="picoA_per_millivolt">
    <unit units="ampere" prefix="pico"/>
    <unit units="millivolt" exponent="-1"/>
  </units>
  
  <units name="picoA_per_millimolar">
    <unit units="ampere" prefix="pico"/>
    <unit units="millimolar" exponent="-1"/>
  </units>
  
  <units name="coulomb_millivolt_per_kelvin_millimole">
    <unit units="coulomb"/>
    <unit units="millivolt"/>
    <unit units="kelvin" exponent="-1"/>
    <unit units="millimole" exponent="-1"/>
  </units>
  
  <units name="coulomb_per_millimole">
    <unit units="coulomb"/>
    <unit units="mole" prefix="milli" exponent="-1"/>
  </units>
  
  <component name="environment">
    <variable units="millisecond" public_interface="out" name="time"/>
  </component>
  
  <component name="membrane">
    <variable units="millivolt" public_interface="out" name="Vm" initial_value="-85.87"/>
    <variable units="coulomb_millivolt_per_kelvin_millimole" public_interface="out" name="R" initial_value="8.3143"/>
    <variable units="kelvin" public_interface="out" name="T" initial_value="310.0"/>
    <variable units="coulomb_per_millimole" public_interface="out" name="F" initial_value="96.4867"/>
    <variable units="picoF" public_interface="out" name="Cm" initial_value="132.0"/>
    <variable units="picoA" public_interface="out" name="i_ext" initial_value="0.0"/>
    
    <variable units="picoA" name="i_tot"/>
    <variable units="picoA" name="i_1"/>
    
    <variable units="millisecond" public_interface="in" name="time"/>
    <variable units="picoA" public_interface="in" name="i_Na"/>
    <variable units="picoA" public_interface="in" name="i_Ca_L"/>
    <variable units="picoA" public_interface="in" name="i_Ca_T"/>
    <variable units="picoA" public_interface="in" name="i_K1"/>
    <variable units="picoA" public_interface="in" name="i_Kr"/>
    <variable units="picoA" public_interface="in" name="i_Ks"/>
    <variable units="picoA" public_interface="in" name="i_to"/>
    <variable units="picoA" public_interface="in" name="i_NaK"/>
    <variable units="picoA" public_interface="in" name="i_NaCa"/>
    <variable units="picoA" public_interface="in" name="i_bNSC"/>
    <variable units="picoA" public_interface="in" name="i_Cab"/>
    <variable units="picoA" public_interface="in" name="i_Kpl"/>
    <variable units="picoA" public_interface="in" name="i_lCa"/>
    <variable units="picoA" public_interface="in" name="i_KATP"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="membrane_voltage_diff_eq">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> Vm </ci>
        </apply>
        <apply>
          <divide/>
          <apply>
            <minus/>
            <apply>
              <plus/>
              <ci> i_tot </ci>
              <ci> i_ext </ci>
            </apply>
          </apply>
          <ci> Cm </ci>
        </apply>
      </apply>
      
      <apply id="i_tot_calculation">
        <eq/>
        <ci> i_tot </ci>
        <apply>
          <plus/>
          <ci> i_Na </ci>
          <ci> i_Ca_L </ci>
          <ci> i_Ca_T </ci>
          <ci> i_K1 </ci>
          <ci> i_Kr </ci>
          <ci> i_Ks </ci>
          <ci> i_to </ci>
          <ci> i_1 </ci>
          <ci> i_NaK </ci>
          <ci> i_NaCa </ci>
        </apply>
      </apply>
      
      <apply id="i_1_calculation">
        <eq/>
        <ci> i_1 </ci>
        <apply>
          <plus/>
          <ci> i_bNSC </ci>
          <ci> i_Cab </ci>
          <ci> i_Kpl </ci>
          <ci> i_lCa </ci>
          <ci> i_KATP </ci>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="external_ion_concentrations">
    <variable units="millimolar" public_interface="out" name="Nao" initial_value="140.0"/>
    <variable units="millimolar" public_interface="out" name="Cao" initial_value="1.8"/>
    <variable units="millimolar" public_interface="out" name="Ko" initial_value="5.4"/>
  </component>
  
  <component name="internal_ion_concentrations">
    <variable units="millimolar" public_interface="out" name="Nai" initial_value="4.772"/>
    <variable units="millimolar" public_interface="out" name="Cai" initial_value="3.943E-4"/>
    <variable units="millimolar" public_interface="out" name="Ki" initial_value="1.425E2"/>
    <variable units="micrometre3" public_interface="out" name="Vi" initial_value="8000.0"/>
    
    <variable units="picoA" name="i_net_Na"/>
    <variable units="picoA" name="i_net_K"/>
    <variable units="picoA" name="i_net_Ca"/>  
    
    <variable units="millisecond" public_interface="in" name="time"/>
    <variable units="coulomb_per_millimole" public_interface="in" name="F"/>
    <variable units="picoA" public_interface="in" name="i_Na_Na"/>
    <variable units="picoA" public_interface="in" name="i_CaL_Na"/>
    <variable units="picoA" public_interface="in" name="i_bNSC_Na"/>
    <variable units="picoA" public_interface="in" name="i_lCa_Na"/>
    <variable units="picoA" public_interface="in" name="i_NaK"/>
    <variable units="picoA" public_interface="in" name="i_K1"/>
    <variable units="picoA" public_interface="in" name="i_Kr"/>
    <variable units="picoA" public_interface="in" name="i_to_K"/>
    <variable units="picoA" public_interface="in" name="i_to_Na"/>
    <variable units="picoA" public_interface="in" name="i_Ks_K"/>
    <variable units="picoA" public_interface="in" name="i_Ks_Na"/>
    <variable units="picoA" public_interface="in" name="i_KATP"/>
    <variable units="picoA" public_interface="in" name="i_Na_K"/>
    <variable units="picoA" public_interface="in" name="i_CaL_K"/>
    <variable units="picoA" public_interface="in" name="i_bNSC_K"/>
    <variable units="picoA" public_interface="in" name="i_lCa_K"/>
    <variable units="picoA" public_interface="in" name="i_Kpl"/>
    <variable units="picoA" public_interface="in" name="i_CaL_Ca"/>
    <variable units="picoA" public_interface="in" name="i_Ca_T"/>
    <variable units="picoA" public_interface="in" name="i_NaCa"/>
    <variable units="picoA" public_interface="in" name="i_bNSC"/>
    <variable units="picoA" public_interface="in" name="i_Cab"/>
    <variable units="picoA" public_interface="in" name="i_ext"/>
    <variable units="picoA" public_interface="in" name="i_RyR"/>
    <variable units="picoA" public_interface="in" name="i_SR_U"/>
    <variable units="picoA" public_interface="in" name="i_SR_L"/>
    <variable units="dimensionless" public_interface="in" name="TT"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="Nai_diff_eq">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> Nai </ci>
        </apply>
        <apply>
          <divide/>
          <apply>
            <minus/>
            <ci> i_net_Na </ci>
          </apply>
          <apply>
            <times/>
            <ci> F </ci>
            <ci> Vi </ci>
          </apply>
        </apply>
      </apply>
      
      <apply id="Ki_diff_eq">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> Ki </ci>
        </apply>
        <apply>
          <divide/>
          <apply>
            <minus/>
            <apply>
              <plus/>
              <ci> i_net_K </ci>
              <ci> i_ext </ci>
            </apply>
          </apply>
          <apply>
            <times/>
            <ci> F </ci>
            <ci> Vi </ci>
          </apply>
        </apply>
      </apply>
      
      <apply id="Cai_diff_eq">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> Cai </ci>
        </apply>
        <apply>
          <plus/>
          <apply>
            <divide/>
            <apply>
              <minus/>
              <apply>
                <plus/>
                <ci> i_net_Ca </ci>
                <ci> i_SR_U </ci>
                <apply>
                  <minus/>
                  <apply>
                    <plus/>
                    <ci> i_RyR </ci>
                    <ci> i_SR_L </ci>
                  </apply>
                </apply>
              </apply>
            </apply>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> 2.0 </cn>
              <ci> F </ci>
              <ci> Vi </ci>
            </apply>
          </apply>
          <ci> TT </ci>
        </apply>
      </apply>
      
      <apply id="i_net_Na_calculation">
        <eq/>
        <ci> i_net_Na </ci>
        <apply>
          <plus/>
          <ci> i_Na_Na </ci>
          <ci> i_Ks_Na </ci>
          <ci> i_to_Na </ci>
          <ci> i_CaL_Na </ci>
          <ci> i_bNSC_Na </ci>
          <ci> i_lCa_Na </ci>
          <apply>
            <times/>
            <cn cellml:units="dimensionless"> 3.0 </cn>
            <ci> i_NaK </ci>
          </apply>
          <apply>
            <times/>
            <cn cellml:units="dimensionless"> 3.0 </cn>
            <ci> i_NaCa </ci>
          </apply>
        </apply>
      </apply>
      
      <apply id="i_net_K_calculation">
        <eq/>
        <ci> i_net_K </ci>
        <apply>
          <plus/>
          <ci> i_K1 </ci>
          <ci> i_Kr </ci>
          <ci> i_to_K </ci>
          <ci> i_Ks_K </ci>
          <ci> i_Na_K </ci>
          <ci> i_CaL_K </ci>
          <ci> i_bNSC_K </ci>
          <ci> i_lCa_K </ci>
          <ci> i_Kpl </ci>
          <apply>
            <times/>
            <cn cellml:units="dimensionless"> -2.0 </cn>
            <ci> i_NaK </ci>
          </apply>
        </apply>
      </apply>
      
      <apply id="i_net_Ca_calculation">
        <eq/>
        <ci> i_net_Ca </ci>
        <apply>
          <plus/>
          <ci> i_CaL_Ca </ci>
          <ci> i_Ca_T </ci>
          <ci> i_Cab </ci>
          <apply>
            <times/>
            <cn cellml:units="dimensionless"> -2.0 </cn>
            <ci> i_NaCa </ci>
          </apply>
        </apply>
      </apply>
    </math>
  </component> 
  
  <component name="constant_field_equations">
    <variable units="millimolar" public_interface="out" name="CF_Na"/>
    <variable units="millimolar" public_interface="out" name="CF_Ca"/>
    <variable units="millimolar" public_interface="out" name="CF_K"/>
    
    <variable units="millimolar" public_interface="in" name="Nai"/>
    <variable units="millimolar" public_interface="in" name="Nao"/>
    <variable units="millimolar" public_interface="in" name="Cai"/>
    <variable units="millimolar" public_interface="in" name="Cao"/>
    <variable units="millimolar" public_interface="in" name="Ki"/>
    <variable units="millimolar" public_interface="in" name="Ko"/>
    <variable units="coulomb_millivolt_per_kelvin_millimole" public_interface="in" name="R"/>
    <variable units="kelvin" public_interface="in" name="T"/>
    <variable units="coulomb_per_millimole" public_interface="in" name="F"/>
    <variable units="millivolt" public_interface="in" name="Vm"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="CF_Na_calculation">
        <eq/>
        <ci> CF_Na </ci>
        <apply>
          <times/>
          <apply>
            <divide/>
            <apply>
              <times/>
              <ci> F </ci>
              <ci> Vm </ci>
            </apply>
            <apply>
              <times/>
              <ci> R </ci>
              <ci> T </ci>
            </apply>
          </apply>
          <apply>
            <divide/>
            <apply>
              <minus/>
              <ci> Nai </ci>
              <apply>
                <times/>
                <ci> Nao </ci>
                <apply>
                  <exp/>
                  <apply>
                    <divide/>
                    <apply>
                      <times/>
                      <cn cellml:units="dimensionless"> -1.0 </cn>
                      <ci> F </ci>
                      <ci> Vm </ci>
                    </apply>
                    <apply>
                      <times/>
                      <ci> R </ci>
                      <ci> T </ci>
                    </apply>
                  </apply>
                </apply>
              </apply>
            </apply>
            <apply>
              <minus/>
              <cn cellml:units="dimensionless"> 1.0 </cn>
              <apply>
                <exp/>
               <apply>
                  <divide/>
                 <apply>
                    <times/>
                   <cn cellml:units="dimensionless"> -1.0 </cn>
                   <ci> F </ci>
                   <ci> Vm </ci>
                 </apply>
                 <apply>
                    <times/>
                   <ci> R </ci>
                   <ci> T </ci>
                 </apply>
               </apply>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply id="CF_Ca_calculation">
        <eq/>
        <ci> CF_Ca </ci>
        <apply>
          <times/>
          <apply>
            <divide/>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> 2.0 </cn>
              <ci> F </ci>
              <ci> Vm </ci>
            </apply>
            <apply>
              <times/>
              <ci> R </ci>
              <ci> T </ci>
            </apply>
          </apply>
          <apply>
            <divide/>
            <apply>
              <minus/>
              <ci> Cai </ci>
              <apply>
                <times/>
                <ci> Cao </ci>
                <apply>
                  <exp/>
                  <apply>
                    <divide/>
                    <apply>
                      <times/>
                      <cn cellml:units="dimensionless"> -2.0 </cn>
                      <ci> F </ci>
                      <ci> Vm </ci>
                    </apply>
                    <apply>
                      <times/>
                      <ci> R </ci>
                      <ci> T </ci>
                    </apply>
                  </apply>
                </apply>
              </apply>
            </apply>
            <apply>
              <minus/>
              <cn cellml:units="dimensionless"> 1.0 </cn>
              <apply>
                <exp/>
               <apply>
                  <divide/>
                 <apply>
                    <times/>
                   <cn cellml:units="dimensionless"> -2.0 </cn>
                   <ci> F </ci>
                   <ci> Vm </ci>
                 </apply>
                 <apply>
                    <times/>
                   <ci> R </ci>
                   <ci> T </ci>
                 </apply>
               </apply>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply id="CF_K_calculation">
        <eq/>
        <ci> CF_K </ci>
        <apply>
          <times/>
          <apply>
            <divide/>
            <apply>
              <times/>
              <ci> F </ci>
              <ci> Vm </ci>
            </apply>
            <apply>
              <times/>
              <ci> R </ci>
              <ci> T </ci>
            </apply>
          </apply>
          <apply>
            <divide/>
            <apply>
              <minus/>
              <ci> Ki </ci>
              <apply>
                <times/>
                <ci> Ko </ci>
                <apply>
                  <exp/>
                  <apply>
                    <divide/>
                    <apply>
                      <times/>
                      <cn cellml:units="dimensionless"> -1.0 </cn>
                      <ci> F </ci>
                      <ci> Vm </ci>
                    </apply>
                    <apply>
                      <times/>
                      <ci> R </ci>
                      <ci> T </ci>
                    </apply>
                  </apply>
                </apply>
              </apply>
            </apply>
            <apply>
              <minus/>
              <cn cellml:units="dimensionless"> 1.0 </cn>
              <apply>
                <exp/>
               <apply>
                  <divide/>
                 <apply>
                    <times/>
                   <cn cellml:units="dimensionless"> -1.0 </cn>
                   <ci> F </ci>
                   <ci> Vm </ci>
                 </apply>
                 <apply>
                    <times/>
                   <ci> R </ci>
                   <ci> T </ci>
                 </apply>
               </apply>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="ATP_production">
    <variable units="millimolar" public_interface="out" name="ATPi" initial_value="4.663"/>
    
    <variable units="millimolar" name="ADP"/>
    <variable units="millimolar" name="ATP_contraction"/>
    <variable units="micromolar" name="troponin" initial_value="70.0"/>
    
    <variable units="millisecond" public_interface="in" name="time"/>
    <variable units="millivolt" public_interface="in" name="Vm"/>
    <variable units="dimensionless" public_interface="in" name="TCa_"/>
    <variable units="coulomb_per_millimole" public_interface="in" name="F"/>
    <variable units="micrometre3" public_interface="in" name="Vi"/>
    <variable units="picoA" public_interface="in" name="i_NaK"/>
    <variable units="picoA" public_interface="in" name="i_SR_U"/>
  
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="ATP_contraction_diff_eq">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> ATP_contraction </ci>
        </apply>    
        <apply>
          <times/>
          <cn cellml:units="dimensionless"> 0.4 </cn>
          <ci> TCa_ </ci>
          <ci> troponin </ci>
        </apply>
      </apply>
      
      <apply id="ATP_synthesis_eq">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> ATPi </ci>
        </apply>
        <apply>
          <minus/>
          <apply>
            <times/>
            <cn cellml:units="dimensionless"> 0.003 </cn>
            <ci> ADP </ci>
          </apply>
          <apply>
            <plus/>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> 0.4 </cn>
              <ci> TCa_ </ci>
              <ci> troponin </ci>
            </apply>
            <apply>
              <divide/>
              <ci> i_NaK </ci>
              <apply>
                <times/>
                <ci> F </ci>
                <ci> Vi </ci>
              </apply>
            </apply>
            <apply>
              <divide/>
              <ci> i_SR_U </ci>
              <apply>
                <times/>
                <cn cellml:units="dimensionless"> 4.0 </cn>
                <ci> F </ci>
                <ci> Vi </ci>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply id="ADP_calculation">
        <eq/>
        <ci> ADP </ci>
        <apply>
          <minus/>
          <cn cellml:units="millimolar"> 5.0 </cn>
          <ci> ATPi </ci>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="sodium_current">
    <variable units="picoA" public_interface="out" name="i_Na"/>
    <variable units="picoA" public_interface="out" name="i_Na_Na"/>
    <variable units="picoA" public_interface="out" name="i_Na_K"/>
    
    <variable units="picoA_per_millimolar" name="P_Na" initial_value="2860.0"/>
    
    <variable units="millisecond" public_interface="in" private_interface="out" name="time"/>
    <variable units="millivolt" public_interface="in" private_interface="out" name="Vm"/>
    <variable units="millimolar" public_interface="in" name="CF_Na"/>
    <variable units="millimolar" public_interface="in" name="CF_K"/>
    
    <variable units="dimensionless" private_interface="in" name="p_AP_Na"/>
    <variable units="dimensionless" private_interface="in" name="y"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">    
      <apply id="i_Na_calculation">
        <eq/>
        <ci> i_Na </ci>
        <apply>
          <plus/>
          <ci> i_Na_Na </ci>
          <ci> i_Na_K </ci>
        </apply>
      </apply>
      
      <apply id="i_Na_Na_calculation">
        <eq/>
        <ci> i_Na_Na </ci>
        <apply>
          <times/>
          <ci> P_Na </ci>
          <ci> CF_Na </ci>
          <ci> p_AP_Na </ci>
          <ci> y </ci>
        </apply>
      </apply>
      
      <apply id="i_Na_K_calculation">
        <eq/>
        <ci> i_Na_K </ci>
        <apply>
          <times/>
          <cn cellml:units="dimensionless"> 0.1 </cn>
          <ci> P_Na </ci>
          <ci> CF_K </ci>
          <ci> p_AP_Na </ci>
          <ci> y </ci>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="sodium_current_voltage_dependent_gate">
    <variable units="dimensionless" public_interface="out" name="p_AP_Na"/> 
    
    <variable units="dimensionless" name="p_RI_Na" initial_value="1.691E-5"/>
    <variable units="dimensionless" name="p_RP_Na" initial_value="3.515E-1"/>
    <variable units="dimensionless" name="p_AI_Na" initial_value="4.068E-1"/>
    <variable units="first_order_rate_constant" name="k_RP_AP"/>
    <variable units="first_order_rate_constant" name="k_AP_RP"/>
    <variable units="first_order_rate_constant" name="k_RI_AI"/>
    <variable units="first_order_rate_constant" name="k_AI_RI"/>
    <variable units="first_order_rate_constant" name="k_AP_AI"/>
    <variable units="first_order_rate_constant" name="k_AI_AP" initial_value="0.0000875"/>
    <variable units="first_order_rate_constant" name="k_RP_RI"/>
    <variable units="first_order_rate_constant" name="k_RI_RP"/>
    
    <variable units="millisecond" public_interface="in" name="time"/>
    <variable units="millivolt" public_interface="in" name="Vm"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="p_AP_Na_diff_eq">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> p_AP_Na </ci>
        </apply>    
        <apply>
          <minus/>
          <apply>
            <plus/>
            <apply>
              <times/>
              <ci> p_RP_Na </ci>
              <ci> k_RP_AP </ci>
            </apply>
            <apply>
              <times/>
              <ci> p_AI_Na </ci>
              <ci> k_AI_AP </ci>
            </apply>
          </apply>
          <apply>
            <plus/>
            <apply>
              <times/>
              <ci> p_AP_Na </ci>
              <ci> k_AP_RP </ci>
            </apply>
            <apply>
              <times/>
              <ci> p_AP_Na </ci>
              <ci> k_AP_AI </ci>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply id="p_RP_Na_diff_eq">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> p_RP_Na </ci>
        </apply>    
        <apply>
          <minus/>
          <apply>
            <plus/>
            <apply>
              <times/>
              <ci> p_AP_Na </ci>
              <ci> k_AP_RP </ci>
            </apply>
            <apply>
              <times/>
              <ci> p_RI_Na </ci>
              <ci> k_RI_RP </ci>
            </apply>
          </apply>
          <apply>
            <plus/>
            <apply>
              <times/>
              <ci> p_RP_Na </ci>
              <ci> k_RP_RI </ci>
            </apply>
            <apply>
              <times/>
              <ci> p_RP_Na </ci>
              <ci> k_RP_AP </ci>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply id="p_AI_Na_diff_eq">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> p_AI_Na </ci>
        </apply>    
        <apply>
          <minus/>
          <apply>
            <plus/>
            <apply>
              <times/>
              <ci> p_RI_Na </ci>
              <ci> k_RI_AI </ci>
            </apply>
            <apply>
              <times/>
              <ci> p_AP_Na </ci>
              <ci> k_AP_AI </ci>
            </apply>
          </apply>
          <apply>
            <plus/>
            <apply>
              <times/>
              <ci> p_AI_Na </ci>
              <ci> k_AI_RI </ci>
            </apply>
            <apply>
              <times/>
              <ci> p_AI_Na </ci>
              <ci> k_AI_AP </ci>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply id="p_RI_Na_diff_eq">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> p_RI_Na </ci>
        </apply>    
        <apply>
          <minus/>
          <apply>
            <plus/>
            <apply>
              <times/>
              <ci> p_RP_Na </ci>
              <ci> k_RP_RI </ci>
            </apply>
            <apply>
              <times/>
              <ci> p_AI_Na </ci>
              <ci> k_AI_RI </ci>
            </apply>
          </apply>
          <apply>
            <plus/>
            <apply>
              <times/>
              <ci> p_RI_Na </ci>
              <ci> k_RI_RP </ci>
            </apply>
            <apply>
              <times/>
              <ci> p_RI_Na </ci>
              <ci> k_RI_AI </ci>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply id="sodium_current_voltage_dependent_gate_k_RP_AP_calculation">
        <eq/>
        <ci> k_RP_AP </ci>
        <apply>
          <divide/>
          <cn cellml:units="first_order_rate_constant"> 1.0 </cn>
          <apply>
            <plus/>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> 0.1027 </cn>
              <apply>
                <exp/>
                <apply>
                  <divide/>
                  <apply>
                    <minus/>
                    <ci> Vm </ci>
                  </apply>
                  <cn cellml:units="millivolt"> 8.0 </cn>
                </apply>
              </apply>
            </apply>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> 0.25 </cn>
              <apply>
                <exp/>
                <apply>
                  <divide/>
                  <apply>
                    <minus/>
                    <ci> Vm </ci>
                  </apply>
                  <cn cellml:units="millivolt"> 50.0 </cn>
                </apply>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply id="sodium_current_voltage_dependent_gate_k_AP_RP_calculation">
        <eq/>
        <ci> k_AP_RP </ci>
        <apply>
          <divide/>
          <cn cellml:units="first_order_rate_constant"> 1.0 </cn>
          <apply>
            <plus/>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> 26.0 </cn>
              <apply>
                <exp/>
                <apply>
                  <divide/>
                  <ci> Vm </ci>
                  <cn cellml:units="millivolt"> 17.0 </cn>
                </apply>
              </apply>
            </apply>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> 0.02 </cn>
              <apply>
                <exp/>
                <apply>
                  <divide/>
                  <ci> Vm </ci>
                  <cn cellml:units="millivolt"> 800.0 </cn>
                </apply>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply id="sodium_current_voltage_dependent_gate_k_RI_AI_calculation">
        <eq/>
        <ci> k_RI_AI </ci>
        <apply>
          <divide/>
          <cn cellml:units="first_order_rate_constant"> 1.0 </cn>
          <apply>
            <plus/>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> 0.0001027 </cn>
              <apply>
                <exp/>
                <apply>
                  <divide/>
                  <apply>
                    <minus/>
                    <ci> Vm </ci>
                  </apply>
                  <cn cellml:units="millivolt"> 8.0 </cn>
                </apply>
              </apply>
            </apply>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> 5.0 </cn>
              <apply>
                <exp/>
                <apply>
                  <divide/>
                  <apply>
                    <minus/>
                    <ci> Vm </ci>
                  </apply>
                  <cn cellml:units="millivolt"> 400.0 </cn>
                </apply>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply id="sodium_current_voltage_dependent_gate_k_AI_RI_calculation">
        <eq/>
        <ci> k_AI_RI </ci>
        <apply>
          <divide/>
          <cn cellml:units="first_order_rate_constant"> 1.0 </cn>
          <apply>
            <plus/>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> 1300.0 </cn>
              <apply>
                <exp/>
                <apply>
                  <divide/>
                  <ci> Vm </ci>
                  <cn cellml:units="millivolt"> 20.0 </cn>
                </apply>
              </apply>
            </apply>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> 0.04 </cn>
              <apply>
                <exp/>
                <apply>
                  <divide/>
                  <ci> Vm </ci>
                  <cn cellml:units="millivolt"> 800.0 </cn>
                </apply>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply id="k_AP_AI_calculation">
        <eq/>
        <ci> k_AP_AI </ci>
        <apply>
          <divide/>
          <cn cellml:units="first_order_rate_constant"> 1.0 </cn>
          <apply>
            <times/>
            <cn cellml:units="dimensionless"> 0.8 </cn>
            <apply>
              <exp/>
              <apply>
                <divide/>
                <apply>
                  <minus/>
                  <ci> Vm </ci>
                </apply>
                <cn cellml:units="millivolt"> 400.0 </cn>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply id="sodium_current_voltage_dependent_gate_k_RP_RI_calculation">
        <eq/>
        <ci> k_RP_RI </ci>
        <apply>
          <divide/>
          <cn cellml:units="first_order_rate_constant"> 0.01 </cn>
          <apply>
            <plus/>
            <cn cellml:units="first_order_rate_constant"> 1.0 </cn>
            <apply>
              <divide/>
              <apply>
                <times/>
                <ci> k_AI_AP </ci>
                <ci> k_AP_RP </ci>
                <ci> k_RI_AI </ci>
              </apply>
              <apply>
                <times/>
                <ci> k_AP_AI </ci>
                <ci> k_RP_AP </ci>
                <ci> k_AI_RI </ci>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply id="sodium_current_voltage_dependent_gate_k_RI_RP_calculation">
        <eq/>
        <ci> k_RI_RP </ci>
        <apply>
          <minus/>
          <cn cellml:units="first_order_rate_constant"> 0.01 </cn>
          <ci> k_RP_RI </ci>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="sodium_current_ultra_slow_gate">
    <variable units="dimensionless" public_interface="out" name="y" initial_value="5.810E-1"/> 

    <variable units="first_order_rate_constant" name="alpha_y"/>
    <variable units="first_order_rate_constant" name="beta_y"/>
    
    <variable units="millivolt" public_interface="in" name="Vm"/>
    <variable units="millisecond" public_interface="in" name="time"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="sodium_current_ultra_slow_gate_alpha_y_calculation">
        <eq/>
        <ci> alpha_y </ci>
        <apply>
          <divide/>
          <cn cellml:units="first_order_rate_constant"> 1.0 </cn>
          <apply>
            <plus/>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> 9000000000.0 </cn>
              <apply>
                <exp/>
                <apply>
                  <divide/>
                  <ci> Vm </ci>
                  <cn cellml:units="millivolt"> 5.0 </cn>
                </apply>
              </apply>
            </apply>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> 8000.0 </cn>
              <apply>
                <exp/>
                <apply>
                  <divide/>
                  <ci> Vm </ci>
                  <cn cellml:units="millivolt"> 100.0 </cn>
                </apply>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
    
      <apply id="sodium_current_ultra_slow_gate_beta_y_calculation">
        <eq/>
        <ci> beta_y </ci>
        <apply>
          <divide/>
          <cn cellml:units="first_order_rate_constant"> 1.0 </cn>
          <apply>
            <plus/>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> 0.014 </cn>
              <apply>
                <exp/>
                <apply>
                  <divide/>
                  <apply>
                    <minus/>
                    <ci> Vm </ci>
                  </apply>
                  <cn cellml:units="millivolt"> 5.0 </cn>
                </apply>
              </apply>
            </apply>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> 4000.0 </cn>
              <apply>
                <exp/>
                <apply>
                  <divide/>
                  <apply>
                    <minus/>
                    <ci> Vm </ci>
                  </apply>
                  <cn cellml:units="millivolt"> 100.0 </cn>
                </apply>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
    
      <apply id="sodium_current_ultra_slow_gate_dy_dt">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> y </ci>
        </apply>
        <apply>
          <minus/>
          <apply>
            <times/>
            <ci> alpha_y </ci>
            <apply>
              <minus/>
              <cn cellml:units="dimensionless"> 1.0 </cn>
              <ci> y </ci>
            </apply>
          </apply>
          <apply>
            <times/>
            <ci> beta_y </ci>
            <ci> y </ci>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="L_type_Ca_channel">
    <variable units="picoA" public_interface="out" name="i_Ca_L"/>
    <variable units="picoA" public_interface="out" name="i_CaL_Na"/>
    <variable units="picoA" public_interface="out" name="i_CaL_Ca"/>
    <variable units="picoA" public_interface="out" name="i_CaL_K"/>
    <variable units="dimensionless" public_interface="out" name="p_open_CaL"/>
    <variable units="picoA" public_interface="out" private_interface="in" name="iCaL"/>
    
    <variable units="picoA_per_millimolar" name="P_CaL" initial_value="8712.0"/>
     
    <variable units="millisecond" public_interface="in" private_interface="out" name="time"/>
    <variable units="millivolt" public_interface="in" private_interface="out" name="Vm"/>
    <variable units="millimolar" public_interface="in" name="CF_Na"/>
    <variable units="millimolar" public_interface="in" private_interface="out" name="CF_Ca"/>
    <variable units="millimolar" public_interface="in" private_interface="out" name="Cai"/>
    <variable units="millimolar" public_interface="in" name="CF_K"/>
    <variable units="millimolar" public_interface="in" name="ATPi"/>
    
    <variable units="dimensionless" private_interface="in" name="p_AP_CaL"/>
    <variable units="dimensionless" private_interface="in" name="p_U"/>
    <variable units="dimensionless" private_interface="in" name="p_UCa"/>
    <variable units="dimensionless" private_interface="in" name="y"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">    
      <apply id="i_Ca_L_calculation">
        <eq/>
        <ci> i_Ca_L </ci>
        <apply>
          <plus/>
          <ci> i_CaL_Na </ci>
          <ci> i_CaL_Ca </ci>
          <ci> i_CaL_K </ci>
        </apply>
      </apply>
      
      <apply id="i_CaL_Ca_calculation">
        <eq/>
        <ci> i_CaL_Ca </ci>
        <apply>
          <times/>
          <ci> P_CaL </ci>
          <ci> CF_Ca </ci>
          <ci> p_open_CaL </ci>
        </apply>
      </apply>
      
      <apply id="i_CaL_Na_calculation">
        <eq/>
        <ci> i_CaL_Na </ci>
        <apply>
          <times/>
          <cn cellml:units="dimensionless"> 0.0000185 </cn>
          <ci> P_CaL </ci>
          <ci> CF_Na </ci>
          <ci> p_open_CaL </ci>
        </apply>
      </apply>
      
      <apply id="i_CaL_K_calculation">
        <eq/>
        <ci> i_CaL_K </ci>
        <apply>
          <times/>
          <cn cellml:units="dimensionless"> 0.000365 </cn>
          <ci> P_CaL </ci>
          <ci> CF_K </ci>
          <ci> p_open_CaL </ci>
        </apply>
      </apply>
      
      <apply id="p_open_CaL_calculation">
        <eq/>
        <ci> p_open_CaL </ci>
        <apply>
          <times/>
          <ci> p_AP_CaL </ci>
          <apply>
            <plus/>
            <ci> p_U </ci>
            <ci> p_UCa </ci>
          </apply>
          <apply>
            <divide/>
            <ci> y </ci>
            <apply>
              <plus/>
              <cn cellml:units="dimensionless"> 1.0 </cn>
              <apply>
                <power/>
                <apply>
                  <divide/>
                  <cn cellml:units="millimolar"> 1.4 </cn>
                  <ci> ATPi </ci>
                </apply>
                <cn cellml:units="dimensionless"> 3.0 </cn>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="L_type_Ca_channel_voltage_dependent_gate">
    <variable units="dimensionless" public_interface="out" name="p_AP_CaL"/> 
    
    <variable units="dimensionless" name="p_RI_CaL" initial_value="1.621E-6"/>
    <variable units="dimensionless" name="p_RP_CaL" initial_value="9.966E-1"/>
    <variable units="dimensionless" name="p_AI_CaL" initial_value="9.784E-4"/>
    <variable units="first_order_rate_constant" name="k_RP_AP"/>
    <variable units="first_order_rate_constant" name="k_AP_RP"/>
    <variable units="first_order_rate_constant" name="k_RI_AI"/>
    <variable units="first_order_rate_constant" name="k_AI_RI"/>
    <variable units="first_order_rate_constant" name="k_AP_AI" initial_value="0.004"/>
    <variable units="first_order_rate_constant" name="k_AI_AP" initial_value="0.001"/>
    <variable units="first_order_rate_constant" name="k_RP_RI"/>
    <variable units="first_order_rate_constant" name="k_RI_RP"/>
    
    <variable units="millisecond" public_interface="in" name="time"/>
    <variable units="millivolt" public_interface="in" name="Vm"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="p_AP_CaL_diff_eq">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> p_AP_CaL </ci>
        </apply>    
        <apply>
          <minus/>
          <apply>
            <plus/>
            <apply>
              <times/>
              <ci> p_RP_CaL </ci>
              <ci> k_RP_AP </ci>
            </apply>
            <apply>
              <times/>
              <ci> p_AI_CaL </ci>
              <ci> k_AI_AP </ci>
            </apply>
          </apply>
          <apply>
            <plus/>
            <apply>
              <times/>
              <ci> p_AP_CaL </ci>
              <ci> k_AP_RP </ci>
            </apply>
            <apply>
              <times/>
              <ci> p_AP_CaL </ci>
              <ci> k_AP_AI </ci>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply id="p_RP_CaL_diff_eq">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> p_RP_CaL </ci>
        </apply>    
        <apply>
          <minus/>
          <apply>
            <plus/>
            <apply>
              <times/>
              <ci> p_AP_CaL </ci>
              <ci> k_AP_RP </ci>
            </apply>
            <apply>
              <times/>
              <ci> p_RI_CaL </ci>
              <ci> k_RI_RP </ci>
            </apply>
          </apply>
          <apply>
            <plus/>
            <apply>
              <times/>
              <ci> p_RP_CaL </ci>
              <ci> k_RP_RI </ci>
            </apply>
            <apply>
              <times/>
              <ci> p_RP_CaL </ci>
              <ci> k_RP_AP </ci>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply id="p_AI_CaL_diff_eq">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> p_AI_CaL </ci>
        </apply>    
        <apply>
          <minus/>
          <apply>
            <plus/>
            <apply>
              <times/>
              <ci> p_RI_CaL </ci>
              <ci> k_RI_AI </ci>
            </apply>
            <apply>
              <times/>
              <ci> p_AP_CaL </ci>
              <ci> k_AP_AI </ci>
            </apply>
          </apply>
          <apply>
            <plus/>
            <apply>
              <times/>
              <ci> p_AI_CaL </ci>
              <ci> k_AI_RI </ci>
            </apply>
            <apply>
              <times/>
              <ci> p_AI_CaL </ci>
              <ci> k_AI_AP </ci>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply id="p_RI_CaL_diff_eq">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> p_RI_CaL </ci>
        </apply>    
        <apply>
          <minus/>
          <apply>
            <plus/>
            <apply>
              <times/>
              <ci> p_RP_CaL </ci>
              <ci> k_RP_RI </ci>
            </apply>
            <apply>
              <times/>
              <ci> p_AI_CaL </ci>
              <ci> k_AI_RI </ci>
            </apply>
          </apply>
          <apply>
            <plus/>
            <apply>
              <times/>
              <ci> p_RI_CaL </ci>
              <ci> k_RI_RP </ci>
            </apply>
            <apply>
              <times/>
              <ci> p_RI_CaL </ci>
              <ci> k_RI_AI </ci>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply id="L_type_Ca_channel_voltage_dependent_gate_k_RP_AP_calculation">
        <eq/>
        <ci> k_RP_AP </ci>
        <apply>
          <divide/>
          <cn cellml:units="first_order_rate_constant"> 1.0 </cn>
          <apply>
            <plus/>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> 0.27 </cn>
              <apply>
                <exp/>
                <apply>
                  <divide/>
                  <apply>
                    <minus/>
                    <ci> Vm </ci>
                  </apply>
                  <cn cellml:units="millivolt"> 5.9 </cn>
                </apply>
              </apply>
            </apply>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> 1.5 </cn>
              <apply>
                <exp/>
                <apply>
                  <divide/>
                  <apply>
                    <minus/>
                    <ci> Vm </ci>
                  </apply>
                  <cn cellml:units="millivolt"> 65.0 </cn>
                </apply>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply id="L_type_Ca_channel_voltage_dependent_gate_k_AP_RP_calculation">
        <eq/>
        <ci> k_AP_RP </ci>
        <apply>
          <divide/>
          <cn cellml:units="first_order_rate_constant"> 1.0 </cn>
          <apply>
            <plus/>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> 480.0 </cn>
              <apply>
                <exp/>
                <apply>
                  <divide/>
                  <ci> Vm </ci>
                  <cn cellml:units="millivolt"> 7.0 </cn>
                </apply>
              </apply>
            </apply>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> 2.2 </cn>
              <apply>
                <exp/>
                <apply>
                  <divide/>
                  <ci> Vm </ci>
                  <cn cellml:units="millivolt"> 65.0 </cn>
                </apply>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply id="L_type_Ca_channel_voltage_dependent_gate_k_RI_AI_calculation">
        <eq/>
        <ci> k_RI_AI </ci>
        <apply>
          <divide/>
          <cn cellml:units="first_order_rate_constant"> 1.0 </cn>
          <apply>
            <plus/>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> 0.0018 </cn>
              <apply>
                <exp/>
                <apply>
                  <divide/>
                  <apply>
                    <minus/>
                    <ci> Vm </ci>
                  </apply>
                  <cn cellml:units="millivolt"> 7.4 </cn>
                </apply>
              </apply>
            </apply>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> 2.0 </cn>
              <apply>
                <exp/>
                <apply>
                  <divide/>
                  <apply>
                    <minus/>
                    <ci> Vm </ci>
                  </apply>
                  <cn cellml:units="millivolt"> 100.0 </cn>
                </apply>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      
            <apply id="L_type_Ca_channel_voltage_dependent_gate_k_AI_RI_calculation">
        <eq/>
        <ci> k_AI_RI </ci>
        <apply>
          <divide/>
          <cn cellml:units="first_order_rate_constant"> 1.0 </cn>
          <apply>
            <plus/>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> 2200000.0 </cn>
              <apply>
                <exp/>
                <apply>
                  <divide/>
                  <ci> Vm </ci>
                  <cn cellml:units="millivolt"> 7.4 </cn>
                </apply>
              </apply>
            </apply>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> 11.0 </cn>
              <apply>
                <exp/>
                <apply>
                  <divide/>
                  <ci> Vm </ci>
                  <cn cellml:units="millivolt"> 100.0 </cn>
                </apply>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply id="L_type_Ca_channel_voltage_dependent_gate_k_RP_RI_calculation">
        <eq/>
        <ci> k_RP_RI </ci>
        <apply>
          <divide/>
          <cn cellml:units="first_order_rate_constant"> 0.04 </cn>
          <apply>
            <plus/>
            <cn cellml:units="first_order_rate_constant"> 1.0 </cn>
            <apply>
              <divide/>
              <apply>
                <times/>
                <ci> k_AI_AP </ci>
                <ci> k_AP_RP </ci>
                <ci> k_RI_AI </ci>
              </apply>
              <apply>
                <times/>
                <ci> k_AP_AI </ci>
                <ci> k_RP_AP </ci>
                <ci> k_AI_RI </ci>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply id="L_type_Ca_channel_voltage_dependent_gate_k_RI_RP_calculation">
        <eq/>
        <ci> k_RI_RP </ci>
        <apply>
          <minus/>
          <cn cellml:units="first_order_rate_constant"> 0.04 </cn>
          <ci> k_RP_RI </ci>
        </apply>
      </apply>
    </math> 
  </component>
  
  <component name="L_type_Ca_channel_Ca_dependent_gate">
    <variable units="dimensionless" public_interface="out" name="p_U" initial_value="1.715E-1"/>
    <variable units="dimensionless" public_interface="out" name="p_UCa" initial_value="6.118E-5"/> 
    <variable units="picoA" public_interface="out" name="iCaL"/>
    
    <variable units="millimolar" name="Cacm"/> 
    <variable units="dimensionless" name="p_CCa"/>
    <variable units="dimensionless" name="p_C" initial_value="4.228E-1"/>
    <variable units="first_order_rate_constant" name="k_CCa_UCa" initial_value="0.0003"/>
    <variable units="first_order_rate_constant" name="k_UCa_CCa" initial_value="0.35"/>
    <variable units="first_order_rate_constant" name="k_C_U" initial_value="0.143"/>
    <variable units="first_order_rate_constant" name="k_U_C" initial_value="0.35"/>
    <variable units="first_order_rate_constant" name="k_UCa_U"/>
    <variable units="second_order_rate_constant" name="k_U_UCa" initial_value="6.954"/>
    <variable units="first_order_rate_constant" name="k_CCa_C" initial_value="0.0042"/>
    <variable units="second_order_rate_constant" name="k_C_CCa" initial_value="6.954"/>
    
    <variable units="millisecond" public_interface="in" name="time"/>
    <variable units="millivolt" public_interface="in" name="Vm"/>
    <variable units="millimolar" public_interface="in" name="Cai"/>
    <variable units="millimolar" public_interface="in" name="CF_Ca"/>
    <variable units="dimensionless" public_interface="in" name="p_AP_CaL"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="iCaL_calculation">
        <eq/>
        <ci> iCaL </ci>
        <apply>
          <times/>
          <cn cellml:units="picoA_per_millimolar"> 0.0676 </cn>
          <ci> CF_Ca </ci>
        </apply>
      </apply>
        
      <apply id="Cacm_calculation">
        <eq/>
        <ci> Cacm </ci> 
        <apply>
          <minus/>
          <ci> Cai </ci>
          <apply>
            <times/>
            <cn cellml:units="millimolar_per_picoA"> 0.3 </cn>
            <ci> iCaL </ci>
          </apply>
        </apply>
      </apply>
     
      <apply id="p_U_diff_eq">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> p_U </ci>
        </apply>    
        <apply>
          <minus/>
          <apply>
            <plus/>
            <apply>
              <times/>
              <ci> p_C </ci>
              <ci> k_C_U </ci>
            </apply>
            <apply>
              <times/>
              <ci> p_UCa </ci>
              <ci> k_UCa_U </ci>
            </apply>
          </apply>
          <apply>
            <plus/>
            <apply>
              <times/>
              <ci> p_U </ci>
              <apply>
                <times/>
                <ci> k_U_UCa </ci> 
                <apply>
                  <plus/>
                  <apply>
                    <times/>
                    <ci> Cacm </ci>
                    <ci> p_AP_CaL </ci>
                  </apply>
                  <apply>
                    <times/>
                    <ci> Cai </ci>
                    <apply>
                      <minus/>
                      <cn cellml:units="dimensionless"> 1.0 </cn>
                      <ci> p_AP_CaL </ci>
                    </apply>
                  </apply>
                </apply>
              </apply>
            </apply>
            <apply>
              <times/>
              <ci> p_U </ci>
              <ci> k_U_C </ci>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply id="p_UCa_diff_eq">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> p_UCa </ci>
        </apply>   
        <apply>
          <minus/>
          <apply>
            <plus/>
            <apply>
              <times/>
              <ci> p_U </ci>
              <apply>
                <times/>
                <ci> k_U_UCa </ci> 
                <apply>
                  <plus/>
                  <apply>
                    <times/>
                    <ci> Cacm </ci>
                    <ci> p_AP_CaL </ci>
                  </apply>
                  <apply>
                    <times/>
                    <ci> Cai </ci>
                    <apply>
                      <minus/>
                      <cn cellml:units="dimensionless"> 1.0 </cn>
                      <ci> p_AP_CaL </ci>
                    </apply>
                  </apply>
                </apply>
              </apply>
            </apply>
            <apply>
              <times/>
              <ci> p_CCa </ci>
              <ci> k_CCa_UCa </ci>
            </apply>
          </apply>
          <apply>
            <plus/>
            <apply>
              <times/>
              <ci> p_UCa </ci>
              <ci> k_UCa_CCa </ci>
            </apply>
            <apply>
              <times/>
              <ci> p_UCa </ci>
              <ci> k_UCa_U </ci>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply id="p_C_diff_eq">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> p_C </ci>
        </apply>   
        <apply>
          <minus/>
          <apply>
            <plus/>
            <apply>
              <times/>
              <ci> p_CCa </ci>
              <ci> k_CCa_C </ci>
            </apply> 
            <apply>
              <times/>
              <ci> p_U </ci>
              <ci> k_U_C </ci>
            </apply>
          </apply>
          <apply>
            <plus/>
            <apply>
              <times/>
              <ci> p_C </ci>
              <ci> k_C_U </ci>
            </apply> 
            <apply>
              <times/>
              <ci> p_C </ci>
              <ci> k_C_U </ci>
              <ci> Cacm </ci>
              <ci> p_AP_CaL </ci>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply id="p_CCa_calculation">
        <eq/>
        <ci> p_CCa </ci>
        <apply>
          <minus/>
          <cn cellml:units="dimensionless"> 1.0 </cn>
          <apply>
            <plus/>
            <ci> p_C </ci>
            <ci> p_U </ci>
            <ci> p_UCa </ci>
          </apply>
        </apply>
      </apply>
      
      <apply id="k_UCa_U_calculation">
        <eq/>
        <ci> k_CCa_C </ci>
        <apply>
          <divide/>
          <apply>
            <times/>
            <ci> k_CCa_C </ci>
            <ci> k_C_U </ci>
            <ci> k_U_UCa </ci>
            <ci> k_UCa_CCa </ci>
          </apply>
          <apply>
            <times/>
            <ci> k_U_C </ci>
            <ci> k_C_CCa </ci>
            <ci> k_CCa_UCa </ci>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="L_type_Ca_channel_ultra_slow_gate">
    <variable units="dimensionless" public_interface="out" name="y" initial_value="9.984E-1"/> 

    <variable units="first_order_rate_constant" name="alpha_y"/>
    <variable units="first_order_rate_constant" name="beta_y"/>
    
    <variable units="millivolt" public_interface="in" name="Vm"/>
    <variable units="millisecond" public_interface="in" name="time"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="L_type_Ca_channel_ultra_slow_gate_alpha_y_calculation">
        <eq/>
        <ci> alpha_y </ci>
        <apply>
          <divide/>
          <cn cellml:units="first_order_rate_constant"> 1.0 </cn>
          <apply>
            <plus/>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> 250000.0 </cn>
              <apply>
                <exp/>
                <apply>
                  <divide/>
                  <ci> Vm </ci>
                  <cn cellml:units="millivolt"> 9.0 </cn>
                </apply>
              </apply>
            </apply>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> 58.0 </cn>
              <apply>
                <exp/>
                <apply>
                  <divide/>
                  <ci> Vm </ci>
                  <cn cellml:units="millivolt"> 65.0 </cn>
                </apply>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
    
      <apply id="L_type_Ca_channel_ultra_slow_gate_beta_y_calculation">
        <eq/>
        <ci> beta_y </ci>
        <apply>
          <divide/>
          <cn cellml:units="first_order_rate_constant"> 1.0 </cn>
          <apply>
            <plus/>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> 1800.0 </cn>
              <apply>
                <exp/>
                <apply>
                  <divide/>
                  <apply>
                    <minus/>
                    <ci> Vm </ci>
                  </apply>
                  <cn cellml:units="millivolt"> 14.0 </cn>
                </apply>
              </apply>
            </apply>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> 66.0 </cn>
              <apply>
                <exp/>
                <apply>
                  <divide/>
                  <apply>
                    <minus/>
                    <ci> Vm </ci>
                  </apply>
                  <cn cellml:units="millivolt"> 65.0 </cn>
                </apply>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
    
      <apply id="L_type_Ca_channel_ultra_slow_gate_dy_dt">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> y </ci>
        </apply>
        <apply>
          <minus/>
          <apply>
            <times/>
            <ci> alpha_y </ci>
            <apply>
              <minus/>
              <cn cellml:units="dimensionless"> 1.0 </cn>
              <ci> y </ci>
            </apply>
          </apply>
          <apply>
            <times/>
            <ci> beta_y </ci>
            <ci> y </ci>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
            
  <component name="T_type_Ca_channel">
    <variable units="picoA" public_interface="out" name="i_Ca_T"/>
    
    <variable units="millisecond" public_interface="in" private_interface="out" name="time"/>
    <variable units="millivolt" public_interface="in" private_interface="out" name="Vm"/>
    <variable units="millimolar" public_interface="in" name="CF_Ca"/>
  
    <variable units="dimensionless" private_interface="in" name="y1"/>
    <variable units="dimensionless" private_interface="in" name="y2"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">    
      <apply id="i_Ca_T_calculation">
        <eq/>
        <ci> i_Ca_T </ci>
        <apply>
          <times/>
          <cn cellml:units="dimensionless"> 612.0 </cn>
          <ci> CF_Ca </ci>
          <ci> y1 </ci>
          <ci> y2 </ci>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="T_type_Ca_channel_y1_gate">
    <variable units="dimensionless" public_interface="out" name="y1" initial_value="1.714E-5"/> 

    <variable units="first_order_rate_constant" name="alpha_y1"/>
    <variable units="first_order_rate_constant" name="beta_y1"/>
    
    <variable units="millivolt" public_interface="in" name="Vm"/>
    <variable units="millisecond" public_interface="in" name="time"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="T_type_Ca_channel_y1_gate_alpha_y1_calculation">
        <eq/>
        <ci> alpha_y1 </ci>
        <apply>
          <divide/>
          <cn cellml:units="first_order_rate_constant"> 1.0 </cn>
          <apply>
            <plus/>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> 0.019 </cn>
              <apply>
                <exp/>
                <apply>
                  <divide/>
                  <apply>
                    <minus/>
                    <ci> Vm </ci>
                  </apply>
                  <cn cellml:units="millivolt"> 5.6 </cn>
                </apply>
              </apply>
            </apply>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> 0.82 </cn>
              <apply>
                <exp/>
                <apply>
                  <divide/>
                  <apply>
                    <minus/>
                    <ci> Vm </ci>
                  </apply>
                  <cn cellml:units="millivolt"> 250.0 </cn>
                </apply>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
    
      <apply id="T_type_Ca_channel_y1_gate_beta_y1_calculation">
        <eq/>
        <ci> beta_y1 </ci>
        <apply>
          <divide/>
          <cn cellml:units="first_order_rate_constant"> 1.0 </cn>
          <apply>
            <plus/>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> 40.0 </cn>
              <apply>
                <exp/>
                <apply>
                  <divide/>
                  <ci> Vm </ci>
                  <cn cellml:units="millivolt"> 6.3 </cn>
                </apply>
              </apply>
            </apply>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> 1.5 </cn>
              <apply>
                <exp/>
                <apply>
                  <divide/>
                  <ci> Vm </ci>
                  <cn cellml:units="millivolt"> 10000.0 </cn>
                </apply>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
    
      <apply id="T_type_Ca_channel_y1_gate_dy1_dt">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> y1 </ci>
        </apply>
        <apply>
          <minus/>
          <apply>
            <times/>
            <ci> alpha_y1 </ci>
            <apply>
              <minus/>
              <cn cellml:units="dimensionless"> 1.0 </cn>
              <ci> y1 </ci>
            </apply>
          </apply>
          <apply>
            <times/>
            <ci> beta_y1 </ci>
            <ci> y1 </ci>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="T_type_Ca_channel_y2_gate">
    <variable units="dimensionless" public_interface="out" name="y2" initial_value="8.567E-1"/> 

    <variable units="first_order_rate_constant" name="alpha_y2"/>
    <variable units="first_order_rate_constant" name="beta_y2"/>
    
    <variable units="millivolt" public_interface="in" name="Vm"/>
    <variable units="millisecond" public_interface="in" name="time"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="T_type_Ca_channel_y2_gate_alpha_y2_calculation">
        <eq/>
        <ci> alpha_y2 </ci>
        <apply>
          <divide/>
          <cn cellml:units="first_order_rate_constant"> 1.0 </cn>
          <apply>
            <plus/>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> 62000.0 </cn>
              <apply>
                <exp/>
                <apply>
                  <divide/>
                  <ci> Vm </ci>
                  <cn cellml:units="millivolt"> 10.1 </cn>
                </apply>
              </apply>
            </apply>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> 30.0 </cn>
              <apply>
                <exp/>
                <apply>
                  <divide/>
                  <ci> Vm </ci>
                  <cn cellml:units="millivolt"> 3000.0 </cn>
                </apply>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
    
      <apply id="T_type_Ca_channel_y2_gate_beta_y2_calculation">
        <eq/>
        <ci> beta_y2 </ci>
        <apply>
          <divide/>
          <cn cellml:units="first_order_rate_constant"> 1.0 </cn>
          <apply>
            <plus/>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> 0.0006 </cn>
              <apply>
                <exp/>
                <apply>
                  <divide/>
                  <apply>
                    <minus/>
                    <ci> Vm </ci>
                  </apply>
                  <cn cellml:units="millivolt"> 6.7 </cn>
                </apply>
              </apply>
            </apply>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> 1.2 </cn>
              <apply>
                <exp/>
                <apply>
                  <divide/>
                  <apply>
                    <minus/>
                    <ci> Vm </ci>
                  </apply>
                  <cn cellml:units="millivolt"> 25.0 </cn>
                </apply>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
    
      <apply id="T_type_Ca_channel_y2_gate_dy2_dt">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> y2 </ci>
        </apply>
        <apply>
          <minus/>
          <apply>
            <times/>
            <ci> alpha_y2 </ci>
            <apply>
              <minus/>
              <cn cellml:units="dimensionless"> 1.0 </cn>
              <ci> y2 </ci>
            </apply>
          </apply>
          <apply>
            <times/>
            <ci> beta_y2 </ci>
            <ci> y2 </ci>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="time_independent_potassium_current">
    <variable units="picoA" public_interface="out" name="i_K1"/>
    <variable units="millivolt" public_interface="out" private_interface="out" name="E_K"/>
    
    <variable units="picoA_per_millivolt" name="g_K1"/>
    <variable units="dimensionless" name="fo"/>
    <variable units="dimensionless" name="fB"/>
    <variable units="first_order_rate_constant" name="mu"/>
    <variable units="first_order_rate_constant" name="lamda"/>
    
    <variable units="millisecond" public_interface="in" private_interface="out" name="time"/>
    <variable units="millivolt" public_interface="in" private_interface="out" name="Vm"/>
    <variable units="millimolar" public_interface="in" name="Ko"/>
    <variable units="millimolar" public_interface="in" name="Ki"/>
    <variable units="coulomb_millivolt_per_kelvin_millimole" public_interface="in" name="R"/>
    <variable units="kelvin" public_interface="in" name="T"/>
    <variable units="coulomb_per_millimole" public_interface="in" name="F"/>
    
    <variable units="dimensionless" private_interface="in" name="y"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="i_K1_calculation">
        <eq/>
        <ci> i_K1 </ci>
        <apply>
          <times/>
          <ci> g_K1 </ci>
          <apply>
            <minus/>     
            <ci> Vm </ci>
            <ci> E_K </ci>
          </apply>
          <apply>
            <plus/>
            <apply>
              <power/>
              <ci> fo </ci>
              <cn cellml:units="dimensionless"> 4.0 </cn>
            </apply>
            <apply>
              <times/>
              <apply>
                <divide/>
                <cn cellml:units="dimensionless"> 8.0 </cn>
                <cn cellml:units="dimensionless"> 3.0 </cn>
              </apply>
              <apply>
                <power/>
                <ci> fo </ci>
                <cn cellml:units="dimensionless"> 3.0 </cn>
              </apply>
              <ci> fB </ci>
            </apply>
            <apply>
              <times/>  
              <apply>
                <divide/>
                <cn cellml:units="dimensionless"> 6.0 </cn>
                <cn cellml:units="dimensionless"> 3.0 </cn>
              </apply>
              <apply>
                <power/>
                <ci> fo </ci>
                <cn cellml:units="dimensionless"> 2.0 </cn>
              </apply>
              <apply>
                <power/>
                <ci> fB </ci>
                <cn cellml:units="dimensionless"> 2.0 </cn>
              </apply>
            </apply>
          </apply>
          <ci> y </ci>  
        </apply>
      </apply>
      
      <apply id="E_K_calculation">
        <eq/>
        <ci> E_K </ci>
        <apply>
          <times/>
          <apply>
            <divide/>
            <apply>
              <times/>
              <ci> R </ci>
              <ci> T </ci>
            </apply>
            <ci> F </ci>
          </apply>
          <apply>
            <ln/>
            <apply>
              <divide/>
              <ci> Ko </ci>
              <ci> Ki </ci>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply id="g_K1_calculation">
        <eq/>
        <ci> g_K1 </ci>
        <apply>
          <times/>
          <cn cellml:units="picoA_per_millivolt"> 151.5 </cn>
          <apply>
            <power/>
            <apply>
              <divide/>
              <ci> Ko </ci>
              <cn cellml:units="millimolar"> 5.4 </cn>
            </apply>
            <cn cellml:units="dimensionless"> 0.4 </cn>
          </apply>
        </apply>
      </apply>
      
      <apply id="fB_calculation">
        <eq/>
        <ci> fB </ci>
        <apply>
          <divide/>
          <ci> mu </ci>
          <apply>
            <plus/>
            <ci> mu </ci>
            <ci> lamda </ci>
          </apply>
        </apply>
      </apply>
      
      <apply id="fo_calculation">
        <eq/>
        <ci> fo </ci>
        <apply>
          <divide/>
          <ci> lamda </ci>
          <apply>
            <plus/>
            <ci> mu </ci>
            <ci> lamda </ci>
          </apply>
        </apply>
      </apply> 
      
      <apply id="mu_calculation">
        <eq/>
        <ci> mu </ci>
        <apply>
          <times/>
          <cn cellml:units="first_order_rate_constant"> 0.75 </cn>
          <apply>
            <divide/>
            <apply>
              <exp/>
              <apply>
                <times/>
                <cn cellml:units="dimensionless"> 0.035 </cn>
                <apply>
                  <minus/>
                  <ci> Vm </ci>
                  <apply>
                    <plus/>
                    <ci> E_K </ci>
                    <cn cellml:units="millivolt"> 10.0 </cn>
                  </apply>
                </apply>
              </apply>
            </apply>
            <apply>
              <plus/>
              <cn cellml:units="dimensionless"> 1.0 </cn>
              <apply>
                <exp/>
                <apply>
                  <times/>
                  <cn cellml:units="dimensionless"> 0.015 </cn>
                  <apply>
                    <minus/>
                    <ci> Vm </ci>
                    <apply>
                      <plus/>
                      <ci> E_K </ci>
                      <cn cellml:units="millivolt"> 140.0 </cn>
                    </apply>
                  </apply>
                </apply>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply id="lamda_calculation">
        <eq/>
        <ci> lamda </ci>
        <apply>
          <times/>
          <cn cellml:units="first_order_rate_constant"> 3.0 </cn>
          <apply>
            <exp/>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> -0.048 </cn>
              <apply>
                <minus/>
                <ci> Vm </ci>
                <apply>
                  <plus/>
                  <ci> E_K </ci>
                  <cn cellml:units="millivolt"> 10.0 </cn>
                </apply>
              </apply>
            </apply>
          </apply>
          <apply>
            <divide/>
            <apply>
              <plus/>
              <cn cellml:units="dimensionless"> 1.0 </cn>
              <apply>
                <exp/>
                <apply>
                  <times/>
                  <cn cellml:units="dimensionless"> 0.064 </cn>
                  <apply>
                    <minus/>
                    <ci> Vm </ci>
                    <apply>
                      <plus/>
                      <ci> E_K </ci>
                      <cn cellml:units="millivolt"> 38.0 </cn>
                    </apply>
                  </apply>
                </apply>
              </apply>
            </apply>
            <apply>
              <plus/>
              <cn cellml:units="dimensionless"> 1.0 </cn>
              <apply>
                <exp/>
                <apply>
                  <times/>
                  <cn cellml:units="dimensionless"> 0.03 </cn>
                  <apply>
                    <minus/>
                    <ci> Vm </ci>
                    <apply>
                      <plus/>
                      <ci> E_K </ci>
                      <cn cellml:units="millivolt"> 70.0 </cn>
                    </apply>
                  </apply>
                </apply>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="time_independent_potassium_current_y_gate">
    <variable units="dimensionless" public_interface="out" name="y" initial_value="6.082E-1"/> 

    <variable units="first_order_rate_constant" name="alpha_y"/>
    <variable units="first_order_rate_constant" name="beta_y"/>
    
    <variable units="millivolt" public_interface="in" name="Vm"/>
    <variable units="millivolt" public_interface="in" name="E_K"/>
    <variable units="millisecond" public_interface="in" name="time"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="time_independent_potassium_current_y_gate_alpha_y_calculation">
        <eq/>
        <ci> alpha_y </ci>
        <apply>
          <divide/>
          <cn cellml:units="first_order_rate_constant"> 1.0 </cn>
          <apply>
            <plus/>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> 8000.0 </cn>
              <apply>
                <exp/>
                <apply>
                  <divide/>
                  <apply>
                    <minus/>
                    <ci> Vm </ci>
                    <apply>
                      <plus/>
                      <ci> E_K </ci>
                      <cn cellml:units="millivolt"> 97.0 </cn>
                    </apply>
                  </apply>
                  <cn cellml:units="millivolt"> 8.5 </cn>
                </apply>
              </apply>
            </apply>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> 7.0 </cn>
              <apply>
                <exp/>
                <apply>
                  <divide/>
                  <apply>
                    <minus/>
                    <ci> Vm </ci>
                    <apply>
                      <plus/>
                      <ci> E_K </ci>
                      <cn cellml:units="millivolt"> 97.0 </cn>
                    </apply>
                  </apply>
                  <cn cellml:units="millivolt"> 300.0 </cn>
                </apply>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
    
      <apply id="time_independent_potassium_current_y_gate_beta_y_calculation">
        <eq/>
        <ci> beta_y </ci>
        <apply>
          <divide/>
          <cn cellml:units="first_order_rate_constant"> 1.0 </cn>
          <apply>
            <plus/>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> 0.00014 </cn>
              <apply>
                <exp/>
                <apply>
                  <divide/>
                  <apply>
                    <minus/>
                    <apply>
                      <minus/>
                      <ci> Vm </ci>
                      <apply>
                        <plus/>
                        <ci> E_K </ci>
                        <cn cellml:units="millivolt"> 97.0 </cn>
                      </apply>
                    </apply>
                  </apply>
                  <cn cellml:units="millivolt"> 9.1 </cn>
                </apply>
              </apply>
            </apply>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> 0.2 </cn>
              <apply>
                <exp/>
                <apply>
                  <divide/>
                  <apply>
                    <minus/>
                    <apply>
                      <minus/>
                      <ci> Vm </ci>
                      <apply>
                        <plus/>
                        <ci> E_K </ci>
                        <cn cellml:units="millivolt"> 97.0 </cn>
                      </apply>
                    </apply>
                  </apply>
                  <cn cellml:units="millivolt"> 500.0 </cn>
                </apply>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
    
      <apply id="time_independent_potassium_current_y_gate_dy_dt">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> y </ci>
        </apply>
        <apply>
          <minus/>
          <apply>
            <times/>
            <ci> alpha_y </ci>
            <apply>
              <minus/>
              <cn cellml:units="dimensionless"> 1.0 </cn>
              <ci> y </ci>
            </apply>
          </apply>
          <apply>
            <times/>
            <ci> beta_y </ci>
            <ci> y </ci>
          </apply>
        </apply>
      </apply>
    </math>
  </component>

  <component name="rapid_time_dependent_potassium_current">
    <variable units="picoA" public_interface="out" name="i_Kr"/>
    
    <variable units="picoA_per_millivolt" name="g_Kr"/>
    
    <variable units="millisecond" public_interface="in" private_interface="out" name="time"/>
    <variable units="millivolt" public_interface="in" private_interface="out" name="Vm"/>
    <variable units="millivolt" public_interface="in" private_interface="out" name="E_K"/>
    <variable units="millimolar" public_interface="in" name="Ko"/>
    <variable units="picoF" public_interface="in" name="Cm"/>
    
    <variable units="dimensionless" private_interface="in" name="y1"/>
    <variable units="dimensionless" private_interface="in" name="y2"/>
    <variable units="dimensionless" private_interface="in" name="y3"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="i_Kr_calculation">
        <eq/>
        <ci> i_Kr </ci>
        <apply>
          <times/>
          <ci> g_Kr </ci>
          <ci> Cm </ci>
          <apply>
            <minus/>     
            <ci> Vm </ci>
            <ci> E_K </ci>
          </apply>
          <apply>
            <plus/>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> 0.6 </cn>
              <ci> y1 </ci>
            </apply>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> 0.4 </cn>
              <ci> y2 </ci>
            </apply>
          </apply>
          <ci> y3 </ci>  
        </apply>
      </apply>
      
      <apply id="g_Kr_calculation">
        <eq/>
        <ci> g_Kr </ci>
        <apply>
          <times/>
          <cn cellml:units="picoA_per_millivolt"> 0.00864 </cn>
          <apply>
            <power/>
            <apply>
              <divide/>
              <ci> Ko </ci>
              <cn cellml:units="millimolar"> 5.4 </cn>
            </apply>
            <cn cellml:units="dimensionless"> 0.2 </cn>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="rapid_time_dependent_potassium_current_y1_gate">
    <variable units="dimensionless" public_interface="out" name="y1" initial_value="1.935E-3"/> 

    <variable units="first_order_rate_constant" name="alpha_y1"/>
    <variable units="first_order_rate_constant" name="beta_y1"/>
    
    <variable units="millivolt" public_interface="in" name="Vm"/>
    <variable units="millivolt" public_interface="in" name="E_K"/>
    <variable units="millisecond" public_interface="in" name="time"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="rapid_time_dependent_potassium_current_y1_gate_alpha_y1_calculation">
        <eq/>
        <ci> alpha_y1 </ci>
        <apply>
          <divide/>
          <cn cellml:units="first_order_rate_constant"> 1.0 </cn>
          <apply>
            <plus/>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> 20.0 </cn>
              <apply>
                <exp/>
                <apply>
                  <divide/>
                  <apply>
                    <minus/>
                    <ci> Vm </ci>
                  </apply>
                  <cn cellml:units="millivolt"> 11.5 </cn>
                </apply>
              </apply>
            </apply>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> 5.0 </cn>
              <apply>
                <exp/>
                <apply>
                  <divide/>
                  <apply>
                    <minus/>
                    <ci> Vm </ci>
                  </apply>
                  <cn cellml:units="millivolt"> 300.0 </cn>
                </apply>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
    
      <apply id="rapid_time_dependent_potassium_current_y1_gate_beta_y1_calculation">
        <eq/>
        <ci> beta_y1 </ci>
        <apply>
          <plus/>
          <apply>
            <divide/>
            <cn cellml:units="first_order_rate_constant"> 1.0 </cn>
            <apply>
              <plus/>
              <apply>
                <times/>
                <cn cellml:units="dimensionless"> 160.0 </cn>
                <apply>
                  <exp/>
                  <apply>
                    <divide/>
                    <ci> Vm </ci>
                    <cn cellml:units="millivolt"> 28.0 </cn>
                  </apply>
                </apply>
              </apply>
              <apply>
                <times/>
                <cn cellml:units="dimensionless"> 200.0 </cn>
                <apply>
                  <exp/>
                  <apply>
                    <divide/>
                    <ci> Vm </ci>
                    <cn cellml:units="millivolt"> 1000.0 </cn>
                  </apply>
                </apply>
              </apply>
            </apply>
          </apply>
          <apply>
            <divide/>
            <cn cellml:units="first_order_rate_constant"> 1.0 </cn>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> 2500.0 </cn>
              <apply>
                <exp/>
                <apply>
                  <divide/>
                  <ci> Vm </ci>
                  <cn cellml:units="millivolt"> 20.0 </cn>
                </apply>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
    
      <apply id="rapid_time_dependent_potassium_current_y1_gate_dy1_dt">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> y1 </ci>
        </apply>
        <apply>
          <minus/>
          <apply>
            <times/>
            <ci> alpha_y1 </ci>
            <apply>
              <minus/>
              <cn cellml:units="dimensionless"> 1.0 </cn>
              <ci> y1 </ci>
            </apply>
          </apply>
          <apply>
            <times/>
            <ci> beta_y1 </ci>
            <ci> y1 </ci>
          </apply>
        </apply>
      </apply>
    </math>
  </component>

  <component name="rapid_time_dependent_potassium_current_y2_gate">
    <variable units="dimensionless" public_interface="out" name="y2" initial_value="2.090E-1"/> 

    <variable units="first_order_rate_constant" name="alpha_y2"/>
    <variable units="first_order_rate_constant" name="beta_y2"/>
    
    <variable units="millivolt" public_interface="in" name="Vm"/>
    <variable units="millivolt" public_interface="in" name="E_K"/>
    <variable units="millisecond" public_interface="in" name="time"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="rapid_time_dependent_potassium_current_y2_gate_alpha_y2_calculation">
        <eq/>
        <ci> alpha_y2 </ci>
        <apply>
          <divide/>
          <cn cellml:units="first_order_rate_constant"> 1.0 </cn>
          <apply>
            <plus/>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> 200.0 </cn>
              <apply>
                <exp/>
                <apply>
                  <divide/>
                  <apply>
                    <minus/>
                    <ci> Vm </ci>
                  </apply>
                  <cn cellml:units="millivolt"> 13.0 </cn>
                </apply>
              </apply>
            </apply>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> 20.0 </cn>
              <apply>
                <exp/>
                <apply>
                  <divide/>
                  <apply>
                    <minus/>
                    <ci> Vm </ci>
                  </apply>
                  <cn cellml:units="millivolt"> 300.0 </cn>
                </apply>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
    
      <apply id="rapid_time_dependent_potassium_current_y2_gate_beta_y2_calculation">
        <eq/>
        <ci> beta_y2 </ci>
        <apply>
          <plus/>
          <apply>
            <divide/>
            <cn cellml:units="first_order_rate_constant"> 1.0 </cn>
            <apply>
              <plus/>
              <apply>
                <times/>
                <cn cellml:units="dimensionless"> 1600.0 </cn>
                <apply>
                  <exp/>
                  <apply>
                    <divide/>
                    <ci> Vm </ci>
                    <cn cellml:units="millivolt"> 28.0 </cn>
                  </apply>
                </apply>
              </apply>
              <apply>
                <times/>
                <cn cellml:units="dimensionless"> 2000.0 </cn>
                <apply>
                  <exp/>
                  <apply>
                    <divide/>
                    <ci> Vm </ci>
                    <cn cellml:units="millivolt"> 1000.0 </cn>
                  </apply>
                </apply>
              </apply>
            </apply>
          </apply>
          <apply>
            <divide/>
            <cn cellml:units="first_order_rate_constant"> 1.0 </cn>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> 10000.0 </cn>
              <apply>
                <exp/>
                <apply>
                  <divide/>
                  <ci> Vm </ci>
                  <cn cellml:units="millivolt"> 20.0 </cn>
                </apply>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
    
      <apply id="rapid_time_dependent_potassium_current_y2_gate_dy2_dt">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> y2 </ci>
        </apply>
        <apply>
          <minus/>
          <apply>
            <times/>
            <ci> alpha_y2 </ci>
            <apply>
              <minus/>
              <cn cellml:units="dimensionless"> 1.0 </cn>
              <ci> y2 </ci>
            </apply>
          </apply>
          <apply>
            <times/>
            <ci> beta_y2 </ci>
            <ci> y2 </ci>
          </apply>
        </apply>
      </apply>
    </math>
  </component>

  <component name="rapid_time_dependent_potassium_current_y3_gate">
    <variable units="dimensionless" public_interface="out" name="y3" initial_value="9.677E-1"/> 

    <variable units="first_order_rate_constant" name="alpha_y3"/>
    <variable units="first_order_rate_constant" name="beta_y3"/>
    
    <variable units="millivolt" public_interface="in" name="Vm"/>
    <variable units="millivolt" public_interface="in" name="E_K"/>
    <variable units="millisecond" public_interface="in" name="time"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="alpha_y3_calculation">
        <eq/>
        <ci> alpha_y3 </ci>
        <apply>
          <divide/>
          <cn cellml:units="first_order_rate_constant"> 1.0 </cn>
          <apply>
            <plus/>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> 10.0 </cn>
              <apply>
                <exp/>
                <apply>
                  <divide/>
                  <ci> Vm </ci>
                  <cn cellml:units="millivolt"> 17.0 </cn>
                </apply>
              </apply>
            </apply>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> 2.5 </cn>
              <apply>
                <exp/>
                <apply>
                  <divide/>
                  <ci> Vm </ci>
                  <cn cellml:units="millivolt"> 300.0 </cn>
                </apply>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
    
      <apply id="beta_y3_calculation">
        <eq/>
        <ci> beta_y3 </ci>
        <apply>
          <divide/>
          <cn cellml:units="first_order_rate_constant"> 1.0 </cn>
          <apply>
            <plus/>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> 0.35 </cn>
              <apply>
                <exp/>
                <apply>
                  <divide/>
                  <apply>
                    <minus/>
                    <ci> Vm </ci>
                  </apply>
                  <cn cellml:units="millivolt"> 17.0 </cn>
                </apply>
              </apply>
            </apply>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> 2.0 </cn>
              <apply>
                <exp/>
                <apply>
                  <divide/>
                  <apply>
                    <minus/>
                    <ci> Vm </ci>
                  </apply>
                  <cn cellml:units="millivolt"> 150.0 </cn>
                </apply>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>  
    
      <apply id="dy3_dt">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> y3 </ci>
        </apply>
        <apply>
          <minus/>
          <apply>
            <times/>
            <ci> alpha_y3 </ci>
            <apply>
              <minus/>
              <cn cellml:units="dimensionless"> 1.0 </cn>
              <ci> y3 </ci>
            </apply>
          </apply>
          <apply>
            <times/>
            <ci> beta_y3 </ci>
            <ci> y3 </ci>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="slow_time_dependent_potassium_current">
    <variable units="picoA" public_interface="out" name="i_Ks"/>
    <variable units="picoA" public_interface="out" name="i_Ks_Na"/>
    <variable units="picoA" public_interface="out" name="i_Ks_K"/>
    
    <variable units="millisecond" public_interface="in" private_interface="out" name="time"/>
    <variable units="millivolt" public_interface="in" private_interface="out" name="Vm"/>
    <variable units="millimolar" public_interface="in" private_interface="out" name="Cai"/>
    <variable units="millimolar" public_interface="in" name="CF_K"/>
    <variable units="millimolar" public_interface="in" name="CF_Na"/>
    
    <variable units="dimensionless" private_interface="in" name="y1"/>
    <variable units="dimensionless" private_interface="in" name="y2"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="i_Ks_calculation">
        <eq/>
        <ci> i_Ks </ci>
        <apply>
          <plus/>
          <ci> i_Ks_Na </ci>
          <ci> i_Ks_K </ci>
        </apply>
      </apply>
      
      <apply id="i_Ks_K_calculation">
        <eq/>
        <ci> i_Ks_K </ci>
        <apply>
          <times/>
          <cn cellml:units="picoA_per_millimolar"> 5.04 </cn>
          <ci> CF_K </ci>
          <apply>
            <power/>     
            <ci> y1 </ci>
            <cn cellml:units="dimensionless"> 2.0 </cn>
          </apply>
          <apply>
            <plus/>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> 0.9 </cn>
              <ci> y2 </ci>
            </apply>
            <cn cellml:units="dimensionless"> 0.1 </cn>
          </apply>
        </apply>
      </apply>
      
      <apply id="i_Ks_Na_calculation">
        <eq/>
        <ci> i_Ks_Na </ci>
        <apply>
          <times/>
          <cn cellml:units="picoA_per_millimolar"> 0.2016 </cn>
          <ci> CF_Na </ci>
          <apply>
            <power/>     
            <ci> y1 </ci>
            <cn cellml:units="dimensionless"> 2.0 </cn>
          </apply>
          <apply>
            <plus/>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> 0.9 </cn>
              <ci> y2 </ci>
            </apply>
            <cn cellml:units="dimensionless"> 0.1 </cn>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="slow_time_dependent_potassium_current_y1_gate">
    <variable units="dimensionless" public_interface="out" name="y1" initial_value="9.996E-2"/> 

    <variable units="first_order_rate_constant" name="alpha_y1"/>
    <variable units="first_order_rate_constant" name="beta_y1"/>
    
    <variable units="millivolt" public_interface="in" name="Vm"/>
    <variable units="millisecond" public_interface="in" name="time"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="slow_time_dependent_potassium_current_y1_gate_alpha_y1_calculation">
        <eq/>
        <ci> alpha_y1 </ci>
        <apply>
          <divide/>
          <cn cellml:units="first_order_rate_constant"> 1.0 </cn>
          <apply>
            <plus/>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> 85.0 </cn>
              <apply>
                <exp/>
                <apply>
                  <divide/>
                  <apply>
                    <minus/>
                    <ci> Vm </ci>
                  </apply>
                  <cn cellml:units="millivolt"> 10.5 </cn>
                </apply>
              </apply>
            </apply>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> 370.0 </cn>
              <apply>
                <exp/>
                <apply>
                  <divide/>
                  <apply>
                    <minus/>
                    <ci> Vm </ci>
                  </apply>
                  <cn cellml:units="millivolt"> 62.0 </cn>
                </apply>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
    
      <apply id="slow_time_dependent_potassium_current_y1_gate_beta_y1_calculation">
        <eq/>
        <ci> beta_y1 </ci>
        <apply>
          <divide/>
          <cn cellml:units="first_order_rate_constant"> 1.0 </cn>
          <apply>
            <plus/>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> 1450.0 </cn>
              <apply>
                <exp/>
                <apply>
                  <divide/>
                  <ci> Vm </ci>
                  <cn cellml:units="millivolt"> 20.0 </cn>
                </apply>
              </apply>
            </apply>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> 260.0 </cn>
              <apply>
                <exp/>
                <apply>
                  <divide/>
                  <ci> Vm </ci>
                  <cn cellml:units="millivolt"> 100.0 </cn>
                </apply>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
     
      <apply id="slow_time_dependent_potassium_current_y1_gate_dy1_dt">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> y1 </ci>
        </apply>
        <apply>
          <minus/>
          <apply>
            <times/>
            <ci> alpha_y1 </ci>
            <apply>
              <minus/>
              <cn cellml:units="dimensionless"> 1.0 </cn>
              <ci> y1 </ci>
            </apply>
          </apply>
          <apply>
            <times/>
            <ci> beta_y1 </ci>
            <ci> y1 </ci>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="slow_time_dependent_potassium_current_y2_gate">
    <variable units="dimensionless" public_interface="out" name="y2" initial_value="9.266E-2"/> 

    <variable units="first_order_rate_constant" name="alpha_y2"/>
    <variable units="first_order_rate_constant" name="beta_y2" initial_value="0.004444"/>
    
    <variable units="millimolar" public_interface="in" name="Cai"/>
    <variable units="millivolt" public_interface="in" name="Vm"/>
    <variable units="millisecond" public_interface="in" name="time"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="slow_time_dependent_potassium_current_y2_gate_alpha_y2_calculation">
        <eq/>
        <ci> alpha_y2 </ci>
        <apply>
          <times/>
          <cn cellml:units="second_order_rate_constant"> 3.7 </cn>
          <ci> Cai </ci>
        </apply>
      </apply>
    
      <apply id="slow_time_dependent_potassium_current_y2_gate_dy2_dt">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> y2 </ci>
        </apply>
        <apply>
          <minus/>
          <apply>
            <times/>
            <ci> alpha_y2 </ci>
            <apply>
              <minus/>
              <cn cellml:units="dimensionless"> 1.0 </cn>
              <ci> y2 </ci>
            </apply>
          </apply>
          <apply>
            <times/>
            <ci> beta_y2 </ci>
            <ci> y2 </ci>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="transient_outward_current">
    <variable units="picoA" public_interface="out" name="i_to"/>
    <variable units="picoA" public_interface="out" name="i_to_Na"/>
    <variable units="picoA" public_interface="out" name="i_to_K"/>
    
    <variable units="millisecond" public_interface="in" private_interface="out" name="time"/>
    <variable units="millivolt" public_interface="in" private_interface="out" name="Vm"/>
    <variable units="millimolar" public_interface="in" name="CF_K"/>
    <variable units="millimolar" public_interface="in" name="CF_Na"/>
    
    <variable units="dimensionless" private_interface="in" name="y1"/>
    <variable units="dimensionless" private_interface="in" name="y2"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="i_to_calculation">
        <eq/>
        <ci> i_to </ci>
        <apply>
          <plus/>
          <ci> i_to_Na </ci>
          <ci> i_to_K </ci>
        </apply>
      </apply>
      
      <apply id="i_to_K_calculation">
        <eq/>
        <ci> i_to_K </ci>
        <apply>
          <times/>
          <cn cellml:units="picoA_per_millimolar"> 0.033 </cn>
          <ci> CF_K </ci>
          <apply>
            <power/>     
            <ci> y1 </ci>
            <cn cellml:units="dimensionless"> 3.0 </cn>
          </apply>
          <ci> y2 </ci>
        </apply>
      </apply>
      
      <apply id="i_to_Na_calculation">
        <eq/>
        <ci> i_to_Na </ci>
        <apply>
          <times/>
          <cn cellml:units="picoA_per_millimolar"> 0.00297 </cn>
          <ci> CF_Na </ci>
          <apply>
            <power/>     
            <ci> y1 </ci>
            <cn cellml:units="dimensionless"> 3.0 </cn>
          </apply>
          <ci> y2 </ci>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="transient_outward_current_y1_gate">
    <variable units="dimensionless" public_interface="out" name="y1" initial_value="7.985E-4"/> 

    <variable units="first_order_rate_constant" name="alpha_y1"/>
    <variable units="first_order_rate_constant" name="beta_y1"/>
    
    <variable units="millivolt" public_interface="in" name="Vm"/>
    <variable units="millisecond" public_interface="in" name="time"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="transient_outward_current_y1_gate_alpha_y1_calculation">
        <eq/>
        <ci> alpha_y1 </ci>
        <apply>
          <divide/>
          <cn cellml:units="first_order_rate_constant"> 1.0 </cn>
          <apply>
            <plus/>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> 11.0 </cn>
              <apply>
                <exp/>
                <apply>
                  <divide/>
                  <apply>
                    <minus/>
                    <ci> Vm </ci>
                  </apply>
                  <cn cellml:units="millivolt"> 28.0 </cn>
                </apply>
              </apply>
            </apply>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> 0.2 </cn>
              <apply>
                <exp/>
                <apply>
                  <divide/>
                  <apply>
                    <minus/>
                    <ci> Vm </ci>
                  </apply>
                  <cn cellml:units="millivolt"> 400.0 </cn>
                </apply>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
    
      <apply id="transient_outward_current_y1_gate_beta_y1_calculation">
        <eq/>
        <ci> beta_y1 </ci>
        <apply>
          <divide/>
          <cn cellml:units="first_order_rate_constant"> 1.0 </cn>
          <apply>
            <plus/>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> 4.4 </cn>
              <apply>
                <exp/>
                <apply>
                  <divide/>
                  <ci> Vm </ci>
                  <cn cellml:units="millivolt"> 16.0 </cn>
                </apply>
              </apply>
            </apply>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> 0.2 </cn>
              <apply>
                <exp/>
                <apply>
                  <divide/>
                  <ci> Vm </ci>
                  <cn cellml:units="millivolt"> 500.0 </cn>
                </apply>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
     
      <apply id="transient_outward_current_y1_gate_dy1_dt">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> y1 </ci>
        </apply>
        <apply>
          <minus/>
          <apply>
            <times/>
            <ci> alpha_y1 </ci>
            <apply>
              <minus/>
              <cn cellml:units="dimensionless"> 1.0 </cn>
              <ci> y1 </ci>
            </apply>
          </apply>
          <apply>
            <times/>
            <ci> beta_y1 </ci>
            <ci> y1 </ci>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="transient_outward_current_y2_gate">
    <variable units="dimensionless" public_interface="out" name="y2" initial_value="9.999E-1"/> 

    <variable units="first_order_rate_constant" name="alpha_y2"/>
    <variable units="first_order_rate_constant" name="beta_y2"/>
    
    <variable units="millivolt" public_interface="in" name="Vm"/>
    <variable units="millisecond" public_interface="in" name="time"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="transient_outward_current_y2_gate_alpha_y2_calculation">
        <eq/>
        <ci> alpha_y2 </ci>
        <apply>
          <divide/>
          <apply>
            <times/>
            <cn cellml:units="first_order_rate_constant"> 0.0038 </cn>
            <apply>
              <exp/>
              <apply>
                <divide/>
                <apply>
                  <minus/>
                  <apply>
                    <plus/>
                    <ci> Vm </ci>
                    <cn cellml:units="millivolt"> 13.5 </cn>
                  </apply>
                </apply>
                <cn cellml:units="millivolt"> 11.3 </cn>
              </apply>
            </apply>
          </apply>
          <apply>
            <plus/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> 0.051335 </cn>
              <apply>
                <exp/>
                <apply>
                  <divide/>
                  <apply>
                    <minus/>
                    <apply>
                      <plus/>
                      <ci> Vm </ci>
                      <cn cellml:units="millivolt"> 13.5 </cn>
                    </apply>
                  </apply>
                  <cn cellml:units="millivolt"> 11.3 </cn>
                </apply>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
    
      <apply id="transient_outward_current_y2_gate_beta_y2_calculation">
        <eq/>
        <ci> beta_y2 </ci>
        <apply>
          <divide/>
          <apply>
            <times/>
            <cn cellml:units="first_order_rate_constant"> 0.0038 </cn>
            <apply>
              <exp/>
              <apply>
                <divide/>
                <apply>
                  <plus/>
                  <ci> Vm </ci>
                  <cn cellml:units="millivolt"> 13.5 </cn>
                </apply>
                <cn cellml:units="millivolt"> 11.3 </cn>
              </apply>
            </apply>
          </apply>
          <apply>
            <plus/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> 0.067083 </cn>
              <apply>
                <exp/>
                <apply>
                  <divide/>
                  <apply>
                    <plus/>
                    <ci> Vm </ci>
                    <cn cellml:units="millivolt"> 13.5 </cn>
                  </apply>
                  <cn cellml:units="millivolt"> 11.3 </cn>
                </apply>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply id="transient_outward_current_y2_gate_dy2_dt">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> y2 </ci>
        </apply>
        <apply>
          <minus/>
          <apply>
            <times/>
            <ci> alpha_y2 </ci>
            <apply>
              <minus/>
              <cn cellml:units="dimensionless"> 1.0 </cn>
              <ci> y2 </ci>
            </apply>
          </apply>
          <apply>
            <times/>
            <ci> beta_y2 </ci>
            <ci> y2 </ci>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="background_NSC_current">
    <variable units="picoA" public_interface="out" name="i_bNSC"/>
    <variable units="picoA" public_interface="out" name="i_bNSC_Na"/>
    <variable units="picoA" public_interface="out" name="i_bNSC_K"/>
    
    <variable units="picoA_per_millimolar" name="P_bNSC" initial_value="0.385"/>
    
    <variable units="millisecond" public_interface="in" name="time"/>
    <variable units="millimolar" public_interface="in" name="CF_K"/>
    <variable units="millimolar" public_interface="in" name="CF_Na"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="i_bNSC_calculation">
        <eq/>
        <ci> i_bNSC </ci>
        <apply>
          <plus/>
          <ci> i_bNSC_K </ci>
          <ci> i_bNSC_Na </ci>
        </apply>
      </apply>
      
      <apply id="i_bNSC_K_calculation">
        <eq/>
        <ci> i_bNSC_K </ci>
        <apply>
          <times/>
          <ci> P_bNSC </ci>
          <ci> CF_K </ci>
        </apply> 
      </apply>
      
      <apply id="i_bNSC_Na_calculation">
        <eq/>
        <ci> i_bNSC_Na </ci>
        <apply>
          <times/>
          <cn cellml:units="picoA_per_millimolar"> 0.4 </cn>
          <ci> P_bNSC </ci>
          <ci> CF_Na </ci>
        </apply> 
      </apply>
    </math>
  </component>
  
  <component name="background_Kpl_current">
    <variable units="picoA" public_interface="out" name="i_Kpl"/>
    
    <variable units="picoA_per_millimolar" name="P_Kpl"/>
    
    <variable units="millisecond" public_interface="in" name="time"/>
    <variable units="millivolt" public_interface="in" name="Vm"/>
    <variable units="millimolar" public_interface="in" name="CF_K"/>
    <variable units="millimolar" public_interface="in" name="Ko"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="i_Kpl_calculation">
        <eq/>
        <ci> i_Kpl </ci>
        <apply>
          <times/>
          <ci> P_Kpl </ci>
          <ci> CF_K </ci>
          <apply>
            <divide/>
            <apply>
              <plus/>
              <ci> Vm </ci>
              <cn cellml:units="millivolt"> 3.0 </cn>
            </apply>
            <apply>
              <minus/>
              <cn cellml:units="dimensionless"> 1.0 </cn>
              <apply>
                <exp/>
                <apply>
                  <divide/>
                  <apply>
                    <minus/>
                    <apply>
                      <plus/>
                      <ci> Vm </ci>
                      <cn cellml:units="millivolt"> 3.0 </cn>
                    </apply>  
                  </apply>
                  <cn cellml:units="millivolt"> 13.0 </cn>
                </apply>
              </apply>
            </apply>
          </apply> 
        </apply> 
      </apply>
      
      <apply id="P_Kpl_calculation">
        <eq/>
        <ci> P_Kpl </ci>
        <apply>
          <times/>
          <cn cellml:units="picoA_per_millimolar"> 0.00011 </cn>
          <apply>
            <power/>
            <apply>
              <divide/>
              <ci> Ko </ci>
              <cn cellml:units="millimolar"> 5.4 </cn>
            </apply>
            <cn cellml:units="dimensionless"> 0.16 </cn>
          </apply>
        </apply> 
      </apply>
    </math>
  </component>
  
  <component name="background_lCa_current">
    <variable units="picoA" public_interface="out" name="i_lCa"/>
    <variable units="picoA" public_interface="out" name="i_lCa_Na"/>
    <variable units="picoA" public_interface="out" name="i_lCa_K"/>
    
    <variable units="picoA_per_millimolar" name="P_lCa" initial_value="0.11"/>
    <variable units="dimensionless" name="p_open"/>
    
    <variable units="millisecond" public_interface="in" name="time"/>
    <variable units="millimolar" public_interface="in" name="CF_K"/>
    <variable units="millimolar" public_interface="in" name="CF_Na"/>
    <variable units="millimolar" public_interface="in" name="Cai"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="i_lCa_calculation">
        <eq/>
        <ci> i_lCa </ci>
        <apply>
          <plus/>
          <ci> i_lCa_K </ci>
          <ci> i_lCa_Na </ci>
        </apply>
      </apply>
      
      <apply id="i_lCa_K_calculation">
        <eq/>
        <ci> i_lCa_K </ci>
        <apply>
          <times/>
          <ci> P_lCa </ci>
          <ci> CF_K </ci>
          <ci> p_open </ci>
        </apply> 
      </apply>
      
      <apply id="i_lCa_Na_calculation">
        <eq/>
        <ci> i_lCa_Na </ci>
        <apply>
          <times/>
          <ci> P_lCa </ci>
          <ci> CF_Na </ci>
          <ci> p_open </ci>
        </apply> 
      </apply>
      
      <apply id="background_lCa_current_p_open_calculation">
        <eq/>
        <ci> p_open </ci>
        <apply>
          <divide/>
          <cn cellml:units="dimensionless"> 1.0 </cn>
          <apply>
            <plus/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
              <power/>
              <apply>
                <divide/>
                <cn cellml:units="millimolar"> 0.0012 </cn>
                <ci> Cai </ci>
              </apply>
              <cn cellml:units="dimensionless"> 3.0 </cn>
            </apply>
          </apply>
        </apply>
      </apply> 
    </math>
  </component>
  
  <component name="background_KATP_current">
    <variable units="picoA" public_interface="out" name="i_KATP"/>
    
    <variable units="dimensionless" name="N" initial_value="2333.0"/>
    <variable units="dimensionless" name="p_open"/>
    <variable units="dimensionless" name="gamma"/>
    
    <variable units="millisecond" public_interface="in" name="time"/>
    <variable units="millivolt" public_interface="in" name="Vm"/>
    <variable units="millivolt" public_interface="in" name="E_K"/>
    <variable units="millimolar" public_interface="in" name="Ko"/>
    <variable units="millimolar" public_interface="in" name="ATPi"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="i_KATP_calculation">
        <eq/>
        <ci> i_KATP </ci>
        <apply>
          <times/>
          <ci> N </ci>
          <ci> gamma </ci>
          <apply>
            <minus/>
            <ci> Vm </ci>
            <ci> E_K </ci>
          </apply> 
          <ci> p_open </ci>
        </apply> 
      </apply>
      
      <apply id="gamma_calculation">
        <eq/>
        <ci> gamma </ci>
        <apply>
          <times/>
          <cn cellml:units="dimensionless"> 0.0236 </cn>
          <apply>
            <power/>
            <ci> Ko </ci>
            <cn cellml:units="dimensionless"> 0.24 </cn>
          </apply>
        </apply> 
      </apply>
      
      <apply id="background_KATP_current_p_open_calculation">
        <eq/>
        <ci> p_open </ci>
        <apply>
          <divide/>
          <cn cellml:units="dimensionless"> 0.8 </cn>
          <apply>
            <plus/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
              <power/>
              <apply>
                <divide/>
                <ci> ATPi </ci>
                <cn cellml:units="millimolar"> 0.1 </cn>
              </apply>
              <cn cellml:units="dimensionless"> 2.0 </cn>
            </apply>
          </apply>
        </apply>
      </apply> 
    </math>
  </component>
  
  <component name="background_Cab_current">
    <variable units="picoA" public_interface="out" name="i_Cab"/>
    
    <variable units="picoA_per_millimolar" name="P_Cab" initial_value="0.04"/>
    
    <variable units="millisecond" public_interface="in" name="time"/>
    <variable units="millimolar" public_interface="in" name="CF_Ca"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="i_Cab_calculation">
        <eq/>
        <ci> i_Cab </ci>
        <apply>
          <times/>
          <ci> P_Cab </ci>
          <ci> CF_Ca </ci>
        </apply> 
      </apply>
    </math>
  </component>
  
  <component name="sodium_calcium_exchanger">
    <variable units="picoA" public_interface="out" name="i_NaCa"/>
    
    <variable units="dimensionless" private_interface="out" name="p_E2Na"/>
    <variable units="dimensionless" private_interface="out" name="p_E1Na"/>
    <variable units="dimensionless" private_interface="out" name="p_E1Ca"/>
    <variable units="dimensionless" private_interface="out" name="p_E2Ca"/>
    <variable units="first_order_rate_constant" private_interface="out" name="k1"/>
    <variable units="first_order_rate_constant" private_interface="out" name="k2"/>
    <variable units="first_order_rate_constant" private_interface="out" name="k3" initial_value="1.0"/>
    <variable units="first_order_rate_constant" private_interface="out" name="k4" initial_value="1.0"/>
    <variable units="millimolar" name="Km_Nai" initial_value="8.75"/>
    <variable units="millimolar" name="Km_Nao" initial_value="87.5"/>
    <variable units="millimolar" name="Km_Cai" initial_value="0.00138"/>
    <variable units="millimolar" name="Km_Cao" initial_value="1.38"/>
    
    <variable units="millimolar" public_interface="in" name="Nai"/>
    <variable units="millimolar" public_interface="in" name="Nao"/>
    <variable units="millimolar" public_interface="in" name="Cai"/>  
    <variable units="millimolar" public_interface="in" name="Cao"/>
    <variable units="coulomb_millivolt_per_kelvin_millimole" public_interface="in" name="R"/>
    <variable units="coulomb_per_millimole" public_interface="in" name="F"/>
    <variable units="kelvin" public_interface="in" name="T"/>
    <variable units="millivolt" public_interface="in" name="Vm"/>
    <variable units="picoF" public_interface="in" name="Cm"/>
    <variable units="millisecond" public_interface="in" private_interface="out" name="time"/>
    
    <variable units="dimensionless" private_interface="in" name="y"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="i_NaCa_calculation">
        <eq/>
        <ci> i_NaCa </ci>
        <apply>
          <times/>
          <cn cellml:units="dimensionless"> 6.81 </cn>
          <ci> Cm </ci>
          <apply>
            <minus/>
            <apply>
              <times/>
              <ci> k1 </ci>
              <ci> p_E1Na </ci>
              <ci> y </ci>
            </apply>
            <apply>
              <times/>
              <ci> k2 </ci>
              <ci> p_E2Na </ci>
              <apply>
                <minus/>
                <cn cellml:units="dimensionless"> 1.0 </cn>
                <ci> y </ci>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply id="sodium_calcium_exchanger_p_E1Na_calculation">
        <eq/>
        <ci> p_E1Na </ci>
        <apply>
          <divide/>
          <cn cellml:units="dimensionless"> 1.0 </cn>
          <apply>
            <plus/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
              <times/>
              <apply>
                <power/>
                <apply>
                  <divide/>
                  <ci> Km_Nai </ci>
                  <ci> Nai </ci>
                </apply>
                <cn cellml:units="dimensionless"> 3.0 </cn>
              </apply>
              <apply>
                <plus/>
                <cn cellml:units="dimensionless"> 1.0 </cn>
                <apply>
                  <divide/>
                  <ci> Cai </ci>
                  <ci> Km_Cai </ci>
                </apply>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply> 
      
      <apply id="sodium_calcium_exchanger_p_E2Na_calculation">
        <eq/>
        <ci> p_E2Na </ci>
        <apply>
          <divide/>
          <cn cellml:units="dimensionless"> 1.0 </cn>
          <apply>
            <plus/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
              <times/>
              <apply>
                <power/>
                <apply>
                  <divide/>
                  <ci> Km_Nao </ci>
                  <ci> Nao </ci>
                </apply>
                <cn cellml:units="dimensionless"> 3.0 </cn>
              </apply>
              <apply>
                <plus/>
                <cn cellml:units="dimensionless"> 1.0 </cn>
                <apply>
                  <divide/>
                  <ci> Cao </ci>
                  <ci> Km_Cao </ci>
                </apply>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply id="sodium_calcium_exchanger_p_E1Ca_calculation">
        <eq/>
        <ci> p_E1Ca </ci>
        <apply>
          <divide/>
          <cn cellml:units="dimensionless"> 1.0 </cn>
          <apply>
            <plus/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
              <times/>
              <apply>
                <divide/>
                <ci> Km_Cai </ci>
                <ci> Cai </ci>
              </apply>
              <apply>
                <plus/>
                <cn cellml:units="dimensionless"> 1.0 </cn>
                <apply>
                  <power/>
                  <apply>
                    <divide/>
                    <ci> Nai </ci>
                    <ci> Km_Nai </ci>
                  </apply>
                  <cn cellml:units="dimensionless"> 3.0 </cn>
                </apply>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>  
      
      <apply id="sodium_calcium_exchanger_p_E2Ca_calculation">
        <eq/>
        <ci> p_E2Ca </ci>
        <apply>
          <divide/>
          <cn cellml:units="dimensionless"> 1.0 </cn>
          <apply>
            <plus/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
              <times/>
              <apply>
                <divide/>
                <ci> Km_Cao </ci>
                <ci> Cao </ci>
              </apply>
              <apply>
                <plus/>
                <cn cellml:units="dimensionless"> 1.0 </cn>
                <apply>
                  <power/>
                  <apply>
                    <divide/>
                    <ci> Nao </ci>
                    <ci> Km_Nao </ci>
                  </apply>
                  <cn cellml:units="dimensionless"> 3.0 </cn>
                </apply>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply id="sodium_calcium_exchanger_k1_calculation">
        <eq/>
        <ci> k1 </ci>
        <apply>
          <exp/>
          <apply>
            <times/>
            <cn cellml:units="dimensionless"> 0.32 </cn>
            <ci> F </ci>
            <apply>
              <divide/>
              <ci> Vm </ci>
              <apply>
                <times/>
                <ci> R </ci>
                <ci> T </ci>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply> 
      
      <apply id="sodium_calcium_exchanger_k2_calculation">
        <eq/>
        <ci> k2 </ci>
        <apply>
          <exp/>
          <apply>
            <times/>
            <apply>
              <minus/>
              <cn cellml:units="dimensionless"> 0.32 </cn>
              <cn cellml:units="dimensionless"> 1.0 </cn>
            </apply>
            <ci> F </ci>
            <apply>
              <divide/>
              <ci> Vm </ci>
              <apply>
                <times/>
                <ci> R </ci>
                <ci> T </ci>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
    </math>
  </component> 
  
  <component name="sodium_calcium_exchanger_y_gate">
    <variable units="dimensionless" public_interface="out" name="y" initial_value="9.901E-1"/> 

    <variable units="first_order_rate_constant" name="alpha_y"/>
    <variable units="first_order_rate_constant" name="beta_y"/>
    
    <variable units="millisecond" public_interface="in" name="time"/>
    <variable units="dimensionless" public_interface="in" name="p_E1Na"/>
    <variable units="dimensionless" public_interface="in" name="p_E2Na"/>
    <variable units="dimensionless" public_interface="in" name="p_E1Ca"/>
    <variable units="dimensionless" public_interface="in" name="p_E2Ca"/>
    <variable units="first_order_rate_constant" public_interface="in" name="k1"/>
    <variable units="first_order_rate_constant" public_interface="in" name="k2"/>
    <variable units="first_order_rate_constant" public_interface="in" name="k3"/>
    <variable units="first_order_rate_constant" public_interface="in" name="k4"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="sodium_calcium_exchanger_y_gate_alpha_y_calculation">
        <eq/>
        <ci> alpha_y </ci>
        <apply>
          <plus/>
          <apply>
            <times/>
            <ci> k2 </ci>
            <ci> p_E2Na </ci>
          </apply>
          <apply>
            <times/>
            <ci> k4 </ci>
            <ci> p_E2Ca </ci>
          </apply>
        </apply>
      </apply>
    
      <apply id="sodium_calcium_exchanger_y_gate_beta_y_calculation">
        <eq/>
        <ci> beta_y </ci>
        <apply>
          <plus/>
          <apply>
            <times/>
            <ci> k1 </ci>
            <ci> p_E1Na </ci>
          </apply>
          <apply>
            <times/>
            <ci> k3 </ci>
            <ci> p_E1Ca </ci>
          </apply>
        </apply>
      </apply>
    
      <apply id="sodium_calcium_exchanger_y_gate_dy_dt">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> y </ci>
        </apply>
        <apply>
          <minus/>
          <apply>
            <times/>
            <ci> alpha_y </ci>
            <apply>
              <minus/>
              <cn cellml:units="dimensionless"> 1.0 </cn>
              <ci> y </ci>
            </apply>
          </apply>
          <apply>
            <times/>
            <ci> beta_y </ci>
            <ci> y </ci>
          </apply>
        </apply>
      </apply>
    </math>
  </component>  
             
  <component name="sodium_potassium_pump">
    <variable units="picoA" public_interface="out" name="i_NaK"/>
    
    <variable units="dimensionless" private_interface="out" name="p_E2Na"/>
    <variable units="dimensionless" private_interface="out" name="p_E1Na"/>
    <variable units="dimensionless" private_interface="out" name="p_E1K"/>
    <variable units="dimensionless" private_interface="out" name="p_E2K"/>
    <variable units="first_order_rate_constant" private_interface="out" name="k1"/>
    <variable units="first_order_rate_constant" private_interface="out" name="k2" initial_value="0.04"/>
    <variable units="first_order_rate_constant" private_interface="out" name="k3" initial_value="0.01"/>
    <variable units="first_order_rate_constant" private_interface="out" name="k4" initial_value="0.165"/>
    <variable units="millimolar" name="Km_Nai" initial_value="4.05"/>
    <variable units="millimolar" name="Km_Nao" initial_value="69.8"/>
    <variable units="millimolar" name="Km_Ki" initial_value="32.88"/>
    <variable units="millimolar" name="Km_Ko" initial_value="0.258"/>
    <variable units="millimolar" name="Km_ATP" initial_value="0.094"/>
    <variable units="millimolar" name="Naeff"/>
    
    <variable units="millimolar" public_interface="in" name="ATPi"/>
    <variable units="millimolar" public_interface="in" name="Nai"/>
    <variable units="millimolar" public_interface="in" name="Nao"/>
    <variable units="millimolar" public_interface="in" name="Ki"/>  
    <variable units="millimolar" public_interface="in" name="Ko"/>
    <variable units="coulomb_millivolt_per_kelvin_millimole" public_interface="in" name="R"/>
    <variable units="coulomb_per_millimole" public_interface="in" name="F"/>
    <variable units="kelvin" public_interface="in" name="T"/>
    <variable units="millivolt" public_interface="in" name="Vm"/>
    <variable units="picoF" public_interface="in" name="Cm"/>
    <variable units="millisecond" public_interface="in" private_interface="out" name="time"/>
              
    <variable units="dimensionless" private_interface="in" name="y"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="i_NaK_calculation">
        <eq/>
        <ci> i_NaK </ci>
        <apply>
          <times/>
          <cn cellml:units="dimensionless"> 21.0 </cn>
          <ci> Cm </ci>
          <apply>
            <minus/>
            <apply>
              <times/>
              <ci> k1 </ci>
              <ci> p_E1Na </ci>
              <ci> y </ci>
            </apply>
            <apply>
              <times/>
              <ci> k2 </ci>
              <ci> p_E2Na </ci>
              <apply>
                <minus/>
                <cn cellml:units="dimensionless"> 1.0 </cn>
                <ci> y </ci>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply id="sodium_potassium_pump_p_E1Na_calculation">
        <eq/>
        <ci> p_E1Na </ci>
        <apply>
          <divide/>
          <cn cellml:units="dimensionless"> 1.0 </cn>
          <apply>
            <plus/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
              <times/>
              <apply>
                <power/>
                <apply>
                  <divide/>
                  <ci> Km_Nai </ci>
                  <ci> Nai </ci>
                </apply>
                <cn cellml:units="dimensionless"> 1.06 </cn>
              </apply>
              <apply>
                <plus/>
                <cn cellml:units="dimensionless"> 1.0 </cn>
                <apply>
                  <power/>
                  <apply>
                    <divide/>
                    <ci> Ki </ci>
                    <ci> Km_Ki </ci>
                  </apply>
                  <cn cellml:units="dimensionless"> 1.12 </cn>
                </apply>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply> 
      
      <apply id="sodium_potassium_pump_p_E2Na_calculation">
        <eq/>
        <ci> p_E2Na </ci>
        <apply>
          <divide/>
          <cn cellml:units="dimensionless"> 1.0 </cn>
          <apply>
            <plus/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
              <times/>
              <apply>
                <power/>
                <apply>
                  <divide/>
                  <ci> Km_Nao </ci>
                  <ci> Naeff </ci>
                </apply>
                <cn cellml:units="dimensionless"> 1.06 </cn>
              </apply>
              <apply>
                <plus/>
                <cn cellml:units="dimensionless"> 1.0 </cn>
                <apply>
                  <power/>
                  <apply>
                    <divide/>
                    <ci> Ko </ci>
                    <ci> Km_Ko </ci>
                  </apply>
                  <cn cellml:units="dimensionless"> 1.12 </cn>
                </apply>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply id="p_E1K_calculation">
        <eq/>
        <ci> p_E1K </ci>
        <apply>
          <divide/>
          <cn cellml:units="dimensionless"> 1.0 </cn>
          <apply>
            <plus/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
              <times/>
              <apply>
                <power/>
                <apply>
                  <divide/>
                  <ci> Km_Ki </ci>
                  <ci> Ki </ci>
                </apply>
                <cn cellml:units="dimensionless"> 1.12 </cn>
              </apply>
              <apply>
                <plus/>
                <cn cellml:units="dimensionless"> 1.0 </cn>
                <apply>
                  <power/>
                  <apply>
                    <divide/>
                    <ci> Nai </ci>
                    <ci> Km_Nai </ci>
                  </apply>
                  <cn cellml:units="dimensionless"> 1.06 </cn>
                </apply>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>  
      
      <apply id="p_E2K_calculation">
        <eq/>
        <ci> p_E2K </ci>
        <apply>
          <divide/>
          <cn cellml:units="dimensionless"> 1.0 </cn>
          <apply>
            <plus/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
              <times/>
              <apply>
                <power/>
                <apply>
                  <divide/>
                  <ci> Km_Ko </ci>
                  <ci> Ko </ci>
                </apply>
                <cn cellml:units="dimensionless"> 1.12 </cn>
              </apply>
              <apply>
                <plus/>
                <cn cellml:units="dimensionless"> 1.0 </cn>
                <apply>
                  <power/>
                  <apply>
                    <divide/>
                    <ci> Naeff </ci>
                    <ci> Km_Nao </ci>
                  </apply>
                  <cn cellml:units="dimensionless"> 1.06 </cn>
                </apply>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>  
      
      <apply id="sodium_potassium_pump_k1_calculation">
        <eq/>
        <ci> k1 </ci>
        <apply>
          <times/>
          <cn cellml:units="first_order_rate_constant"> 0.37 </cn>
          <apply>
            <divide/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
              <plus/>
              <cn cellml:units="millimolar"> 1.0 </cn>
              <apply>
                <divide/>
                <ci> Km_ATP </ci>
                <ci> ATPi </ci>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply> 
      
      <apply id="Naeff_calculation">
        <eq/>
        <ci> Naeff </ci>
        <apply>
          <times/>
          <ci> Nao </ci>
          <apply>
            <exp/>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> -0.82 </cn>
              <ci> F </ci>
              <apply>
                <divide/>
                <ci> Vm </ci>
                <apply>
                  <times/>
                  <ci> R </ci>
                  <ci> T </ci>
                </apply>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
    </math>
  </component> 
  
  <component name="sodium_potassium_pump_y_gate">
    <variable units="dimensionless" public_interface="out" name="y" initial_value="5.984E-1"/> 

    <variable units="first_order_rate_constant" name="alpha_y"/>
    <variable units="first_order_rate_constant" name="beta_y"/>
    
    <variable units="millisecond" public_interface="in" name="time"/>
    <variable units="dimensionless" public_interface="in" name="p_E1Na"/>
    <variable units="dimensionless" public_interface="in" name="p_E2Na"/>
    <variable units="dimensionless" public_interface="in" name="p_E1K"/>
    <variable units="dimensionless" public_interface="in" name="p_E2K"/>
    <variable units="first_order_rate_constant" public_interface="in" name="k1"/>
    <variable units="first_order_rate_constant" public_interface="in" name="k2"/>
    <variable units="first_order_rate_constant" public_interface="in" name="k3"/>
    <variable units="first_order_rate_constant" public_interface="in" name="k4"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="sodium_potassium_pump_y_gate_alpha_y_calculation">
        <eq/>
        <ci> alpha_y </ci>
        <apply>
          <plus/>
          <apply>
            <times/>
            <ci> k2 </ci>
            <ci> p_E2Na </ci>
          </apply>
          <apply>
            <times/>
            <ci> k4 </ci>
            <ci> p_E2K </ci>
          </apply>
        </apply>
      </apply>
    
      <apply id="sodium_potassium_pump_y_gate_beta_y_calculation">
        <eq/>
        <ci> beta_y </ci>
        <apply>
          <plus/>
          <apply>
            <times/>
            <ci> k1 </ci>
            <ci> p_E1Na </ci>
          </apply>
          <apply>
            <times/>
            <ci> k3 </ci>
            <ci> p_E1K </ci>
          </apply>
        </apply>
      </apply>
    
      <apply id="sodium_potassium_pump_y_gate_dy_dt">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> y </ci>
        </apply>
        <apply>
          <minus/>
          <apply>
            <times/>
            <ci> alpha_y </ci>
            <apply>
              <minus/>
              <cn cellml:units="dimensionless"> 1.0 </cn>
              <ci> y </ci>
            </apply>
          </apply>
          <apply>
            <times/>
            <ci> beta_y </ci>
            <ci> y </ci>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="SR_calcium_pump">
    <variable units="picoA" public_interface="out" name="i_SR_U"/>
    
    <variable units="dimensionless" private_interface="out" name="p_E2Ca"/>
    <variable units="dimensionless" private_interface="out" name="p_E1Ca"/>
    <variable units="dimensionless" private_interface="out" name="p_E1"/>
    <variable units="dimensionless" private_interface="out" name="p_E2"/>
    <variable units="first_order_rate_constant" private_interface="out" name="k1" initial_value="0.01"/>
    <variable units="first_order_rate_constant" private_interface="out" name="k2" initial_value="1.0"/>
    <variable units="first_order_rate_constant" private_interface="out" name="k3"/>
    <variable units="first_order_rate_constant" private_interface="out" name="k4" initial_value="0.01"/>
    <variable units="millimolar" name="Km_Cai" initial_value="0.08"/>
    <variable units="millimolar" name="Km_Cao" initial_value="0.0008"/>
    <variable units="picoA" name="i_max" initial_value="162500.0"/>
    
    <variable units="millimolar" public_interface="in" name="ATPi"/>
    <variable units="millimolar" public_interface="in" name="Cai"/>
    <variable units="millimolar" public_interface="in" name="Caup"/>
    <variable units="millisecond" public_interface="in" private_interface="out" name="time"/>
              
    <variable units="dimensionless" private_interface="in" name="y"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="i_SR_U_calculation">
        <eq/>
        <ci> i_SR_U </ci>
        <apply>
          <times/>
          <ci> i_max </ci>
          <apply>
            <minus/>
            <apply>
              <times/>
              <ci> k2 </ci>
              <ci> p_E2Ca </ci>
              <apply>
                <minus/>
                <cn cellml:units="dimensionless"> 1.0 </cn>
                <ci> y </ci>
              </apply>
            </apply>
            <apply>
              <times/>
              <ci> k1 </ci>
              <ci> p_E1Ca </ci>
              <ci> y </ci>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply id="SR_calcium_pump_p_E1Ca_calculation">
        <eq/>
        <ci> p_E1Ca </ci>
        <apply>
          <divide/>
          <cn cellml:units="dimensionless"> 1.0 </cn>
          <apply>
            <plus/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
              <divide/>
              <ci> Km_Cai </ci>
              <ci> Caup </ci>
            </apply>
          </apply>
        </apply>
      </apply> 
      
      <apply id="SR_calcium_pump_p_E2Ca_calculation">
        <eq/>
        <ci> p_E2Ca </ci>
        <apply>
          <divide/>
          <cn cellml:units="dimensionless"> 1.0 </cn>
          <apply>
            <plus/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
              <divide/>
              <ci> Km_Cao </ci>
              <ci> Cai </ci>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply id="p_E1_calculation">
        <eq/>
        <ci> p_E1 </ci>
        <apply>
          <minus/>
          <cn cellml:units="dimensionless"> 1.0 </cn>
          <ci> p_E1Ca </ci>
        </apply>
      </apply>  
      
      <apply id="p_E2_calculation">
        <eq/>
        <ci> p_E2 </ci>
        <apply>
          <minus/>
          <cn cellml:units="dimensionless"> 1.0 </cn>
          <ci> p_E2Ca </ci>
        </apply>
      </apply>
      
      <apply id="SR_calcium_pump_k3_calculation">
        <eq/>
        <ci> k3 </ci>
        <apply>
          <divide/>
          <cn cellml:units="first_order_rate_constant"> 1.0 </cn>
          <apply>
            <plus/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
              <divide/>
              <cn cellml:units="millimolar"> 0.1 </cn>
              <ci> ATPi </ci>
            </apply>
          </apply>
        </apply>
      </apply> 
    </math>
  </component> 
  
  <component name="SR_calcium_pump_y_gate">
    <variable units="dimensionless" public_interface="out" name="y" initial_value="4.565E-1"/> 

    <variable units="first_order_rate_constant" name="alpha_y"/>
    <variable units="first_order_rate_constant" name="beta_y"/>
    
    <variable units="millisecond" public_interface="in" name="time"/>
    <variable units="dimensionless" public_interface="in" name="p_E1Ca"/>
    <variable units="dimensionless" public_interface="in" name="p_E2Ca"/>
    <variable units="dimensionless" public_interface="in" name="p_E1"/>
    <variable units="dimensionless" public_interface="in" name="p_E2"/>
    <variable units="first_order_rate_constant" public_interface="in" name="k1"/>
    <variable units="first_order_rate_constant" public_interface="in" name="k2"/>
    <variable units="first_order_rate_constant" public_interface="in" name="k3"/>
    <variable units="first_order_rate_constant" public_interface="in" name="k4"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="SR_calcium_pump_y_gate_alpha_y_calculation">
        <eq/>
        <ci> alpha_y </ci>
        <apply>
          <plus/>
          <apply>
            <times/>
            <ci> k2 </ci>
            <ci> p_E2Ca </ci>
          </apply>
          <apply>
            <times/>
            <ci> k4 </ci>
            <ci> p_E2 </ci>
          </apply>
        </apply>
      </apply>
    
      <apply id="SR_calcium_pump_y_gate_beta_y_calculation">
        <eq/>
        <ci> beta_y </ci>
        <apply>
          <plus/>
          <apply>
            <times/>
            <ci> k1 </ci>
            <ci> p_E1Ca </ci>
          </apply>
          <apply>
            <times/>
            <ci> k3 </ci>
            <ci> p_E1 </ci>
          </apply>
        </apply>
      </apply>
    
      <apply id="SR_calcium_pump_y_gate_dy_dt">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> y </ci>
        </apply>
        <apply>
          <minus/>
          <apply>
            <times/>
            <ci> alpha_y </ci>
            <apply>
              <minus/>
              <cn cellml:units="dimensionless"> 1.0 </cn>
              <ci> y </ci>
            </apply>
          </apply>
          <apply>
            <times/>
            <ci> beta_y </ci>
            <ci> y </ci>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="RyR_channel">
    <variable units="picoA" public_interface="out" name="i_RyR"/>
    
    <variable units="picoA_per_millimolar" name="P_RyR" initial_value="62000.0"/>
    <variable units="first_order_rate_constant" name="k1"/>
    <variable units="first_order_rate_constant" name="k2"/>
    <variable units="first_order_rate_constant" name="k3"/>
    <variable units="first_order_rate_constant" name="k4" initial_value="0.000849"/>
    <variable units="dimensionless" name="p_open_RyR"/>
    <variable units="dimensionless" name="p_close_RyR"/>
    <variable units="dimensionless" name="open_RyR" initial_value="3.330E-4"/>
    <variable units="dimensionless" name="close_RyR" initial_value="1.883E-1"/>
    
    <variable units="dimensionless" public_interface="in" name="p_open_CaL"/>
    <variable units="millimolar" public_interface="in" name="Cai"/>
    <variable units="millimolar" public_interface="in" name="Carel"/>
    <variable units="millisecond" public_interface="in" name="time"/>
    <variable units="picoA" public_interface="in" name="iCaL"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="i_RyR_calculation">
        <eq/>
        <ci> i_RyR </ci>
        <apply>
          <times/>
          <ci> P_RyR </ci>
          <apply>
            <minus/>
            <ci> Carel </ci>
            <ci> Cai </ci>
          </apply>
          <ci> p_open_RyR </ci>
        </apply>
      </apply>
      
      <apply id="p_open_RyR_diff_calculation">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> p_open_RyR </ci>
        </apply>
        <apply>
          <minus/>
          <apply>
            <times/>
            <ci> p_close_RyR </ci>
            <ci> k1 </ci>
          </apply>
          <apply>
            <times/>
            <ci> p_open_RyR </ci>
            <ci> k2 </ci>
          </apply>
        </apply>
      </apply>
      
      <apply id="p_close_RyR_diff_calculation">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> p_close_RyR </ci>
        </apply>
        <apply>
          <minus/>
          <apply>
            <times/>
            <ci> k3 </ci>
            <apply>
              <minus/>
              <cn cellml:units="dimensionless"> 1.0 </cn>
              <apply>
                <plus/>
                <ci> p_open_RyR </ci>
                <ci> p_close_RyR </ci>
              </apply>
            </apply>
          </apply>
          <apply>
            <times/>
            <apply>
              <plus/>
              <ci> k1 </ci>
              <ci> k4 </ci>
            </apply>
            <ci> p_close_RyR </ci>
          </apply>
        </apply>
      </apply>
      
      <apply id="RyR_channel_k1_calculation">
        <eq/>
        <ci> k1 </ci>
        <apply>
          <times/>
          <cn cellml:units="first_order_rate_constant"> 280000 </cn>
          <apply>
            <minus/>
            <apply>
              <power/>
              <ci> Cai </ci>
              <cn cellml:units="dimensionless"> 2.0 </cn>
            </apply>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> 150.0 </cn>
              <ci> iCaL </ci>
              <ci> p_open_CaL </ci>
            </apply>
          </apply>
        </apply>
      </apply> 
      
      <apply id="RyR_channel_k2_calculation">
        <eq/>
        <ci> k2 </ci>
        <apply>
          <divide/>
          <cn cellml:units="first_order_rate_constant"> 0.08 </cn>
          <apply>
            <plus/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
              <divide/>
              <cn cellml:units="millimolar"> 0.36 </cn>
              <ci> Carel </ci>
            </apply>
          </apply>
        </apply>
      </apply>
        
      <apply id="RyR_channel_k3_calculation">
        <eq/>
        <ci> k3 </ci>
        <apply>
          <times/>
          <cn cellml:units="first_order_rate_constant"> 0.000377 </cn>
          <apply>
            <power/>
            <ci> Carel </ci>
            <cn cellml:units="dimensionless"> 2.0 </cn>
          </apply>
        </apply>
      </apply> 
    </math>
  </component>
  
  <component name="SR_T_current">
    <variable units="picoA" public_interface="out" name="i_SR_T"/>
    
    <variable units="picoA_per_millimolar" name="P_SR_T" initial_value="386.0"/>
    
    <variable units="millimolar" public_interface="in" name="Carel"/>
    <variable units="millimolar" public_interface="in" name="Caup"/>
    <variable units="millisecond" public_interface="in" name="time"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="i_SR_T_calculation">
        <eq/>
        <ci> i_SR_T </ci>
        <apply>
          <times/>
          <ci> P_SR_T </ci>
          <apply>
            <minus/>
            <ci> Caup </ci>
            <ci> Carel </ci>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="SR_L_current">
    <variable units="picoA" public_interface="out" name="i_SR_L"/>
    
    <variable units="picoA_per_millimolar" name="P_SR_L" initial_value="459.0"/>
    
    <variable units="millimolar" public_interface="in" name="Cai"/>
    <variable units="millimolar" public_interface="in" name="Caup"/>
    <variable units="millisecond" public_interface="in" name="time"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="i_SR_L_calculation">
        <eq/>
        <ci> i_SR_L </ci>
        <apply>
          <times/>
          <ci> P_SR_L </ci>
          <apply>
            <minus/>
            <ci> Caup </ci>
            <ci> Cai </ci>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="Ca_concentrations_in_SR">
    <variable units="millimolar" public_interface="out" name="Carel" initial_value="9.422"/>
    <variable units="millimolar" public_interface="out" name="Caup" initial_value="2.590"/>
    
    <variable units="micrometre3" name="V_rel" initial_value="160.0"/>
    <variable units="micrometre3" name="V_up" initial_value="400.0"/>
    
    <variable units="millisecond" public_interface="in" name="time"/>
    <variable units="coulomb_per_millimole" public_interface="in" name="F"/>
    <variable units="picoA" public_interface="in" name="i_SR_L"/>
    <variable units="picoA" public_interface="in" name="i_SR_T"/>
    <variable units="picoA" public_interface="in" name="i_RyR"/>
    <variable units="picoA" public_interface="in" name="i_SR_U"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="Carel_diff_eq">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> Carel </ci>
        </apply>
        <apply>
          <divide/>
          <apply>
            <minus/>
            <ci> i_SR_T </ci>
            <ci> i_RyR </ci>
          </apply>
          <apply>
            <times/>
            <cn cellml:units="dimensionless"> 2.0 </cn>
            <ci> F </ci>
            <ci> V_rel </ci>
          </apply>
        </apply>
      </apply>
      
      <apply id="Caup_diff_eq">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> Caup </ci>
        </apply>
        <apply>
          <divide/>
          <apply>
            <minus/>
            <ci> i_SR_U </ci>
            <apply>
              <plus/>
              <ci> i_SR_T </ci>
              <ci> i_SR_L </ci>
            </apply>
          </apply>
          <apply>
            <times/>
            <cn cellml:units="dimensionless"> 2.0 </cn>
            <ci> F </ci>
            <ci> V_up </ci>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="Ca_buffers">
    <variable units="millimolar" name="total_calmodulin" initial_value="0.05"/>
    <variable units="millimolar" name="calmodulin_Km" initial_value="0.00238"/>
    <variable units="millimolar" name="total_calsequestrin" initial_value="10.0"/>
    <variable units="millimolar" name="calsequestrin_Km" initial_value="0.8"/> 
    <variable units="millimolar" name="calmodulin_Ca"/>
    <variable units="millimolar" name="calsequestrin_Ca"/>
    <variable units="millimolar" name="free_calmodulin"/>
    <variable units="millimolar" name="free_calsequestrin"/>
    <variable units="millimolar" name="Cai_total"/>
    <variable units="millimolar" name="Ca_rel"/>
    <variable units="millimolar" name="Ca_rel_total"/>
    
    <variable units="millimolar" public_interface="in" name="Cai"/>
    <variable units="millisecond" public_interface="in" name="time"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="calmodulin_Ca_calculation">
        <eq/>
        <ci> calmodulin_Ca </ci>
        <apply>
          <minus/>
          <ci> total_calmodulin </ci>
          <ci> free_calmodulin </ci>
        </apply>
      </apply>
      
      <apply id="calsequestrin_Ca_calculation">
        <eq/>
        <ci> calsequestrin_Ca </ci>
        <apply>
          <minus/>
          <ci> total_calsequestrin </ci>
          <ci> free_calsequestrin </ci>
        </apply>
      </apply>
      
      <apply id="Cai_total_calculation">
        <eq/>
        <ci> Cai_total </ci>
        <apply>
          <plus/>
          <ci> Cai </ci>
          <ci> calmodulin_Ca </ci>
        </apply>
      </apply>
     
     <!-- 
      <apply id="Cai_calculation"><eq />
        <ci> Cai </ci>
        <apply><divide />
          <ci> Cai_total </ci>
          <apply><plus />
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply><divide />
              <ci> free_calmodulin </ci>
              <ci> calmodulin_Km </ci>
            </apply>
          </apply>
        </apply>
      </apply>
      -->
      
      <apply id="free_calsequestrin_calculation">
        <eq/>
        <ci> free_calsequestrin </ci>
        <apply>
          <divide/>
          <ci> total_calsequestrin </ci>
          <apply>
            <plus/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
              <divide/>
              <ci> Ca_rel </ci>
              <ci> calsequestrin_Km </ci>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply id="free_calmodulin_calculation">
        <eq/>
        <ci> free_calmodulin </ci>
        <apply>
          <divide/>
          <ci> total_calmodulin </ci>
          <apply>
            <plus/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
              <divide/>
              <ci> Cai </ci>
              <ci> calmodulin_Km </ci>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply id="Ca_rel_total_calculation">
        <eq/>
        <ci> Ca_rel_total </ci>
        <apply>
          <plus/>
          <ci> Ca_rel </ci>
          <ci> calsequestrin_Ca </ci>
        </apply>
      </apply>
      
      <apply id="Ca_rel_calculation">
        <eq/>
        <ci> Ca_rel </ci>
        <apply>
          <divide/>
          <ci> Ca_rel_total </ci>
          <apply>
            <plus/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
              <divide/>
              <ci> free_calsequestrin </ci>
              <ci> calsequestrin_Km </ci>
            </apply>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="NL_model">
    <variable units="dimensionless" public_interface="out" name="TCa_" initial_value="1.918E-3"/>
    <variable units="dimensionless" public_interface="out" name="TT"/>
    
    <variable units="dimensionless" name="T_NL"/>
    <variable units="dimensionless" name="T_" initial_value="4.077E-4"/>
    <variable units="dimensionless" name="TCa_eff"/>
    <variable units="dimensionless" name="TCa" initial_value="2.444E-2"/>
    <variable units="flux" name="Qa"/>
    <variable units="flux" name="Qb"/>
    <variable units="flux" name="Qr"/>
    <variable units="flux" name="Qd"/>
    <variable units="flux" name="Qd1"/>
    <variable units="flux" name="Qd2"/>
    <variable units="flux" name="Qrel"/>
    <variable units="flux" name="Qm" initial_value="1600.0"/>
    <variable units="flux" name="Qpump"/>
    <variable units="flux" name="Qpump_rest"/>
    <variable units="micromolar" name="Ca_rest" initial_value="0.04"/>
    <variable units="millisecond" name="t1" initial_value="25.0"/>
    <variable units="flux" name="Kp" initial_value="150.0"/>
    <variable units="micromolar" name="Km" initial_value="0.1"/>
    <variable units="flux" name="Y1" initial_value="39.0"/>
    <variable units="first_order_rate_constant" name="Y2" initial_value="0.0039"/>
    <variable units="first_order_rate_constant" name="Y3" initial_value="0.03"/>
    <variable units="first_order_rate_constant" name="Y4" initial_value="0.12"/>
    <variable units="second_per_micrometre2" name="Yd" initial_value="9.0"/>
    <variable units="first_order_rate_constant" name="Z1" initial_value="0.03"/>
    <variable units="first_order_rate_constant" name="Z2" initial_value="0.0039"/>
    <variable units="flux" name="Z3" initial_value="1560.0"/>
    <variable units="micromolar" name="Tt" initial_value="70.0"/>
    <variable units="first_order_rate_constant" name="B" initial_value="1.2"/>
    <variable units="micrometre" name="hc" initial_value="0.005"/>
    <variable units="micrometre" name="h" initial_value="4.996E-3"/>
    <variable units="micrometre" name="La" initial_value="1.17"/>
    <variable units="micrometre" name="L" initial_value="9.627E-1"/>
    <variable units="per_micrometre2" name="R" initial_value="20.0"/>
    <variable units="mN_per_mm2" name="Fb"/>
    <variable units="mN_per_mm2_micrometre_micromolar" name="A" initial_value="1800.0"/>
    <variable units="micrometre" name="X"/>
    
    <variable units="millisecond" public_interface="in" name="time"/>
    <variable units="millimolar" public_interface="in" name="Cai"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply id="dTCa_dt">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> TCa </ci>
        </apply>
        <apply>
          <minus/>
          <ci> Qb </ci>
          <ci> Qa </ci>
        </apply>
      </apply>
      
      <apply id="dTCa__dt">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> TCa_ </ci>
        </apply>
        <apply>
          <minus/>
          <ci> Qa </ci>
          <apply>
            <plus/>
            <ci> Qr </ci>
            <ci> Qd2 </ci>
          </apply>
        </apply>
      </apply>
      
      <apply id="dT__dt">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> T_ </ci>
        </apply>
        <apply>
          <minus/>
          <ci> Qr </ci>
          <apply>
            <plus/>
            <ci> Qd </ci>
            <ci> Qd1 </ci>
          </apply>
        </apply>
      </apply>
      
      <apply id="Qb_calculation">
        <eq/>
        <ci> Qb </ci>
        <apply>
          <times/>
          <ci> Y1 </ci>
          <ci> Cai </ci>
          <apply>
            <minus/>
            <ci> T_NL </ci>
            <ci> Z1 </ci>
          </apply>
          <ci> TCa </ci>
        </apply>
      </apply>
      
      <apply id="T_NL_calculation">
        <eq/>
        <ci> T_NL </ci>
        <apply>
          <minus/>
          <ci> Tt </ci>
          <apply>
            <plus/>
            <ci> TCa </ci>
            <ci> TCa_ </ci>
            <ci> T_ </ci>
          </apply>
        </apply>
      </apply>
      
      <apply id="Qa_calculation">
        <eq/>
        <ci> Qa </ci>
        <apply>
          <times/>
          <ci> Y2 </ci>
          <apply>
            <minus/>
            <ci> TCa_eff </ci>
            <ci> Z2 </ci>
          </apply>
          <ci> TCa_ </ci>
        </apply>
      </apply>
      
      <apply id="TCa_eff_calculation">
        <eq/>
        <ci> TCa_eff </ci>
        <apply>
          <times/>
          <ci> TCa </ci>
          <apply>
            <exp/>
            <apply>
              <times/>
              <apply>
                <minus/>
                <ci> R </ci>
              </apply>
              <apply>
                <power/>
                <apply>
                  <minus/>
                  <ci> L </ci>
                  <ci> La </ci>
                </apply>
                <cn cellml:units="dimensionless"> 2.0 </cn>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply id="Qr_calculation">
        <eq/>
        <ci> Qr </ci>
        <apply>
          <times/>
          <ci> Y3 </ci>
          <apply>
            <minus/>
            <ci> TCa_ </ci>
            <ci> Z3 </ci>
          </apply>
          <ci> T_ </ci>
          <ci> Cai </ci>
        </apply>
      </apply>
      
      <apply id="Qd_calculation">
        <eq/>
        <ci> Qd </ci>
        <apply>
          <times/>
          <ci> Y4 </ci>
          <ci> T_ </ci>
        </apply>
      </apply>
      
      <apply id="Qd1_calculation">
        <eq/>
        <ci> Qd1 </ci>
        <apply>
          <times/>
          <ci> Yd </ci>
          <apply>
            <power/>
            <apply>
              <times/>
              <ci> B </ci>
              <apply>
                <minus/>
                <ci> h </ci>
                <ci> hc </ci>
              </apply>
            </apply>
            <cn cellml:units="dimensionless"> 2.0 </cn>
          </apply>
          <ci> T_ </ci>
        </apply>
      </apply>
      
      <apply id="Qd2_calculation">
        <eq/>
        <ci> Qd2 </ci>
        <apply>
          <times/>
          <ci> Yd </ci>
          <apply>
            <power/>
            <apply>
              <times/>
              <ci> B </ci>
              <apply>
                <minus/>
                <ci> h </ci>
                <ci> hc </ci>
              </apply>
            </apply>
            <cn cellml:units="dimensionless"> 2.0 </cn>
          </apply>
          <ci> TCa_ </ci>
        </apply>
      </apply>
      
      <apply id="Qpump_calculation">
        <eq/>
        <ci> Qpump </ci>
        <apply>
          <divide/>
          <ci> Kp </ci>
          <apply>
            <plus/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
              <power/>
              <apply>
                <divide/>
                <ci> Km </ci>
                <ci> Cai </ci>
              </apply>
              <cn cellml:units="dimensionless"> 2.0 </cn>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply id="Qpump_rest_calculation">
        <eq/>
        <ci> Qpump_rest </ci>
        <apply>
          <divide/>
          <ci> Kp </ci>
          <apply>
            <plus/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
              <power/>
              <apply>
                <divide/>
                <ci> Km </ci>
                <ci> Ca_rest </ci>
              </apply>
              <cn cellml:units="dimensionless"> 2.0 </cn>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply id="Qrel_calculation">
        <eq/>
        <ci> Qrel </ci>
        <apply>
          <plus/>
          <apply>
            <times/>
            <ci> Qm </ci>
            <apply>
              <power/>
              <apply>
                <divide/>
                <ci> time </ci>
                <ci> t1 </ci>
              </apply>
              <cn cellml:units="dimensionless"> 4.0 </cn>
            </apply>
            <apply>
              <exp/>
              <apply>
                <times/>
                <cn cellml:units="dimensionless"> 4.0 </cn>
                <apply>
                  <minus/>
                  <cn cellml:units="dimensionless"> 1.0 </cn>
                  <apply>
                    <divide/>
                    <ci> time </ci>
                    <ci> t1 </ci>
                  </apply>
                </apply>
              </apply>
            </apply>
          </apply>
          <ci> Qpump_rest </ci>
        </apply>
      </apply> 
      
      <apply id="dX_dt">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> X </ci>
        </apply>
        <apply>
          <times/>
          <ci> B </ci>
          <apply>
            <minus/>
            <ci> h </ci>
            <ci> hc </ci>
          </apply>
        </apply>
      </apply>
      
      <apply id="Fb_calculation">
        <eq/>
        <ci> Fb </ci>
        <apply>
          <times/>
          <ci> A </ci>
          <apply>
            <plus/>
            <ci> TCa_ </ci>
            <ci> T_ </ci>
          </apply>
          <apply>
            <minus/>
            <ci> L </ci>
            <ci> X </ci>
          </apply>
          <ci> TCa_ </ci>
        </apply>
      </apply>
      
      <apply id="TT_calculation">
        <eq/>
        <ci> TT </ci>
        <apply>
          <plus/>
          <ci> Qd2 </ci>
          <ci> Qr </ci>
          <apply>
            <minus/>
            <ci> Qb </ci>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  
  <group>
    <relationship_ref relationship="containment"/>
    <component_ref component="membrane">
      <component_ref component="internal_ion_concentrations"/>
      <component_ref component="constant_field_equations"/>
      <component_ref component="ATP_production"/>
      <component_ref component="sodium_current">
        <component_ref component="sodium_current_voltage_dependent_gate"/>
        <component_ref component="sodium_current_ultra_slow_gate"/>
      </component_ref>
      <component_ref component="L_type_Ca_channel">
        <component_ref component="L_type_Ca_channel_voltage_dependent_gate"/>
        <component_ref component="L_type_Ca_channel_Ca_dependent_gate"/>
        <component_ref component="L_type_Ca_channel_ultra_slow_gate"/>
      </component_ref>
      <component_ref component="T_type_Ca_channel">
        <component_ref component="T_type_Ca_channel_y1_gate"/>
        <component_ref component="T_type_Ca_channel_y2_gate"/>
      </component_ref>
      <component_ref component="rapid_time_dependent_potassium_current">
        <component_ref component="rapid_time_dependent_potassium_current_y1_gate"/>
        <component_ref component="rapid_time_dependent_potassium_current_y2_gate"/>
        <component_ref component="rapid_time_dependent_potassium_current_y3_gate"/>
      </component_ref>
      <component_ref component="slow_time_dependent_potassium_current">
        <component_ref component="slow_time_dependent_potassium_current_y1_gate"/>
        <component_ref component="slow_time_dependent_potassium_current_y2_gate"/>
      </component_ref>
      <component_ref component="transient_outward_current">
        <component_ref component="transient_outward_current_y1_gate"/>
        <component_ref component="transient_outward_current_y2_gate"/>
      </component_ref>
      <component_ref component="background_NSC_current"/>
      <component_ref component="background_Kpl_current"/>
      <component_ref component="background_lCa_current"/>
      <component_ref component="background_Cab_current"/>
      <component_ref component="sodium_calcium_exchanger">
        <component_ref component="sodium_calcium_exchanger_y_gate"/>             
      </component_ref>
      <component_ref component="sodium_potassium_pump">
        <component_ref component="sodium_potassium_pump_y_gate"/>
      </component_ref>
      <component_ref component="SR_calcium_pump">
        <component_ref component="SR_calcium_pump_y_gate"/>
      </component_ref>
      <component_ref component="Ca_buffers"/>
      <component_ref component="RyR_channel"/>
      <component_ref component="SR_T_current"/>
      <component_ref component="SR_L_current"/>
      <component_ref component="Ca_concentrations_in_SR"/>
    </component_ref>
  </group>
  
  <group>
    <relationship_ref relationship="encapsulation"/>
    <component_ref component="sodium_current">
      <component_ref component="sodium_current_voltage_dependent_gate"/>
      <component_ref component="sodium_current_ultra_slow_gate"/>
    </component_ref>
    <component_ref component="L_type_Ca_channel">
      <component_ref component="L_type_Ca_channel_voltage_dependent_gate"/>
      <component_ref component="L_type_Ca_channel_Ca_dependent_gate"/>
      <component_ref component="L_type_Ca_channel_ultra_slow_gate"/>
    </component_ref>
    <component_ref component="T_type_Ca_channel">
      <component_ref component="T_type_Ca_channel_y1_gate"/>
      <component_ref component="T_type_Ca_channel_y2_gate"/>
    </component_ref>
    <component_ref component="time_independent_potassium_current">
      <component_ref component="time_independent_potassium_current_y_gate"/>
    </component_ref>
    <component_ref component="rapid_time_dependent_potassium_current">
      <component_ref component="rapid_time_dependent_potassium_current_y1_gate"/>
      <component_ref component="rapid_time_dependent_potassium_current_y2_gate"/>
      <component_ref component="rapid_time_dependent_potassium_current_y3_gate"/>
    </component_ref>
    <component_ref component="slow_time_dependent_potassium_current">
      <component_ref component="slow_time_dependent_potassium_current_y1_gate"/>
      <component_ref component="slow_time_dependent_potassium_current_y2_gate"/>
    </component_ref>
    <component_ref component="transient_outward_current">
      <component_ref component="transient_outward_current_y1_gate"/>
      <component_ref component="transient_outward_current_y2_gate"/>
    </component_ref>
    <component_ref component="sodium_calcium_exchanger">
      <component_ref component="sodium_calcium_exchanger_y_gate"/>             
    </component_ref>
    <component_ref component="sodium_potassium_pump">
      <component_ref component="sodium_potassium_pump_y_gate"/>
    </component_ref>
    <component_ref component="SR_calcium_pump">
      <component_ref component="SR_calcium_pump_y_gate"/>
    </component_ref>
  </group>
  
  <connection>
    <map_components component_2="environment" component_1="membrane"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="environment" component_1="internal_ion_concentrations"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="environment" component_1="ATP_production"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>

  <connection>
    <map_components component_2="environment" component_1="sodium_current"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="environment" component_1="L_type_Ca_channel"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>

  <connection>
    <map_components component_2="environment" component_1="T_type_Ca_channel"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="environment" component_1="SR_T_current"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="environment" component_1="SR_L_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="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="transient_outward_current"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="environment" component_1="background_NSC_current"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="environment" component_1="background_Kpl_current"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="environment" component_1="background_lCa_current"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="environment" component_1="background_KATP_current"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="environment" component_1="background_Cab_current"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
   
  <connection>
    <map_components component_2="environment" component_1="sodium_calcium_exchanger"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="environment" component_1="sodium_potassium_pump"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="environment" component_1="SR_calcium_pump"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="environment" component_1="RyR_channel"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="environment" component_1="Ca_concentrations_in_SR"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="environment" component_1="Ca_buffers"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="environment" component_1="NL_model"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="internal_ion_concentrations" component_1="membrane"/>
    <map_variables variable_2="F" variable_1="F"/>
    <map_variables variable_2="i_ext" variable_1="i_ext"/>
  </connection>
  
  <connection>
    <map_components component_2="constant_field_equations" component_1="membrane"/>
    <map_variables variable_2="Vm" variable_1="Vm"/>
    <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="ATP_production" component_1="membrane"/>
    <map_variables variable_2="Vm" variable_1="Vm"/>
    <map_variables variable_2="F" variable_1="F"/>
  </connection>
  
  <connection>
    <map_components component_2="sodium_current" component_1="membrane"/>
    <map_variables variable_2="Vm" variable_1="Vm"/>
    <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="Vm" variable_1="Vm"/>
    <map_variables variable_2="i_Ca_L" variable_1="i_Ca_L"/>
  </connection>
  
  <connection>
    <map_components component_2="T_type_Ca_channel" component_1="membrane"/>
    <map_variables variable_2="Vm" variable_1="Vm"/>
    <map_variables variable_2="i_Ca_T" variable_1="i_Ca_T"/>
  </connection>
  
  <connection>
    <map_components component_2="time_independent_potassium_current" component_1="membrane"/>
    <map_variables variable_2="Vm" variable_1="Vm"/>
    <map_variables variable_2="R" variable_1="R"/>
    <map_variables variable_2="T" variable_1="T"/>
    <map_variables variable_2="F" variable_1="F"/>
    <map_variables variable_2="i_K1" variable_1="i_K1"/>
  </connection>
  
  <connection>
    <map_components component_2="rapid_time_dependent_potassium_current" component_1="membrane"/>
    <map_variables variable_2="Vm" variable_1="Vm"/>
    <map_variables variable_2="Cm" variable_1="Cm"/>
    <map_variables variable_2="i_Kr" variable_1="i_Kr"/>
  </connection>
  
  <connection>
    <map_components component_2="slow_time_dependent_potassium_current" component_1="membrane"/>
    <map_variables variable_2="Vm" variable_1="Vm"/>
    <map_variables variable_2="i_Ks" variable_1="i_Ks"/>
  </connection>
  
  <connection>
    <map_components component_2="transient_outward_current" component_1="membrane"/>
    <map_variables variable_2="Vm" variable_1="Vm"/>
    <map_variables variable_2="i_to" variable_1="i_to"/>
  </connection>
  
  <connection>
    <map_components component_2="background_NSC_current" component_1="membrane"/>
    <map_variables variable_2="i_bNSC" variable_1="i_bNSC"/>
  </connection>
  
  <connection>
    <map_components component_2="background_Kpl_current" component_1="membrane"/>
    <map_variables variable_2="i_Kpl" variable_1="i_Kpl"/>
    <map_variables variable_2="Vm" variable_1="Vm"/>
  </connection>
  
  <connection>
    <map_components component_2="background_lCa_current" component_1="membrane"/>
    <map_variables variable_2="i_lCa" variable_1="i_lCa"/>
  </connection>
  
  <connection>
    <map_components component_2="background_KATP_current" component_1="membrane"/>
    <map_variables variable_2="i_KATP" variable_1="i_KATP"/>
    <map_variables variable_2="Vm" variable_1="Vm"/>
  </connection>
  
  <connection>
    <map_components component_2="background_Cab_current" component_1="membrane"/>
    <map_variables variable_2="i_Cab" variable_1="i_Cab"/>
  </connection>
  
  <connection>
    <map_components component_2="sodium_calcium_exchanger" component_1="membrane"/>
    <map_variables variable_2="Vm" variable_1="Vm"/>
    <map_variables variable_2="Cm" variable_1="Cm"/>
    <map_variables variable_2="R" variable_1="R"/>
    <map_variables variable_2="T" variable_1="T"/>
    <map_variables variable_2="F" variable_1="F"/>
    <map_variables variable_2="i_NaCa" variable_1="i_NaCa"/>
  </connection>
  
  <connection>
    <map_components component_2="sodium_potassium_pump" component_1="membrane"/>
    <map_variables variable_2="Vm" variable_1="Vm"/>
    <map_variables variable_2="Cm" variable_1="Cm"/>
    <map_variables variable_2="R" variable_1="R"/>
    <map_variables variable_2="T" variable_1="T"/>
    <map_variables variable_2="F" variable_1="F"/>
    <map_variables variable_2="i_NaK" variable_1="i_NaK"/>
  </connection>
  
  <connection>
    <map_components component_2="Ca_concentrations_in_SR" component_1="membrane"/>
    <map_variables variable_2="F" variable_1="F"/>
  </connection>
  
  <connection>
    <map_components component_2="constant_field_equations" component_1="internal_ion_concentrations"/>
    <map_variables variable_2="Nai" variable_1="Nai"/>
    <map_variables variable_2="Ki" variable_1="Ki"/>
    <map_variables variable_2="Cai" variable_1="Cai"/>
  </connection>
  
  <connection>
    <map_components component_2="NL_model" component_1="internal_ion_concentrations"/>
    <map_variables variable_2="Cai" variable_1="Cai"/>
    <map_variables variable_2="TT" variable_1="TT"/>
  </connection>
  
  <connection>
    <map_components component_2="constant_field_equations" component_1="external_ion_concentrations"/>
    <map_variables variable_2="Nao" variable_1="Nao"/>
    <map_variables variable_2="Ko" variable_1="Ko"/>
    <map_variables variable_2="Cao" variable_1="Cao"/>
  </connection>
  
  <connection>
    <map_components component_2="time_independent_potassium_current" component_1="external_ion_concentrations"/>
    <map_variables variable_2="Ko" variable_1="Ko"/>
  </connection>
  
  <connection>
    <map_components component_2="rapid_time_dependent_potassium_current" component_1="external_ion_concentrations"/>
    <map_variables variable_2="Ko" variable_1="Ko"/>
  </connection>
  
  <connection>
    <map_components component_2="background_Kpl_current" component_1="external_ion_concentrations"/>
    <map_variables variable_2="Ko" variable_1="Ko"/>
  </connection>
  
  <connection>
    <map_components component_2="background_KATP_current" component_1="external_ion_concentrations"/>
    <map_variables variable_2="Ko" variable_1="Ko"/>
  </connection>
  
  <connection>
    <map_components component_2="background_lCa_current" component_1="internal_ion_concentrations"/>
    <map_variables variable_2="Cai" variable_1="Cai"/>
    <map_variables variable_2="i_lCa_Na" variable_1="i_lCa_Na"/>
    <map_variables variable_2="i_lCa_K" variable_1="i_lCa_K"/>
  </connection>
  
  <connection>
    <map_components component_2="sodium_calcium_exchanger" component_1="internal_ion_concentrations"/>
    <map_variables variable_2="Nai" variable_1="Nai"/>
    <map_variables variable_2="Cai" variable_1="Cai"/>
    <map_variables variable_2="i_NaCa" variable_1="i_NaCa"/>
  </connection>
  
  <connection>
    <map_components component_2="sodium_calcium_exchanger" component_1="external_ion_concentrations"/>
    <map_variables variable_2="Nao" variable_1="Nao"/>
    <map_variables variable_2="Cao" variable_1="Cao"/>
  </connection>
  
  <connection>
    <map_components component_2="sodium_potassium_pump" component_1="internal_ion_concentrations"/>
    <map_variables variable_2="Nai" variable_1="Nai"/>
    <map_variables variable_2="Ki" variable_1="Ki"/>
    <map_variables variable_2="i_NaK" variable_1="i_NaK"/>
  </connection>
  
  <connection>
    <map_components component_2="sodium_potassium_pump" component_1="external_ion_concentrations"/>
    <map_variables variable_2="Nao" variable_1="Nao"/>
    <map_variables variable_2="Ko" variable_1="Ko"/>
  </connection>
  
  <connection>
    <map_components component_2="SR_calcium_pump" component_1="internal_ion_concentrations"/>
    <map_variables variable_2="Cai" variable_1="Cai"/>
    <map_variables variable_2="i_SR_U" variable_1="i_SR_U"/>
  </connection>
  
  <connection>
    <map_components component_2="Ca_buffers" component_1="internal_ion_concentrations"/>
    <map_variables variable_2="Cai" variable_1="Cai"/>
  </connection>
  
  <connection>
    <map_components component_2="RyR_channel" component_1="internal_ion_concentrations"/>
    <map_variables variable_2="i_RyR" variable_1="i_RyR"/>
    <map_variables variable_2="Cai" variable_1="Cai"/>
  </connection>
  
  <connection>
    <map_components component_2="SR_L_current" component_1="internal_ion_concentrations"/>
    <map_variables variable_2="Cai" variable_1="Cai"/>
    <map_variables variable_2="i_SR_L" variable_1="i_SR_L"/>
  </connection>
  
  <connection>
    <map_components component_2="sodium_current" component_1="internal_ion_concentrations"/>
    <map_variables variable_2="i_Na_Na" variable_1="i_Na_Na"/>
    <map_variables variable_2="i_Na_K" variable_1="i_Na_K"/>
  </connection>
  
  <connection>
    <map_components component_2="L_type_Ca_channel" component_1="internal_ion_concentrations"/>
    <map_variables variable_2="i_CaL_Na" variable_1="i_CaL_Na"/>
    <map_variables variable_2="i_CaL_Ca" variable_1="i_CaL_Ca"/>
    <map_variables variable_2="i_CaL_K" variable_1="i_CaL_K"/>
    <map_variables variable_2="Cai" variable_1="Cai"/>
  </connection>
  
  <connection>
    <map_components component_2="T_type_Ca_channel" component_1="internal_ion_concentrations"/>
    <map_variables variable_2="i_Ca_T" variable_1="i_Ca_T"/>
  </connection>
  
  <connection>
    <map_components component_2="time_independent_potassium_current" component_1="internal_ion_concentrations"/>
    <map_variables variable_2="i_K1" variable_1="i_K1"/>
    <map_variables variable_2="Ki" variable_1="Ki"/>
  </connection>
  
  <connection>
    <map_components component_2="rapid_time_dependent_potassium_current" component_1="internal_ion_concentrations"/>
    <map_variables variable_2="i_Kr" variable_1="i_Kr"/>
  </connection>
  
  <connection>
    <map_components component_2="slow_time_dependent_potassium_current" component_1="internal_ion_concentrations"/>
    <map_variables variable_2="i_Ks_Na" variable_1="i_Ks_Na"/>
    <map_variables variable_2="i_Ks_K" variable_1="i_Ks_K"/>
    <map_variables variable_2="Cai" variable_1="Cai"/>
  </connection>
  
  <connection>
    <map_components component_2="transient_outward_current" component_1="internal_ion_concentrations"/>
    <map_variables variable_2="i_to_Na" variable_1="i_to_Na"/>
    <map_variables variable_2="i_to_K" variable_1="i_to_K"/>
  </connection>
  
  <connection>
    <map_components component_2="background_NSC_current" component_1="internal_ion_concentrations"/>
    <map_variables variable_2="i_bNSC" variable_1="i_bNSC"/>
    <map_variables variable_2="i_bNSC_Na" variable_1="i_bNSC_Na"/>
    <map_variables variable_2="i_bNSC_K" variable_1="i_bNSC_K"/>
  </connection>
  
  <connection>
    <map_components component_2="background_Kpl_current" component_1="internal_ion_concentrations"/>
    <map_variables variable_2="i_Kpl" variable_1="i_Kpl"/>
  </connection>
  
  <connection>
    <map_components component_2="background_KATP_current" component_1="internal_ion_concentrations"/>
    <map_variables variable_2="i_KATP" variable_1="i_KATP"/>
  </connection>
  
  <connection>
    <map_components component_2="background_Cab_current" component_1="internal_ion_concentrations"/>
    <map_variables variable_2="i_Cab" variable_1="i_Cab"/>
  </connection>
  
  <connection>
    <map_components component_2="sodium_current" component_1="constant_field_equations"/>
    <map_variables variable_2="CF_Na" variable_1="CF_Na"/>
    <map_variables variable_2="CF_K" variable_1="CF_K"/>
  </connection>
  
  <connection>
    <map_components component_2="L_type_Ca_channel" component_1="constant_field_equations"/>
    <map_variables variable_2="CF_Ca" variable_1="CF_Ca"/>
    <map_variables variable_2="CF_Na" variable_1="CF_Na"/>
    <map_variables variable_2="CF_K" variable_1="CF_K"/>
  </connection>
  
  <connection>
    <map_components component_2="T_type_Ca_channel" component_1="constant_field_equations"/>
    <map_variables variable_2="CF_Ca" variable_1="CF_Ca"/>
  </connection>
  
  <connection>
    <map_components component_2="slow_time_dependent_potassium_current" component_1="constant_field_equations"/>
    <map_variables variable_2="CF_Na" variable_1="CF_Na"/>
    <map_variables variable_2="CF_K" variable_1="CF_K"/>
  </connection>
  
  <connection>
    <map_components component_2="transient_outward_current" component_1="constant_field_equations"/>
    <map_variables variable_2="CF_Na" variable_1="CF_Na"/>
    <map_variables variable_2="CF_K" variable_1="CF_K"/>
  </connection>
  
  <connection>
    <map_components component_2="background_NSC_current" component_1="constant_field_equations"/>
    <map_variables variable_2="CF_Na" variable_1="CF_Na"/>
    <map_variables variable_2="CF_K" variable_1="CF_K"/>
  </connection>
  
  <connection>
    <map_components component_2="background_Kpl_current" component_1="constant_field_equations"/>
    <map_variables variable_2="CF_K" variable_1="CF_K"/>
  </connection>
  
  <connection>
    <map_components component_2="background_lCa_current" component_1="constant_field_equations"/>
    <map_variables variable_2="CF_Na" variable_1="CF_Na"/>
    <map_variables variable_2="CF_K" variable_1="CF_K"/>
  </connection>
  
  <connection>
    <map_components component_2="background_Cab_current" component_1="constant_field_equations"/>
    <map_variables variable_2="CF_Ca" variable_1="CF_Ca"/>
  </connection>
  
  <connection>
    <map_components component_2="NL_model" component_1="ATP_production"/>
    <map_variables variable_2="TCa_" variable_1="TCa_"/>
  </connection>
  
  <connection>
    <map_components component_2="L_type_Ca_channel" component_1="ATP_production"/>
    <map_variables variable_2="ATPi" variable_1="ATPi"/>
  </connection>
  
  <connection>
    <map_components component_2="background_KATP_current" component_1="ATP_production"/>
    <map_variables variable_2="ATPi" variable_1="ATPi"/>
  </connection>
  
  <connection>
    <map_components component_2="sodium_potassium_pump" component_1="ATP_production"/>
    <map_variables variable_2="i_NaK" variable_1="i_NaK"/>
    <map_variables variable_2="ATPi" variable_1="ATPi"/>
  </connection>
  
  <connection>
    <map_components component_2="SR_calcium_pump" component_1="ATP_production"/>
    <map_variables variable_2="ATPi" variable_1="ATPi"/>
    <map_variables variable_2="i_SR_U" variable_1="i_SR_U"/>
  </connection>
  
  <connection>
    <map_components component_2="internal_ion_concentrations" component_1="ATP_production"/>
    <map_variables variable_2="Vi" variable_1="Vi"/>
  </connection>
  
  <connection>
    <map_components component_2="rapid_time_dependent_potassium_current" component_1="time_independent_potassium_current"/>
    <map_variables variable_2="E_K" variable_1="E_K"/>
  </connection>
  
  <connection>
    <map_components component_2="background_KATP_current" component_1="time_independent_potassium_current"/>
    <map_variables variable_2="E_K" variable_1="E_K"/>
  </connection>
  
  <connection>
    <map_components component_2="SR_calcium_pump" component_1="Ca_concentrations_in_SR"/>
    <map_variables variable_2="Caup" variable_1="Caup"/>
    <map_variables variable_2="i_SR_U" variable_1="i_SR_U"/>
  </connection>
  
  <connection>
    <map_components component_2="RyR_channel" component_1="L_type_Ca_channel"/>
    <map_variables variable_2="p_open_CaL" variable_1="p_open_CaL"/>
    <map_variables variable_2="iCaL" variable_1="iCaL"/>
  </connection>
  
  <connection>
    <map_components component_2="RyR_channel" component_1="Ca_concentrations_in_SR"/>
    <map_variables variable_2="Carel" variable_1="Carel"/>
    <map_variables variable_2="i_RyR" variable_1="i_RyR"/>
  </connection>
  
  <connection>
    <map_components component_2="SR_T_current" component_1="Ca_concentrations_in_SR"/>
    <map_variables variable_2="Caup" variable_1="Caup"/>
    <map_variables variable_2="Carel" variable_1="Carel"/>
    <map_variables variable_2="i_SR_T" variable_1="i_SR_T"/>
  </connection>
  
  <connection>
    <map_components component_2="SR_L_current" component_1="Ca_concentrations_in_SR"/>
    <map_variables variable_2="Caup" variable_1="Caup"/>
    <map_variables variable_2="i_SR_L" variable_1="i_SR_L"/>
  </connection>
  
  <connection>
    <map_components component_2="sodium_current_voltage_dependent_gate" component_1="sodium_current"/>
    <map_variables variable_2="p_AP_Na" variable_1="p_AP_Na"/>
    <map_variables variable_2="time" variable_1="time"/>
    <map_variables variable_2="Vm" variable_1="Vm"/>
  </connection>
  
  <connection>
    <map_components component_2="sodium_current_ultra_slow_gate" component_1="sodium_current"/>
    <map_variables variable_2="y" variable_1="y"/>
    <map_variables variable_2="time" variable_1="time"/>
    <map_variables variable_2="Vm" variable_1="Vm"/>
  </connection>
  
  <connection>
    <map_components component_2="L_type_Ca_channel_voltage_dependent_gate" component_1="L_type_Ca_channel"/>
    <map_variables variable_2="p_AP_CaL" variable_1="p_AP_CaL"/>
    <map_variables variable_2="time" variable_1="time"/>
    <map_variables variable_2="Vm" variable_1="Vm"/>
  </connection>
  
  <connection>
    <map_components component_2="L_type_Ca_channel_voltage_dependent_gate" component_1="L_type_Ca_channel_Ca_dependent_gate"/>
    <map_variables variable_2="p_AP_CaL" variable_1="p_AP_CaL"/>
  </connection>
  
  <connection>
    <map_components component_2="L_type_Ca_channel_Ca_dependent_gate" component_1="L_type_Ca_channel"/>
    <map_variables variable_2="p_U" variable_1="p_U"/>
    <map_variables variable_2="p_UCa" variable_1="p_UCa"/>
    <map_variables variable_2="Cai" variable_1="Cai"/>
    <map_variables variable_2="CF_Ca" variable_1="CF_Ca"/>
    <map_variables variable_2="time" variable_1="time"/>
    <map_variables variable_2="Vm" variable_1="Vm"/>
    <map_variables variable_2="iCaL" variable_1="iCaL"/>
  </connection>
  
  <connection>
    <map_components component_2="L_type_Ca_channel_ultra_slow_gate" component_1="L_type_Ca_channel"/>
    <map_variables variable_2="y" variable_1="y"/>
    <map_variables variable_2="time" variable_1="time"/>
    <map_variables variable_2="Vm" variable_1="Vm"/>
  </connection>
  
  <connection>
    <map_components component_2="T_type_Ca_channel_y1_gate" component_1="T_type_Ca_channel"/>
    <map_variables variable_2="y1" variable_1="y1"/>
    <map_variables variable_2="time" variable_1="time"/>
    <map_variables variable_2="Vm" variable_1="Vm"/>
  </connection>
  
  <connection>
    <map_components component_2="T_type_Ca_channel_y2_gate" component_1="T_type_Ca_channel"/>
    <map_variables variable_2="y2" variable_1="y2"/>
    <map_variables variable_2="time" variable_1="time"/>
    <map_variables variable_2="Vm" variable_1="Vm"/>
  </connection>
  
  <connection>
    <map_components component_2="time_independent_potassium_current_y_gate" component_1="time_independent_potassium_current"/>
    <map_variables variable_2="y" variable_1="y"/>
    <map_variables variable_2="E_K" variable_1="E_K"/>
    <map_variables variable_2="time" variable_1="time"/>
    <map_variables variable_2="Vm" variable_1="Vm"/>
  </connection>
  
  <connection>
    <map_components component_2="rapid_time_dependent_potassium_current_y1_gate" component_1="rapid_time_dependent_potassium_current"/>
    <map_variables variable_2="y1" variable_1="y1"/>
    <map_variables variable_2="E_K" variable_1="E_K"/>
    <map_variables variable_2="time" variable_1="time"/>
    <map_variables variable_2="Vm" variable_1="Vm"/>
  </connection>
  
  <connection>
    <map_components component_2="rapid_time_dependent_potassium_current_y2_gate" component_1="rapid_time_dependent_potassium_current"/>
    <map_variables variable_2="y2" variable_1="y2"/>
    <map_variables variable_2="E_K" variable_1="E_K"/>
    <map_variables variable_2="time" variable_1="time"/>
    <map_variables variable_2="Vm" variable_1="Vm"/>
  </connection>
  
  <connection>
    <map_components component_2="rapid_time_dependent_potassium_current_y3_gate" component_1="rapid_time_dependent_potassium_current"/>
    <map_variables variable_2="y3" variable_1="y3"/>
    <map_variables variable_2="E_K" variable_1="E_K"/>
    <map_variables variable_2="time" variable_1="time"/>
    <map_variables variable_2="Vm" variable_1="Vm"/>
  </connection>
  
  <connection>
    <map_components component_2="slow_time_dependent_potassium_current_y1_gate" component_1="slow_time_dependent_potassium_current"/>
    <map_variables variable_2="y1" variable_1="y1"/>
    <map_variables variable_2="time" variable_1="time"/>
    <map_variables variable_2="Vm" variable_1="Vm"/>
  </connection>
  
  <connection>
    <map_components component_2="slow_time_dependent_potassium_current_y2_gate" component_1="slow_time_dependent_potassium_current"/>
    <map_variables variable_2="y2" variable_1="y2"/>
    <map_variables variable_2="Cai" variable_1="Cai"/>
    <map_variables variable_2="time" variable_1="time"/>
    <map_variables variable_2="Vm" variable_1="Vm"/>
  </connection>
  
  <connection>
    <map_components component_2="transient_outward_current_y1_gate" component_1="transient_outward_current"/>
    <map_variables variable_2="y1" variable_1="y1"/>
    <map_variables variable_2="time" variable_1="time"/>
    <map_variables variable_2="Vm" variable_1="Vm"/>
  </connection>
  
  <connection>
    <map_components component_2="transient_outward_current_y2_gate" component_1="transient_outward_current"/>
    <map_variables variable_2="y2" variable_1="y2"/>
    <map_variables variable_2="time" variable_1="time"/>
    <map_variables variable_2="Vm" variable_1="Vm"/>
  </connection>
  
  <connection>
    <map_components component_2="sodium_calcium_exchanger_y_gate" component_1="sodium_calcium_exchanger"/>
    <map_variables variable_2="y" variable_1="y"/>
    <map_variables variable_2="p_E1Na" variable_1="p_E1Na"/>
    <map_variables variable_2="p_E2Na" variable_1="p_E2Na"/>
    <map_variables variable_2="p_E1Ca" variable_1="p_E1Ca"/>
    <map_variables variable_2="p_E2Ca" variable_1="p_E2Ca"/>
    <map_variables variable_2="time" variable_1="time"/>
    <map_variables variable_2="k1" variable_1="k1"/>
    <map_variables variable_2="k2" variable_1="k2"/>
    <map_variables variable_2="k3" variable_1="k3"/>
    <map_variables variable_2="k4" variable_1="k4"/>
  </connection>
  
  <connection>
    <map_components component_2="sodium_potassium_pump_y_gate" component_1="sodium_potassium_pump"/>
    <map_variables variable_2="y" variable_1="y"/>
    <map_variables variable_2="p_E1Na" variable_1="p_E1Na"/>
    <map_variables variable_2="p_E2Na" variable_1="p_E2Na"/>
    <map_variables variable_2="p_E1K" variable_1="p_E1K"/>
    <map_variables variable_2="p_E2K" variable_1="p_E2K"/>
    <map_variables variable_2="time" variable_1="time"/>
    <map_variables variable_2="k1" variable_1="k1"/>
    <map_variables variable_2="k2" variable_1="k2"/>
    <map_variables variable_2="k3" variable_1="k3"/>
    <map_variables variable_2="k4" variable_1="k4"/>
  </connection>
  
  <connection>
    <map_components component_2="SR_calcium_pump_y_gate" component_1="SR_calcium_pump"/>
    <map_variables variable_2="y" variable_1="y"/>
    <map_variables variable_2="p_E1" variable_1="p_E1"/>
    <map_variables variable_2="p_E2" variable_1="p_E2"/>
    <map_variables variable_2="p_E1Ca" variable_1="p_E1Ca"/>
    <map_variables variable_2="p_E2Ca" variable_1="p_E2Ca"/>
    <map_variables variable_2="time" variable_1="time"/>
    <map_variables variable_2="k1" variable_1="k1"/>
    <map_variables variable_2="k2" variable_1="k2"/>
    <map_variables variable_2="k3" variable_1="k3"/>
    <map_variables variable_2="k4" variable_1="k4"/>
  </connection>
  







<rdf:RDF>
  <rdf:Bag rdf:about="rdf:#60ac114e-be61-41bd-8b87-33230a869cc0">
    <rdf:li>ventricular myocyte</rdf:li>
    <rdf:li>cardiac</rdf:li>
    <rdf:li>electrophysiology</rdf:li>
  </rdf:Bag>
  <rdf:Seq rdf:about="rdf:#0c4ff849-c9c6-4539-bdfb-73495089fcf5">
    <rdf:li rdf:resource="rdf:#0ed772da-ba9c-4a5f-b07e-a037294ef256"/>
    <rdf:li rdf:resource="rdf:#00b46e5d-2f26-4543-8611-02616bcf9722"/>
    <rdf:li rdf:resource="rdf:#afb46809-bf15-4d1d-a975-756811c6077e"/>
    <rdf:li rdf:resource="rdf:#1293204e-5823-4251-98cc-256b7a68ccdb"/>
    <rdf:li rdf:resource="rdf:#98fd97ed-b640-4c55-998a-1eec6a3c0d8d"/>
  </rdf:Seq>
  <rdf:Description rdf:about="rdf:#15c32719-d5e8-4f74-939e-08c88ed58868">
    <dcterms:modified rdf:resource="rdf:#6c46dc05-d3fd-4afb-b7b2-6d7124f2eac8"/>
    <rdf:value>
          Made MathML id's unique
        </rdf:value>
    <cmeta:modifier rdf:resource="rdf:#92a3e687-33c4-456f-83ff-d9ccb1f7d321"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#a1451727-99c9-41ac-9b67-fc9cebf7fbf5">
    <vCard:Given>Kyoichi</vCard:Given>
    <vCard:Family>Ono</vCard:Family>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#98fd97ed-b640-4c55-998a-1eec6a3c0d8d">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#756f8690-ff84-4833-80cf-0206e8b1daa5"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#c16905ba-1488-4f6a-8faa-9d53697c626c">
    <vCard:Given>Catherine</vCard:Given>
    <vCard:Family>Lloyd</vCard:Family>
    <vCard:Other>May</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#0ed772da-ba9c-4a5f-b07e-a037294ef256">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#b7d7f517-f8d3-4a2b-83c7-23cc2d106716"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#afb46809-bf15-4d1d-a975-756811c6077e">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#27e6b5e3-be81-41ef-aa8e-f635e7c476d5"/>
  </rdf:Description>
  <rdf:Description rdf:about="">
    <dc:publisher>The University of Auckland, Bioengineering Institute</dc:publisher>
    <cmeta:comment rdf:resource="rdf:#ec5790a7-3d96-4ec8-a673-da7359c1e998"/>
    <dcterms:created rdf:resource="rdf:#38bc1627-0ba4-442c-937d-df85fd34176b"/>
    <dc:creator rdf:resource="rdf:#b653b501-be49-4358-924e-c66c178d2cbe"/>
    <cmeta:modification rdf:resource="rdf:#15c32719-d5e8-4f74-939e-08c88ed58868"/>
    <cmeta:modification rdf:resource="rdf:#3793f56d-fa1e-4008-8cdc-527fe735224f"/>
    <cmeta:modification rdf:resource="rdf:#44c5746e-636a-4d5b-964f-55c9aa58bbca"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#1e26060c-0cee-461b-a890-156089dad9ab">
    <dcterms:W3CDTF>2003-01-01</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#b7d7f517-f8d3-4a2b-83c7-23cc2d106716">
    <vCard:Given>Satoshi</vCard:Given>
    <vCard:Family>Matsuoka</vCard:Family>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#575ce958-0e05-407e-a12c-2ce9b5585a91">
    <vCard:FN/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#32eae4e1-7462-4298-818c-fe85f1b71cca">
    <dcterms:W3CDTF>2003-08-13</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#00b46e5d-2f26-4543-8611-02616bcf9722">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#28cb479f-cb15-42d7-84bd-b1384f31a86f"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#4a8781d4-7c40-4d5f-be41-6cc0eba99b8f">
    <dc:title>Japanese Journal of Physiology</dc:title>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#46464609-0036-4898-b312-312a7bc5c16d">
    <vCard:Orgname>The University of Auckland</vCard:Orgname>
    <vCard:Orgunit>The Bioengineering Institute</vCard:Orgunit>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#98511fdd-4de6-4da9-92f8-1ebd9f76936e">
    <vCard:Given>Peter</vCard:Given>
    <vCard:Family>Villiger</vCard:Family>
    <vCard:Other>J</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#1293204e-5823-4251-98cc-256b7a68ccdb">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#a1451727-99c9-41ac-9b67-fc9cebf7fbf5"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#33a6ccaa-97cc-4eeb-9244-883d49c84ba6">
    <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:#db7fae06-d248-4556-a44f-cdfbf0bde0dd">
    <vCard:Given>Catherine</vCard:Given>
    <vCard:Family>Lloyd</vCard:Family>
    <vCard:Other>May</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#44c5746e-636a-4d5b-964f-55c9aa58bbca">
    <dcterms:modified rdf:resource="rdf:#bbd8a98e-5627-448b-ae59-1515b5739e12"/>
    <rdf:value>
          Altered some equations based on the advice of the model author.
        </rdf:value>
    <cmeta:modifier rdf:resource="rdf:#633b1fc9-af50-4457-99e8-b0690c8a099a"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#27e6b5e3-be81-41ef-aa8e-f635e7c476d5">
    <vCard:Given>Shinobu</vCard:Given>
    <vCard:Family>Kuratomi</vCard:Family>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#dec6c971-2ca4-41b9-b985-925fff28a0b5">
    <bqs:Pubmed_id/>
    <bqs:JournalArticle rdf:resource="rdf:#7607a50a-e00a-4607-8707-bf504069e255"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#7607a50a-e00a-4607-8707-bf504069e255">
    <dc:creator rdf:resource="rdf:#0c4ff849-c9c6-4539-bdfb-73495089fcf5"/>
    <dc:title>Role of individual ionic current systems in ventricular cells hypothesized by a model study</dc:title>
    <bqs:volume>53(2)</bqs:volume>
    <bqs:first_page>105</bqs:first_page>
    <bqs:Journal rdf:resource="rdf:#4a8781d4-7c40-4d5f-be41-6cc0eba99b8f"/>
    <dcterms:issued rdf:resource="rdf:#1e26060c-0cee-461b-a890-156089dad9ab"/>
    <bqs:last_page>123</bqs:last_page>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#28cb479f-cb15-42d7-84bd-b1384f31a86f">
    <vCard:Given>Nobuaki</vCard:Given>
    <vCard:Family>Sarai</vCard:Family>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#633b1fc9-af50-4457-99e8-b0690c8a099a">
    <vCard:N rdf:resource="rdf:#5fdc73c5-c24b-49e3-8ffd-787b970e9b2f"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#3793f56d-fa1e-4008-8cdc-527fe735224f">
    <dcterms:modified rdf:resource="rdf:#32eae4e1-7462-4298-818c-fe85f1b71cca"/>
    <rdf:value>
          Altered some equations based on the advice of the model author.
        </rdf:value>
    <cmeta:modifier rdf:resource="rdf:#0fd26e19-0f27-4686-aa29-532bccda64d2"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#38bc1627-0ba4-442c-937d-df85fd34176b">
    <dcterms:W3CDTF>2003-03-14T00:00:00+00:00</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#92a3e687-33c4-456f-83ff-d9ccb1f7d321">
    <vCard:N rdf:resource="rdf:#98511fdd-4de6-4da9-92f8-1ebd9f76936e"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#37ccda2a-cabf-4494-a4fc-d422bf5d6d19">
    <bqs:subject_type>keyword</bqs:subject_type>
    <rdf:value rdf:resource="rdf:#60ac114e-be61-41bd-8b87-33230a869cc0"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#ec5790a7-3d96-4ec8-a673-da7359c1e998">
    <dc:creator rdf:resource="rdf:#575ce958-0e05-407e-a12c-2ce9b5585a91"/>
    <rdf:value/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#6c46dc05-d3fd-4afb-b7b2-6d7124f2eac8">
    <dcterms:W3CDTF>2005-04-20</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#756f8690-ff84-4833-80cf-0206e8b1daa5">
    <vCard:Given>Akinori</vCard:Given>
    <vCard:Family>Noma</vCard:Family>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#b653b501-be49-4358-924e-c66c178d2cbe">
    <vCard:ORG rdf:resource="rdf:#46464609-0036-4898-b312-312a7bc5c16d"/>
    <vCard:EMAIL rdf:resource="rdf:#33a6ccaa-97cc-4eeb-9244-883d49c84ba6"/>
    <vCard:N rdf:resource="rdf:#c16905ba-1488-4f6a-8faa-9d53697c626c"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#e571738f-9960-41e7-9cd9-23f253aebe64">
    <dc:subject rdf:resource="rdf:#37ccda2a-cabf-4494-a4fc-d422bf5d6d19"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#bbd8a98e-5627-448b-ae59-1515b5739e12">
    <dcterms:W3CDTF>2003-07-21</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#0fd26e19-0f27-4686-aa29-532bccda64d2">
    <vCard:N rdf:resource="rdf:#db7fae06-d248-4556-a44f-cdfbf0bde0dd"/>
  </rdf:Description>
  <rdf:Description rdf:about="#matsuoka_sarai_kuratomi_ono_noma_2003_version02">
    <dc:title>
        Matsuoka et al's 2003 mathematical model of the roles of individual ion 
        currents in ventricular cells. 
      </dc:title>
    <cmeta:bio_entity>Ventricular Myocyte</cmeta:bio_entity>
    <cmeta:comment rdf:resource="rdf:#6781a22b-0b67-4806-8472-65e58a40a842"/>
    <bqs:reference rdf:resource="rdf:#dec6c971-2ca4-41b9-b985-925fff28a0b5"/>
    <bqs:reference rdf:resource="rdf:#e571738f-9960-41e7-9cd9-23f253aebe64"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#5fdc73c5-c24b-49e3-8ffd-787b970e9b2f">
    <vCard:Given>Catherine</vCard:Given>
    <vCard:Family>Lloyd</vCard:Family>
    <vCard:Other>May</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#6781a22b-0b67-4806-8472-65e58a40a842">
    <dc:creator rdf:resource="rdf:#403d973e-5f4b-407e-8e00-3ceab4ccc817"/>
    <rdf:value>This is the CellML description of Matsuoka et al's 2003 mathematical model of the roles of individual ion currents in ventricular cells.</rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#403d973e-5f4b-407e-8e00-3ceab4ccc817">
    <vCard:FN>Catherine Lloyd</vCard:FN>
  </rdf:Description>
</rdf:RDF>
</model>