Location: Shannon, Wang, Puglisi, Weber, Bers, 2004 @ 2771759c3eb7 / shannon_wang_puglisi_weber_bers_2004_a.cellml

Author:
pmr2.import <nobody@models.cellml.org>
Date:
2008-02-25 00:26:43+13:00
Desc:
committing version01 of shannon_wang_puglisi_weber_bers_2004
Permanent Source URI:
https://staging.physiomeproject.org/workspace/shannon_wang_puglisi_weber_bers_2004/rawfile/2771759c3eb707fe241d991d33f4430194364780/shannon_wang_puglisi_weber_bers_2004_a.cellml

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

CREATED :  6th November 2004

LAST MODIFIED : 6th November 2004

AUTHOR :  Catherine Lloyd
          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/1/02 CellML Metadata 1.0 Specification.

DESCRIPTION :  This file contains a CellML description of Shannon et al.'s 2004 mathematical model of integrated calcium dynamics within the ventricular myocyte.

CHANGES:  
  
--><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="shannon_wang_puglisi_weber_bers_2004_version01" name="shannon_wang_puglisi_weber_bers_2004_version01">
<documentation xmlns="http://cellml.org/tmp-documentation">
<article>
  <articleinfo>
  <title>Modelling The Integrated Ca Dynamics Of The Ventricular Myocyte</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>
Cardiac myocyte excitation-contraction coupling (ECC) is a complex process which involves the integration of a wide range of systems, complexes and proteins.  Experimentally studying the interactions between several of these processes is difficult, however, mathematical modelling can be used to mimic laboratory experiments, and make the process of studying them much easier.  There are many examples of published mathematical models which describe physiological processes that occur in cardiac myocytes.  In parallel with the increased sophistication of experimental techniques, and the generation of more data, mathematical models have become more accurate at describing biological phenomena. 
</para>

<para>
Using new experimental data, Shannon <emphasis>et al.</emphasis> have recently developed a new mathematical model of rabbit cardiac myocyte Ca and Na homeostasis.  In developing this model they had five aims:
</para>
<itemizedlist>
  <listitem>
            <para>1) To track ion influx and efflux such that the model reaches a steady state solution with a realistic balance of Ca fluxes;</para>
          </listitem>
  <listitem>
            <para>2) To include a component which describes calcium-induced calcium release;</para>
          </listitem>
  <listitem>
            <para>3) To use physiologically realistic parameters that are based on experimental measurements;</para>
          </listitem>
  <listitem>
            <para>4) To be able to accurately simulate basic physiological processes when the components of the model are combined; and</para>
          </listitem>
  <listitem>
            <para>5) To make reasonable compromises to make this model computationally efficient, and possible to solve on a desktop computer.</para>
          </listitem>
</itemizedlist>  

<para>
Shannon <emphasis>et al.</emphasis> conclude that this current model is more robust than many previously existing cardiac myocyte models, and it is able to accurately reproduce physiolgical processes by incorporating parameters which have been obtained experimentally.
</para>

<para>
The model has been described here in CellML (the raw CellML description of the Shannon <emphasis>et al.</emphasis> 2004 model can be downloaded in various formats as described in <xref linkend="sec_download_this_model"/>).     
</para>

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

<para>
<ulink url="http://www.biophysj.org/cgi/content/abstract/87/5/3351">A Mathematical Treatment of Integrated Ca Dynamics within the Ventricular Myocyte</ulink>, Thomas R. Shannon, Fei Wang, Jose Puglisi, Christopher Webber, and Donald M. Bers, 2004, <ulink url="http://www.biophysj.org/">
            <emphasis>Biophysical Journal</emphasis>
          </ulink>, 87, 3351-3371.  (<ulink url="http://www.biophysj.org/cgi/content/full/87/5/3351">Full text (HTML)</ulink> and <ulink url="http://www.biophysj.org/cgi/reprint/87/5/3351.pdf">PDF</ulink> versions of the article are available to subscribers on the <emphasis>Biophysical Journal</emphasis> website.)  <ulink url="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?cmd=Retrieve&amp;db=PubMed&amp;list_uids=15347581&amp;dopt=Abstract">PubMed ID: 15347581</ulink>
</para>


</sect1>
</article>
</documentation>
  
  
  <!--
    Below, we define some additional units for association with variables and
    mitochondrial_membrane within the model. 
  -->
  
  <units name="millivolt">
    <unit units="volt" prefix="milli"/>
  </units>
  
  <units name="per_millivolt">
    <unit units="volt" prefix="milli" exponent="-1"/>
  </units>
  
  <units name="per_millivolt_millisecond">
    <unit units="volt" prefix="milli" exponent="-1"/>
    <unit units="second" prefix="milli" exponent="-1"/>
  </units>
  
  <units name="millisecond">
    <unit units="second" 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="micromole">
    <unit units="mole" prefix="micro"/>
  </units>
  
  <units name="milliS_per_microF">
    <unit units="siemens" prefix="milli"/>
    <unit units="farad" prefix="micro" exponent="-1"/>
  </units>
  
  <units name="microF_per_cm2">
    <unit units="farad" prefix="micro"/>
    <unit units="metre" prefix="centi" exponent="-2"/>
  </units>
  
  <units name="A_per_F">
    <unit units="ampere"/>
    <unit units="farad" exponent="-1"/>
  </units>
  
  <units name="micromolar_per_second">
    <unit units="micromolar" exponent="-1"/>
    <unit units="second" exponent="-1"/>
  </units>
  
  <units name="flux">
    <unit units="millimolar"/>
    <unit units="millisecond" 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="second" exponent="-1"/>
  </units>
  
  <units name="third_order_rate_constant">
    <unit units="millimolar" exponent="-2"/>
    <unit units="millisecond" exponent="-1"/>
  </units>
  
  <units name="micrometre">
    <unit units="metre" prefix="micro"/>
  </units>
  
  <units name="nanometre">
    <unit units="metre" prefix="nano"/>
  </units>
  
  <units name="cm2">
    <unit units="metre" prefix="centi" exponent="2"/>
  </units>
  
  <units name="cm2_per_second">
    <unit units="metre" prefix="centi" exponent="2"/>
    <unit units="second" exponent="-1"/>
  </units>
  
  <units name="cm_per_second">
    <unit units="metre" prefix="centi"/>
    <unit units="second" exponent="-1"/>
  </units>
  
  <units name="micromole_per_litre">
    <unit units="mole" prefix="micro"/>
    <unit units="litre" exponent="-1"/>
  </units>

  <units name="joule_per_mole_kelvin">
    <unit units="joule"/>
    <unit units="mole" exponent="-1"/>
    <unit units="kelvin" exponent="-1"/>
  </units>
  
  <units name="coulomb_per_mole">
    <unit units="coulomb"/>
    <unit units="mole" exponent="-1"/>
  </units>
  
   
  <component name="environment">
    <variable units="second" public_interface="out" name="time"/>
  </component>
  
  <component name="membrane">
    <variable units="millivolt" public_interface="out" name="V"/>
    <variable units="microF_per_cm2" name="Cm" initial_value="2.0"/>
    
    <variable units="second" public_interface="in" name="time"/>
    <variable units="A_per_F" public_interface="in" name="i_Na"/>
    <variable units="A_per_F" public_interface="in" name="i_NaBk"/>
    <variable units="A_per_F" public_interface="in" name="i_NaK"/>
    <variable units="A_per_F" public_interface="in" name="i_Kr"/>
    <variable units="A_per_F" public_interface="in" name="i_Ks"/>
    <variable units="A_per_F" public_interface="in" name="i_tos"/>
    <variable units="A_per_F" public_interface="in" name="i_tof"/>
    <variable units="A_per_F" public_interface="in" name="i_K1"/>
    <variable units="A_per_F" public_interface="in" name="i_NCX"/>
    <variable units="A_per_F" public_interface="in" name="i_ClCa"/>
    <variable units="A_per_F" public_interface="in" name="i_ClBk"/>
    <variable units="A_per_F" public_interface="in" name="i_CaL"/>
    <variable units="A_per_F" public_interface="in" name="i_CaBk"/>
    <variable units="A_per_F" public_interface="in" name="i_SLCaP"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="membrane_voltage_diff_eq">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> V </ci>
        </apply>
        <apply>
          <minus/>
          <apply>
            <divide/>
            <apply>
              <plus/>
              <ci> i_Na </ci>
              <ci> i_NaBk </ci>
              <ci> i_NaK </ci>
              <ci> i_Kr </ci>
              <ci> i_Ks </ci>
              <ci> i_tos </ci>
              <ci> i_tof </ci>
              <ci> i_K1 </ci>
              <ci> i_NCX </ci>
              <ci> i_ClCa </ci>
              <ci> i_ClBk </ci>
              <ci> i_CaL </ci>
              <ci> i_CaBk </ci>
              <ci> i_SLCaP </ci>
            </apply>
            <ci> Cm </ci>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="reversal_potentials">
    <variable units="millivolt" public_interface="out" name="E_Na"/>
    <variable units="millivolt" public_interface="out" name="E_Ca"/>
    <variable units="millivolt" public_interface="out" name="E_K"/>
    <variable units="millivolt" public_interface="out" name="E_Cl"/>
    
    <variable units="joule_per_mole_kelvin" public_interface="in" name="R"/>
    <variable units="kelvin" public_interface="in" name="T"/>
    <variable units="coulomb_per_mole" public_interface="in" name="F"/>
    <variable units="millimolar" public_interface="in" name="Nai"/>
    <variable units="millimolar" public_interface="in" name="Nao"/>
    <variable units="millimolar" public_interface="in" name="Cai"/>
    <variable units="millimolar" public_interface="in" name="Cao"/>
    <variable units="millimolar" public_interface="in" name="Ki"/>
    <variable units="millimolar" public_interface="in" name="Ko"/>
    <variable units="millimolar" public_interface="in" name="Cli"/>
    <variable units="millimolar" public_interface="in" name="Clo"/>
   
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="E_Na_calculation">
        <eq/>
        <ci> E_Na </ci>
        <apply>
          <times/>
          <apply>
            <divide/>
            <apply>
              <times/>
              <ci> R </ci>
              <ci> T </ci>
            </apply>
            <ci> F </ci>
          </apply>
          <apply>
            <ln/>                
            <apply>
              <divide/>
              <ci> Nao </ci>
              <ci> Nai </ci>
            </apply>
          </apply>
        </apply>
      </apply> 
      
      <apply id="E_Ca_calculation">
        <eq/>
        <ci> E_Ca </ci>
        <apply>
          <times/>
          <apply>
            <divide/>
            <apply>
              <times/>
              <ci> R </ci>
              <ci> T </ci>
            </apply>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> 2.0 </cn>
              <ci> F </ci>
            </apply>
          </apply>
          <apply>
            <ln/>                
            <apply>
              <divide/>
              <ci> Cao </ci>
              <ci> Cai </ci>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply id="E_K_calculation">
        <eq/>
        <ci> E_K </ci>
        <apply>
          <times/>
          <apply>
            <divide/>
            <apply>
              <times/>
              <ci> R </ci>
              <ci> T </ci>
            </apply>
            <ci> F </ci>
          </apply>
          <apply>
            <ln/>                
            <apply>
              <divide/>
              <ci> Ko </ci>
              <ci> Ki </ci>
            </apply>
          </apply>
        </apply>
      </apply> 
      
      <apply id="E_Cl_calculation">
        <eq/>
        <ci> E_Cl </ci>
        <apply>
          <times/>
          <apply>
            <divide/>
            <apply>
              <times/>
              <ci> R </ci>
              <ci> T </ci>
            </apply>
            <ci> F </ci>
          </apply>
          <apply>
            <ln/>                
            <apply>
              <divide/>
              <ci> Clo </ci>
              <ci> Cli </ci>
            </apply>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="ion_diffusion">
    <variable units="micromole_per_litre" public_interface="out" name="J_Na_jct_SL"/>
    <variable units="micromole_per_litre" public_interface="out" name="J_Na_SL_cytosol"/>
    <variable units="micromole_per_litre" public_interface="out" name="J_Ca_jct_SL"/>
    <variable units="micromole_per_litre" public_interface="out" name="J_Ca_SL_cytosol"/>
    
    <variable units="cm2_per_second" name="D_Ca_jct_SL" initial_value="1.64E-6"/>
    <variable units="cm2_per_second" name="D_Ca_SL_cytosol" initial_value="1.22E-6"/>
    <variable units="cm2_per_second" name="D_Na_jct_SL" initial_value="1.09E-5"/>
    <variable units="cm2_per_second" name="D_Na_SL_cytosol" initial_value="1.79E-5"/>
    <variable units="micromole" name="x_jct_SL" initial_value="0.5"/>
    <variable units="micromole" name="x_SL_cytosol" initial_value="0.45"/>
    <variable units="cm2" name="A_jct_SL" initial_value="3.01E-6"/>
    <variable units="cm2" name="A_SL_cytosol" initial_value="1.30E-4"/> 
    
    <variable units="millimolar" public_interface="in" name="Na_SL"/>
    <variable units="millimolar" public_interface="in" name="Na_jct"/>
    <variable units="millimolar" public_interface="in" name="Nai"/>
    <variable units="millimolar" public_interface="in" name="Ca_SL"/>
    <variable units="millimolar" public_interface="in" name="Cai"/>
    <variable units="millimolar" public_interface="in" name="Ca_jct"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="J_Na_jct_SL_calculation">
        <eq/>
        <ci> J_Na_jct_SL </ci>
        <apply>
          <times/>
          <ci> D_Na_jct_SL </ci>
          <ci> A_jct_SL </ci>
          <apply>
            <divide/>
            <apply>
              <minus/>
              <ci> Na_jct </ci>
              <ci> Na_SL </ci>
            </apply>
            <ci> x_jct_SL </ci>
          </apply>
        </apply>
      </apply>
      
      <apply id="J_Na_SL_cytosol_calculation">
        <eq/>
        <ci> J_Na_SL_cytosol </ci>
        <apply>
          <times/>
          <ci> D_Na_SL_cytosol </ci>
          <ci> A_SL_cytosol </ci>
          <apply>
            <divide/>
            <apply>
              <minus/>
              <ci> Na_SL </ci>
              <ci> Nai </ci>
            </apply>
            <ci> x_SL_cytosol </ci>
          </apply>
        </apply>
      </apply> 
      
      <apply id="J_Ca_jct_SL_calculation">
        <eq/>
        <ci> J_Ca_jct_SL </ci>
        <apply>
          <times/>
          <ci> D_Ca_jct_SL </ci>
          <ci> A_jct_SL </ci>
          <apply>
            <divide/>
            <apply>
              <minus/>
              <ci> Ca_jct </ci>
              <ci> Ca_SL </ci>
            </apply>
            <ci> x_jct_SL </ci>
          </apply>
        </apply>
      </apply>
      
      <apply id="J_Ca_SL_cytosol_calculation">
        <eq/>
        <ci> J_Ca_SL_cytosol </ci>
        <apply>
          <times/>
          <ci> D_Ca_SL_cytosol </ci>
          <ci> A_SL_cytosol </ci>
          <apply>
            <divide/>
            <apply>
              <minus/>
              <ci> Ca_SL </ci>
              <ci> Cai </ci>
            </apply>
            <ci> x_SL_cytosol </ci>
          </apply>
        </apply>
      </apply> 
    </math>
  </component>
  
  <component name="sodium_buffering">
    <variable units="millimolar" public_interface="out" name="Na_SL"/>
    <variable units="millimolar" public_interface="out" name="Na_jct"/>
    <variable units="millimolar" public_interface="out" name="Nai"/>
    
    <variable units="millimolar" name="Na_SL_L"/>
    <variable units="millimolar" name="Na_jct_L"/>
    <variable units="millimolar" name="Nai_L"/>
    <variable units="millimolar" name="Bmax_SL" initial_value="1.65"/>
    <variable units="millimolar" name="Bmax_jct" initial_value="3.7"/>
    <variable units="millimolar" name="Bmax_cytosol" initial_value="0.0"/>
    <variable units="second_order_rate_constant" name="kon" initial_value="0.1"/>
    <variable units="first_order_rate_constant" name="koff" initial_value="1.0"/>
    
    <variable units="second" public_interface="in" name="time"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="dNa_SL_L_dt">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> Na_SL_L </ci>
        </apply>
        <apply>
          <minus/>
          <apply>
            <times/>
            <ci> kon </ci>
            <ci> Na_SL </ci>
            <apply>
              <minus/>
              <ci> Bmax_SL </ci>
              <ci> Na_SL_L </ci>
            </apply>
          </apply>
          <apply>
            <times/>
            <ci> koff </ci>
            <ci> Na_SL_L </ci>
          </apply>
        </apply>
      </apply>
      
      <apply id="dNa_jct_L_dt">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> Na_jct_L </ci>
        </apply>
        <apply>
          <minus/>
          <apply>
            <times/>
            <ci> kon </ci>
            <ci> Na_jct </ci>
            <apply>
              <minus/>
              <ci> Bmax_jct </ci>
              <ci> Na_jct_L </ci>
            </apply>
          </apply>
          <apply>
            <times/>
            <ci> koff </ci>
            <ci> Na_jct_L </ci>
          </apply>
        </apply>
      </apply>
      
      <apply id="dNai_L_dt">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> Nai_L </ci>
        </apply>
        <apply>
          <minus/>
          <apply>
            <times/>
            <ci> kon </ci>
            <ci> Nai </ci>
            <apply>
              <minus/>
              <ci> Bmax_cytosol </ci>
              <ci> Nai_L </ci>
            </apply>
          </apply>
          <apply>
            <times/>
            <ci> koff </ci>
            <ci> Nai_L </ci>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="fast_sodium_current">
    <variable units="A_per_F" public_interface="out" name="i_Na"/>
    
    <variable units="milliS_per_microF" name="G_INa" initial_value="16.0"/>
    <variable units="dimensionless" name="Fx_Na_jct" initial_value="0.11"/>
    
    <variable units="second" public_interface="in" private_interface="out" name="time"/>
    <variable units="millivolt" public_interface="in" private_interface="out" name="V"/>
    <variable units="millivolt" public_interface="in" name="E_Na"/>
    
    <variable units="dimensionless" private_interface="in" name="m"/>
    <variable units="dimensionless" private_interface="in" name="h"/>
    <variable units="dimensionless" private_interface="in" name="j"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="i_Na_calculation">
        <eq/>
        <ci> i_Na </ci>
        <apply>
          <times/>
          <ci> Fx_Na_jct </ci>
          <ci> G_INa </ci>
          <apply>
            <power/>
            <ci> m </ci>
            <cn cellml:units="dimensionless"> 3.0 </cn>
          </apply>
          <ci> h </ci>
          <ci> j </ci>
          <apply>
            <minus/>
            <ci> V </ci>
            <ci> E_Na </ci>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="fast_sodium_current_h_gate">
    <variable units="dimensionless" public_interface="out" name="h"/>
        
    <variable units="first_order_rate_constant" name="alpha_h"/>
    <variable units="first_order_rate_constant" name="beta_h"/>
    <variable units="second" name="tau_h"/>
    <variable units="dimensionless" name="h_infinity"/>
    
    <variable units="millivolt" public_interface="in" name="V"/>
    <variable units="second" public_interface="in" name="time"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="alpha_h_calculation">
        <eq/>
        <ci> alpha_h </ci>
        <piecewise>
          <piece>
            <apply>
              <times/>
              <cn cellml:units="first_order_rate_constant"> 0.135 </cn>
              <apply>
                <exp/>
                <apply>
                  <divide/>
                  <apply>
                    <plus/>
                     <cn cellml:units="millivolt"> 80.0 </cn>
                     <ci> V </ci>
                   </apply>
                  <cn cellml:units="millivolt"> -6.8 </cn>
                </apply>
              </apply>
            </apply>
            <apply>
              <lt/>
              <ci> V </ci> 
              <cn cellml:units="millivolt"> -40.0 </cn>
            </apply>
          </piece>
          <otherwise>
            <cn cellml:units="first_order_rate_constant"> 0.0 </cn>
          </otherwise>
        </piecewise>
      </apply>
       
      <apply id="beta_h_calculation">
        <eq/>
        <ci> beta_h </ci>
        <piecewise>
          <piece>
            <apply>
              <plus/>
              <apply>
                <times/>
                <cn cellml:units="first_order_rate_constant"> 3.56 </cn>
                <apply>
                  <exp/>
                  <apply>
                    <times/>
                    <cn cellml:units="millivolt"> 0.079 </cn>
                    <ci> V </ci>
                  </apply>
                </apply>
              </apply>
              <apply>
                <times/>
                <cn cellml:units="first_order_rate_constant"> 3.1E5 </cn>
                <apply>
                  <exp/>
                  <apply>
                    <times/>
                    <cn cellml:units="per_millivolt"> 0.35 </cn>
                    <ci> V </ci>
                  </apply>
                </apply>
              </apply>
            </apply>
            <apply>
              <lt/>
              <ci> V </ci>
              <cn cellml:units="millivolt"> -40.0 </cn>
            </apply>
          </piece>
          <otherwise>
            <apply>
              <divide/>
              <cn cellml:units="dimensionless"> 1.0 </cn>
              <apply>
                <times/>
                <cn cellml:units="millisecond"> 0.13 </cn>
                <apply>
                  <plus/>
                  <cn cellml:units="dimensionless"> 1.0 </cn>
                  <apply>
                    <exp/>
                    <apply>
                      <divide/>
                      <apply>
                        <plus/>
                        <ci> V </ci>
                        <cn cellml:units="millivolt"> 10.66 </cn>
                      </apply>
                      <cn cellml:units="millivolt"> -11.1 </cn>
                    </apply>
                  </apply>
                </apply>
              </apply>
            </apply>
          </otherwise>
        </piecewise>
      </apply>
      
      <apply id="tau_h_calculation">
        <eq/>
        <ci> tau_h </ci>
        <apply>
          <divide/>
          <cn cellml:units="dimensionless"> 1.0 </cn>
          <apply>
            <plus/>
            <ci> alpha_h </ci>
            <ci> beta_h </ci>
          </apply>
        </apply>
      </apply>
      
      <apply id="h_infinity_calculation">
        <eq/>
        <ci> h_infinity </ci>
        <apply>
          <divide/>
          <ci> alpha_h </ci>
          <apply>
            <plus/>
            <ci> alpha_h </ci>
            <ci> beta_h </ci>
          </apply>
        </apply>
      </apply>
      
      <apply id="dh_dt">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> h </ci>
        </apply>
        <apply>
          <divide/>
          <apply>
            <minus/>
            <ci> h_infinity </ci>
            <ci> h </ci>
          </apply>
          <ci> tau_h </ci>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="fast_sodium_current_j_gate">
    <variable units="dimensionless" public_interface="out" name="j"/>
        
    <variable units="first_order_rate_constant" name="alpha_j"/>
    <variable units="first_order_rate_constant" name="beta_j"/>
    <variable units="second" name="tau_j"/>
    <variable units="dimensionless" name="j_infinity"/>
    
    <variable units="millivolt" public_interface="in" name="V"/>
    <variable units="second" public_interface="in" name="time"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="alpha_j_calculation">
        <eq/>
        <ci> alpha_j </ci>
        <piecewise>
          <piece>
            <apply>
              <divide/>
              <apply>
                <times/>
                <apply>
                  <minus/>
                  <apply>
                    <times/>
                    <cn cellml:units="per_millivolt_millisecond"> -1.2714E5
  </cn>
                    <apply>
                      <exp/>
                      <apply>
                        <times/> 
                        <cn cellml:units="per_millivolt"> 0.2444 </cn>
                        <ci> V </ci>
                      </apply>
                    </apply>
                  </apply>
                  <apply>
                    <times/>
                    <cn cellml:units="per_millivolt_millisecond"> 3.474E-5 </cn>
                    <apply>
                      <exp/>
                      <apply>
                        <times/>
                        <cn cellml:units="per_millivolt"> -0.04391 </cn>
                        <ci> V </ci>
                      </apply>
                    </apply>
                  </apply>
                </apply>
                <apply>
                  <plus/>
                  <ci> V </ci>
                  <cn cellml:units="millivolt"> 37.78 </cn>
                </apply>
              </apply>
              <apply>
                <plus/>
                <cn cellml:units="dimensionless"> 1.0 </cn>
                <apply>
                  <exp/>
                  <apply>
                    <times/>
                    <cn cellml:units="per_millivolt"> 0.311 </cn>
                    <apply>
                      <plus/>
                      <ci> V </ci>
                      <cn cellml:units="millivolt"> 79.23 </cn>
                    </apply>
                  </apply>
                </apply>
              </apply>
            </apply>
            <apply>
              <lt/>
              <ci> V </ci>
              <cn cellml:units="millivolt"> -40.0 </cn>
            </apply>
          </piece>
          <otherwise>
            <cn cellml:units="first_order_rate_constant"> 0.0 </cn>
          </otherwise>
        </piecewise>
      </apply>
       
      <apply id="beta_j_calculation">
        <eq/>
        <ci> beta_j </ci>
        <piecewise>
          <piece>
            <apply>
              <divide/>
              <apply>
                <times/>
                <cn cellml:units="first_order_rate_constant"> 0.1212 </cn>
                <apply>
                  <exp/>
                  <apply>
                    <times/>
                    <cn cellml:units="per_millivolt"> -0.01052 </cn>
                    <ci> V </ci>
                  </apply>
                </apply>
              </apply>
              <apply>
                <plus/>
                <cn cellml:units="dimensionless"> 1.0 </cn>
                <apply>
                  <exp/>
                  <apply>
                    <times/>
                    <cn cellml:units="per_millivolt"> -0.1378 </cn>
                    <apply>
                      <plus/>       
                      <ci> V </ci>
                      <cn cellml:units="millivolt"> 40.14 </cn>
                    </apply>
                  </apply>
                </apply>
              </apply>
            </apply>
            <apply>
              <lt/>
              <ci> V </ci>
              <cn cellml:units="millivolt"> -40.0 </cn>
            </apply>
          </piece>
          <otherwise>
            <apply>
              <divide/>
              <apply>
                <times/>
                <cn cellml:units="first_order_rate_constant"> 0.3 </cn>
                <apply>
                  <exp/>
                  <apply>
                    <times/>
                    <cn cellml:units="per_millivolt"> -2.535E-7 </cn>
                    <ci> V </ci>                                                
                  </apply>
                </apply>
              </apply>
              <apply>
                <plus/>
                <cn cellml:units="dimensionless"> 1.0 </cn>
                <apply>
                  <exp/>
                  <apply>
                    <times/>
                    <cn cellml:units="per_millivolt"> -0.1 </cn>
                    <apply>
                      <plus/>
                      <ci> V </ci>
                      <cn cellml:units="millivolt"> 32.0 </cn>
                    </apply>
                  </apply>
                </apply>
              </apply>
            </apply>                                                            
          </otherwise> 
        </piecewise>
      </apply>
      
      <apply id="tau_j_calculation">
        <eq/>
        <ci> tau_j </ci>
        <apply>
          <divide/>
          <cn cellml:units="dimensionless"> 1.0 </cn>
          <apply>
            <plus/>
            <ci> alpha_j </ci>
            <ci> beta_j </ci>
          </apply>
        </apply>
      </apply>
      
      <apply id="j_infinity_calculation">
        <eq/>
        <ci> j_infinity </ci>
        <apply>
          <divide/>
          <ci> alpha_j </ci>
          <apply>
            <plus/>
            <ci> alpha_j </ci>
            <ci> beta_j </ci>
          </apply>
        </apply>
      </apply>
      
      <apply id="dj_dt">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> j </ci>
        </apply>
        <apply>
          <divide/>
          <apply>
            <minus/>
            <ci> j_infinity </ci>
            <ci> j </ci>
          </apply>
          <ci> tau_j </ci>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="fast_sodium_current_m_gate">
    <variable units="dimensionless" public_interface="out" name="m"/>
        
    <variable units="first_order_rate_constant" name="alpha_m"/>
    <variable units="first_order_rate_constant" name="beta_m"/>
    <variable units="second" name="tau_m"/>
    <variable units="dimensionless" name="m_infinity"/>
    
    <variable units="millivolt" public_interface="in" name="V"/>
    <variable units="second" public_interface="in" name="time"/>
   
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="alpha_m_calculation">
        <eq/>
        <ci> alpha_m </ci>
        <apply>
          <divide/>
          <apply>
            <times/>
            <cn cellml:units="first_order_rate_constant"> 0.32 </cn>
            <apply>
              <plus/>
              <ci> V </ci>
              <cn cellml:units="millivolt"> 47.13 </cn>
            </apply>
          </apply>
          <apply>
            <minus/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
              <exp/>
              <apply>
                <times/>
                <cn cellml:units="dimensionless"> -0.1 </cn>
                <apply>
                  <plus/>
                  <ci> V </ci>
                  <cn cellml:units="millivolt"> 47.13 </cn>
                </apply>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply id="beta_m_calculation">
        <eq/>
        <ci> beta_m </ci>
        <apply>
          <times/>
          <cn cellml:units="first_order_rate_constant"> 0.8 </cn>
          <apply>
            <exp/>
            <apply>
              <divide/>
              <apply>
                <minus/>
                <ci> V </ci>
              </apply>
              <cn cellml:units="millivolt"> 11.0 </cn>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply id="tau_m_calculation">
        <eq/>
        <ci> tau_m </ci>
        <apply>
          <times/>
          <ci> alpha_m </ci>
          <ci> beta_m </ci>
        </apply>
      </apply>
      
      <apply id="m_infinity_calculation">
        <eq/>
        <ci> m_infinity </ci>
        <apply>
          <divide/>
          <ci> alpha_m </ci>
          <apply>
            <plus/>
            <ci> alpha_m </ci>
            <ci> beta_m </ci>
          </apply>
        </apply>
      </apply>
      
      <apply id="dm_dt">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> m </ci>
        </apply>
        <apply>
          <divide/>
          <apply>
            <minus/>
            <ci> m_infinity </ci>
            <ci> m </ci>
          </apply>
          <ci> tau_m </ci>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="Na_leak_current">
    <variable units="A_per_F" public_interface="out" name="i_NaBk"/>
    
    <variable units="milliS_per_microF" name="G_NaBk" initial_value="0.297E-3"/>
    <variable units="dimensionless" name="Fx_NaBk_jct" initial_value="0.11"/>
    
    <variable units="millivolt" public_interface="in" name="V"/>
    <variable units="millivolt" public_interface="in" name="E_Na"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="i_NaBk_calculation">
        <eq/>
        <ci> i_NaBk </ci>
        <apply>
          <times/>
          <ci> Fx_NaBk_jct </ci>
          <ci> G_NaBk </ci>
          <apply>
            <minus/>
            <ci> V </ci>
            <ci> E_Na </ci>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="Na_K_pump">
    <variable units="A_per_F" public_interface="out" name="i_NaK"/>
    
    <variable units="A_per_F" name="f_NaK"/>
    <variable units="dimensionless" name="Q10_NaK" initial_value="1.63"/>
    <variable units="dimensionless" name="Q_NaK"/>
    <variable units="dimensionless" name="H_NaK" initial_value="4.0"/>
    <variable units="dimensionless" name="Q10_Km_Nai" initial_value="1.49"/>
    <variable units="dimensionless" name="Q_Km_Nai"/>
    <variable units="millimolar" name="Km_Nai" initial_value="11.0"/>
    <variable units="millimolar" name="Km_Ko" initial_value="1.5"/>
    <variable units="A_per_F" name="I_NaK" initial_value="1.91"/>
    <variable units="dimensionless" name="sigma"/>
    <variable units="dimensionless" name="Fx_NaK_jct" initial_value="0.11"/>
    
    <variable units="millivolt" public_interface="in" name="V"/>
    <variable units="joule_per_mole_kelvin" public_interface="in" name="R"/>
    <variable units="kelvin" public_interface="in" name="T"/>
    <variable units="coulomb_per_mole" public_interface="in" name="F"/>  
    <variable units="millimolar" public_interface="in" name="Nai"/>
    <variable units="millimolar" public_interface="in" name="Nao"/>
    <variable units="millimolar" public_interface="in" name="Ko"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="sigma_calculation">
        <eq/>
        <ci> sigma </ci>
        <apply>
          <divide/>
          <apply>
            <minus/>
            <apply>
              <exp/>
              <apply>
                <divide/>
                <ci> Nao </ci>
                <cn cellml:units="millimolar"> 67.3 </cn>
              </apply>
            </apply>
            <cn cellml:units="dimensionless"> 1.0 </cn>
          </apply>
          <cn cellml:units="millimolar"> 7.0 </cn>
        </apply>
      </apply>
      
      <apply id="f_NaK_calculation">
        <eq/>
        <ci> f_NaK </ci>
        <apply>
          <divide/>
          <cn cellml:units="dimensionless"> 1.0 </cn>
          <apply>
            <plus/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> 0.1245 </cn>
              <apply>
                <exp/>
                <apply>
                  <times/>
                  <cn cellml:units="dimensionless"> -0.1 </cn>
                  <apply>
                    <divide/>
                    <apply>
                      <times/>
                      <ci> V </ci>
                      <ci> F </ci>
                    </apply>
                    <apply>
                      <times/>
                      <ci> R </ci>
                      <ci> T </ci>
                    </apply>
                  </apply>
                </apply>
              </apply>
            </apply>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> 0.0365 </cn>
              <ci> sigma </ci>
              <apply>
                <exp/>
                <apply>
                  <divide/>
                  <apply>
                    <times/>
                    <apply>
                      <minus/>
                      <ci> V </ci>
                    </apply>
                    <ci> F </ci>
                  </apply>
                  <apply>
                    <times/>
                    <ci> R </ci>
                    <ci> T </ci>
                  </apply>
                </apply>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
     
      <apply id="i_NaK_calculation">
        <eq/>
        <ci> i_NaK </ci>
        <apply>
          <divide/>
          <apply>
            <times/>
            <ci> Fx_NaK_jct </ci>
            <ci> Q_NaK </ci>
            <ci> I_NaK </ci>
            <ci> f_NaK </ci>
            <ci> Ko </ci>
          </apply>
          <apply>
            <plus/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
              <times/>
              <apply>
                <divide/>
                <apply>
                  <times/>
                  <ci> Q_Km_Nai </ci>
                  <ci> Km_Nai </ci>
                </apply>
                <ci> Nai </ci>
              </apply>
              <ci> H_NaK </ci>
              <apply>
                <plus/>
                <ci> Ko </ci>
                <ci> Km_Ko </ci>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply id="Q_NaK_calculation">
        <eq/>
        <ci> Q_NaK </ci>
        <apply>
          <power/>
          <ci> Q10_NaK </ci>
          <apply>
            <divide/>
            <apply>
              <minus/>
              <ci> T </ci>
              <cn cellml:units="kelvin"> 310.3 </cn>
            </apply>
            <cn cellml:units="kelvin"> 10.0 </cn>
          </apply>
        </apply>
      </apply>
      
      <apply id="Q_Km_Nai_calculation">
        <eq/>
        <ci> Q_Km_Nai </ci>
        <apply>
          <power/>
          <ci> Q10_Km_Nai </ci>
          <apply>
            <divide/>
            <apply>
              <minus/>
              <ci> T </ci>
              <cn cellml:units="kelvin"> 310.3 </cn>
            </apply>
            <cn cellml:units="kelvin"> 10.0 </cn>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="rapid_delayed_rectifier_current"> 
    <variable units="A_per_F" public_interface="out" name="i_Kr"/>          
    <variable units="milliS_per_microF" name="G_IKr"/>
    
    <variable units="second" public_interface="in" private_interface="out" name="time"/>
    <variable units="millivolt" public_interface="in" private_interface="out" name="V"/>   
    <variable units="millivolt" public_interface="in" name="E_K"/>
    <variable units="millimolar" public_interface="in" name="Ko"/>            
    
    <variable units="dimensionless" private_interface="in" name="Xr"/>
    <variable units="dimensionless" private_interface="in" name="Rr"/>
      
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="i_Kr_calculation">
        <eq/>
        <ci> i_Kr </ci>
        <apply>
          <times/>
          <ci> G_IKr </ci>
          <ci> Xr </ci>
          <ci> Rr </ci>
          <apply>
            <minus/>
            <ci> V </ci>
            <ci> E_K </ci>
          </apply>
        </apply>
      </apply> 
      
      <apply id="G_IKr_calculation">
        <eq/>
        <ci> G_IKr </ci>
        <apply>
          <times/>
          <cn cellml:units="milliS_per_microF"> 0.03 </cn>
          <apply>
            <power/>
            <apply>
              <divide/>
              <ci> Ko </ci>
              <cn cellml:units="millimolar"> 5.4 </cn>
            </apply>
            <cn cellml:units="dimensionless"> 0.5 </cn>
          </apply>
        </apply>
      </apply>  
    </math>
  </component>
  
  <component name="rapid_delayed_rectifier_current_Xr_gate">
    <variable units="dimensionless" public_interface="out" name="Xr"/>
   
    <variable units="dimensionless" name="Xr_infinity"/>
    <variable units="second" name="tau_Xr"/>
    
    <variable units="millivolt" public_interface="in" name="V"/>
    <variable units="second" public_interface="in" name="time"/>
   
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="Xr_diff_eq">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> Xr </ci>
        </apply>
        <apply>
          <divide/>
          <apply>
            <minus/>
            <ci> Xr_infinity </ci>
            <ci> Xr </ci>
          </apply>
          <ci> tau_Xr </ci>
        </apply>
      </apply>
      
      <apply id="Xr_infinity_calculation">
        <eq/>
        <ci> Xr_infinity </ci>
        <apply>
          <divide/>
          <cn cellml:units="dimensionless"> 1.0 </cn>
          <apply>
            <plus/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
              <exp/>
              <apply>
                <minus/>
                <apply>
                  <divide/>
                  <apply>
                    <plus/>
                    <cn cellml:units="millivolt"> 50.0 </cn>
                    <ci> V </ci>
                  </apply>
                  <cn cellml:units="millivolt"> 7.5 </cn>
                </apply>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply id="tau_Xr_calculation">
        <eq/>
        <ci> tau_Xr </ci>
        <apply>
          <divide/>
          <cn cellml:units="second"> 1.0 </cn>
          <apply>
            <minus/>
            <apply>
              <plus/>
              <apply>
                <divide/>
                <apply>
                  <times/>
                  <cn cellml:units="dimensionless"> 0.00138 </cn>
                  <apply>
                    <plus/>
                    <ci> V </ci>
                    <cn cellml:units="millivolt"> 7.0 </cn>
                  </apply>
                </apply>
                <apply>
                  <minus/>
                  <cn cellml:units="dimensionless"> 1.0 </cn>
                  <apply>
                    <exp/>
                    <apply>
                      <times/>
                      <cn cellml:units="dimensionless"> -0.123 </cn>
                      <apply>
                        <plus/>
                        <ci> V </ci>
                        <cn cellml:units="millivolt"> 7.0 </cn>
                      </apply>
                    </apply>
                  </apply>
                </apply>
              </apply>
              <apply>
                <divide/>
                <apply>
                  <times/>
                  <cn cellml:units="dimensionless"> 0.00061 </cn>
                  <apply>
                    <plus/>
                    <ci> V </ci>
                    <cn cellml:units="millivolt"> 10.0 </cn>
                  </apply>
                </apply>
                <apply>
                  <exp/>
                  <apply>
                    <times/>
                    <cn cellml:units="dimensionless"> 0.145 </cn>
                    <apply>
                      <plus/>
                      <ci> V </ci>
                      <cn cellml:units="millivolt"> 7.0 </cn>
                    </apply>
                  </apply>
                </apply>
              </apply>
            </apply>
            <cn cellml:units="dimensionless"> 1.0 </cn>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="rapid_delayed_rectifier_current_Rr_gate">
    <variable units="dimensionless" public_interface="out" name="Rr"/>
   
    <variable units="millivolt" public_interface="in" name="V"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="Rr_calculation">
        <eq/>
        <ci> Rr </ci>
        <apply>
          <divide/>
          <cn cellml:units="dimensionless"> 1.0 </cn>
          <apply>
            <plus/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
              <exp/>
              <apply>
                <divide/>
                <apply>
                  <plus/>
                  <cn cellml:units="millivolt"> 33.0 </cn>
                  <ci> V </ci>
                </apply>
                <cn cellml:units="millivolt"> 22.4 </cn>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="slow_delayed_rectifier_current"> 
    <variable units="A_per_F" public_interface="out" name="i_Ks"/>
              
    <variable units="milliS_per_microF" name="G_Ks" initial_value="0.245"/>
    <variable units="dimensionless" name="Fx_Ks_jct" initial_value="0.11"/>
    <variable units="millivolt" name="E_Ks"/>
    <variable units="dimensionless" name="pKNa" initial_value="0.01833"/>
    <variable units="dimensionless" name="pCai"/>
    
    <variable units="joule_per_mole_kelvin" public_interface="in" name="R"/>
    <variable units="kelvin" public_interface="in" name="T"/>
    <variable units="coulomb_per_mole" public_interface="in" name="F"/>
    <variable units="millimolar" public_interface="in" name="Na_jct"/>
    <variable units="millimolar" public_interface="in" name="Nao"/>
    <variable units="millimolar" public_interface="in" name="Ki"/>
    <variable units="millimolar" public_interface="in" name="Cai"/>
    <variable units="millimolar" public_interface="in" name="Ko"/>
    <variable units="millisecond" public_interface="in" private_interface="out" name="time"/>
    <variable units="millivolt" public_interface="in" private_interface="out" name="V"/>   
    
    <variable units="dimensionless" private_interface="in" name="Xs"/>
      
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="pCai_calculation">
        <eq/>
        <ci> pCai </ci>
        <apply>
          <minus/>
          <apply>
            <log/>
            <apply>
              <plus/>
              <ci> Cai </ci>
              <cn cellml:units="millimolar"> 3.0 </cn>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply id="E_Ks_calculation">
        <eq/>
        <ci> E_Ks </ci>
        <apply>
          <times/>
          <apply>
            <divide/>
            <apply>
              <times/>
              <ci> R </ci>
              <ci> T </ci>
            </apply>
            <ci> F </ci>
          </apply>
          <apply>
            <ln/>                
            <apply>
              <divide/>
              <apply>
                <plus/>
                <ci> Ko </ci>
                <apply>
                  <times/>
                  <ci> pKNa </ci>
                  <ci> Nao </ci>
                </apply>
              </apply>
              <apply>
                <plus/>
                <ci> Ki </ci>
                <apply>
                  <times/>
                  <ci> pKNa </ci>
                  <ci> Na_jct </ci>
                </apply>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply> 
      
      <apply id="i_Ks_calculation">
        <eq/>
        <ci> i_Ks </ci>
        <apply>
          <times/>
          <ci> G_Ks </ci>
          <apply>
            <power/>
            <ci> Xs </ci>
            <cn cellml:units="dimensionless"> 2.0 </cn>
          </apply>
          <apply>
            <minus/>
            <ci> V </ci>
            <ci> E_Ks </ci>
          </apply>
        </apply>
      </apply>  
    </math>
  </component>
  
  <component name="slow_delayed_rectifier_current_Xs_gate">
    <variable units="dimensionless" public_interface="out" name="Xs"/>
   
    <variable units="dimensionless" name="Xs_infinity"/>
    <variable units="second" name="tau_Xs"/>
    
    <variable units="millivolt" public_interface="in" name="V"/>
    <variable units="second" public_interface="in" name="time"/>
   
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="Xs_diff_eq">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> Xs </ci>
        </apply>
        <apply>
          <divide/>
          <apply>
            <minus/>
            <ci> Xs_infinity </ci>
            <ci> Xs </ci>
          </apply>
          <ci> tau_Xs </ci>
        </apply>
      </apply>
      
      <apply id="Xs_infinity_calculation">
        <eq/>
        <ci> Xs_infinity </ci>
        <apply>
          <divide/>
          <cn cellml:units="dimensionless"> 1.0 </cn>
          <apply>
            <plus/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
              <exp/>
              <apply>
                <minus/>
                <apply>
                  <divide/>
                  <apply>
                    <minus/>
                    <ci> V </ci>
                    <cn cellml:units="millivolt"> 1.5 </cn>
                  </apply>
                  <cn cellml:units="millivolt"> 16.7 </cn>
                </apply>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply id="tau_Xs_calculation">
        <eq/>
        <ci> tau_Xs </ci>
        <apply>
          <divide/>
          <cn cellml:units="second"> 1.0 </cn>
          <apply>
            <plus/>
            <apply>
              <divide/>
              <apply>
                <times/>
                <cn cellml:units="dimensionless"> 7.19E-5 </cn>
                <apply>
                  <plus/>
                  <ci> V </ci>
                  <cn cellml:units="millivolt"> 30.0 </cn>
                </apply>
              </apply>
              <apply>
                <minus/>
                <cn cellml:units="dimensionless"> 1.0 </cn>
                <apply>
                  <exp/>
                  <apply>
                    <times/>
                    <cn cellml:units="dimensionless"> -0.148 </cn>
                    <apply>
                      <plus/>
                      <ci> V </ci>
                      <cn cellml:units="millivolt"> 30.0 </cn>
                    </apply>
                  </apply>
                </apply>
              </apply>
            </apply>
            <apply>
              <divide/>
              <apply>
                <times/>
                <cn cellml:units="dimensionless"> 1.31E-4 </cn>
                <apply>
                  <plus/>
                  <ci> V </ci>
                  <cn cellml:units="millivolt"> 30.0 </cn>
                </apply>
              </apply>
              <apply>
                <plus/>
                <cn cellml:units="dimensionless"> -1.0 </cn>
                <apply>
                  <exp/>
                  <apply>
                    <times/>
                    <cn cellml:units="dimensionless"> 0.0687 </cn>
                    <apply>
                      <plus/>
                      <ci> V </ci>
                      <cn cellml:units="millivolt"> 30.0 </cn>
                    </apply>
                  </apply>
                </apply>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="slow_transient_outward_current">
    <variable units="A_per_F" public_interface="out" name="i_tos"/>
    
    <variable units="milliS_per_microF" name="G_tos" initial_value="0.2"/>
    
    <variable units="millivolt" public_interface="in" private_interface="out" name="V"/>
    <variable units="second" public_interface="in" private_interface="out" name="time"/>
    <variable units="millivolt" public_interface="in" name="E_K"/>
    
    <variable units="dimensionless" private_interface="in" name="Y_tos"/>
    <variable units="dimensionless" private_interface="in" name="X_tos"/>
    <variable units="dimensionless" private_interface="in" name="Rs_infinity"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="i_tos_calculation">
        <eq/>
        <ci> i_tos </ci>
        <apply>
          <times/>
          <ci> G_tos </ci>
          <ci> X_tos </ci>
          <apply>
            <plus/>
            <ci> Y_tos </ci>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> 0.5 </cn>
              <ci> Rs_infinity </ci>
            </apply>
          </apply>
          <apply>
            <minus/>
            <ci> V </ci>
            <ci> E_K </ci>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="slow_transient_outward_current_X_gate">
    <variable units="dimensionless" public_interface="out" name="X_tos"/>
    
    <variable units="dimensionless" name="X_tos_infinity"/>
    <variable units="second" name="tau_X_tos"/>
  
    <variable units="millivolt" public_interface="in" name="V"/>
    <variable units="second" public_interface="in" name="time"/>
   
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="X_tos_diff_eq">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> X_tos </ci>
        </apply>
        <apply>
          <divide/>
          <apply>
            <minus/>
            <ci> X_tos_infinity </ci>
            <ci> X_tos </ci>
          </apply>
          <ci> tau_X_tos </ci>
        </apply>
      </apply>
      
      <apply id="X_tos_infinity_calculation">
        <eq/>
        <ci> X_tos_infinity </ci>
        <apply>
          <divide/>
          <cn cellml:units="dimensionless"> 1.0 </cn>
          <apply>
            <plus/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
              <exp/>
              <apply>
                <divide/>
                <apply>
                  <minus/>
                  <apply>
                    <plus/>
                    <cn cellml:units="millivolt"> 3.0 </cn>
                    <ci> V </ci>
                  </apply>
                </apply>
                <cn cellml:units="millivolt"> 15.0 </cn>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply id="tau_X_tos_calculation">
        <eq/>
        <ci> tau_X_tos </ci>
        <apply>
          <plus/>
          <apply>
            <divide/>
            <cn cellml:units="second"> 9.0 </cn>
            <apply>
              <plus/>
              <cn cellml:units="dimensionless"> 1.0 </cn>
              <apply>
                <exp/>
                <apply>
                  <divide/>
                  <apply>
                    <plus/>
                    <ci> V </ci>
                    <cn cellml:units="millivolt"> 3.0 </cn>
                  </apply>
                  <cn cellml:units="millivolt"> 15.0 </cn>
                </apply>
              </apply>
            </apply>
          </apply>
          <cn cellml:units="second"> 0.5 </cn>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="slow_transient_outward_current_Y_gate">
    <variable units="dimensionless" public_interface="out" name="Y_tos"/>
    
    <variable units="dimensionless" name="Y_tos_infinity"/>
    <variable units="second" name="tau_Y_tos"/>
  
    <variable units="millivolt" public_interface="in" name="V"/>
    <variable units="second" public_interface="in" name="time"/>
   
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="Y_tos_diff_eq">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> Y_tos </ci>
        </apply>
        <apply>
          <divide/>
          <apply>
            <minus/>
            <ci> Y_tos_infinity </ci>
            <ci> Y_tos </ci>
          </apply>
          <ci> tau_Y_tos </ci>
        </apply>
      </apply>
      
      <apply id="Y_tos_infinity_calculation">
        <eq/>
        <ci> Y_tos_infinity </ci>
        <apply>
          <divide/>
          <cn cellml:units="dimensionless"> 1.0 </cn>
          <apply>
            <plus/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
              <exp/>
              <apply>
                <divide/>
                <apply>
                  <plus/>
                  <cn cellml:units="millivolt"> 33.5 </cn>
                  <ci> V </ci>
                </apply>
                <cn cellml:units="millivolt"> 10.0 </cn>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply id="tau_Y_tos_calculation">
        <eq/>
        <ci> tau_Y_tos </ci>
        <apply>
          <plus/>
          <apply>
            <divide/>
            <cn cellml:units="second"> 3000.0 </cn>
            <apply>
              <plus/>
              <cn cellml:units="dimensionless"> 1.0 </cn>
              <apply>
                <exp/>
                <apply>
                  <divide/>
                  <apply>
                    <plus/>
                    <ci> V </ci>
                    <cn cellml:units="millivolt"> 60.0 </cn>
                  </apply>
                  <cn cellml:units="millivolt"> 10.0 </cn>
                </apply>
              </apply>
            </apply>
          </apply>
          <cn cellml:units="second"> 30.0 </cn>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="slow_transient_outward_current_R_gate">
    <variable units="dimensionless" public_interface="out" name="Rs_infinity"/>
    
    <variable units="millivolt" public_interface="in" name="V"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="Rs_infinity_calculation">
        <eq/>
        <ci> Rs_infinity </ci>
        <apply>
          <divide/>
          <cn cellml:units="dimensionless"> 1.0 </cn>
          <apply>
            <plus/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
              <exp/>
              <apply>
                <divide/>
                <apply>
                  <plus/>
                  <cn cellml:units="millivolt"> 33.5 </cn>
                  <ci> V </ci>
                </apply>
                <cn cellml:units="millivolt"> 10.0 </cn>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
     
  <component name="fast_transient_outward_current">
    <variable units="A_per_F" public_interface="out" name="i_tof"/>
    
    <variable units="milliS_per_microF" name="G_tof" initial_value="0.06"/>
    
    <variable units="millivolt" public_interface="in" private_interface="out" name="V"/>
    <variable units="second" public_interface="in" private_interface="out" name="time"/>
    <variable units="millivolt" public_interface="in" name="E_K"/>
    
    <variable units="dimensionless" private_interface="in" name="Y_tof"/>
    <variable units="dimensionless" private_interface="in" name="X_tof"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="i_tof_calculation">
        <eq/>
        <ci> i_tof </ci>
        <apply>
          <times/>
          <ci> G_tof </ci>
          <ci> X_tof </ci>
          <ci> Y_tof </ci>
          <apply>
            <minus/>
            <ci> V </ci>
            <ci> E_K </ci>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="fast_transient_outward_current_X_gate">
    <variable units="dimensionless" public_interface="out" name="X_tof"/>
    
    <variable units="dimensionless" name="X_tof_infinity"/>
    <variable units="second" name="tau_X_tof"/>
  
    <variable units="millivolt" public_interface="in" name="V"/>
    <variable units="second" public_interface="in" name="time"/>
   
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="X_tof_diff_eq">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> X_tof </ci>
        </apply>
        <apply>
          <divide/>
          <apply>
            <minus/>
            <ci> X_tof_infinity </ci>
            <ci> X_tof </ci>
          </apply>
          <ci> tau_X_tof </ci>
        </apply>
      </apply>
      
      <apply id="X_tof_infinity_calculation">
        <eq/>
        <ci> X_tof_infinity </ci>
        <apply>
          <divide/>
          <cn cellml:units="dimensionless"> 1.0 </cn>
          <apply>
            <plus/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
              <exp/>
              <apply>
                <divide/>
                <apply>
                  <minus/>
                  <apply>
                    <plus/>
                    <cn cellml:units="millivolt"> 3.0 </cn>
                    <ci> V </ci>
                  </apply>
                </apply>
                <cn cellml:units="millivolt"> 15.0 </cn>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply id="tau_X_tof_calculation">
        <eq/>
        <ci> tau_X_tof </ci>
        <apply>
          <plus/>
          <apply>
            <times/>
            <cn cellml:units="second"> 3.5 </cn>
            <apply>
              <exp/>
              <apply>
                <times/>
                <cn cellml:units="dimensionless"> -1.0 </cn>
                <apply>
                  <divide/>
                  <ci> V </ci>
                  <cn cellml:units="millivolt"> 30.0 </cn>
                </apply>
                <apply>
                  <divide/>
                  <ci> V </ci>
                  <cn cellml:units="millivolt"> 30.0 </cn>
                </apply>
              </apply>
            </apply>
          </apply>
          <cn cellml:units="second"> 1.5 </cn>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="fast_transient_outward_current_Y_gate">
    <variable units="dimensionless" public_interface="out" name="Y_tof"/>
    
    <variable units="dimensionless" name="Y_tof_infinity"/>
    <variable units="second" name="tau_Y_tof"/>
  
    <variable units="millivolt" public_interface="in" name="V"/>
    <variable units="second" public_interface="in" name="time"/>
   
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="Y_tof_diff_eq">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> Y_tof </ci>
        </apply>
        <apply>
          <divide/>
          <apply>
            <minus/>
            <ci> Y_tof_infinity </ci>
            <ci> Y_tof </ci>
          </apply>
          <ci> tau_Y_tof </ci>
        </apply>
      </apply>
      
      <apply id="Y_tof_infinity_calculation">
        <eq/>
        <ci> Y_tof_infinity </ci>
        <apply>
          <divide/>
          <cn cellml:units="dimensionless"> 1.0 </cn>
          <apply>
            <plus/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
              <exp/>
              <apply>
                <divide/>
                <apply>
                  <plus/>
                  <cn cellml:units="millivolt"> 33.5 </cn>
                  <ci> V </ci>
                </apply>
                <cn cellml:units="millivolt"> 10.0 </cn>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply id="tau_Y_tof_calculation">
        <eq/>
        <ci> tau_Y_tof </ci>
        <apply>
          <plus/>
          <apply>
            <divide/>
            <cn cellml:units="second"> 20.0 </cn>
            <apply>
              <plus/>
              <cn cellml:units="dimensionless"> 1.0 </cn>
              <apply>
                <exp/>
                <apply>
                  <divide/>
                  <apply>
                    <plus/>
                    <ci> V </ci>
                    <cn cellml:units="millivolt"> 33.5 </cn>
                  </apply>
                  <cn cellml:units="millivolt"> 10.0 </cn>
                </apply>
              </apply>
            </apply>
          </apply>
          <cn cellml:units="second"> 20.0 </cn>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="inward_rectifier_potassium_current"> 
    <variable units="A_per_F" public_interface="out" name="i_K1"/>          
    <variable units="milliS_per_microF" name="G_K1"/>
    
    <variable units="millivolt" public_interface="in" private_interface="out" name="V"/>   
    <variable units="millivolt" public_interface="in" private_interface="out" name="E_K"/>
    <variable units="millimolar" public_interface="in" name="Ko"/>            
    
    <variable units="dimensionless" private_interface="in" name="K1_infinity"/>
      
    <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>
          <ci> K1_infinity </ci>
          <apply>
            <minus/>
            <ci> V </ci>
            <ci> E_K </ci>
          </apply>
        </apply>
      </apply> 
      
      <apply id="G_K1_calculation">
        <eq/>
        <ci> G_K1 </ci>
        <apply>
          <times/>
          <cn cellml:units="milliS_per_microF"> 0.9 </cn>
          <apply>
            <power/>
            <apply>
              <divide/>
              <ci> Ko </ci>
              <cn cellml:units="millimolar"> 5.4 </cn>
            </apply>
            <cn cellml:units="dimensionless"> 0.5 </cn>
          </apply>
        </apply>
      </apply>  
    </math>
  </component>
  
  <component name="inward_rectifier_potassium_current_K1_gate">
    <variable units="dimensionless" public_interface="out" name="K1_infinity"/>
   
    <variable units="first_order_rate_constant" name="alpha_K1"/>
    <variable units="first_order_rate_constant" name="beta_K1"/>
    
    <variable units="millivolt" public_interface="in" name="E_K"/>
    <variable units="millivolt" public_interface="in" name="V"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="K1_infinity_calculation">
        <eq/>
        <ci> K1_infinity </ci>
        <apply>
          <divide/>
          <ci> alpha_K1 </ci>
          <apply>
            <plus/>
            <ci> alpha_K1 </ci>
            <ci> beta_K1 </ci>
          </apply>
        </apply>
      </apply>
      
      <apply id="alpha_K1_calculation">
        <eq/>
        <ci> alpha_K1 </ci>
        <apply>
          <divide/>
          <cn cellml:units="first_order_rate_constant"> 1.02 </cn>
          <apply>
            <plus/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
              <exp/>
              <apply>
                <times/>
                <cn cellml:units="dimensionless"> 0.2385 </cn>
                <apply>
                  <minus/>
                  <ci> V </ci>
                  <apply>
                    <plus/>
                    <ci> E_K </ci>
                    <cn cellml:units="millivolt"> 59.215 </cn>
                  </apply>
                </apply>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply id="beta_K1_calculation">
        <eq/>
        <ci> beta_K1 </ci>
        <apply>
          <divide/>
          <apply>
            <plus/>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> 0.49124 </cn>
              <apply>
                <exp/>
                <apply>
                  <times/>
                  <cn cellml:units="dimensionless"> 0.08032 </cn>
                  <apply>
                    <plus/>
                    <apply>
                      <minus/>
                      <ci> V </ci>
                      <ci> E_K </ci>
                    </apply>
                    <cn cellml:units="millivolt"> 5.476 </cn>
                  </apply>
                </apply>
              </apply>
            </apply>
            <apply>
              <exp/>
              <apply>
                <times/>
                <cn cellml:units="dimensionless"> 0.06175 </cn>
                <apply>
                  <minus/>
                  <ci> V </ci>
                  <apply>
                    <plus/>
                    <ci> E_K </ci>
                    <cn cellml:units="millivolt"> 594.31 </cn>
                  </apply>
                </apply>
              </apply>
            </apply>
          </apply>
          <apply>
            <plus/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
              <exp/>
              <apply>
                <times/>
                <cn cellml:units="dimensionless"> -0.5143 </cn>
                <apply>
                  <plus/>
                  <apply>
                    <minus/>
                    <ci> V </ci>
                    <ci> E_K </ci>
                  </apply>
                  <cn cellml:units="millivolt"> 4.753 </cn>
                </apply>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="calcium_dependent_chloride_current"> 
    <variable units="A_per_F" public_interface="out" name="i_ClCa"/>          
    
    <variable units="milliS_per_microF" name="G_Cl" initial_value="0.109625"/>
    <variable units="micromolar" name="Kd_ClCa" initial_value="100.0"/>
    <variable units="dimensionless" name="Fx_Cl_jct" initial_value="0.11"/>
    
    <variable units="millivolt" public_interface="in" name="V"/>   
    <variable units="millivolt" public_interface="in" name="E_Cl"/>
    <variable units="millimolar" public_interface="in" name="Ca_jct"/>          
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="i_ClCa_calculation">
        <eq/>
        <ci> i_ClCa </ci>
        <apply>
          <times/>
          <ci> Fx_Cl_jct </ci>
          <apply>
            <divide/>
            <apply>
              <times/>
              <ci> G_Cl </ci>
              <apply>
                <minus/>
                <ci> V </ci>
                <ci> E_Cl </ci>
              </apply>
            </apply>
            <apply>
              <plus/>
              <cn cellml:units="dimensionless"> 1.0 </cn>
              <apply>
                <divide/>
                <ci> Kd_ClCa </ci>
                <ci> Ca_jct </ci>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply> 
    </math>
  </component>
  
  <component name="background_chloride_current"> 
    <variable units="A_per_F" public_interface="out" name="i_ClBk"/>          
    
    <variable units="milliS_per_microF" name="G_ClBk" initial_value="0.009"/>
    
    <variable units="millivolt" public_interface="in" name="V"/>   
    <variable units="millivolt" public_interface="in" name="E_Cl"/>
     
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="i_ClBk_calculation">
        <eq/>
        <ci> i_ClBk </ci>
        <apply>
          <times/>
          <ci> G_ClBk </ci>
          <apply>
            <minus/>
            <ci> V </ci>
            <ci> E_Cl </ci>
          </apply>
        </apply>
      </apply> 
    </math>
  </component>
  
  <component name="calcium_buffering">
    <variable units="millimolar" public_interface="out" name="Ca_SL"/>
    <variable units="millimolar" public_interface="out" name="Ca_SR"/>
    <variable units="millimolar" public_interface="out" name="Ca_jct"/>
    <variable units="millimolar" public_interface="out" name="Cai"/>
    
    <variable units="micromolar" name="TroponinC"/>
    <variable units="micromolar" name="TroponinC_Ca"/>
    <variable units="micromolar" name="TroponinC_Mg"/>
    <variable units="micromolar" name="Calmodulin"/>
    <variable units="micromolar" name="Myosin_Ca"/>
    <variable units="micromolar" name="Myosin_Mg"/>
    <variable units="micromolar" name="SR"/>
    <variable units="micromolar" name="SL_SL"/>
    <variable units="micromolar" name="SL_jct"/>
    <variable units="micromolar" name="SLHigh_SL"/>
    <variable units="micromolar" name="SLHigh_jct"/>
    <variable units="micromolar" name="Indo1_Cytosol"/>
    <variable units="micromolar" name="Indo1_SL"/>
    <variable units="micromolar" name="Indo1_jct"/>
    <variable units="micromolar" name="Fluo3_Cytosol"/>
    <variable units="micromolar" name="Fluo3_SL"/>
    <variable units="micromolar" name="Fluo3_jct"/>
    <variable units="micromolar" name="Calsequestrin"/>
    <variable units="micromolar" name="TroponinC_L"/>
    <variable units="micromolar" name="TroponinC_Ca_L"/>
    <variable units="micromolar" name="TroponinC_Mg_L"/>
    <variable units="micromolar" name="Calmodulin_L"/>
    <variable units="micromolar" name="Myosin_Ca_L"/>
    <variable units="micromolar" name="Myosin_Mg_L"/>
    <variable units="micromolar" name="SR_L"/>
    <variable units="micromolar" name="SL_SL_L"/>
    <variable units="micromolar" name="SL_jct_L"/>
    <variable units="micromolar" name="SLHigh_SL_L"/>
    <variable units="micromolar" name="SLHigh_jct_L"/>
    <variable units="micromolar" name="Indo1_Cytosol_L"/>
    <variable units="micromolar" name="Indo1_SL_L"/>
    <variable units="micromolar" name="Indo1_jct_L"/>
    <variable units="micromolar" name="Fluo3_Cytosol_L"/>
    <variable units="micromolar" name="Fluo3_SL_L"/>
    <variable units="micromolar" name="Fluo3_jct_L"/>
    <variable units="micromolar" name="Calsequestrin_L"/>
    <variable units="micromolar" name="Bmax_TroponinC" initial_value="70.0"/>
    <variable units="micromolar" name="Bmax_TroponinC_Ca" initial_value="140.0"/>
    <variable units="micromolar" name="Bmax_TroponinC_Mg" initial_value="140.0"/>
    <variable units="micromolar" name="Bmax_Calmodulin" initial_value="24.0"/>
    <variable units="micromolar" name="Bmax_Myosin_Ca" initial_value="140.0"/>
    <variable units="micromolar" name="Bmax_Myosin_Mg" initial_value="140.0"/>
    <variable units="micromolar" name="Bmax_SR" initial_value="19.0"/>
    <variable units="micromolar" name="Bmax_SL_SL" initial_value="37.4"/>
    <variable units="micromolar" name="Bmax_SL_jct" initial_value="4.6"/>
    <variable units="micromolar" name="Bmax_SLHigh_SL" initial_value="13.4"/>
    <variable units="micromolar" name="Bmax_SLHigh_jct" initial_value="1.65"/>
    <variable units="micromolar" name="Bmax_Indo1_Cytosol" initial_value="25.0"/>
    <variable units="micromolar" name="Bmax_Indo1_SL" initial_value="0.77"/>
    <variable units="micromolar" name="Bmax_Indo1_jct" initial_value="0.02"/>
    <variable units="micromolar" name="Bmax_Fluo3_Cytosol" initial_value="25.0"/>
    <variable units="micromolar" name="Bmax_Fluo3_SL" initial_value="0.77"/>
    <variable units="micromolar" name="Bmax_Fluo3_jct" initial_value="0.02"/>
    <variable units="micromolar" name="Bmax_Calsequestrin" initial_value="140.0"/>
    <variable units="second_order_rate_constant" name="kon_TroponinC" initial_value="32.7"/>
    <variable units="second_order_rate_constant" name="kon_TroponinC_Ca" initial_value="2.37"/>
    <variable units="second_order_rate_constant" name="kon_TroponinC_Mg" initial_value="0.003"/>
    <variable units="second_order_rate_constant" name="kon_Calmodulin" initial_value="34.0"/>
    <variable units="second_order_rate_constant" name="kon_Myosin_Ca" initial_value="13.8"/>
    <variable units="second_order_rate_constant" name="kon_Myosin_Mg" initial_value="0.0157"/>
    <variable units="second_order_rate_constant" name="kon_SR" initial_value="100.0"/>
    <variable units="second_order_rate_constant" name="kon_SL_SL" initial_value="100.0"/>
    <variable units="second_order_rate_constant" name="kon_SL_jct" initial_value="100.0"/>
    <variable units="second_order_rate_constant" name="kon_SLHigh_SL" initial_value="100.0"/>
    <variable units="second_order_rate_constant" name="kon_SLHigh_jct" initial_value="100.0"/>
    <variable units="second_order_rate_constant" name="kon_Indo1_Cytosol" initial_value="100.0"/>
    <variable units="second_order_rate_constant" name="kon_Indo1_SL" initial_value="100.0"/>
    <variable units="second_order_rate_constant" name="kon_Indo1_jct" initial_value="100.0"/>
    <variable units="second_order_rate_constant" name="kon_Fluo3_Cytosol" initial_value="100.0"/>
    <variable units="second_order_rate_constant" name="kon_Fluo3_SL" initial_value="100.0"/>
    <variable units="second_order_rate_constant" name="kon_Fluo3_jct" initial_value="100.0"/>
    <variable units="second_order_rate_constant" name="kon_Calsequestrin" initial_value="100.0"/>
    <variable units="first_order_rate_constant" name="koff_TroponinC" initial_value="19.6"/>
    <variable units="first_order_rate_constant" name="koff_TroponinC_Ca" initial_value="0.032"/>
    <variable units="first_order_rate_constant" name="koff_TroponinC_Mg" initial_value="3.33"/>
    <variable units="first_order_rate_constant" name="koff_Calmodulin" initial_value="238.0"/>
    <variable units="first_order_rate_constant" name="koff_Myosin_Ca" initial_value="0.46"/>
    <variable units="first_order_rate_constant" name="koff_Myosin_Mg" initial_value="0.057"/>
    <variable units="first_order_rate_constant" name="koff_SR" initial_value="60.0"/>
    <variable units="first_order_rate_constant" name="koff_SL_SL" initial_value="1300.0"/>
    <variable units="first_order_rate_constant" name="koff_SL_jct" initial_value="1300.0"/>
    <variable units="first_order_rate_constant" name="koff_SLHigh_SL" initial_value="30.0"/>
    <variable units="first_order_rate_constant" name="koff_SLHigh_jct" initial_value="30.0"/>
    <variable units="first_order_rate_constant" name="koff_Indo1_Cytosol" initial_value="60.0"/>
    <variable units="first_order_rate_constant" name="koff_Indo1_SL" initial_value="60.0"/>
    <variable units="first_order_rate_constant" name="koff_Indo1_jct" initial_value="60.0"/>
    <variable units="first_order_rate_constant" name="koff_Fluo3_Cytosol" initial_value="110.0"/>
    <variable units="first_order_rate_constant" name="koff_Fluo3_SL" initial_value="110.0"/>
    <variable units="first_order_rate_constant" name="koff_Fluo3_jct" initial_value="110.0"/>
    <variable units="first_order_rate_constant" name="koff_Calsequestrin" initial_value="65000.0"/>
    
    <variable units="second" public_interface="in" name="time"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="dTroponinC_L_dt">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> TroponinC_L </ci>
        </apply>
        <apply>
          <minus/>
          <apply>
            <times/>
            <ci> kon_TroponinC </ci>
            <ci> TroponinC </ci>
            <apply>
              <minus/>
              <ci> Bmax_TroponinC </ci>
              <ci> TroponinC_L </ci>
            </apply>
          </apply>
          <apply>
            <times/>
            <ci> koff_TroponinC </ci>
            <ci> TroponinC_L </ci>
          </apply>
        </apply>
      </apply>
      
      <apply id="dTroponinC_Ca_L_dt">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> TroponinC_Ca_L </ci>
        </apply>
        <apply>
          <minus/>
          <apply>
            <times/>
            <ci> kon_TroponinC_Ca </ci>
            <ci> TroponinC_Ca </ci>
            <apply>
              <minus/>
              <ci> Bmax_TroponinC_Ca </ci>
              <ci> TroponinC_Ca_L </ci>
            </apply>
          </apply>
          <apply>
            <times/>
            <ci> koff_TroponinC_Ca </ci>
            <ci> TroponinC_Ca_L </ci>
          </apply>
        </apply>
      </apply>
      
      <apply id="dTroponinC_Mg_L_dt">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> TroponinC_Mg_L </ci>
        </apply>
        <apply>
          <minus/>
          <apply>
            <times/>
            <ci> kon_TroponinC_Mg </ci>
            <ci> TroponinC_Mg </ci>
            <apply>
              <minus/>
              <ci> Bmax_TroponinC_Mg </ci>
              <ci> TroponinC_Mg_L </ci>
            </apply>
          </apply>
          <apply>
            <times/>
            <ci> koff_TroponinC_Mg </ci>
            <ci> TroponinC_Mg_L </ci>
          </apply>
        </apply>
      </apply>
      
      <apply id="dCalmodulin_L_dt">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> Calmodulin_L </ci>
        </apply>
        <apply>
          <minus/>
          <apply>
            <times/>
            <ci> kon_Calmodulin </ci>
            <ci> Calmodulin </ci>
            <apply>
              <minus/>
              <ci> Bmax_Calmodulin </ci>
              <ci> Calmodulin_L </ci>
            </apply>
          </apply>
          <apply>
            <times/>
            <ci> koff_Calmodulin </ci>
            <ci> Calmodulin_L </ci>
          </apply>
        </apply>
      </apply>
      
      <apply id="dMyosin_Ca_L_dt">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> Myosin_Ca_L </ci>
        </apply>
        <apply>
          <minus/>
          <apply>
            <times/>
            <ci> kon_Myosin_Ca </ci>
            <ci> Myosin_Ca </ci>
            <apply>
              <minus/>
              <ci> Bmax_Myosin_Ca </ci>
              <ci> Myosin_Ca_L </ci>
            </apply>
          </apply>
          <apply>
            <times/>
            <ci> koff_Myosin_Ca </ci>
            <ci> Myosin_Ca_L </ci>
          </apply>
        </apply>
      </apply>
      
      <apply id="dMyosin_Mg_L_dt">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> Myosin_Mg_L </ci>
        </apply>
        <apply>
          <minus/>
          <apply>
            <times/>
            <ci> kon_Myosin_Mg </ci>
            <ci> Myosin_Mg </ci>
            <apply>
              <minus/>
              <ci> Bmax_Myosin_Mg </ci>
              <ci> Myosin_Mg_L </ci>
            </apply>
          </apply>
          <apply>
            <times/>
            <ci> koff_Myosin_Mg </ci>
            <ci> Myosin_Mg_L </ci>
          </apply>
        </apply>
      </apply>
      
      <apply id="dSR_L_dt">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> SR_L </ci>
        </apply>
        <apply>
          <minus/>
          <apply>
            <times/>
            <ci> kon_SR </ci>
            <ci> SR </ci>
            <apply>
              <minus/>
              <ci> Bmax_SR </ci>
              <ci> SR_L </ci>
            </apply>
          </apply>
          <apply>
            <times/>
            <ci> koff_SR </ci>
            <ci> SR_L </ci>
          </apply>
        </apply>
      </apply>
      
      <apply id="dSL_SL_L_dt">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> SL_SL_L </ci>
        </apply>
        <apply>
          <minus/>
          <apply>
            <times/>
            <ci> kon_SL_SL </ci>
            <ci> SL_SL </ci>
            <apply>
              <minus/>
              <ci> Bmax_SL_SL </ci>
              <ci> SL_SL_L </ci>
            </apply>
          </apply>
          <apply>
            <times/>
            <ci> koff_SL_SL </ci>
            <ci> SL_SL_L </ci>
          </apply>
        </apply>
      </apply>
      
      <apply id="dSL_jct_L_dt">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> SL_jct_L </ci>
        </apply>
        <apply>
          <minus/>
          <apply>
            <times/>
            <ci> kon_SL_jct </ci>
            <ci> SL_jct </ci>
            <apply>
              <minus/>
              <ci> Bmax_SL_jct </ci>
              <ci> SL_jct_L </ci>
            </apply>
          </apply>
          <apply>
            <times/>
            <ci> koff_SL_jct </ci>
            <ci> SL_jct_L </ci>
          </apply>
        </apply>
      </apply>
      
      <apply id="dSLHigh_SL_L_dt">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> SLHigh_SL_L </ci>
        </apply>
        <apply>
          <minus/>
          <apply>
            <times/>
            <ci> kon_SLHigh_SL </ci>
            <ci> SLHigh_SL </ci>
            <apply>
              <minus/>
              <ci> Bmax_SLHigh_SL </ci>
              <ci> SLHigh_SL_L </ci>
            </apply>
          </apply>
          <apply>
            <times/>
            <ci> koff_SLHigh_SL </ci>
            <ci> SLHigh_SL_L </ci>
          </apply>
        </apply>
      </apply>
      
      <apply id="dSLHigh_jct_L_dt">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> SLHigh_jct_L </ci>
        </apply>
        <apply>
          <minus/>
          <apply>
            <times/>
            <ci> kon_SLHigh_jct </ci>
            <ci> SLHigh_jct </ci>
            <apply>
              <minus/>
              <ci> Bmax_SLHigh_jct </ci>
              <ci> SLHigh_jct_L </ci>
            </apply>
          </apply>
          <apply>
            <times/>
            <ci> koff_SLHigh_jct </ci>
            <ci> SLHigh_jct_L </ci>
          </apply>
        </apply>
      </apply>
      
      <apply id="dIndo1_Cytosol_L_dt">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> Indo1_Cytosol_L </ci>
        </apply>
        <apply>
          <minus/>
          <apply>
            <times/>
            <ci> kon_Indo1_Cytosol </ci>
            <ci> Indo1_Cytosol </ci>
            <apply>
              <minus/>
              <ci> Bmax_Indo1_Cytosol </ci>
              <ci> Indo1_Cytosol_L </ci>
            </apply>
          </apply>
          <apply>
            <times/>
            <ci> koff_Indo1_Cytosol </ci>
            <ci> Indo1_Cytosol_L </ci>
          </apply>
        </apply>
      </apply>
      
      <apply id="dIndo1_SL_L_dt">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> Indo1_SL_L </ci>
        </apply>
        <apply>
          <minus/>
          <apply>
            <times/>
            <ci> kon_Indo1_SL </ci>
            <ci> Indo1_SL </ci>
            <apply>
              <minus/>
              <ci> Bmax_Indo1_SL </ci>
              <ci> Indo1_SL_L </ci>
            </apply>
          </apply>
          <apply>
            <times/>
            <ci> koff_Indo1_SL </ci>
            <ci> Indo1_SL_L </ci>
          </apply>
        </apply>
      </apply>
      
      <apply id="dIndo1_jct_L_dt">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> Indo1_jct_L </ci>
        </apply>
        <apply>
          <minus/>
          <apply>
            <times/>
            <ci> kon_Indo1_jct </ci>
            <ci> Indo1_jct </ci>
            <apply>
              <minus/>
              <ci> Bmax_Indo1_jct </ci>
              <ci> Indo1_jct_L </ci>
            </apply>
          </apply>
          <apply>
            <times/>
            <ci> koff_Indo1_jct </ci>
            <ci> Indo1_jct_L </ci>
          </apply>
        </apply>
      </apply>
      
      <apply id="dFluo3_Cytosol_L_dt">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> Fluo3_Cytosol_L </ci>
        </apply>
        <apply>
          <minus/>
          <apply>
            <times/>
            <ci> kon_Fluo3_Cytosol </ci>
            <ci> Fluo3_Cytosol </ci>
            <apply>
              <minus/>
              <ci> Bmax_Fluo3_Cytosol </ci>
              <ci> Fluo3_Cytosol_L </ci>
            </apply>
          </apply>
          <apply>
            <times/>
            <ci> koff_Fluo3_Cytosol </ci>
            <ci> Fluo3_Cytosol_L </ci>
          </apply>
        </apply>
      </apply>
      
      <apply id="dFluo3_SL_L_dt">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> Fluo3_SL_L </ci>
        </apply>
        <apply>
          <minus/>
          <apply>
            <times/>
            <ci> kon_Fluo3_SL </ci>
            <ci> Fluo3_SL </ci>
            <apply>
              <minus/>
              <ci> Bmax_Fluo3_SL </ci>
              <ci> Fluo3_SL_L </ci>
            </apply>
          </apply>
          <apply>
            <times/>
            <ci> koff_Fluo3_SL </ci>
            <ci> Fluo3_SL_L </ci>
          </apply>
        </apply>
      </apply>
      
      <apply id="dFluo3_jct_L_dt">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> Fluo3_jct_L </ci>
        </apply>
        <apply>
          <minus/>
          <apply>
            <times/>
            <ci> kon_Fluo3_jct </ci>
            <ci> Fluo3_jct </ci>
            <apply>
              <minus/>
              <ci> Bmax_Fluo3_jct </ci>
              <ci> Fluo3_jct_L </ci>
            </apply>
          </apply>
          <apply>
            <times/>
            <ci> koff_Fluo3_jct </ci>
            <ci> Fluo3_jct_L </ci>
          </apply>
        </apply>
      </apply>
      
      <apply id="dCalsequestrin_L_dt">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> Calsequestrin_L </ci>
        </apply>
        <apply>
          <minus/>
          <apply>
            <times/>
            <ci> kon_Calsequestrin </ci>
            <ci> Calsequestrin </ci>
            <apply>
              <minus/>
              <ci> Bmax_Calsequestrin </ci>
              <ci> Calsequestrin_L </ci>
            </apply>
          </apply>
          <apply>
            <times/>
            <ci> koff_Calsequestrin </ci>
            <ci> Calsequestrin_L </ci>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="L_type_calcium_current">
    <variable units="A_per_F" public_interface="out" name="i_CaL"/>
    
    <variable units="A_per_F" name="IP_Ca"/>
    <variable units="A_per_F" name="IP_Na"/>
    <variable units="A_per_F" name="IP_K"/>
    <variable units="A_per_F" name="ICa"/>
    <variable units="A_per_F" name="INa"/>
    <variable units="A_per_F" name="IK"/>
    <variable units="cm_per_second" name="PCa" initial_value="5.4E-4"/>
    <variable units="cm_per_second" name="PNa" initial_value="1.5E-8"/>
    <variable units="cm_per_second" name="PK" initial_value="2.7E-7"/>
    <variable units="dimensionless" name="Fx_CaL" initial_value="0.9"/>
    <variable units="dimensionless" name="gamma_Cai" initial_value="0.341"/>
    <variable units="dimensionless" name="gamma_Cao" initial_value="0.341"/>
    <variable units="dimensionless" name="gamma_Nai" initial_value="0.75"/>
    <variable units="dimensionless" name="gamma_Nao" initial_value="0.75"/>
    <variable units="dimensionless" name="gamma_Ki" initial_value="0.75"/>
    <variable units="dimensionless" name="gamma_Ko" initial_value="0.75"/>
    <variable units="dimensionless" name="Q10_CaL" initial_value="1.8"/>
    <variable units="dimensionless" name="Q_CaL"/>
    
    <variable units="millivolt" public_interface="in" private_interface="out" name="V"/>
    <variable units="second" public_interface="in" private_interface="out" name="time"/>
    <variable units="millimolar" public_interface="in" private_interface="out" name="Cai"/>
    <variable units="millimolar" public_interface="in" name="Na_SL"/>
    <variable units="millimolar" public_interface="in" name="Nao"/>
    <variable units="millimolar" public_interface="in" name="Ca_SL"/>
    <variable units="millimolar" public_interface="in" name="Cao"/>
    <variable units="millimolar" public_interface="in" name="K_SL"/>
    <variable units="millimolar" public_interface="in" name="Ko"/>
    <variable units="joule_per_mole_kelvin" public_interface="in" name="R"/>
    <variable units="kelvin" public_interface="in" name="T"/>
    <variable units="coulomb_per_mole" public_interface="in" name="F"/> 
    
    <variable units="dimensionless" private_interface="in" name="d"/>
    <variable units="dimensionless" private_interface="in" name="f"/>
    <variable units="dimensionless" private_interface="in" name="fCa"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="i_CaL_calculation">
        <eq/>
        <ci> i_CaL </ci>
        <apply>
          <times/>
          <ci> Fx_CaL </ci>
          <apply>
            <plus/>
            <ci> IP_Ca </ci>
            <ci> IP_Na </ci>
            <ci> IP_K </ci>
          </apply>
        </apply>
      </apply>
      
      <apply id="Q_CaL_calculation">
        <eq/>
        <ci> Q_CaL </ci>
        <apply>
          <power/>
          <ci> Q10_CaL </ci>
          <apply>
            <divide/>
            <apply>
              <minus/>
              <ci> T </ci>
              <cn cellml:units="kelvin"> 310.3 </cn>
            </apply>
            <cn cellml:units="kelvin"> 10.0 </cn>
          </apply>
        </apply>
      </apply>
      
      <apply id="IP_Ca_calculation">
        <eq/>
        <ci> IP_Ca </ci>
        <apply>
          <times/>
          <ci> d </ci>
          <ci> fCa </ci>
          <ci> Q_CaL </ci>
          <ci> ICa </ci>
        </apply>
      </apply>
      
      <apply id="IP_Na_calculation">
        <eq/>
        <ci> IP_Na </ci>
        <apply>
          <times/>
          <ci> d </ci>
          <ci> fCa </ci>
          <ci> Q_CaL </ci>
          <ci> INa </ci>
        </apply>
      </apply>
      
      <apply id="IP_K_calculation">
        <eq/>
        <ci> IP_K </ci>
        <apply>
          <times/>
          <ci> d </ci>
          <ci> fCa </ci>
          <ci> Q_CaL </ci>
          <ci> IK </ci>
        </apply>
      </apply>
      
      <apply id="ICa_calculation">
        <eq/>
        <ci> ICa </ci>
        <apply>
          <times/>
          <ci> PCa </ci>
          <apply>
            <power/>
            <cn cellml:units="dimensionless"> 2.0 </cn>
            <cn cellml:units="dimensionless"> 2.0 </cn>
          </apply>
          <ci> V </ci>
          <apply>
            <divide/>
            <apply>
              <power/>
              <ci> F </ci>
              <cn cellml:units="dimensionless"> 2.0 </cn>
            </apply>
            <apply>
              <times/>
              <ci> R </ci>
              <ci> T </ci>
            </apply>
          </apply>
          <apply>
            <divide/>
            <apply>
              <minus/>
              <apply>
                <times/>
                <ci> gamma_Cai </ci>
                <ci> Ca_SL </ci>
                <apply>
                  <exp/>
                  <apply>
                    <divide/>
                    <apply>
                      <times/>
                      <cn cellml:units="dimensionless"> 2.0 </cn>
                      <ci> V </ci>
                      <ci> F </ci>
                    </apply>
                    <apply>
                      <times/>
                      <ci> R </ci>
                      <ci> T </ci>
                    </apply>
                  </apply>
                </apply>
              </apply>
              <apply>
                <times/>
                <ci> gamma_Cao </ci>
                <ci> Cao </ci>
              </apply>
            </apply>
            <apply>
              <minus/>
              <apply>
                <exp/>
                <apply>
                  <divide/>
                  <apply>
                    <times/>
                    <cn cellml:units="dimensionless"> 2.0 </cn>
                    <ci> V </ci>
                    <ci> F </ci>
                  </apply>
                  <apply>
                    <times/>
                    <ci> R </ci>
                    <ci> T </ci>
                  </apply>
                </apply>
              </apply>
              <cn cellml:units="dimensionless"> 1.0 </cn>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply id="INa_calculation">
        <eq/>
        <ci> INa </ci>
        <apply>
          <times/>
          <ci> PNa </ci>
          <apply>
            <power/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <cn cellml:units="dimensionless"> 2.0 </cn>
          </apply>
          <ci> V </ci>
          <apply>
            <divide/>
            <apply>
              <power/>
              <ci> F </ci>
              <cn cellml:units="dimensionless"> 2.0 </cn>
            </apply>
            <apply>
              <times/>
              <ci> R </ci>
              <ci> T </ci>
            </apply>
          </apply>
          <apply>
            <divide/>
            <apply>
              <minus/>
              <apply>
                <times/>
                <ci> gamma_Nai </ci>
                <ci> Na_SL </ci>
                <apply>
                  <exp/>
                  <apply>
                    <divide/>
                    <apply>
                      <times/>
                      <cn cellml:units="dimensionless"> 1.0 </cn>
                      <ci> V </ci>
                      <ci> F </ci>
                    </apply>
                    <apply>
                      <times/>
                      <ci> R </ci>
                      <ci> T </ci>
                    </apply>
                  </apply>
                </apply>
              </apply>
              <apply>
                <times/>
                <ci> gamma_Nao </ci>
                <ci> Nao </ci>
              </apply>
            </apply>
            <apply>
              <minus/>
              <apply>
                <exp/>
                <apply>
                  <divide/>
                  <apply>
                    <times/>
                    <cn cellml:units="dimensionless"> 1.0 </cn>
                    <ci> V </ci>
                    <ci> F </ci>
                  </apply>
                  <apply>
                    <times/>
                    <ci> R </ci>
                    <ci> T </ci>
                  </apply>
                </apply>
              </apply>
              <cn cellml:units="dimensionless"> 1.0 </cn>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply id="IK_calculation">
        <eq/>
        <ci> IK </ci>
        <apply>
          <times/>
          <ci> PK </ci>
          <apply>
            <power/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <cn cellml:units="dimensionless"> 2.0 </cn>
          </apply>
          <ci> V </ci>
          <apply>
            <divide/>
            <apply>
              <power/>
              <ci> F </ci>
              <cn cellml:units="dimensionless"> 2.0 </cn>
            </apply>
            <apply>
              <times/>
              <ci> R </ci>
              <ci> T </ci>
            </apply>
          </apply>
          <apply>
            <divide/>
            <apply>
              <minus/>
              <apply>
                <times/>
                <ci> gamma_Ki </ci>
                <ci> K_SL </ci>
                <apply>
                  <exp/>
                  <apply>
                    <divide/>
                    <apply>
                      <times/>
                      <cn cellml:units="dimensionless"> 1.0 </cn>
                      <ci> V </ci>
                      <ci> F </ci>
                    </apply>
                    <apply>
                      <times/>
                      <ci> R </ci>
                      <ci> T </ci>
                    </apply>
                  </apply>
                </apply>
              </apply>
              <apply>
                <times/>
                <ci> gamma_Ko </ci>
                <ci> Ko </ci>
              </apply>
            </apply>
            <apply>
              <minus/>
              <apply>
                <exp/>
                <apply>
                  <divide/>
                  <apply>
                    <times/>
                    <cn cellml:units="dimensionless"> 1.0 </cn>
                    <ci> V </ci>
                    <ci> F </ci>
                  </apply>
                  <apply>
                    <times/>
                    <ci> R </ci>
                    <ci> T </ci>
                  </apply>
                </apply>
              </apply>
              <cn cellml:units="dimensionless"> 1.0 </cn>
            </apply>
          </apply>
        </apply>
      </apply>
      
      
    </math>
  </component>
  
  <component name="L_type_calcium_current_d_gate">
    <variable units="dimensionless" public_interface="out" name="d"/>
    
    <variable units="dimensionless" name="d_infinity"/>
    <variable units="second" name="tau_d"/>
    <variable units="first_order_rate_constant" name="alpha_d"/>
    <variable units="first_order_rate_constant" name="beta_d"/>
  
    <variable units="millivolt" public_interface="in" name="V"/>
    <variable units="second" public_interface="in" name="time"/>
   
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="d_diff_eq">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> d </ci>
        </apply>
        <apply>
          <divide/>
          <apply>
            <minus/>
            <ci> d_infinity </ci>
            <ci> d </ci>
          </apply>
          <ci> tau_d </ci>
        </apply>
      </apply>
      
      <apply id="d_infinity_calculation">
        <eq/>
        <ci> d_infinity </ci>
        <apply>
          <divide/>
          <cn cellml:units="dimensionless"> 1.0 </cn>
          <apply>
            <plus/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
              <exp/>
              <apply>
                <divide/>
                <apply>
                  <minus/>
                  <apply>
                    <plus/>
                    <cn cellml:units="millivolt"> 14.5 </cn>
                    <ci> V </ci>
                  </apply>
                </apply>
                <cn cellml:units="millivolt"> 6.0 </cn>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply id="tau_d_calculation">
        <eq/>
        <ci> tau_d </ci>
        <apply>
          <times/>
          <ci> d_infinity </ci>
          <apply>
            <divide/>
            <apply>
              <minus/>
              <cn cellml:units="dimensionless"> 1.0 </cn>
              <apply>
                <exp/>
                <apply>
                  <divide/>
                  <apply>
                    <minus/>
                    <apply>
                      <plus/>
                      <cn cellml:units="millivolt"> 14.5 </cn>
                      <ci> V </ci>
                    </apply>
                  </apply>
                  <cn cellml:units="millivolt"> 6.0 </cn>
                </apply>
              </apply>
            </apply>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> 0.035 </cn>
              <apply>
                <plus/>
                <cn cellml:units="millivolt"> 14.5 </cn>
                <ci> V </ci>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply id="alpha_d_calculation">
        <eq/>
        <ci> alpha_d </ci>
        <apply>
          <divide/>
          <ci> d_infinity </ci>
          <ci> tau_d </ci>
        </apply>
      </apply>
      
      <apply id="beta_d_calculation">
        <eq/>
        <ci> beta_d </ci>
        <apply>
          <divide/>
          <apply>
            <minus/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <ci> d_infinity </ci>
          </apply>
          <ci> tau_d </ci>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="L_type_calcium_current_f_gate">
    <variable units="dimensionless" public_interface="out" name="f"/>
    
    <variable units="dimensionless" name="f_infinity"/>
    <variable units="second" name="tau_f"/>
    <variable units="first_order_rate_constant" name="alpha_f"/>
    <variable units="first_order_rate_constant" name="beta_f"/>
  
    <variable units="millivolt" public_interface="in" name="V"/>
    <variable units="second" public_interface="in" name="time"/>
   
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="f_diff_eq">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> f </ci>
        </apply>
        <apply>
          <divide/>
          <apply>
            <minus/>
            <ci> f_infinity </ci>
            <ci> f </ci>
          </apply>
          <ci> tau_f </ci>
        </apply>
      </apply>
      
      <apply id="f_infinity_calculation">
        <eq/>
        <ci> f_infinity </ci>
        <apply>
          <plus/>
          <apply>
            <divide/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
              <plus/>
              <cn cellml:units="dimensionless"> 1.0 </cn>
              <apply>
                <exp/>
                <apply>
                  <divide/>
                  <apply>
                    <plus/>
                    <cn cellml:units="millivolt"> 35.06 </cn>
                    <ci> V </ci>
                  </apply>
                  <cn cellml:units="millivolt"> 3.6 </cn>
                </apply>
              </apply>
            </apply>
          </apply>
          <apply>
            <divide/>
            <cn cellml:units="dimensionless"> 0.6 </cn>
            <apply>
              <plus/>
              <cn cellml:units="dimensionless"> 1.0 </cn>
              <apply>
                <exp/>
                <apply>
                  <divide/>
                  <apply>
                    <minus/>
                    <cn cellml:units="millivolt"> 50.0 </cn>
                    <ci> V </ci>
                  </apply>
                  <cn cellml:units="millivolt"> 20.0 </cn>
                </apply>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply id="tau_f_calculation">
        <eq/>
        <ci> tau_f </ci>
        <apply>
          <divide/>
          <cn cellml:units="second"> 1.0 </cn>
          <apply>
            <plus/>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> 0.0197 </cn>
              <apply>
                <exp/>
                <apply>
                  <minus/>
                  <apply>
                    <power/>
                    <apply>
                      <times/>
                      <cn cellml:units="dimensionless"> 0.0337 </cn>
                      <apply>
                        <plus/>
                        <cn cellml:units="millivolt"> 14.5 </cn>
                        <ci> V </ci>
                      </apply>
                    </apply>
                    <cn cellml:units="dimensionless"> 2.0 </cn>
                  </apply>
                </apply>
              </apply>
            </apply>
            <cn cellml:units="dimensionless"> 0.02 </cn>
          </apply>
        </apply>
      </apply>
      
      <apply id="alpha_f_calculation">
        <eq/>
        <ci> alpha_f </ci>
        <apply>
          <divide/>
          <ci> f_infinity </ci>
          <ci> tau_f </ci>
        </apply>
      </apply>
      
      <apply id="beta_f_calculation">
        <eq/>
        <ci> beta_f </ci>
        <apply>
          <divide/>
          <apply>
            <minus/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <ci> f_infinity </ci>
          </apply>
          <ci> tau_f </ci>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="L_type_calcium_current_fCa_gate">
    <variable units="dimensionless" public_interface="out" name="fCa"/>
    
    <variable units="dimensionless" name="fCaB"/>
    
    <variable units="millimolar" public_interface="in" name="Cai"/>
    <variable units="second" public_interface="in" name="time"/>
   
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="fCa_calculation">
        <eq/>
        <ci> fCa </ci>
        <apply>
          <minus/>
          <cn cellml:units="dimensionless"> 1.0 </cn>
          <ci> fCaB </ci>
        </apply>
      </apply>
      
      <apply id="fCaB_diff_eq">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> fCaB </ci>
        </apply>
        <apply>
          <minus/>
          <apply>
            <times/>
            <cn cellml:units="dimensionless"> 1.7 </cn>
            <ci> Cai </ci>
            <apply>
              <minus/>
              <cn cellml:units="dimensionless"> 1.0 </cn>
              <ci> fCaB </ci>
            </apply>
          </apply>
          <apply>
            <times/>
            <cn cellml:units="dimensionless"> 11.9 </cn>
            <ci> fCaB </ci>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="SL_Ca_leak"> 
    <variable units="A_per_F" public_interface="out" name="i_CaBk"/>          
    
    <variable units="milliS_per_microF" name="G_CaBk" initial_value="0.0002513"/>
    <variable units="dimensionless" name="Fx_CaBk_jct" initial_value="0.11"/>
    
    <variable units="millivolt" public_interface="in" name="V"/>   
    <variable units="millivolt" public_interface="in" name="E_Ca"/>
      
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="i_CaBk_calculation">
        <eq/>
        <ci> i_CaBk </ci>
        <apply>
          <times/>
          <ci> Fx_CaBk_jct </ci>
          <ci> G_CaBk </ci>
          <apply>
            <minus/>
            <ci> V </ci>
            <ci> E_Ca </ci>
          </apply>
        </apply>
      </apply>  
    </math>
  </component>
  
  <component name="Na_Ca_exchange">
    <variable units="A_per_F" public_interface="out" name="i_NCX"/>
    
    <variable units="A_per_F" name="V_max" initial_value="9.0"/>
    <variable units="dimensionless" name="Q10_NCX" initial_value="1.57"/>
    <variable units="dimensionless" name="Q_NCX"/>
    <variable units="millimolar" name="K_mNai" initial_value="12.29"/>
    <variable units="millimolar" name="K_mCao" initial_value="1.3"/>
    <variable units="millimolar" name="K_mNao" initial_value="87.5"/>
    <variable units="micromolar" name="K_mCai" initial_value="3.59"/>
    <variable units="dimensionless" name="Ka"/>
    <variable units="micromolar" name="Kd_act" initial_value="0.256"/>
    <variable units="dimensionless" name="ksat" initial_value="0.27"/>
    <variable units="dimensionless" name="eta" initial_value="0.35"/>
    <variable units="dimensionless" name="HNa" initial_value="3.0"/>
    
    <variable units="millivolt" public_interface="in" name="V"/>
    <variable units="joule_per_mole_kelvin" public_interface="in" name="R"/>
    <variable units="kelvin" public_interface="in" name="T"/>
    <variable units="coulomb_per_mole" public_interface="in" name="F"/>
    <variable units="millimolar" public_interface="in" name="Nai"/>
    <variable units="millimolar" public_interface="in" name="Nao"/>
    <variable units="millimolar" public_interface="in" name="Cai"/>
    <variable units="millimolar" public_interface="in" name="Cao"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="Na_Ca_exchange">
        <eq/>
        <ci> i_NCX </ci>
        <apply>
          <divide/>
          <apply>
            <minus/>
            <apply>
              <times/>
              <ci> Ka </ci>
              <ci> Q_NCX </ci>
              <ci> V_max </ci>
              <apply>
                <exp/>
                <apply>
                  <divide/>
                  <apply>
                    <times/>
                    <ci> eta </ci>
                    <ci> V </ci>
                    <ci> F </ci>
                  </apply>
                  <apply>
                    <times/>
                    <ci> R </ci>
                    <ci> T </ci>
                  </apply>
                </apply>
              </apply>
              <apply>
                <power/>
                <ci> Nai </ci>
                <ci> HNa </ci>
              </apply>
              <ci> Cao </ci>
            </apply>
            <apply>
              <times/>
              <apply>
                <exp/>
                <apply>
                  <divide/>
                  <apply>
                    <times/>
                    <apply>
                      <minus/>
                      <ci> eta </ci>
                      <cn cellml:units="dimensionless"> 1.0 </cn>
                    </apply>
                    <ci> V </ci>
                    <ci> F </ci>
                  </apply>
                  <apply>
                    <times/>
                    <ci> R </ci>
                    <ci> T </ci>
                  </apply>
                </apply>
              </apply>
              <apply>
                <power/>
                <ci> Nao </ci>
                <ci> HNa </ci>
              </apply>
              <ci> Cai </ci>
            </apply>
          </apply>
          <apply>
            <times/>
            <apply>
              <plus/>
              <apply>
                <times/>
                <ci> K_mCai </ci>
                <apply>
                  <power/>
                  <ci> Nao </ci>
                  <ci> HNa </ci>
                </apply>
                <apply>
                  <plus/>
                  <cn cellml:units="dimensionless"> 1.0 </cn>
                  <apply>
                    <power/>
                    <apply>
                      <divide/>
                      <ci> Nai </ci>
                      <ci> K_mNai </ci>
                    </apply>
                    <ci> HNa </ci>
                  </apply>
                </apply>
              </apply>
              <apply>
                <times/>
                <apply>
                  <power/>
                  <ci> K_mNao </ci>
                  <ci> HNa </ci>
                </apply>
                <ci> Cai </ci>
                <apply>
                  <plus/>
                  <cn cellml:units="dimensionless"> 1.0 </cn>
                  <apply>
                    <divide/>
                    <ci> Cai </ci>
                    <ci> K_mCai </ci>
                  </apply>
                </apply>
              </apply>
              <apply>
                <times/>
                <ci> K_mCao </ci>
                <apply>
                  <power/>
                  <ci> Nai </ci>
                  <ci> HNa </ci>
                </apply>
              </apply>
              <apply>
                <times/>
                <apply>
                  <power/>
                  <ci> Nai </ci>
                  <ci> HNa </ci>
                </apply>
                <ci> Cai </ci>
              </apply>
              <apply>
                <times/>
                <apply>
                  <power/>
                  <ci> Nao </ci>
                  <ci> HNa </ci>
                </apply>
                <ci> Cai </ci>
              </apply>
            </apply>
            <apply>
              <plus/>
              <cn cellml:units="dimensionless"> 1.0 </cn>
              <apply>
                <times/>
                <ci> ksat </ci>
                <apply>
                  <exp/>
                  <apply>
                    <times/>
                    <apply>
                      <minus/>
                      <ci> eta </ci>
                      <cn cellml:units="dimensionless"> 1.0 </cn>
                    </apply>
                    <apply>
                      <divide/>
                      <apply>
                        <times/>
                        <ci> V </ci>
                        <ci> F </ci>
                      </apply>
                      <apply>
                        <times/>
                        <ci> R </ci>
                        <ci> T </ci>
                      </apply>
                    </apply>
                  </apply>
                </apply>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply id="Q_NCX_calculation">
        <eq/>
        <ci> Q_NCX </ci>
        <apply>
          <power/>
          <ci> Q10_NCX </ci>
          <apply>
            <divide/>
            <apply>
              <minus/>
              <ci> T </ci>
              <cn cellml:units="kelvin"> 310.3 </cn>
            </apply>
            <cn cellml:units="kelvin"> 10.0 </cn>
          </apply>
        </apply>
      </apply>
     
      <apply id="Ka_calculation">
        <eq/>
        <ci> Ka </ci>
        <apply>
          <divide/>
          <cn cellml:units="dimensionless"> 1.0 </cn>
          <apply>
            <plus/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
              <power/>
              <apply>
                <divide/>
                <ci> Kd_act </ci>
                <ci> Cai </ci>
              </apply>
              <cn cellml:units="dimensionless"> 3.0 </cn>
            </apply>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="SL_Ca_pump">
    <variable units="A_per_F" public_interface="out" name="i_SLCaP"/>
    
    <variable units="micromolar_per_second" name="V_max" initial_value="2.2"/>
    <variable units="dimensionless" name="Fx_SLCaP_jct" initial_value="0.11"/>
    <variable units="dimensionless" name="Q10_SLCaP" initial_value="2.35"/>
    <variable units="dimensionless" name="Q_SLCaP"/>
    <variable units="micromolar" name="Km" initial_value="0.5"/>
    <variable units="dimensionless" name="H" initial_value="1.6"/>
    
    <variable units="millimolar" public_interface="in" name="Cai"/>
    <variable units="kelvin" public_interface="in" name="T"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="i_SLCaP_calculation">
        <eq/>
        <ci> i_SLCaP </ci>
        <apply>
          <divide/>
          <apply>
            <times/>
            <ci> Fx_SLCaP_jct </ci>
            <ci> Q_SLCaP </ci>
            <ci> V_max </ci>
          </apply>
          <apply>
            <plus/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
              <power/>
              <apply>
                <divide/>
                <ci> Km </ci>
                <ci> Cai </ci>
              </apply>
              <ci> H </ci>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply id="Q_SLCaP_calculation">
        <eq/>
        <ci> Q_SLCaP </ci>
        <apply>
          <power/>
          <ci> Q10_SLCaP </ci>
          <apply>
            <divide/>
            <apply>
              <minus/>
              <ci> T </ci>
              <cn cellml:units="kelvin"> 310.3 </cn>
            </apply>
            <cn cellml:units="kelvin"> 10.0 </cn>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="SR_Ca_pump">
    <variable units="flux" public_interface="out" name="J_SRCaP"/>
    
    <variable units="micromolar_per_second" name="V_max" initial_value="286.0"/>
    <variable units="dimensionless" name="Fx_SLCaP_jct" initial_value="0.11"/>
    <variable units="dimensionless" name="Q10_SRCaP" initial_value="2.6"/>
    <variable units="dimensionless" name="Q_SRCaP"/>
    <variable units="micromolar" name="Kmf" initial_value="0.246"/>
    <variable units="millimolar" name="Kmr" initial_value="1.7"/>
    <variable units="dimensionless" name="H" initial_value="1.787"/>
    
    <variable units="millimolar" public_interface="in" name="Cai"/>
    <variable units="millimolar" public_interface="in" name="Ca_SR"/>
    <variable units="kelvin" public_interface="in" name="T"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="J_SRCaP_calculation">
        <eq/>
        <ci> J_SRCaP </ci>
        <apply>
          <times/>
          <ci> Q_SRCaP </ci>
          <apply>
            <divide/>
            <apply>
              <minus/>
              <apply>
                <times/>
                <ci> V_max </ci>
                <apply>
                  <power/>
                  <apply>
                    <divide/>
                    <ci> Cai </ci>
                    <ci> Kmf </ci>
                  </apply>
                  <ci> H </ci>
                </apply>
              </apply>
              <apply>
                <times/>
                <ci> V_max </ci>
                <apply>
                  <power/>
                  <apply>
                    <divide/>
                    <ci> Ca_SR </ci>
                    <ci> Kmr </ci>
                  </apply>
                  <ci> H </ci>
                </apply>
              </apply>
            </apply>
            <apply>
              <plus/>
              <cn cellml:units="dimensionless"> 1.0 </cn>
              <apply>
                <power/>
                <apply>
                  <divide/>
                  <ci> Cai </ci>
                  <ci> Kmf </ci>
                </apply>
                <ci> H </ci>
              </apply>
              <apply>
                <power/>
                <apply>
                  <divide/>
                  <ci> Ca_SR </ci>
                  <ci> Kmr </ci>
                </apply>
                <ci> H </ci>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply id="Q_SRCaP_calculation">
        <eq/>
        <ci> Q_SRCaP </ci>
        <apply>
          <power/>
          <ci> Q10_SRCaP </ci>
          <apply>
            <divide/>
            <apply>
              <minus/>
              <ci> T </ci>
              <cn cellml:units="kelvin"> 310.3 </cn>
            </apply>
            <cn cellml:units="kelvin"> 10.0 </cn>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="SR_Ca_release">
    <variable units="flux" public_interface="out" name="J_SRCaRel"/>
    
    <variable units="dimensionless" name="Max_SR" initial_value="15.0"/>
    <variable units="dimensionless" name="Min_SR" initial_value="1.0"/>
    <variable units="millimolar" name="EC50_SR" initial_value="0.45"/>
    <variable units="dimensionless" name="kCaSR"/>
    <variable units="dimensionless" name="koSRCa"/>
    <variable units="dimensionless" name="kiSRCa"/>
    <variable units="dimensionless" name="R"/>
    <variable units="dimensionless" name="I"/>
    <variable units="dimensionless" name="RI"/>
    <variable units="dimensionless" name="O"/>
    <variable units="first_order_rate_constant" name="ks" initial_value="25.0"/>
    <variable units="third_order_rate_constant" name="koCa" initial_value="10.0"/>
    <variable units="first_order_rate_constant" name="kom" initial_value="0.06"/>
    <variable units="second_order_rate_constant" name="kiCa" initial_value="0.5"/>
    <variable units="first_order_rate_constant" name="kim" initial_value="0.005"/>
    <variable units="dimensionless" name="HSR" initial_value="2.5"/>
    
    <variable units="millimolar" public_interface="in" name="Ca_jct"/>
    <variable units="millimolar" public_interface="in" name="Ca_SR"/>
    <variable units="second" public_interface="in" name="time"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="J_SRCaRel_calculation">
        <eq/>
        <ci> J_SRCaRel </ci>
        <apply>
          <times/>
          <ci> ks </ci>
          <ci> O </ci>
          <apply>
            <minus/>
            <ci> Ca_SR </ci>
            <ci> Ca_jct </ci>
          </apply>
        </apply>
      </apply>
      
      <apply id="kCaSR_calculation">
        <eq/>
        <ci> kCaSR </ci>
        <apply>
          <minus/>
          <ci> Max_SR </ci>
          <apply>
            <divide/>
            <apply>
              <minus/>
              <ci> Max_SR </ci>
              <ci> Min_SR </ci>
            </apply>
            <apply>
              <plus/>
              <cn cellml:units="dimensionless"> 1.0 </cn>
              <apply>
                <power/>
                <apply>
                  <divide/>
                  <ci> EC50_SR </ci>
                  <ci> Ca_SR </ci>
                </apply>
                <ci> HSR </ci>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply id="koSRCa_calculation">
        <eq/>
        <ci> koSRCa </ci>
        <apply>
          <divide/>
          <ci> koCa </ci>
          <ci> kCaSR </ci>
        </apply>
      </apply>
      
      <apply id="kiSRCa_calculation">
        <eq/>
        <ci> kiSRCa </ci>
        <apply>
          <times/>
          <ci> kiCa </ci>
          <ci> kCaSR </ci>
        </apply>
      </apply>
      
      <apply id="R_diff_eq">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> R </ci>
        </apply>
        <apply>
          <minus/>
          <apply>
            <minus/>
            <apply>
              <times/>
              <ci> kim </ci>
              <ci> RI </ci>
            </apply>
            <apply>
              <times/>
              <ci> kiSRCa </ci>
              <ci> Ca_jct </ci>
              <ci> R </ci>
            </apply>
          </apply>
          <apply>
            <minus/>
            <apply>
              <times/>
              <ci> koSRCa </ci>
              <apply>
                <power/>
                <ci> Ca_jct </ci>
                <cn cellml:units="dimensionless"> 2.0 </cn>
              </apply>
              <ci> R </ci>
            </apply>
            <apply>
              <times/>
              <ci> kom </ci>
              <ci> O </ci>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply id="O_diff_eq">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> O </ci>
        </apply>
        <apply>
          <minus/>
          <apply>
            <minus/>
            <apply>
              <times/>
              <ci> koSRCa </ci>
              <apply>
                <power/>
                <ci> Ca_jct </ci>
                <cn cellml:units="dimensionless"> 2.0 </cn>
              </apply>
              <ci> R </ci>
            </apply>
            <apply>
              <times/>
              <ci> kom </ci>
              <ci> O </ci>
            </apply>
          </apply>
          <apply>
            <minus/>
            <apply>
              <times/>
              <ci> kiSRCa </ci>
              <ci> Ca_jct </ci>
              <ci> O </ci>
            </apply>
            <apply>
              <times/>
              <ci> kim </ci>
              <ci> I </ci>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply id="I_diff_eq">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> I </ci>
        </apply>
        <apply>
          <minus/>
          <apply>
            <minus/>
            <apply>
              <times/>
              <ci> kiSRCa </ci>
              <ci> Ca_jct </ci>
              <ci> O </ci>
            </apply>
            <apply>
              <times/>
              <ci> kim </ci>
              <ci> I </ci>
            </apply>
          </apply>
          <apply>
            <minus/>
            <apply>
              <times/>
              <ci> kom </ci>
              <ci> I </ci>
            </apply>
            <apply>
              <times/>
              <ci> koSRCa </ci>
              <apply>
                <power/>
                <ci> Ca_jct </ci>
                <cn cellml:units="dimensionless"> 2.0 </cn>
              </apply>
              <ci> RI </ci>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply id="RI_diff_eq">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> RI </ci>
        </apply>
        <apply>
          <minus/>
          <apply>
            <minus/>
            <apply>
              <times/>
              <ci> kom </ci>
              <ci> I </ci>
            </apply>
            <apply>
              <times/>
              <ci> koSRCa </ci>
              <apply>
                <power/>
                <ci> Ca_jct </ci>
                <cn cellml:units="dimensionless"> 2.0 </cn>
              </apply>
              <ci> RI </ci>
            </apply>
          </apply>
          <apply>
            <minus/>
            <apply>
              <times/>
              <ci> kim </ci>
              <ci> RI </ci>
            </apply>
            <apply>
              <times/>
              <ci> kiSRCa </ci>
              <ci> Ca_jct </ci>
              <ci> R </ci>
            </apply>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="passive_SR_Ca_leak">
    <variable units="flux" public_interface="out" name="J_SRLeak"/>
    
    <variable units="first_order_rate_constant" name="KSRleak" initial_value="5.348E-6"/>
    
    <variable units="millimolar" public_interface="in" name="Ca_jct"/>
    <variable units="millimolar" public_interface="in" name="Ca_SR"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="J_SRLeak_calculation">
        <eq/>
        <ci> J_SRLeak </ci>
        <apply>
          <times/>
          <ci> KSRleak </ci>
          <apply>
            <minus/>
            <ci> Ca_SR </ci>
            <ci> Ca_jct </ci>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="model_parameters">
    <variable units="millimolar" public_interface="out" name="Ko" initial_value="5.4"/>
    <variable units="millimolar" public_interface="out" name="K_SL"/>
    <variable units="millimolar" public_interface="out" name="Nao" initial_value="145.0"/>
    <variable units="millimolar" public_interface="out" name="Cao" initial_value="1.8"/>
    <variable units="millimolar" public_interface="out" name="Clo" initial_value="150.0"/>
    <variable units="millimolar" public_interface="out" name="Ki" initial_value="135.0"/>
    <variable units="millimolar" public_interface="out" name="Mgi" initial_value="1.0"/>
    <variable units="millimolar" public_interface="out" name="Cli" initial_value="15.0"/>
    <variable units="joule_per_mole_kelvin" 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_mole" public_interface="out" name="F" initial_value="96.4867"/>
  </component>
  
  
  <group>
    <relationship_ref relationship="containment"/>
    <component_ref component="membrane">
      <component_ref component="fast_sodium_current">
        <component_ref component="fast_sodium_current_m_gate"/>
        <component_ref component="fast_sodium_current_h_gate"/>
        <component_ref component="fast_sodium_current_j_gate"/>
      </component_ref>
      <component_ref component="L_type_calcium_current">
        <component_ref component="L_type_calcium_current_d_gate"/>
        <component_ref component="L_type_calcium_current_f_gate"/>
        <component_ref component="L_type_calcium_current_fCa_gate"/>
      </component_ref>
      <component_ref component="slow_delayed_rectifier_current">
        <component_ref component="slow_delayed_rectifier_current_Xs_gate"/>
      </component_ref>
      <component_ref component="rapid_delayed_rectifier_current">
        <component_ref component="rapid_delayed_rectifier_current_Xr_gate"/>
        <component_ref component="rapid_delayed_rectifier_current_Rr_gate"/>
      </component_ref>
      <component_ref component="fast_transient_outward_current">
        <component_ref component="fast_transient_outward_current_X_gate"/>
        <component_ref component="fast_transient_outward_current_Y_gate"/>
      </component_ref>
      <component_ref component="slow_transient_outward_current">
        <component_ref component="slow_transient_outward_current_X_gate"/>
        <component_ref component="slow_transient_outward_current_Y_gate"/>
        <component_ref component="slow_transient_outward_current_R_gate"/>
      </component_ref>
      <component_ref component="inward_rectifier_potassium_current">
        <component_ref component="inward_rectifier_potassium_current_K1_gate"/>
      </component_ref>
      <component_ref component="reversal_potentials"/>
      <component_ref component="ion_diffusion"/>
      <component_ref component="sodium_buffering"/>
      <component_ref component="calcium_dependent_chloride_current"/>
      <component_ref component="background_chloride_current"/>
      <component_ref component="Na_Ca_exchange"/>
      <component_ref component="SL_Ca_pump"/>
      <component_ref component="SR_Ca_pump"/>
      <component_ref component="SR_Ca_release"/>
      <component_ref component="passive_SR_Ca_leak"/>
      <component_ref component="Na_leak_current"/>
      <component_ref component="Na_K_pump"/>
      <component_ref component="SL_Ca_leak"/>
    </component_ref>
  </group>
  
  <group>
    <relationship_ref relationship="encapsulation"/>
    <component_ref component="fast_sodium_current">
      <component_ref component="fast_sodium_current_m_gate"/>
      <component_ref component="fast_sodium_current_h_gate"/>
      <component_ref component="fast_sodium_current_j_gate"/>
    </component_ref>
    <component_ref component="L_type_calcium_current">
      <component_ref component="L_type_calcium_current_d_gate"/>
      <component_ref component="L_type_calcium_current_f_gate"/>
      <component_ref component="L_type_calcium_current_fCa_gate"/>
    </component_ref>
    <component_ref component="slow_delayed_rectifier_current">
      <component_ref component="slow_delayed_rectifier_current_Xs_gate"/>
    </component_ref>
    <component_ref component="rapid_delayed_rectifier_current">
      <component_ref component="rapid_delayed_rectifier_current_Xr_gate"/>
      <component_ref component="rapid_delayed_rectifier_current_Rr_gate"/>
    </component_ref>
    <component_ref component="fast_transient_outward_current">
      <component_ref component="fast_transient_outward_current_X_gate"/>
      <component_ref component="fast_transient_outward_current_Y_gate"/>
    </component_ref>
    <component_ref component="slow_transient_outward_current">
      <component_ref component="slow_transient_outward_current_X_gate"/>
      <component_ref component="slow_transient_outward_current_Y_gate"/>
      <component_ref component="slow_transient_outward_current_R_gate"/>
    </component_ref>
    <component_ref component="inward_rectifier_potassium_current">
      <component_ref component="inward_rectifier_potassium_current_K1_gate"/>
    </component_ref>
  </group>
  
  
  <connection>
    <map_components component_2="environment" component_1="membrane"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="environment" component_1="sodium_buffering"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="environment" component_1="fast_sodium_current"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="environment" component_1="slow_transient_outward_current"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="environment" component_1="fast_transient_outward_current"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="environment" component_1="slow_delayed_rectifier_current"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="environment" component_1="rapid_delayed_rectifier_current"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="environment" component_1="SR_Ca_release"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="environment" component_1="calcium_buffering"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="environment" component_1="L_type_calcium_current"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="fast_sodium_current" component_1="membrane"/>
    <map_variables variable_2="V" variable_1="V"/>
    <map_variables variable_2="i_Na" variable_1="i_Na"/>
  </connection>
  
  <connection>
    <map_components component_2="L_type_calcium_current" component_1="membrane"/>
    <map_variables variable_2="V" variable_1="V"/>
    <map_variables variable_2="i_CaL" variable_1="i_CaL"/>
  </connection>
  
  <connection>
    <map_components component_2="slow_transient_outward_current" component_1="membrane"/>
    <map_variables variable_2="V" variable_1="V"/>
    <map_variables variable_2="i_tos" variable_1="i_tos"/>
  </connection>
  
  <connection>
    <map_components component_2="fast_transient_outward_current" component_1="membrane"/>
    <map_variables variable_2="V" variable_1="V"/>
    <map_variables variable_2="i_tof" variable_1="i_tof"/>
  </connection>
  
  <connection>
    <map_components component_2="slow_delayed_rectifier_current" component_1="membrane"/>
    <map_variables variable_2="V" variable_1="V"/>
    <map_variables variable_2="i_Ks" variable_1="i_Ks"/>
  </connection>
  
  <connection>
    <map_components component_2="rapid_delayed_rectifier_current" component_1="membrane"/>
    <map_variables variable_2="V" variable_1="V"/>
    <map_variables variable_2="i_Kr" variable_1="i_Kr"/>
  </connection>
  
  <connection>
    <map_components component_2="inward_rectifier_potassium_current" component_1="membrane"/>
    <map_variables variable_2="V" variable_1="V"/>
    <map_variables variable_2="i_K1" variable_1="i_K1"/>
  </connection>
  
  <connection>
    <map_components component_2="Na_Ca_exchange" component_1="membrane"/>
    <map_variables variable_2="V" variable_1="V"/>
    <map_variables variable_2="i_NCX" variable_1="i_NCX"/>
  </connection>
  
  <connection>
    <map_components component_2="Na_K_pump" component_1="membrane"/>
    <map_variables variable_2="V" variable_1="V"/>
    <map_variables variable_2="i_NaK" variable_1="i_NaK"/>
  </connection>
  
  <connection>
    <map_components component_2="Na_leak_current" component_1="membrane"/>
    <map_variables variable_2="V" variable_1="V"/>
    <map_variables variable_2="i_NaBk" variable_1="i_NaBk"/>
  </connection>
  
  <connection>
    <map_components component_2="SL_Ca_leak" component_1="membrane"/>
    <map_variables variable_2="V" variable_1="V"/>
    <map_variables variable_2="i_CaBk" variable_1="i_CaBk"/>
  </connection>
  
  <connection>
    <map_components component_2="calcium_dependent_chloride_current" component_1="membrane"/>
    <map_variables variable_2="V" variable_1="V"/>
    <map_variables variable_2="i_ClCa" variable_1="i_ClCa"/>
  </connection>
  
  <connection>
    <map_components component_2="background_chloride_current" component_1="membrane"/>
    <map_variables variable_2="V" variable_1="V"/>
    <map_variables variable_2="i_ClBk" variable_1="i_ClBk"/>
  </connection>
  
  <connection>
    <map_components component_2="SL_Ca_pump" component_1="membrane"/>
    <map_variables variable_2="i_SLCaP" variable_1="i_SLCaP"/>
  </connection>
  
  <connection>
    <map_components component_2="fast_sodium_current" component_1="reversal_potentials"/>
    <map_variables variable_2="E_Na" variable_1="E_Na"/>
  </connection>
  
  <connection>
    <map_components component_2="slow_transient_outward_current" component_1="reversal_potentials"/>
    <map_variables variable_2="E_K" variable_1="E_K"/>
  </connection>
  
  <connection>
    <map_components component_2="fast_transient_outward_current" component_1="reversal_potentials"/>
    <map_variables variable_2="E_K" variable_1="E_K"/>
  </connection>
  
  <connection>
    <map_components component_2="rapid_delayed_rectifier_current" component_1="reversal_potentials"/>
    <map_variables variable_2="E_K" variable_1="E_K"/>
  </connection>
  
  <connection>
    <map_components component_2="inward_rectifier_potassium_current" component_1="reversal_potentials"/>
    <map_variables variable_2="E_K" variable_1="E_K"/>
  </connection>
  
  <connection>
    <map_components component_2="Na_leak_current" component_1="reversal_potentials"/>
    <map_variables variable_2="E_Na" variable_1="E_Na"/>
  </connection>
  
  <connection>
    <map_components component_2="SL_Ca_leak" component_1="reversal_potentials"/>
    <map_variables variable_2="E_Ca" variable_1="E_Ca"/>
  </connection>
  
  <connection>
    <map_components component_2="calcium_dependent_chloride_current" component_1="reversal_potentials"/>
    <map_variables variable_2="E_Cl" variable_1="E_Cl"/>
  </connection>
  
  <connection>
    <map_components component_2="background_chloride_current" component_1="reversal_potentials"/>
    <map_variables variable_2="E_Cl" variable_1="E_Cl"/>
  </connection>
  
  <connection>
    <map_components component_2="sodium_buffering" component_1="reversal_potentials"/>
    <map_variables variable_2="Nai" variable_1="Nai"/>
  </connection>
  
  <connection>
    <map_components component_2="calcium_buffering" component_1="reversal_potentials"/>
    <map_variables variable_2="Cai" variable_1="Cai"/>
  </connection>
  
  <connection>
    <map_components component_2="model_parameters" component_1="reversal_potentials"/>
    <map_variables variable_2="Ki" variable_1="Ki"/>
    <map_variables variable_2="Cli" variable_1="Cli"/>
    <map_variables variable_2="Clo" variable_1="Clo"/>
    <map_variables variable_2="Nao" variable_1="Nao"/>
    <map_variables variable_2="Cao" variable_1="Cao"/>
    <map_variables variable_2="Ko" variable_1="Ko"/>
    <map_variables variable_2="R" variable_1="R"/>
    <map_variables variable_2="T" variable_1="T"/>
    <map_variables variable_2="F" variable_1="F"/>
  </connection>
  
  <connection>
    <map_components component_2="L_type_calcium_current" component_1="model_parameters"/>
    <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="Cao" variable_1="Cao"/>
    <map_variables variable_2="Nao" variable_1="Nao"/>
    <map_variables variable_2="Ko" variable_1="Ko"/>
    <map_variables variable_2="K_SL" variable_1="K_SL"/>
  </connection>
  
  <connection>
    <map_components component_2="Na_Ca_exchange" component_1="model_parameters"/>
    <map_variables variable_2="Nao" variable_1="Nao"/>
    <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="Cao" variable_1="Cao"/>
  </connection>
  
  <connection>
    <map_components component_2="Na_K_pump" component_1="model_parameters"/>
    <map_variables variable_2="Ko" variable_1="Ko"/>
    <map_variables variable_2="Nao" variable_1="Nao"/>
    <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="rapid_delayed_rectifier_current" component_1="model_parameters"/>
    <map_variables variable_2="Ko" variable_1="Ko"/>
  </connection>
  
  <connection>
    <map_components component_2="slow_delayed_rectifier_current" component_1="model_parameters"/>
    <map_variables variable_2="Ko" variable_1="Ko"/>
    <map_variables variable_2="Ki" variable_1="Ki"/>
    <map_variables variable_2="Nao" variable_1="Nao"/>
    <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="inward_rectifier_potassium_current" component_1="model_parameters"/>
    <map_variables variable_2="Ko" variable_1="Ko"/>
  </connection>
  
  <connection>
    <map_components component_2="SL_Ca_pump" component_1="model_parameters"/>
    <map_variables variable_2="T" variable_1="T"/>
  </connection>
  
  <connection>
    <map_components component_2="SR_Ca_pump" component_1="model_parameters"/>
    <map_variables variable_2="T" variable_1="T"/>
  </connection>
  
  <connection>
    <map_components component_2="fast_sodium_current_m_gate" component_1="fast_sodium_current"/>
    <map_variables variable_2="m" variable_1="m"/>
    <map_variables variable_2="time" variable_1="time"/>
    <map_variables variable_2="V" variable_1="V"/>
  </connection>
  
  <connection>
    <map_components component_2="fast_sodium_current_h_gate" component_1="fast_sodium_current"/>
    <map_variables variable_2="h" variable_1="h"/>
    <map_variables variable_2="time" variable_1="time"/>
    <map_variables variable_2="V" variable_1="V"/>
  </connection>
  
  <connection>
    <map_components component_2="fast_sodium_current_j_gate" component_1="fast_sodium_current"/>
    <map_variables variable_2="j" variable_1="j"/>
    <map_variables variable_2="time" variable_1="time"/>
    <map_variables variable_2="V" variable_1="V"/>
  </connection>
  
  <connection>
    <map_components component_2="L_type_calcium_current_d_gate" component_1="L_type_calcium_current"/>
    <map_variables variable_2="d" variable_1="d"/>
    <map_variables variable_2="time" variable_1="time"/>
    <map_variables variable_2="V" variable_1="V"/>
  </connection>
  
  <connection>
    <map_components component_2="L_type_calcium_current_f_gate" component_1="L_type_calcium_current"/>
    <map_variables variable_2="f" variable_1="f"/>
    <map_variables variable_2="time" variable_1="time"/>
    <map_variables variable_2="V" variable_1="V"/>
  </connection>
  
  <connection>
    <map_components component_2="L_type_calcium_current_fCa_gate" component_1="L_type_calcium_current"/>
    <map_variables variable_2="fCa" variable_1="fCa"/>
    <map_variables variable_2="time" variable_1="time"/>
    <map_variables variable_2="Cai" variable_1="Cai"/>
  </connection>
  
  <connection>
    <map_components component_2="slow_transient_outward_current_X_gate" component_1="slow_transient_outward_current"/>
    <map_variables variable_2="X_tos" variable_1="X_tos"/>
    <map_variables variable_2="time" variable_1="time"/>
    <map_variables variable_2="V" variable_1="V"/>
  </connection>
  
  <connection>
    <map_components component_2="slow_transient_outward_current_Y_gate" component_1="slow_transient_outward_current"/>
    <map_variables variable_2="Y_tos" variable_1="Y_tos"/>
    <map_variables variable_2="time" variable_1="time"/>
    <map_variables variable_2="V" variable_1="V"/>
  </connection>
  
  <connection>
    <map_components component_2="slow_transient_outward_current_R_gate" component_1="slow_transient_outward_current"/>
    <map_variables variable_2="Rs_infinity" variable_1="Rs_infinity"/>
    <map_variables variable_2="V" variable_1="V"/>
  </connection>
  
  <connection>
    <map_components component_2="fast_transient_outward_current_X_gate" component_1="fast_transient_outward_current"/>
    <map_variables variable_2="X_tof" variable_1="X_tof"/>
    <map_variables variable_2="time" variable_1="time"/>
    <map_variables variable_2="V" variable_1="V"/>
  </connection>
  
  <connection>
    <map_components component_2="fast_transient_outward_current_Y_gate" component_1="fast_transient_outward_current"/>
    <map_variables variable_2="Y_tof" variable_1="Y_tof"/>
    <map_variables variable_2="time" variable_1="time"/>
    <map_variables variable_2="V" variable_1="V"/>
  </connection>
  
  <connection>
    <map_components component_2="slow_delayed_rectifier_current_Xs_gate" component_1="slow_delayed_rectifier_current"/>
    <map_variables variable_2="Xs" variable_1="Xs"/>
    <map_variables variable_2="time" variable_1="time"/>
    <map_variables variable_2="V" variable_1="V"/>
  </connection>
  
  <connection>
    <map_components component_2="rapid_delayed_rectifier_current_Xr_gate" component_1="rapid_delayed_rectifier_current"/>
    <map_variables variable_2="Xr" variable_1="Xr"/>
    <map_variables variable_2="time" variable_1="time"/>
    <map_variables variable_2="V" variable_1="V"/>
  </connection>
 
  <connection>
    <map_components component_2="rapid_delayed_rectifier_current_Rr_gate" component_1="rapid_delayed_rectifier_current"/>
    <map_variables variable_2="Rr" variable_1="Rr"/>
    <map_variables variable_2="V" variable_1="V"/>
  </connection>
  
  <connection>
    <map_components component_2="inward_rectifier_potassium_current_K1_gate" component_1="inward_rectifier_potassium_current"/>
    <map_variables variable_2="K1_infinity" variable_1="K1_infinity"/>
    <map_variables variable_2="V" variable_1="V"/>
    <map_variables variable_2="E_K" variable_1="E_K"/>
  </connection>
  
  <connection>
    <map_components component_2="sodium_buffering" component_1="ion_diffusion"/>
    <map_variables variable_2="Nai" variable_1="Nai"/>
    <map_variables variable_2="Na_jct" variable_1="Na_jct"/>
    <map_variables variable_2="Na_SL" variable_1="Na_SL"/>
  </connection>
  
  <connection>
    <map_components component_2="calcium_buffering" component_1="ion_diffusion"/>
    <map_variables variable_2="Cai" variable_1="Cai"/>
    <map_variables variable_2="Ca_jct" variable_1="Ca_jct"/>
    <map_variables variable_2="Ca_SL" variable_1="Ca_SL"/>
  </connection>
  
  <connection>
    <map_components component_2="calcium_buffering" component_1="calcium_dependent_chloride_current"/>
    <map_variables variable_2="Ca_jct" variable_1="Ca_jct"/>
  </connection>
  
  <connection>
    <map_components component_2="sodium_buffering" component_1="Na_K_pump"/>
    <map_variables variable_2="Nai" variable_1="Nai"/>
  </connection>
  
  <connection>
    <map_components component_2="calcium_buffering" component_1="slow_delayed_rectifier_current"/>
    <map_variables variable_2="Cai" variable_1="Cai"/>
  </connection>
  
  <connection>
    <map_components component_2="sodium_buffering" component_1="slow_delayed_rectifier_current"/>
    <map_variables variable_2="Na_jct" variable_1="Na_jct"/>
  </connection>
  
  <connection>
    <map_components component_2="calcium_buffering" component_1="L_type_calcium_current"/>
    <map_variables variable_2="Cai" variable_1="Cai"/>
    <map_variables variable_2="Ca_SL" variable_1="Ca_SL"/>
  </connection>
  
  <connection>
    <map_components component_2="sodium_buffering" component_1="L_type_calcium_current"/>
    <map_variables variable_2="Na_SL" variable_1="Na_SL"/>
  </connection>
  
  <connection>
    <map_components component_2="calcium_buffering" component_1="Na_Ca_exchange"/>
    <map_variables variable_2="Cai" variable_1="Cai"/>
  </connection>
  
  <connection>
    <map_components component_2="sodium_buffering" component_1="Na_Ca_exchange"/>
    <map_variables variable_2="Nai" variable_1="Nai"/>
  </connection>
  
  <connection>
    <map_components component_2="calcium_buffering" component_1="SR_Ca_pump"/>
    <map_variables variable_2="Cai" variable_1="Cai"/>
    <map_variables variable_2="Ca_SR" variable_1="Ca_SR"/>
  </connection>
  
  <connection>
    <map_components component_2="calcium_buffering" component_1="SL_Ca_pump"/>
    <map_variables variable_2="Cai" variable_1="Cai"/>
  </connection>
  
  <connection>
    <map_components component_2="calcium_buffering" component_1="SR_Ca_release"/>
    <map_variables variable_2="Ca_jct" variable_1="Ca_jct"/>
    <map_variables variable_2="Ca_SR" variable_1="Ca_SR"/>
  </connection>
  
  <connection>
    <map_components component_2="calcium_buffering" component_1="passive_SR_Ca_leak"/>
    <map_variables variable_2="Ca_jct" variable_1="Ca_jct"/>
    <map_variables variable_2="Ca_SR" variable_1="Ca_SR"/>
  </connection>
  



<rdf:RDF>
  <rdf:Bag rdf:about="rdf:#ea586037-0c5e-48e1-862a-23761cdca0f6">
    <rdf:li>calcium dynamics</rdf:li>
    <rdf:li>cardiac</rdf:li>
    <rdf:li>electrophysiology</rdf:li>
  </rdf:Bag>
  <rdf:Seq rdf:about="rdf:#29ed2752-78ea-468c-b924-b48407226c27">
    <rdf:li rdf:resource="rdf:#f26769e8-2e4d-480f-a1ca-173663ff4dd1"/>
    <rdf:li rdf:resource="rdf:#61af7017-bf22-4953-bb1e-77588988edfb"/>
    <rdf:li rdf:resource="rdf:#b7f1bf7b-58c8-4822-aeb5-f6a5a68bc182"/>
    <rdf:li rdf:resource="rdf:#c5026964-af2c-4773-ab3a-05828f472754"/>
    <rdf:li rdf:resource="rdf:#24afc4ca-6586-4e74-9945-b7ed5a7693aa"/>
  </rdf:Seq>
  <rdf:Description rdf:about="rdf:#3f565957-7611-44ae-8cea-3e8aee7a8ff5">
    <bqs:subject_type>keyword</bqs:subject_type>
    <rdf:value rdf:resource="rdf:#ea586037-0c5e-48e1-862a-23761cdca0f6"/>
  </rdf:Description>
  <rdf:Description rdf:about="">
    <dc:publisher>The University of Auckland, Bioengineering Institute</dc:publisher>
    <cmeta:comment rdf:resource="rdf:#0edaa7a9-211f-42b6-b147-bdc2e8c27850"/>
    <dc:title>
        Shannon et al.'s 2004 mathematical model of integrated calcium dynamics within the ventricular myocyte.
      </dc:title>
    <dcterms:created rdf:resource="rdf:#47077df2-adac-4f58-87b4-e25cbd70150b"/>
    <dc:creator rdf:resource="rdf:#9a7690b3-9535-40da-a231-36fb21604a53"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#9a7690b3-9535-40da-a231-36fb21604a53">
    <vCard:ORG rdf:resource="rdf:#ef712018-ef2f-496f-b24e-97b98ce2bc70"/>
    <vCard:EMAIL rdf:resource="rdf:#83c55a3c-a8f7-47a1-9844-30d2840f5808"/>
    <vCard:N rdf:resource="rdf:#7492848e-9caf-467b-b0df-16dd479a7d7f"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#83c55a3c-a8f7-47a1-9844-30d2840f5808">
    <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:#f26769e8-2e4d-480f-a1ca-173663ff4dd1">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#92f04414-3a5b-40ac-b31b-5e178333944c"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#2906e531-503f-4423-8c39-7ac4541081af">
    <vCard:Given>Donald</vCard:Given>
    <vCard:Family>Bers</vCard:Family>
    <vCard:Other>M</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#3c4761bc-e9d3-4b42-b2bd-658a80824d7e">
    <vCard:Given>Fei</vCard:Given>
    <vCard:Family>Wang</vCard:Family>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#47077df2-adac-4f58-87b4-e25cbd70150b">
    <dcterms:W3CDTF>2004-11-06T00:00:00+00:00</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#9c31686d-5971-4344-9735-ef70fcd15483">
    <dc:subject rdf:resource="rdf:#3f565957-7611-44ae-8cea-3e8aee7a8ff5"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#99dcff6d-6670-4d16-a593-29c17c64df80">
    <vCard:FN>Catherine Lloyd</vCard:FN>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#7492848e-9caf-467b-b0df-16dd479a7d7f">
    <vCard:Given>Catherine</vCard:Given>
    <vCard:Family>Lloyd</vCard:Family>
    <vCard:Other>May</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#9c784ead-e791-4746-8d4c-ed856e1ab0a3">
    <dc:title>Biophysical Journal</dc:title>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#7bad8a81-a80b-4431-8808-5ce1f476354f">
    <vCard:FN/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#ef712018-ef2f-496f-b24e-97b98ce2bc70">
    <vCard:Orgname>The University of Auckland</vCard:Orgname>
    <vCard:Orgunit>The Bioengineering Institute</vCard:Orgunit>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#f01b9baf-55e9-4a38-a20b-e7d14792d0b0">
    <dc:creator rdf:resource="rdf:#99dcff6d-6670-4d16-a593-29c17c64df80"/>
    <rdf:value>This is the CellML description of Shannon et al.'s 2004 mathematical model of integrated calcium dynamics within the ventricular myocyte.</rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#2d83f27a-bdff-4ef4-8d23-600997f4cf02">
    <dc:creator rdf:resource="rdf:#29ed2752-78ea-468c-b924-b48407226c27"/>
    <dc:title>A Mathematical Treatment of Integrated Ca Dynamics within the &#13;
            Ventricular Myocyte</dc:title>
    <bqs:volume>87</bqs:volume>
    <bqs:first_page>3351</bqs:first_page>
    <bqs:Journal rdf:resource="rdf:#9c784ead-e791-4746-8d4c-ed856e1ab0a3"/>
    <dcterms:issued rdf:resource="rdf:#3cf257e7-49c3-43bc-ada2-e6be4d3d9014"/>
    <bqs:last_page>3371</bqs:last_page>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#61af7017-bf22-4953-bb1e-77588988edfb">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#3c4761bc-e9d3-4b42-b2bd-658a80824d7e"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#c5026964-af2c-4773-ab3a-05828f472754">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#ecb45b6b-2533-469c-978c-0e98e30b7138"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#92f04414-3a5b-40ac-b31b-5e178333944c">
    <vCard:Given>Thomas</vCard:Given>
    <vCard:Family>Shannon</vCard:Family>
    <vCard:Other>R</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#e0c04a1a-a056-4fb3-b24a-627a35d34a59">
    <bqs:Pubmed_id>15347581</bqs:Pubmed_id>
    <bqs:JournalArticle rdf:resource="rdf:#2d83f27a-bdff-4ef4-8d23-600997f4cf02"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#0edaa7a9-211f-42b6-b147-bdc2e8c27850">
    <dc:creator rdf:resource="rdf:#7bad8a81-a80b-4431-8808-5ce1f476354f"/>
    <rdf:value/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#b7f1bf7b-58c8-4822-aeb5-f6a5a68bc182">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#8fb664ff-3cda-4faf-8497-03f912a87a28"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#8fb664ff-3cda-4faf-8497-03f912a87a28">
    <vCard:Given>Jose</vCard:Given>
    <vCard:Family>Puglisi</vCard:Family>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#3cf257e7-49c3-43bc-ada2-e6be4d3d9014">
    <dcterms:W3CDTF>2004-11</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#24afc4ca-6586-4e74-9945-b7ed5a7693aa">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#2906e531-503f-4423-8c39-7ac4541081af"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#ecb45b6b-2533-469c-978c-0e98e30b7138">
    <vCard:Given>Christopher</vCard:Given>
    <vCard:Family>Weber</vCard:Family>
  </rdf:Description>
  <rdf:Description rdf:about="#shannon_wang_puglisi_weber_bers_2004_version01">
    <bqs:reference rdf:resource="rdf:#9c31686d-5971-4344-9735-ef70fcd15483"/>
    <bqs:reference rdf:resource="rdf:#e0c04a1a-a056-4fb3-b24a-627a35d34a59"/>
    <cmeta:comment rdf:resource="rdf:#f01b9baf-55e9-4a38-a20b-e7d14792d0b0"/>
    <cmeta:bio_entity>ventricular myocyte</cmeta:bio_entity>
  </rdf:Description>
</rdf:RDF>
</model>