Location: Noble, 2000 @ a8221217219a / noble_2000_b.cellml

Author:
pmr2.import <nobody@models.cellml.org>
Date:
2008-11-24 20:21:26+13:00
Desc:
committing version01 of noble_2000
Permanent Source URI:
https://staging.physiomeproject.org/workspace/noble_2000/rawfile/a8221217219a10681c57d8faffc134c01a27962d/noble_2000_b.cellml

<?xml version='1.0' encoding='utf-8'?>
<!--
This CellML file was generated on 22/09/2008 at 11:46:10 using:

COR (0.9.31.1033)
Copyright 2002-2008 Dr Alan Garny
http://COR.physiol.ox.ac.uk/ - COR@physiol.ox.ac.uk

CellML 1.0 was used to generate this model
http://www.CellML.org/
--><model xmlns="http://www.cellml.org/cellml/1.0#" xmlns:cellml="http://www.cellml.org/cellml/1.0#" xmlns:cmeta="http://www.cellml.org/metadata/1.0#" name="noble_model_2000" cmeta:id="noble_model_2000">

<documentation xmlns="http://cellml.org/tmp-documentation">
		<article>
			<articleinfo>
				<title>A model of the guinea-pig ventricular cardiac myocyte incorporating a transverse-axial tubular system</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 CellML model is known to run in both COR and PCEnv to recreate the original results.  The units have been checked and they are consistent. This particular version of the model has multiple stimuli between 200s and 600s.         </para>
			</section>
			<sect1 id="sec_structure">
				<title>Model Structure</title>
				<para>
This CellML model is a description of the Noble 2000 (N2K) guinea pig ventricular cell model. N2K is an updated version of the model published by Noble et al. in 1998, in which the voltage-sensitive L-type Ca2+ current (ICaL) and the Ca2+ dynamics have been refined.  For more information about this model please refer to the description of the <ulink url="http://www.cellml.org/models/noble_varghese_kohl_noble_1998_version08">original 1998 model</ulink> in the CellML Model Repository.
</para>

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

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

<para>
Paper reference is in which the current model is described:
</para>

<para>
Cellular Open Resource (COR): current status and future directions. Alan Garny, Denis Noble, Peter J. Hunter and Peter Kohl, 2009, <emphasis>Philisophical Transactions of the Royal Society</emphasis> (accepted)
</para>

<informalfigure float="0" id="fig_cell_diagram">
<mediaobject>
  <imageobject>
    <imagedata fileref="noble_1998.png"/>
  </imageobject>
</mediaobject>
<caption>A schematic diagram describing the current flows across the cell membrane that are captured in the Noble'98 model.</caption>
</informalfigure>
				
			</sect1>
		</article>
	</documentation>

   <units name="micrometre">
      <unit units="metre" prefix="micro"/>
   </units>
   <units name="per_second">
      <unit units="second" exponent="-1"/>
   </units>
   <units name="millivolt">
      <unit units="volt" prefix="milli"/>
   </units>
   <units name="per_millivolt">
      <unit units="volt" prefix="milli" exponent="-1"/>
   </units>
   <units name="per_millivolt_second">
      <unit units="millivolt" exponent="-1"/>
      <unit units="second" exponent="-1"/>
   </units>
   <units name="microS">
      <unit units="siemens" prefix="micro"/>
   </units>
   <units name="microF">
      <unit units="farad" prefix="micro"/>
   </units>
   <units name="nanoA">
      <unit units="ampere" prefix="nano"/>
   </units>
   <units name="mA_nA">
      <unit units="ampere" prefix="milli"/>
      <unit units="nanoA" exponent="-1"/>
   </units>
   <units name="nanoA_per_millimolar">
      <unit units="ampere" prefix="nano"/>
      <unit units="millimolar" exponent="-1"/>
   </units>
   <units name="millimolar">
      <unit units="mole" prefix="milli"/>
      <unit units="litre" exponent="-1"/>
   </units>
   <units name="millimolar4">
      <unit units="millimolar" exponent="4"/>
   </units>
   <units name="millimolar_per_second">
      <unit units="millimolar"/>
      <unit units="second" exponent="-1"/>
   </units>
   <units name="per_millimolar_second">
      <unit units="millimolar" exponent="-1"/>
      <unit units="second" exponent="-1"/>
   </units>
   <units name="joule_per_kilomole_kelvin">
      <unit units="joule"/>
      <unit units="mole" prefix="kilo" exponent="-1"/>
      <unit units="kelvin" exponent="-1"/>
   </units>
   <units name="coulomb_per_mole">
      <unit units="coulomb"/>
      <unit units="mole" exponent="-1"/>
   </units>
   <units name="micrometre3">
      <unit units="metre" prefix="micro" exponent="3"/>
   </units>
   <units name="litre_micrometre3">
      <unit units="litre"/>
      <unit units="micrometre3" exponent="-1"/>
   </units>
   <component name="environment">
      <variable name="time" units="second" public_interface="out"/>
   </component>
   <component name="membrane">
      <variable name="V" units="millivolt" initial_value="-92.8915042" public_interface="out"/>
      <variable name="R" units="joule_per_kilomole_kelvin" initial_value="8314.472" public_interface="out"/>
      <variable name="T" units="kelvin" initial_value="310" public_interface="out"/>
      <variable name="F" units="coulomb_per_mole" initial_value="96485.3415" public_interface="out"/>
      <variable name="Cm" units="microF" initial_value="9.5e-5"/>
      <variable name="time" units="second" public_interface="in"/>
      <variable name="i_K1" units="nanoA" public_interface="in"/>
      <variable name="i_to" units="nanoA" public_interface="in"/>
      <variable name="i_Kr" units="nanoA" public_interface="in"/>
      <variable name="i_Ks" units="nanoA" public_interface="in"/>
      <variable name="i_Ca_L_K_cyt" units="nanoA" public_interface="in"/>
      <variable name="i_Ca_L_K_ds" units="nanoA" public_interface="in"/>
      <variable name="i_NaK" units="nanoA" public_interface="in"/>
      <variable name="i_Na" units="nanoA" public_interface="in"/>
      <variable name="i_b_Na" units="nanoA" public_interface="in"/>
      <variable name="i_p_Na" units="nanoA" public_interface="in"/>
      <variable name="i_Ca_L_Na_cyt" units="nanoA" public_interface="in"/>
      <variable name="i_Ca_L_Na_ds" units="nanoA" public_interface="in"/>
      <variable name="i_NaCa_cyt" units="nanoA" public_interface="in"/>
      <variable name="i_NaCa_ds" units="nanoA" public_interface="in"/>
      <variable name="i_Ca_L_Ca_cyt" units="nanoA" public_interface="in"/>
      <variable name="i_Ca_L_Ca_ds" units="nanoA" public_interface="in"/>
      <variable name="i_b_Ca" units="nanoA" public_interface="in"/>
      <variable name="i_Stim" units="nanoA"/>
      <variable name="stim_start" units="second" initial_value="200"/>
      <variable name="stim_end" units="second" initial_value="600"/>
      <variable name="stim_period" units="second" initial_value="1"/>
      <variable name="stim_duration" units="second" initial_value="0.002"/>
      <variable name="stim_amplitude" units="nanoA" initial_value="-6"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_Stim</ci>
            <piecewise>
               <piece>
                  <ci>stim_amplitude</ci>
                  <apply>
                     <and/>
                     <apply>
                        <geq/>
                        <ci>time</ci>
                        <ci>stim_start</ci>
                     </apply>
                     <apply>
                        <leq/>
                        <ci>time</ci>
                        <ci>stim_end</ci>
                     </apply>
                     <apply>
                        <leq/>
                        <apply>
                           <minus/>
                           <apply>
                              <minus/>
                              <ci>time</ci>
                              <ci>stim_start</ci>
                           </apply>
                           <apply>
                              <times/>
                              <apply>
                                 <floor/>
                                 <apply>
                                    <divide/>
                                    <apply>
                                       <minus/>
                                       <ci>time</ci>
                                       <ci>stim_start</ci>
                                    </apply>
                                    <ci>stim_period</ci>
                                 </apply>
                              </apply>
                              <ci>stim_period</ci>
                           </apply>
                        </apply>
                        <ci>stim_duration</ci>
                     </apply>
                  </apply>
               </piece>
               <otherwise>
                  <cn cellml:units="nanoA">0</cn>
               </otherwise>
            </piecewise>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>V</ci>
            </apply>
            <apply>
               <times/>
               <apply>
                  <divide/>
                  <apply>
                     <minus/>
                     <cn cellml:units="dimensionless">1</cn>
                  </apply>
                  <ci>Cm</ci>
               </apply>
               <apply>
                  <plus/>
                  <ci>i_Stim</ci>
                  <ci>i_K1</ci>
                  <ci>i_to</ci>
                  <ci>i_Kr</ci>
                  <ci>i_Ks</ci>
                  <ci>i_NaK</ci>
                  <ci>i_Na</ci>
                  <ci>i_b_Na</ci>
                  <ci>i_p_Na</ci>
                  <ci>i_Ca_L_Na_cyt</ci>
                  <ci>i_Ca_L_Na_ds</ci>
                  <ci>i_NaCa_cyt</ci>
                  <ci>i_NaCa_ds</ci>
                  <ci>i_Ca_L_Ca_cyt</ci>
                  <ci>i_Ca_L_Ca_ds</ci>
                  <ci>i_Ca_L_K_cyt</ci>
                  <ci>i_Ca_L_K_ds</ci>
                  <ci>i_b_Ca</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="reversal_potentials">
      <variable name="E_Na" units="millivolt" public_interface="out"/>
      <variable name="E_K" units="millivolt" public_interface="out"/>
      <variable name="E_Ks" units="millivolt" public_interface="out"/>
      <variable name="E_Ca" units="millivolt" public_interface="out"/>
      <variable name="E_mh" units="millivolt" public_interface="out"/>
      <variable name="P_kna" units="dimensionless" initial_value="0.03"/>
      <variable name="K_o" units="millimolar" public_interface="in"/>
      <variable name="Na_o" units="millimolar" public_interface="in"/>
      <variable name="K_i" units="millimolar" public_interface="in"/>
      <variable name="Na_i" units="millimolar" public_interface="in"/>
      <variable name="Ca_o" units="millimolar" public_interface="in"/>
      <variable name="Ca_i" units="millimolar" public_interface="in"/>
      <variable name="R" units="joule_per_kilomole_kelvin" public_interface="in"/>
      <variable name="F" units="coulomb_per_mole" public_interface="in"/>
      <variable name="T" units="kelvin" public_interface="in"/>
      <variable name="time" units="second" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>E_Na</ci>
            <apply>
               <times/>
               <apply>
                  <divide/>
                  <apply>
                     <times/>
                     <ci>R</ci>
                     <ci>T</ci>
                  </apply>
                  <ci>F</ci>
               </apply>
               <apply>
                  <ln/>
                  <apply>
                     <divide/>
                     <ci>Na_o</ci>
                     <ci>Na_i</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>E_K</ci>
            <apply>
               <times/>
               <apply>
                  <divide/>
                  <apply>
                     <times/>
                     <ci>R</ci>
                     <ci>T</ci>
                  </apply>
                  <ci>F</ci>
               </apply>
               <apply>
                  <ln/>
                  <apply>
                     <divide/>
                     <ci>K_o</ci>
                     <ci>K_i</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>E_Ks</ci>
            <apply>
               <times/>
               <apply>
                  <divide/>
                  <apply>
                     <times/>
                     <ci>R</ci>
                     <ci>T</ci>
                  </apply>
                  <ci>F</ci>
               </apply>
               <apply>
                  <ln/>
                  <apply>
                     <divide/>
                     <apply>
                        <plus/>
                        <ci>K_o</ci>
                        <apply>
                           <times/>
                           <ci>P_kna</ci>
                           <ci>Na_o</ci>
                        </apply>
                     </apply>
                     <apply>
                        <plus/>
                        <ci>K_i</ci>
                        <apply>
                           <times/>
                           <ci>P_kna</ci>
                           <ci>Na_i</ci>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>E_Ca</ci>
            <apply>
               <times/>
               <apply>
                  <divide/>
                  <apply>
                     <times/>
                     <cn cellml:units="dimensionless">0.5</cn>
                     <ci>R</ci>
                     <ci>T</ci>
                  </apply>
                  <ci>F</ci>
               </apply>
               <apply>
                  <ln/>
                  <apply>
                     <divide/>
                     <ci>Ca_o</ci>
                     <ci>Ca_i</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>E_mh</ci>
            <apply>
               <times/>
               <apply>
                  <divide/>
                  <apply>
                     <times/>
                     <ci>R</ci>
                     <ci>T</ci>
                  </apply>
                  <ci>F</ci>
               </apply>
               <apply>
                  <ln/>
                  <apply>
                     <divide/>
                     <apply>
                        <plus/>
                        <ci>Na_o</ci>
                        <apply>
                           <times/>
                           <cn cellml:units="dimensionless">0.12</cn>
                           <ci>K_o</ci>
                        </apply>
                     </apply>
                     <apply>
                        <plus/>
                        <ci>Na_i</ci>
                        <apply>
                           <times/>
                           <cn cellml:units="dimensionless">0.12</cn>
                           <ci>K_i</ci>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="time_independent_potassium_current">
      <variable name="i_K1" units="nanoA" public_interface="out"/>
      <variable name="K_mk1" units="millimolar" initial_value="10" public_interface="out"/>
      <variable name="g_K1" units="microS" initial_value="0.5"/>
      <variable name="K_o" units="millimolar" public_interface="in"/>
      <variable name="V" units="millivolt" public_interface="in"/>
      <variable name="E_K" units="millivolt" public_interface="in"/>
      <variable name="R" units="joule_per_kilomole_kelvin" public_interface="in"/>
      <variable name="F" units="coulomb_per_mole" public_interface="in"/>
      <variable name="T" units="kelvin" public_interface="in"/>
      <variable name="time" units="second" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_K1</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <ci>g_K1</ci>
                        <ci>K_o</ci>
                     </apply>
                     <apply>
                        <plus/>
                        <ci>K_o</ci>
                        <ci>K_mk1</ci>
                     </apply>
                  </apply>
                  <apply>
                     <minus/>
                     <ci>V</ci>
                     <ci>E_K</ci>
                  </apply>
               </apply>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <times/>
                           <apply>
                              <minus/>
                              <apply>
                                 <minus/>
                                 <ci>V</ci>
                                 <ci>E_K</ci>
                              </apply>
                              <cn cellml:units="millivolt">10</cn>
                           </apply>
                           <ci>F</ci>
                           <cn cellml:units="dimensionless">1.25</cn>
                        </apply>
                        <apply>
                           <times/>
                           <ci>R</ci>
                           <ci>T</ci>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="rapid_delayed_rectifier_potassium_current">
      <variable name="i_Kr" units="nanoA" public_interface="out"/>
      <variable name="g_Kr1" units="microS" initial_value="0.0028"/>
      <variable name="g_Kr2" units="microS" initial_value="0.0017"/>
      <variable name="E_K" units="millivolt" public_interface="in"/>
      <variable name="time" units="second" public_interface="in" private_interface="out"/>
      <variable name="V" units="millivolt" public_interface="in" private_interface="out"/>
      <variable name="xr1" units="dimensionless" private_interface="in"/>
      <variable name="xr2" units="dimensionless" private_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_Kr</ci>
            <apply>
               <times/>
               <apply>
                  <divide/>
                  <apply>
                     <times/>
                     <apply>
                        <plus/>
                        <apply>
                           <times/>
                           <ci>g_Kr1</ci>
                           <ci>xr1</ci>
                        </apply>
                        <apply>
                           <times/>
                           <ci>g_Kr2</ci>
                           <ci>xr2</ci>
                        </apply>
                     </apply>
                     <cn cellml:units="dimensionless">1</cn>
                  </apply>
                  <apply>
                     <plus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <plus/>
                              <ci>V</ci>
                              <cn cellml:units="millivolt">9</cn>
                           </apply>
                           <cn cellml:units="millivolt">22.4</cn>
                        </apply>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>E_K</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="rapid_delayed_rectifier_potassium_current_xr1_gate">
      <variable name="xr1" units="dimensionless" initial_value="1.02e-5" public_interface="out"/>
      <variable name="V" units="millivolt" public_interface="in"/>
      <variable name="alpha_xr1" units="per_second"/>
      <variable name="beta_xr1" units="per_second"/>
      <variable name="time" units="second" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>alpha_xr1</ci>
            <apply>
               <divide/>
               <cn cellml:units="per_second">50</cn>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <minus/>
                           <apply>
                              <minus/>
                              <ci>V</ci>
                              <cn cellml:units="millivolt">5</cn>
                           </apply>
                        </apply>
                        <cn cellml:units="millivolt">9</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>beta_xr1</ci>
            <apply>
               <times/>
               <cn cellml:units="per_second">0.05</cn>
               <apply>
                  <exp/>
                  <apply>
                     <divide/>
                     <apply>
                        <minus/>
                        <apply>
                           <minus/>
                           <ci>V</ci>
                           <cn cellml:units="millivolt">20</cn>
                        </apply>
                     </apply>
                     <cn cellml:units="millivolt">15</cn>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>xr1</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>alpha_xr1</ci>
                  <apply>
                     <minus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <ci>xr1</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>beta_xr1</ci>
                  <ci>xr1</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="rapid_delayed_rectifier_potassium_current_xr2_gate">
      <variable name="xr2" units="dimensionless" initial_value="2e-7" public_interface="out"/>
      <variable name="V" units="millivolt" public_interface="in"/>
      <variable name="alpha_xr2" units="per_second"/>
      <variable name="beta_xr2" units="per_second"/>
      <variable name="time" units="second" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>alpha_xr2</ci>
            <apply>
               <divide/>
               <cn cellml:units="per_second">50</cn>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <minus/>
                           <apply>
                              <minus/>
                              <ci>V</ci>
                              <cn cellml:units="millivolt">5</cn>
                           </apply>
                        </apply>
                        <cn cellml:units="millivolt">9</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>beta_xr2</ci>
            <apply>
               <times/>
               <cn cellml:units="per_second">0.4</cn>
               <apply>
                  <exp/>
                  <apply>
                     <minus/>
                     <apply>
                        <power/>
                        <apply>
                           <divide/>
                           <apply>
                              <plus/>
                              <ci>V</ci>
                              <cn cellml:units="millivolt">30</cn>
                           </apply>
                           <cn cellml:units="millivolt">30</cn>
                        </apply>
                        <cn cellml:units="dimensionless">3</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>xr2</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>alpha_xr2</ci>
                  <apply>
                     <minus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <ci>xr2</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>beta_xr2</ci>
                  <ci>xr2</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="slow_delayed_rectifier_potassium_current">
      <variable name="i_Ks" units="nanoA" public_interface="out"/>
      <variable name="g_Ks" units="microS" initial_value="0.0016"/>
      <variable name="E_Ks" units="millivolt" public_interface="in"/>
      <variable name="time" units="second" public_interface="in" private_interface="out"/>
      <variable name="V" units="millivolt" public_interface="in" private_interface="out"/>
      <variable name="xs" units="dimensionless" private_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_Ks</ci>
            <apply>
               <times/>
               <ci>g_Ks</ci>
               <apply>
                  <power/>
                  <ci>xs</ci>
                  <cn cellml:units="dimensionless">2</cn>
               </apply>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>E_Ks</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="slow_delayed_rectifier_potassium_current_xs_gate">
      <variable name="xs" units="dimensionless" initial_value="0.0006469" public_interface="out"/>
      <variable name="V" units="millivolt" public_interface="in"/>
      <variable name="alpha_xs" units="per_second"/>
      <variable name="beta_xs" units="per_second"/>
      <variable name="time" units="second" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>alpha_xs</ci>
            <apply>
               <divide/>
               <cn cellml:units="per_second">14</cn>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <minus/>
                           <apply>
                              <minus/>
                              <ci>V</ci>
                              <cn cellml:units="millivolt">40</cn>
                           </apply>
                        </apply>
                        <cn cellml:units="millivolt">9</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>beta_xs</ci>
            <apply>
               <times/>
               <cn cellml:units="per_second">1</cn>
               <apply>
                  <exp/>
                  <apply>
                     <divide/>
                     <apply>
                        <minus/>
                        <ci>V</ci>
                     </apply>
                     <cn cellml:units="millivolt">45</cn>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>xs</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>alpha_xs</ci>
                  <apply>
                     <minus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <ci>xs</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>beta_xs</ci>
                  <ci>xs</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="fast_sodium_current">
      <variable name="i_Na" units="nanoA" public_interface="out"/>
      <variable name="g_Na" units="microS" initial_value="0.5"/>
      <variable name="E_mh" units="millivolt" public_interface="in"/>
      <variable name="time" units="second" public_interface="in" private_interface="out"/>
      <variable name="V" units="millivolt" public_interface="in" private_interface="out"/>
      <variable name="m" units="dimensionless" private_interface="in"/>
      <variable name="h" units="dimensionless" private_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_Na</ci>
            <apply>
               <times/>
               <ci>g_Na</ci>
               <apply>
                  <power/>
                  <ci>m</ci>
                  <cn cellml:units="dimensionless">3</cn>
               </apply>
               <ci>h</ci>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>E_mh</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="fast_sodium_current_m_gate">
      <variable name="m" units="dimensionless" initial_value="0.0016111" public_interface="out"/>
      <variable name="alpha_m" units="per_second"/>
      <variable name="beta_m" units="per_second"/>
      <variable name="V" units="millivolt" public_interface="in"/>
      <variable name="time" units="second" public_interface="in"/>
      <variable name="delta_m" units="millivolt" initial_value="1e-5"/>
      <variable name="E0_m" units="millivolt"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>E0_m</ci>
            <apply>
               <plus/>
               <ci>V</ci>
               <cn cellml:units="millivolt">41</cn>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>alpha_m</ci>
            <piecewise>
               <piece>
                  <cn cellml:units="per_second">2000</cn>
                  <apply>
                     <lt/>
                     <apply>
                        <abs/>
                        <ci>E0_m</ci>
                     </apply>
                     <ci>delta_m</ci>
                  </apply>
               </piece>
               <otherwise>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <cn cellml:units="per_millivolt_second">200</cn>
                        <ci>E0_m</ci>
                     </apply>
                     <apply>
                        <minus/>
                        <cn cellml:units="dimensionless">1</cn>
                        <apply>
                           <exp/>
                           <apply>
                              <times/>
                              <apply>
                                 <minus/>
                                 <cn cellml:units="per_millivolt">0.1</cn>
                              </apply>
                              <ci>E0_m</ci>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </otherwise>
            </piecewise>
         </apply>
         <apply>
            <eq/>
            <ci>beta_m</ci>
            <apply>
               <times/>
               <cn cellml:units="per_second">8000</cn>
               <apply>
                  <exp/>
                  <apply>
                     <times/>
                     <apply>
                        <minus/>
                        <cn cellml:units="per_millivolt">0.056</cn>
                     </apply>
                     <apply>
                        <plus/>
                        <ci>V</ci>
                        <cn cellml:units="millivolt">66</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>m</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>alpha_m</ci>
                  <apply>
                     <minus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <ci>m</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>beta_m</ci>
                  <ci>m</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="fast_sodium_current_h_gate">
      <variable name="h" units="dimensionless" initial_value="0.9944559" public_interface="out"/>
      <variable name="alpha_h" units="per_second"/>
      <variable name="beta_h" units="per_second"/>
      <variable name="shift_h" units="millivolt" initial_value="0"/>
      <variable name="V" units="millivolt" public_interface="in"/>
      <variable name="time" units="second" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>alpha_h</ci>
            <apply>
               <times/>
               <cn cellml:units="per_second">20</cn>
               <apply>
                  <exp/>
                  <apply>
                     <times/>
                     <apply>
                        <minus/>
                        <cn cellml:units="per_millivolt">0.125</cn>
                     </apply>
                     <apply>
                        <minus/>
                        <apply>
                           <plus/>
                           <ci>V</ci>
                           <cn cellml:units="millivolt">75</cn>
                        </apply>
                        <ci>shift_h</ci>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>beta_h</ci>
            <apply>
               <divide/>
               <cn cellml:units="per_second">2000</cn>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <times/>
                     <cn cellml:units="dimensionless">320</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <times/>
                           <apply>
                              <minus/>
                              <cn cellml:units="per_millivolt">0.1</cn>
                           </apply>
                           <apply>
                              <minus/>
                              <apply>
                                 <plus/>
                                 <ci>V</ci>
                                 <cn cellml:units="millivolt">75</cn>
                              </apply>
                              <ci>shift_h</ci>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>h</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>alpha_h</ci>
                  <apply>
                     <minus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <ci>h</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>beta_h</ci>
                  <ci>h</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="persistent_sodium_current">
      <variable name="i_p_Na" units="nanoA" public_interface="out"/>
      <variable name="g_pna" units="microS" initial_value="0.0027"/>
      <variable name="E_Na" units="millivolt" public_interface="in"/>
      <variable name="time" units="second" public_interface="in"/>
      <variable name="V" units="millivolt" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_p_Na</ci>
            <apply>
               <times/>
               <apply>
                  <divide/>
                  <apply>
                     <times/>
                     <ci>g_pna</ci>
                     <cn cellml:units="dimensionless">1</cn>
                  </apply>
                  <apply>
                     <plus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <minus/>
                              <apply>
                                 <plus/>
                                 <ci>V</ci>
                                 <cn cellml:units="millivolt">52</cn>
                              </apply>
                           </apply>
                           <cn cellml:units="millivolt">8</cn>
                        </apply>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>E_Na</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="sodium_background_current">
      <variable name="i_b_Na" units="nanoA" public_interface="out"/>
      <variable name="g_bna" units="microS" initial_value="0.0006"/>
      <variable name="Na_i" units="millimolar" public_interface="in"/>
      <variable name="E_Na" units="millivolt" public_interface="in"/>
      <variable name="time" units="second" public_interface="in"/>
      <variable name="V" units="millivolt" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_b_Na</ci>
            <apply>
               <times/>
               <ci>g_bna</ci>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>E_Na</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="L_type_Ca_channel">
      <variable name="i_Ca_L_Ca_cyt" units="nanoA" public_interface="out"/>
      <variable name="i_Ca_L_K_cyt" units="nanoA" public_interface="out"/>
      <variable name="i_Ca_L_Na_cyt" units="nanoA" public_interface="out"/>
      <variable name="i_Ca_L_Ca_ds" units="nanoA" public_interface="out"/>
      <variable name="i_Ca_L_K_ds" units="nanoA" public_interface="out"/>
      <variable name="i_Ca_L_Na_ds" units="nanoA" public_interface="out"/>
      <variable name="i_Ca_L" units="nanoA" public_interface="out"/>
      <variable name="P_Ca_L" units="nanoA_per_millimolar" initial_value="1"/>
      <variable name="P_CaK" units="dimensionless" initial_value="0.002"/>
      <variable name="P_CaNa" units="dimensionless" initial_value="0.01"/>
      <variable name="Ca_o" units="millimolar" public_interface="in"/>
      <variable name="Ca_i" units="millimolar" public_interface="in" private_interface="out"/>
      <variable name="Ca_ds" units="millimolar" public_interface="in" private_interface="out"/>
      <variable name="K_o" units="millimolar" public_interface="in"/>
      <variable name="K_i" units="millimolar" public_interface="in"/>
      <variable name="Na_o" units="millimolar" public_interface="in"/>
      <variable name="Na_i" units="millimolar" public_interface="in"/>
      <variable name="R" units="joule_per_kilomole_kelvin" public_interface="in"/>
      <variable name="F" units="coulomb_per_mole" public_interface="in"/>
      <variable name="T" units="kelvin" public_interface="in"/>
      <variable name="time" units="second" public_interface="in" private_interface="out"/>
      <variable name="V" units="millivolt" public_interface="in" private_interface="out"/>
      <variable name="d" units="dimensionless" private_interface="in"/>
      <variable name="CaChoncyt" units="dimensionless" private_interface="in"/>
      <variable name="CaChonds" units="dimensionless" private_interface="in"/>
      <variable name="KCaChoff" units="millimolar" initial_value="0.01" private_interface="out"/>
      <variable name="Kmdsinact" units="millimolar" initial_value="0.001" private_interface="out"/>
      <variable name="FrICa" units="dimensionless" initial_value="0.7"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_Ca_L_Ca_cyt</ci>
            <apply>
               <times/>
               <apply>
                  <divide/>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <apply>
                           <minus/>
                           <cn cellml:units="dimensionless">1</cn>
                           <ci>FrICa</ci>
                        </apply>
                        <cn cellml:units="dimensionless">4</cn>
                        <ci>P_Ca_L</ci>
                        <ci>d</ci>
                        <ci>CaChoncyt</ci>
                        <apply>
                           <minus/>
                           <ci>V</ci>
                           <cn cellml:units="millivolt">50</cn>
                        </apply>
                        <ci>F</ci>
                     </apply>
                     <apply>
                        <times/>
                        <ci>R</ci>
                        <ci>T</ci>
                     </apply>
                  </apply>
                  <apply>
                     <minus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <times/>
                              <apply>
                                 <minus/>
                                 <apply>
                                    <minus/>
                                    <ci>V</ci>
                                    <cn cellml:units="millivolt">50</cn>
                                 </apply>
                              </apply>
                              <ci>F</ci>
                              <cn cellml:units="dimensionless">2</cn>
                           </apply>
                           <apply>
                              <times/>
                              <ci>R</ci>
                              <ci>T</ci>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <minus/>
                  <apply>
                     <times/>
                     <ci>Ca_i</ci>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <times/>
                              <cn cellml:units="millivolt">100</cn>
                              <ci>F</ci>
                           </apply>
                           <apply>
                              <times/>
                              <ci>R</ci>
                              <ci>T</ci>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
                  <apply>
                     <times/>
                     <ci>Ca_o</ci>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <times/>
                              <apply>
                                 <minus/>
                                 <apply>
                                    <minus/>
                                    <ci>V</ci>
                                    <cn cellml:units="millivolt">50</cn>
                                 </apply>
                              </apply>
                              <ci>F</ci>
                              <cn cellml:units="dimensionless">2</cn>
                           </apply>
                           <apply>
                              <times/>
                              <ci>R</ci>
                              <ci>T</ci>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>i_Ca_L_K_cyt</ci>
            <apply>
               <times/>
               <apply>
                  <divide/>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <apply>
                           <minus/>
                           <cn cellml:units="dimensionless">1</cn>
                           <ci>FrICa</ci>
                        </apply>
                        <ci>P_CaK</ci>
                        <ci>P_Ca_L</ci>
                        <ci>d</ci>
                        <ci>CaChoncyt</ci>
                        <apply>
                           <minus/>
                           <ci>V</ci>
                           <cn cellml:units="millivolt">50</cn>
                        </apply>
                        <ci>F</ci>
                     </apply>
                     <apply>
                        <times/>
                        <ci>R</ci>
                        <ci>T</ci>
                     </apply>
                  </apply>
                  <apply>
                     <minus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <times/>
                              <apply>
                                 <minus/>
                                 <apply>
                                    <minus/>
                                    <ci>V</ci>
                                    <cn cellml:units="millivolt">50</cn>
                                 </apply>
                              </apply>
                              <ci>F</ci>
                           </apply>
                           <apply>
                              <times/>
                              <ci>R</ci>
                              <ci>T</ci>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <minus/>
                  <apply>
                     <times/>
                     <ci>K_i</ci>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <times/>
                              <cn cellml:units="millivolt">50</cn>
                              <ci>F</ci>
                           </apply>
                           <apply>
                              <times/>
                              <ci>R</ci>
                              <ci>T</ci>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
                  <apply>
                     <times/>
                     <ci>K_o</ci>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <times/>
                              <apply>
                                 <minus/>
                                 <apply>
                                    <minus/>
                                    <ci>V</ci>
                                    <cn cellml:units="millivolt">50</cn>
                                 </apply>
                              </apply>
                              <ci>F</ci>
                           </apply>
                           <apply>
                              <times/>
                              <ci>R</ci>
                              <ci>T</ci>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>i_Ca_L_Na_cyt</ci>
            <apply>
               <times/>
               <apply>
                  <divide/>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <apply>
                           <minus/>
                           <cn cellml:units="dimensionless">1</cn>
                           <ci>FrICa</ci>
                        </apply>
                        <ci>P_CaNa</ci>
                        <ci>P_Ca_L</ci>
                        <ci>d</ci>
                        <ci>CaChoncyt</ci>
                        <apply>
                           <minus/>
                           <ci>V</ci>
                           <cn cellml:units="millivolt">50</cn>
                        </apply>
                        <ci>F</ci>
                     </apply>
                     <apply>
                        <times/>
                        <ci>R</ci>
                        <ci>T</ci>
                     </apply>
                  </apply>
                  <apply>
                     <minus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <times/>
                              <apply>
                                 <minus/>
                                 <apply>
                                    <minus/>
                                    <ci>V</ci>
                                    <cn cellml:units="millivolt">50</cn>
                                 </apply>
                              </apply>
                              <ci>F</ci>
                           </apply>
                           <apply>
                              <times/>
                              <ci>R</ci>
                              <ci>T</ci>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <minus/>
                  <apply>
                     <times/>
                     <ci>Na_i</ci>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <times/>
                              <cn cellml:units="millivolt">50</cn>
                              <ci>F</ci>
                           </apply>
                           <apply>
                              <times/>
                              <ci>R</ci>
                              <ci>T</ci>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
                  <apply>
                     <times/>
                     <ci>Na_o</ci>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <times/>
                              <apply>
                                 <minus/>
                                 <apply>
                                    <minus/>
                                    <ci>V</ci>
                                    <cn cellml:units="millivolt">50</cn>
                                 </apply>
                              </apply>
                              <ci>F</ci>
                           </apply>
                           <apply>
                              <times/>
                              <ci>R</ci>
                              <ci>T</ci>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>i_Ca_L_Ca_ds</ci>
            <apply>
               <times/>
               <apply>
                  <divide/>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <ci>FrICa</ci>
                        <cn cellml:units="dimensionless">4</cn>
                        <ci>P_Ca_L</ci>
                        <ci>d</ci>
                        <ci>CaChonds</ci>
                        <apply>
                           <minus/>
                           <ci>V</ci>
                           <cn cellml:units="millivolt">50</cn>
                        </apply>
                        <ci>F</ci>
                     </apply>
                     <apply>
                        <times/>
                        <ci>R</ci>
                        <ci>T</ci>
                     </apply>
                  </apply>
                  <apply>
                     <minus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <times/>
                              <apply>
                                 <minus/>
                                 <apply>
                                    <minus/>
                                    <ci>V</ci>
                                    <cn cellml:units="millivolt">50</cn>
                                 </apply>
                              </apply>
                              <ci>F</ci>
                              <cn cellml:units="dimensionless">2</cn>
                           </apply>
                           <apply>
                              <times/>
                              <ci>R</ci>
                              <ci>T</ci>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <minus/>
                  <apply>
                     <times/>
                     <ci>Ca_i</ci>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <times/>
                              <cn cellml:units="millivolt">100</cn>
                              <ci>F</ci>
                           </apply>
                           <apply>
                              <times/>
                              <ci>R</ci>
                              <ci>T</ci>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
                  <apply>
                     <times/>
                     <ci>Ca_o</ci>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <times/>
                              <apply>
                                 <minus/>
                                 <apply>
                                    <minus/>
                                    <ci>V</ci>
                                    <cn cellml:units="millivolt">50</cn>
                                 </apply>
                              </apply>
                              <ci>F</ci>
                              <cn cellml:units="dimensionless">2</cn>
                           </apply>
                           <apply>
                              <times/>
                              <ci>R</ci>
                              <ci>T</ci>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>i_Ca_L_K_ds</ci>
            <apply>
               <times/>
               <apply>
                  <divide/>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <ci>FrICa</ci>
                        <ci>P_CaK</ci>
                        <ci>P_Ca_L</ci>
                        <ci>d</ci>
                        <ci>CaChonds</ci>
                        <apply>
                           <minus/>
                           <ci>V</ci>
                           <cn cellml:units="millivolt">50</cn>
                        </apply>
                        <ci>F</ci>
                     </apply>
                     <apply>
                        <times/>
                        <ci>R</ci>
                        <ci>T</ci>
                     </apply>
                  </apply>
                  <apply>
                     <minus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <times/>
                              <apply>
                                 <minus/>
                                 <apply>
                                    <minus/>
                                    <ci>V</ci>
                                    <cn cellml:units="millivolt">50</cn>
                                 </apply>
                              </apply>
                              <ci>F</ci>
                           </apply>
                           <apply>
                              <times/>
                              <ci>R</ci>
                              <ci>T</ci>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <minus/>
                  <apply>
                     <times/>
                     <ci>K_i</ci>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <times/>
                              <cn cellml:units="millivolt">50</cn>
                              <ci>F</ci>
                           </apply>
                           <apply>
                              <times/>
                              <ci>R</ci>
                              <ci>T</ci>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
                  <apply>
                     <times/>
                     <ci>K_o</ci>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <times/>
                              <apply>
                                 <minus/>
                                 <apply>
                                    <minus/>
                                    <ci>V</ci>
                                    <cn cellml:units="millivolt">50</cn>
                                 </apply>
                              </apply>
                              <ci>F</ci>
                           </apply>
                           <apply>
                              <times/>
                              <ci>R</ci>
                              <ci>T</ci>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>i_Ca_L_Na_ds</ci>
            <apply>
               <times/>
               <apply>
                  <divide/>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <ci>FrICa</ci>
                        <ci>P_CaNa</ci>
                        <ci>P_Ca_L</ci>
                        <ci>d</ci>
                        <ci>CaChonds</ci>
                        <apply>
                           <minus/>
                           <ci>V</ci>
                           <cn cellml:units="millivolt">50</cn>
                        </apply>
                        <ci>F</ci>
                     </apply>
                     <apply>
                        <times/>
                        <ci>R</ci>
                        <ci>T</ci>
                     </apply>
                  </apply>
                  <apply>
                     <minus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <times/>
                              <apply>
                                 <minus/>
                                 <apply>
                                    <minus/>
                                    <ci>V</ci>
                                    <cn cellml:units="millivolt">50</cn>
                                 </apply>
                              </apply>
                              <ci>F</ci>
                           </apply>
                           <apply>
                              <times/>
                              <ci>R</ci>
                              <ci>T</ci>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <minus/>
                  <apply>
                     <times/>
                     <ci>Na_i</ci>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <times/>
                              <cn cellml:units="millivolt">50</cn>
                              <ci>F</ci>
                           </apply>
                           <apply>
                              <times/>
                              <ci>R</ci>
                              <ci>T</ci>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
                  <apply>
                     <times/>
                     <ci>Na_o</ci>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <times/>
                              <apply>
                                 <minus/>
                                 <apply>
                                    <minus/>
                                    <ci>V</ci>
                                    <cn cellml:units="millivolt">50</cn>
                                 </apply>
                              </apply>
                              <ci>F</ci>
                           </apply>
                           <apply>
                              <times/>
                              <ci>R</ci>
                              <ci>T</ci>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>i_Ca_L</ci>
            <apply>
               <plus/>
               <ci>i_Ca_L_Ca_cyt</ci>
               <ci>i_Ca_L_K_cyt</ci>
               <ci>i_Ca_L_Na_cyt</ci>
               <ci>i_Ca_L_Ca_ds</ci>
               <ci>i_Ca_L_K_ds</ci>
               <ci>i_Ca_L_Na_ds</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="L_type_Ca_channel_d_gate">
      <variable name="d" units="dimensionless" initial_value="0" public_interface="out"/>
      <variable name="alpha_d" units="per_second"/>
      <variable name="beta_d" units="per_second"/>
      <variable name="V" units="millivolt" public_interface="in"/>
      <variable name="time" units="second" public_interface="in"/>
      <variable name="E0_d" units="millivolt"/>
      <variable name="speed_d" units="dimensionless" initial_value="3"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>E0_d</ci>
            <apply>
               <minus/>
               <apply>
                  <plus/>
                  <ci>V</ci>
                  <cn cellml:units="millivolt">24</cn>
               </apply>
               <cn cellml:units="millivolt">5</cn>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>alpha_d</ci>
            <piecewise>
               <piece>
                  <cn cellml:units="per_second">120</cn>
                  <apply>
                     <lt/>
                     <apply>
                        <abs/>
                        <ci>E0_d</ci>
                     </apply>
                     <cn cellml:units="millivolt">0.0001</cn>
                  </apply>
               </piece>
               <otherwise>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <cn cellml:units="per_millivolt_second">30</cn>
                        <ci>E0_d</ci>
                     </apply>
                     <apply>
                        <minus/>
                        <cn cellml:units="dimensionless">1</cn>
                        <apply>
                           <exp/>
                           <apply>
                              <divide/>
                              <apply>
                                 <minus/>
                                 <ci>E0_d</ci>
                              </apply>
                              <cn cellml:units="millivolt">4</cn>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </otherwise>
            </piecewise>
         </apply>
         <apply>
            <eq/>
            <ci>beta_d</ci>
            <piecewise>
               <piece>
                  <cn cellml:units="per_second">120</cn>
                  <apply>
                     <lt/>
                     <apply>
                        <abs/>
                        <ci>E0_d</ci>
                     </apply>
                     <cn cellml:units="millivolt">0.0001</cn>
                  </apply>
               </piece>
               <otherwise>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <cn cellml:units="per_millivolt_second">12</cn>
                        <ci>E0_d</ci>
                     </apply>
                     <apply>
                        <minus/>
                        <apply>
                           <exp/>
                           <apply>
                              <divide/>
                              <ci>E0_d</ci>
                              <cn cellml:units="millivolt">10</cn>
                           </apply>
                        </apply>
                        <cn cellml:units="dimensionless">1</cn>
                     </apply>
                  </apply>
               </otherwise>
            </piecewise>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>d</ci>
            </apply>
            <apply>
               <times/>
               <ci>speed_d</ci>
               <apply>
                  <minus/>
                  <apply>
                     <times/>
                     <ci>alpha_d</ci>
                     <apply>
                        <minus/>
                        <cn cellml:units="dimensionless">1</cn>
                        <ci>d</ci>
                     </apply>
                  </apply>
                  <apply>
                     <times/>
                     <ci>beta_d</ci>
                     <ci>d</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="L_type_Ca_channel_f_gate">
      <variable name="f" units="dimensionless" initial_value="0"/>
      <variable name="CaChoncyt" units="dimensionless" public_interface="out"/>
      <variable name="CaChonds" units="dimensionless" public_interface="out"/>
      <variable name="CaChoffcyt" units="dimensionless"/>
      <variable name="CaChoffds" units="dimensionless"/>
      <variable name="alpha_f" units="per_second"/>
      <variable name="beta_f" units="per_second"/>
      <variable name="V" units="millivolt" public_interface="in"/>
      <variable name="Ca_i" units="millimolar" public_interface="in"/>
      <variable name="Ca_ds" units="millimolar" public_interface="in"/>
      <variable name="KCaChoff" units="millimolar" public_interface="in"/>
      <variable name="Kmdsinact" units="millimolar" public_interface="in"/>
      <variable name="time" units="second" public_interface="in"/>
      <variable name="speed_f" units="dimensionless" initial_value="0.3"/>
      <variable name="delta_f" units="millivolt" initial_value="0.0001"/>
      <variable name="E0_f" units="millivolt"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>E0_f</ci>
            <apply>
               <plus/>
               <ci>V</ci>
               <cn cellml:units="millivolt">34</cn>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>alpha_f</ci>
            <piecewise>
               <piece>
                  <cn cellml:units="per_second">25</cn>
                  <apply>
                     <lt/>
                     <apply>
                        <abs/>
                        <ci>E0_f</ci>
                     </apply>
                     <ci>delta_f</ci>
                  </apply>
               </piece>
               <otherwise>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <cn cellml:units="per_millivolt_second">6.25</cn>
                        <ci>E0_f</ci>
                     </apply>
                     <apply>
                        <minus/>
                        <apply>
                           <exp/>
                           <apply>
                              <divide/>
                              <ci>E0_f</ci>
                              <cn cellml:units="millivolt">4</cn>
                           </apply>
                        </apply>
                        <cn cellml:units="dimensionless">1</cn>
                     </apply>
                  </apply>
               </otherwise>
            </piecewise>
         </apply>
         <apply>
            <eq/>
            <ci>beta_f</ci>
            <apply>
               <divide/>
               <cn cellml:units="per_second">12</cn>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <times/>
                           <apply>
                              <minus/>
                              <cn cellml:units="dimensionless">1</cn>
                           </apply>
                           <apply>
                              <plus/>
                              <ci>V</ci>
                              <cn cellml:units="millivolt">34</cn>
                           </apply>
                        </apply>
                        <cn cellml:units="millivolt">4</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>CaChoffcyt</ci>
            <apply>
               <divide/>
               <ci>Ca_i</ci>
               <apply>
                  <plus/>
                  <ci>KCaChoff</ci>
                  <ci>Ca_i</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>f</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <apply>
                     <plus/>
                     <apply>
                        <times/>
                        <cn cellml:units="dimensionless">120</cn>
                        <apply>
                           <minus/>
                           <cn cellml:units="dimensionless">1</cn>
                           <ci>f</ci>
                        </apply>
                        <ci>CaChoffcyt</ci>
                     </apply>
                     <apply>
                        <times/>
                        <apply>
                           <minus/>
                           <cn cellml:units="dimensionless">1</cn>
                           <ci>f</ci>
                        </apply>
                        <apply>
                           <minus/>
                           <cn cellml:units="dimensionless">1</cn>
                           <ci>CaChoffcyt</ci>
                        </apply>
                     </apply>
                  </apply>
                  <ci>speed_f</ci>
                  <ci>beta_f</ci>
               </apply>
               <apply>
                  <times/>
                  <ci>alpha_f</ci>
                  <ci>f</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>CaChoncyt</ci>
            <apply>
               <times/>
               <apply>
                  <minus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <ci>f</ci>
               </apply>
               <apply>
                  <minus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <ci>CaChoffcyt</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>CaChoffds</ci>
            <apply>
               <divide/>
               <ci>Ca_ds</ci>
               <apply>
                  <plus/>
                  <ci>Kmdsinact</ci>
                  <ci>Ca_ds</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>CaChonds</ci>
            <apply>
               <times/>
               <apply>
                  <minus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <ci>f</ci>
               </apply>
               <apply>
                  <minus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <ci>CaChoffds</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="calcium_background_current">
      <variable name="i_b_Ca" units="nanoA" public_interface="out"/>
      <variable name="g_bca" units="microS" initial_value="0.00025"/>
      <variable name="E_Ca" units="millivolt" public_interface="in"/>
      <variable name="time" units="second" public_interface="in"/>
      <variable name="V" units="millivolt" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_b_Ca</ci>
            <apply>
               <times/>
               <ci>g_bca</ci>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>E_Ca</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="transient_outward_current">
      <variable name="i_to" units="nanoA" public_interface="out"/>
      <variable name="g_to" units="microS" initial_value="0.005"/>
      <variable name="g_tos" units="dimensionless" initial_value="0"/>
      <variable name="E_K" units="millivolt" public_interface="in"/>
      <variable name="V" units="millivolt" public_interface="in" private_interface="out"/>
      <variable name="time" units="second" public_interface="in" private_interface="out"/>
      <variable name="s" units="dimensionless" private_interface="in"/>
      <variable name="r" units="dimensionless" private_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_to</ci>
            <apply>
               <times/>
               <ci>g_to</ci>
               <apply>
                  <plus/>
                  <ci>g_tos</ci>
                  <apply>
                     <times/>
                     <ci>s</ci>
                     <apply>
                        <minus/>
                        <cn cellml:units="dimensionless">1</cn>
                        <ci>g_tos</ci>
                     </apply>
                  </apply>
               </apply>
               <ci>r</ci>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>E_K</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="transient_outward_current_s_gate">
      <variable name="s" units="dimensionless" initial_value="0.9948645" public_interface="out"/>
      <variable name="alpha_s" units="per_second"/>
      <variable name="beta_s" units="per_second"/>
      <variable name="V" units="millivolt" public_interface="in"/>
      <variable name="time" units="second" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>alpha_s</ci>
            <apply>
               <times/>
               <cn cellml:units="per_second">0.033</cn>
               <apply>
                  <exp/>
                  <apply>
                     <divide/>
                     <apply>
                        <minus/>
                        <ci>V</ci>
                     </apply>
                     <cn cellml:units="millivolt">17</cn>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>beta_s</ci>
            <apply>
               <divide/>
               <cn cellml:units="per_second">33</cn>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <times/>
                        <apply>
                           <minus/>
                           <cn cellml:units="per_millivolt">0.125</cn>
                        </apply>
                        <apply>
                           <plus/>
                           <ci>V</ci>
                           <cn cellml:units="millivolt">10</cn>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>s</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>alpha_s</ci>
                  <apply>
                     <minus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <ci>s</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>beta_s</ci>
                  <ci>s</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="transient_outward_current_r_gate">
      <variable name="r" units="dimensionless" initial_value="0" public_interface="out"/>
      <variable name="V" units="millivolt" public_interface="in"/>
      <variable name="time" units="second" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>r</ci>
            </apply>
            <apply>
               <times/>
               <cn cellml:units="per_second">333</cn>
               <apply>
                  <minus/>
                  <apply>
                     <divide/>
                     <cn cellml:units="dimensionless">1</cn>
                     <apply>
                        <plus/>
                        <cn cellml:units="dimensionless">1</cn>
                        <apply>
                           <exp/>
                           <apply>
                              <divide/>
                              <apply>
                                 <minus/>
                                 <apply>
                                    <plus/>
                                    <ci>V</ci>
                                    <cn cellml:units="millivolt">4</cn>
                                 </apply>
                              </apply>
                              <cn cellml:units="millivolt">5</cn>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
                  <ci>r</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="sodium_potassium_pump">
      <variable name="i_NaK" units="nanoA" public_interface="out"/>
      <variable name="i_NaK_max" units="nanoA" initial_value="0.7"/>
      <variable name="K_mK" units="millimolar" initial_value="1"/>
      <variable name="K_mNa" units="millimolar" initial_value="40"/>
      <variable name="K_o" units="millimolar" public_interface="in"/>
      <variable name="Na_i" units="millimolar" public_interface="in"/>
      <variable name="time" units="second" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_NaK</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <ci>i_NaK_max</ci>
                        <ci>K_o</ci>
                     </apply>
                     <apply>
                        <plus/>
                        <ci>K_mK</ci>
                        <ci>K_o</ci>
                     </apply>
                  </apply>
                  <ci>Na_i</ci>
               </apply>
               <apply>
                  <plus/>
                  <ci>K_mNa</ci>
                  <ci>Na_i</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="sodium_calcium_exchanger">
      <variable name="i_NaCa" units="nanoA" public_interface="out"/>
      <variable name="i_NaCa_cyt" units="nanoA" public_interface="out"/>
      <variable name="i_NaCa_ds" units="nanoA" public_interface="out"/>
      <variable name="k_NaCa" units="nanoA" initial_value="0.0005"/>
      <variable name="n_NaCa" units="dimensionless" initial_value="3"/>
      <variable name="d_NaCa" units="dimensionless" initial_value="0"/>
      <variable name="gamma" units="dimensionless" initial_value="0.2"/>
      <variable name="FRiNaCa" units="dimensionless" initial_value="0.001"/>
      <variable name="Na_i" units="millimolar" public_interface="in"/>
      <variable name="Na_o" units="millimolar" public_interface="in"/>
      <variable name="Ca_i" units="millimolar" public_interface="in"/>
      <variable name="Ca_ds" units="millimolar" public_interface="in"/>
      <variable name="Ca_o" units="millimolar" public_interface="in"/>
      <variable name="R" units="joule_per_kilomole_kelvin" public_interface="in"/>
      <variable name="F" units="coulomb_per_mole" public_interface="in"/>
      <variable name="T" units="kelvin" public_interface="in"/>
      <variable name="V" units="millivolt" public_interface="in"/>
      <variable name="time" units="second" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_NaCa_cyt</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <apply>
                     <minus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <ci>FRiNaCa</ci>
                  </apply>
                  <ci>k_NaCa</ci>
                  <apply>
                     <minus/>
                     <apply>
                        <times/>
                        <apply>
                           <exp/>
                           <apply>
                              <divide/>
                              <apply>
                                 <times/>
                                 <ci>gamma</ci>
                                 <apply>
                                    <minus/>
                                    <ci>n_NaCa</ci>
                                    <cn cellml:units="dimensionless">2</cn>
                                 </apply>
                                 <ci>V</ci>
                                 <ci>F</ci>
                              </apply>
                              <apply>
                                 <times/>
                                 <ci>R</ci>
                                 <ci>T</ci>
                              </apply>
                           </apply>
                        </apply>
                        <apply>
                           <power/>
                           <ci>Na_i</ci>
                           <ci>n_NaCa</ci>
                        </apply>
                        <ci>Ca_o</ci>
                     </apply>
                     <apply>
                        <times/>
                        <apply>
                           <exp/>
                           <apply>
                              <divide/>
                              <apply>
                                 <times/>
                                 <apply>
                                    <minus/>
                                    <ci>gamma</ci>
                                    <cn cellml:units="dimensionless">1</cn>
                                 </apply>
                                 <apply>
                                    <minus/>
                                    <ci>n_NaCa</ci>
                                    <cn cellml:units="dimensionless">2</cn>
                                 </apply>
                                 <ci>V</ci>
                                 <ci>F</ci>
                              </apply>
                              <apply>
                                 <times/>
                                 <ci>R</ci>
                                 <ci>T</ci>
                              </apply>
                           </apply>
                        </apply>
                        <apply>
                           <power/>
                           <ci>Na_o</ci>
                           <ci>n_NaCa</ci>
                        </apply>
                        <ci>Ca_i</ci>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <apply>
                     <plus/>
                     <cn cellml:units="millimolar4">1</cn>
                     <apply>
                        <times/>
                        <ci>d_NaCa</ci>
                        <apply>
                           <plus/>
                           <apply>
                              <times/>
                              <ci>Ca_i</ci>
                              <apply>
                                 <power/>
                                 <ci>Na_o</ci>
                                 <ci>n_NaCa</ci>
                              </apply>
                           </apply>
                           <apply>
                              <times/>
                              <ci>Ca_o</ci>
                              <apply>
                                 <power/>
                                 <ci>Na_i</ci>
                                 <ci>n_NaCa</ci>
                              </apply>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
                  <apply>
                     <plus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <apply>
                        <divide/>
                        <ci>Ca_i</ci>
                        <cn cellml:units="millimolar">0.0069</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>i_NaCa_ds</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <ci>FRiNaCa</ci>
                  <ci>k_NaCa</ci>
                  <apply>
                     <minus/>
                     <apply>
                        <times/>
                        <apply>
                           <exp/>
                           <apply>
                              <divide/>
                              <apply>
                                 <times/>
                                 <ci>gamma</ci>
                                 <apply>
                                    <minus/>
                                    <ci>n_NaCa</ci>
                                    <cn cellml:units="dimensionless">2</cn>
                                 </apply>
                                 <ci>V</ci>
                                 <ci>F</ci>
                              </apply>
                              <apply>
                                 <times/>
                                 <ci>R</ci>
                                 <ci>T</ci>
                              </apply>
                           </apply>
                        </apply>
                        <apply>
                           <power/>
                           <ci>Na_i</ci>
                           <ci>n_NaCa</ci>
                        </apply>
                        <ci>Ca_o</ci>
                     </apply>
                     <apply>
                        <times/>
                        <apply>
                           <exp/>
                           <apply>
                              <divide/>
                              <apply>
                                 <times/>
                                 <apply>
                                    <minus/>
                                    <ci>gamma</ci>
                                    <cn cellml:units="dimensionless">1</cn>
                                 </apply>
                                 <apply>
                                    <minus/>
                                    <ci>n_NaCa</ci>
                                    <cn cellml:units="dimensionless">2</cn>
                                 </apply>
                                 <ci>V</ci>
                                 <ci>F</ci>
                              </apply>
                              <apply>
                                 <times/>
                                 <ci>R</ci>
                                 <ci>T</ci>
                              </apply>
                           </apply>
                        </apply>
                        <apply>
                           <power/>
                           <ci>Na_o</ci>
                           <ci>n_NaCa</ci>
                        </apply>
                        <ci>Ca_ds</ci>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <apply>
                     <plus/>
                     <cn cellml:units="millimolar4">1</cn>
                     <apply>
                        <times/>
                        <ci>d_NaCa</ci>
                        <apply>
                           <plus/>
                           <apply>
                              <times/>
                              <ci>Ca_ds</ci>
                              <apply>
                                 <power/>
                                 <ci>Na_o</ci>
                                 <ci>n_NaCa</ci>
                              </apply>
                           </apply>
                           <apply>
                              <times/>
                              <ci>Ca_o</ci>
                              <apply>
                                 <power/>
                                 <ci>Na_i</ci>
                                 <ci>n_NaCa</ci>
                              </apply>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
                  <apply>
                     <plus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <apply>
                        <divide/>
                        <ci>Ca_ds</ci>
                        <cn cellml:units="millimolar">0.0069</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>i_NaCa</ci>
            <apply>
               <plus/>
               <ci>i_NaCa_cyt</ci>
               <ci>i_NaCa_ds</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="sarcoplasmic_reticulum_calcium_pump">
      <variable name="i_up" units="millimolar_per_second" public_interface="out"/>
      <variable name="K_1" units="dimensionless"/>
      <variable name="K_2" units="millimolar"/>
      <variable name="K_cyca" units="millimolar" initial_value="0.0003"/>
      <variable name="K_xcs" units="dimensionless" initial_value="0.4"/>
      <variable name="K_srca" units="millimolar" initial_value="0.5"/>
      <variable name="alpha_up" units="millimolar_per_second" initial_value="0.4"/>
      <variable name="beta_up" units="millimolar_per_second" initial_value="0.03"/>
      <variable name="Ca_i" units="millimolar" public_interface="in"/>
      <variable name="Ca_up" units="millimolar" public_interface="in"/>
      <variable name="time" units="second" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>K_1</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <ci>K_cyca</ci>
                  <ci>K_xcs</ci>
               </apply>
               <ci>K_srca</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>K_2</ci>
            <apply>
               <plus/>
               <ci>Ca_i</ci>
               <apply>
                  <times/>
                  <ci>Ca_up</ci>
                  <ci>K_1</ci>
               </apply>
               <apply>
                  <times/>
                  <ci>K_cyca</ci>
                  <ci>K_xcs</ci>
               </apply>
               <ci>K_cyca</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>i_up</ci>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <apply>
                     <divide/>
                     <ci>Ca_i</ci>
                     <ci>K_2</ci>
                  </apply>
                  <ci>alpha_up</ci>
               </apply>
               <apply>
                  <times/>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <ci>Ca_up</ci>
                        <ci>K_1</ci>
                     </apply>
                     <ci>K_2</ci>
                  </apply>
                  <ci>beta_up</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="calcium_translocation">
      <variable name="i_trans" units="millimolar_per_second" public_interface="out"/>
      <variable name="Ca_rel" units="millimolar" public_interface="in"/>
      <variable name="Ca_up" units="millimolar" public_interface="in"/>
      <variable name="time" units="second" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_trans</ci>
            <apply>
               <times/>
               <cn cellml:units="per_second">50</cn>
               <apply>
                  <minus/>
                  <ci>Ca_up</ci>
                  <ci>Ca_rel</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="calcium_release">
      <variable name="i_rel" units="millimolar_per_second" public_interface="out"/>
      <variable name="VoltDep" units="dimensionless"/>
      <variable name="RegBindSite" units="dimensionless"/>
      <variable name="CaiReg" units="dimensionless"/>
      <variable name="CadsReg" units="dimensionless"/>
      <variable name="ActRate" units="per_second"/>
      <variable name="InactRate" units="per_second"/>
      <variable name="K_leak_rate" units="per_second" initial_value="0.05"/>
      <variable name="K_m_rel" units="per_second" initial_value="250"/>
      <variable name="K_m_Ca_cyt" units="millimolar" initial_value="0.0005"/>
      <variable name="K_m_Ca_ds" units="millimolar" initial_value="0.01"/>
      <variable name="PrecFrac" units="dimensionless"/>
      <variable name="ActFrac" units="dimensionless" initial_value="0.0049039"/>
      <variable name="ProdFrac" units="dimensionless" initial_value="0.6950649"/>
      <variable name="SpeedRel" units="dimensionless"/>
      <variable name="time" units="second" public_interface="in"/>
      <variable name="V" units="millivolt" public_interface="in"/>
      <variable name="Ca_i" units="millimolar" public_interface="in"/>
      <variable name="Ca_ds" units="millimolar" public_interface="in"/>
      <variable name="Ca_rel" units="millimolar" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>VoltDep</ci>
            <apply>
               <exp/>
               <apply>
                  <times/>
                  <cn cellml:units="per_millivolt">0.08</cn>
                  <apply>
                     <minus/>
                     <ci>V</ci>
                     <cn cellml:units="millivolt">40</cn>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>CaiReg</ci>
            <apply>
               <divide/>
               <ci>Ca_i</ci>
               <apply>
                  <plus/>
                  <ci>Ca_i</ci>
                  <ci>K_m_Ca_cyt</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>CadsReg</ci>
            <apply>
               <divide/>
               <ci>Ca_ds</ci>
               <apply>
                  <plus/>
                  <ci>Ca_ds</ci>
                  <ci>K_m_Ca_ds</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>RegBindSite</ci>
            <apply>
               <plus/>
               <ci>CaiReg</ci>
               <apply>
                  <times/>
                  <apply>
                     <minus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <ci>CaiReg</ci>
                  </apply>
                  <ci>CadsReg</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>ActRate</ci>
            <apply>
               <plus/>
               <apply>
                  <times/>
                  <cn cellml:units="per_second">0</cn>
                  <ci>VoltDep</ci>
               </apply>
               <apply>
                  <times/>
                  <cn cellml:units="per_second">90</cn>
                  <apply>
                     <power/>
                     <ci>RegBindSite</ci>
                     <cn cellml:units="dimensionless">2</cn>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>InactRate</ci>
            <apply>
               <plus/>
               <cn cellml:units="per_second">60</cn>
               <apply>
                  <times/>
                  <cn cellml:units="per_second">500</cn>
                  <apply>
                     <power/>
                     <ci>RegBindSite</ci>
                     <cn cellml:units="dimensionless">2</cn>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>SpeedRel</ci>
            <piecewise>
               <piece>
                  <cn cellml:units="dimensionless">5</cn>
                  <apply>
                     <lt/>
                     <ci>V</ci>
                     <apply>
                        <minus/>
                        <cn cellml:units="millivolt">50</cn>
                     </apply>
                  </apply>
               </piece>
               <otherwise>
                  <cn cellml:units="dimensionless">1</cn>
               </otherwise>
            </piecewise>
         </apply>
         <apply>
            <eq/>
            <ci>PrecFrac</ci>
            <apply>
               <minus/>
               <apply>
                  <minus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <ci>ActFrac</ci>
               </apply>
               <ci>ProdFrac</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>ActFrac</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>PrecFrac</ci>
                  <ci>SpeedRel</ci>
                  <ci>ActRate</ci>
               </apply>
               <apply>
                  <times/>
                  <ci>ActFrac</ci>
                  <ci>SpeedRel</ci>
                  <ci>InactRate</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>ProdFrac</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>ActFrac</ci>
                  <ci>SpeedRel</ci>
                  <ci>InactRate</ci>
               </apply>
               <apply>
                  <times/>
                  <ci>SpeedRel</ci>
                  <cn cellml:units="per_second">1</cn>
                  <ci>ProdFrac</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>i_rel</ci>
            <apply>
               <times/>
               <apply>
                  <plus/>
                  <apply>
                     <times/>
                     <apply>
                        <power/>
                        <apply>
                           <divide/>
                           <ci>ActFrac</ci>
                           <apply>
                              <plus/>
                              <ci>ActFrac</ci>
                              <cn cellml:units="dimensionless">0.25</cn>
                           </apply>
                        </apply>
                        <cn cellml:units="dimensionless">2</cn>
                     </apply>
                     <ci>K_m_rel</ci>
                  </apply>
                  <ci>K_leak_rate</ci>
               </apply>
               <ci>Ca_rel</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="extracellular_sodium_concentration">
      <variable name="Na_o" units="millimolar" initial_value="140" public_interface="out"/>
   </component>
   <component name="intracellular_sodium_concentration">
      <variable name="Na_i" units="millimolar" initial_value="7.6713487" public_interface="out"/>
      <variable name="V_i" units="micrometre3" public_interface="in"/>
      <variable name="time" units="second" public_interface="in"/>
      <variable name="F" units="coulomb_per_mole" public_interface="in"/>
      <variable name="i_Na" units="nanoA" public_interface="in"/>
      <variable name="i_b_Na" units="nanoA" public_interface="in"/>
      <variable name="i_p_Na" units="nanoA" public_interface="in"/>
      <variable name="i_Ca_L_Na_cyt" units="nanoA" public_interface="in"/>
      <variable name="i_Ca_L_Na_ds" units="nanoA" public_interface="in"/>
      <variable name="i_NaK" units="nanoA" public_interface="in"/>
      <variable name="i_NaCa_cyt" units="nanoA" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>Na_i</ci>
            </apply>
            <apply>
               <times/>
               <apply>
                  <divide/>
                  <apply>
                     <minus/>
                     <cn cellml:units="mA_nA">1</cn>
                  </apply>
                  <apply>
                     <times/>
                     <cn cellml:units="litre_micrometre3">1</cn>
                     <ci>V_i</ci>
                     <ci>F</ci>
                  </apply>
               </apply>
               <apply>
                  <plus/>
                  <ci>i_Na</ci>
                  <ci>i_p_Na</ci>
                  <ci>i_b_Na</ci>
                  <apply>
                     <times/>
                     <cn cellml:units="dimensionless">3</cn>
                     <ci>i_NaK</ci>
                  </apply>
                  <apply>
                     <times/>
                     <cn cellml:units="dimensionless">3</cn>
                     <ci>i_NaCa_cyt</ci>
                  </apply>
                  <ci>i_Ca_L_Na_cyt</ci>
                  <ci>i_Ca_L_Na_ds</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="extracellular_calcium_concentration">
      <variable name="Ca_o" units="millimolar" initial_value="2" public_interface="out"/>
   </component>
   <component name="extracellular_potassium_concentration">
      <variable name="K_o" units="millimolar" initial_value="4" public_interface="out"/>
   </component>
   <component name="intracellular_potassium_concentration">
      <variable name="K_i" units="millimolar" initial_value="136.1745362" public_interface="out"/>
      <variable name="V_i" units="micrometre3" public_interface="in"/>
      <variable name="i_K1" units="nanoA" public_interface="in"/>
      <variable name="i_Kr" units="nanoA" public_interface="in"/>
      <variable name="i_Ks" units="nanoA" public_interface="in"/>
      <variable name="i_Ca_L_K_cyt" units="nanoA" public_interface="in"/>
      <variable name="i_Ca_L_K_ds" units="nanoA" public_interface="in"/>
      <variable name="i_to" units="nanoA" public_interface="in"/>
      <variable name="i_NaK" units="nanoA" public_interface="in"/>
      <variable name="time" units="second" public_interface="in"/>
      <variable name="F" units="coulomb_per_mole" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>K_i</ci>
            </apply>
            <apply>
               <times/>
               <apply>
                  <divide/>
                  <apply>
                     <minus/>
                     <cn cellml:units="mA_nA">1</cn>
                  </apply>
                  <apply>
                     <times/>
                     <cn cellml:units="litre_micrometre3">1</cn>
                     <ci>V_i</ci>
                     <ci>F</ci>
                  </apply>
               </apply>
               <apply>
                  <minus/>
                  <apply>
                     <plus/>
                     <ci>i_K1</ci>
                     <ci>i_Kr</ci>
                     <ci>i_Ks</ci>
                     <ci>i_Ca_L_K_cyt</ci>
                     <ci>i_Ca_L_K_ds</ci>
                     <ci>i_to</ci>
                  </apply>
                  <apply>
                     <times/>
                     <cn cellml:units="dimensionless">2</cn>
                     <ci>i_NaK</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="intracellular_calcium_concentration">
      <variable name="Ca_i" units="millimolar" initial_value="1.49e-5" public_interface="out"/>
      <variable name="Ca_ds" units="millimolar" initial_value="2.7e-6" public_interface="out"/>
      <variable name="Ca_up" units="millimolar" initial_value="0.36963" public_interface="out"/>
      <variable name="Ca_rel" units="millimolar" initial_value="0.6460487" public_interface="out"/>
      <variable name="Ca_Calmod" units="millimolar" initial_value="0.0005841"/>
      <variable name="Ca_Trop" units="millimolar" initial_value="0.0003732"/>
      <variable name="Calmod" units="millimolar" initial_value="0.02"/>
      <variable name="Trop" units="millimolar" initial_value="0.05"/>
      <variable name="alpha_Calmod" units="per_millimolar_second" initial_value="100000"/>
      <variable name="beta_Calmod" units="per_second" initial_value="50"/>
      <variable name="alpha_Trop" units="per_millimolar_second" initial_value="100000"/>
      <variable name="beta_Trop" units="per_second" initial_value="200"/>
      <variable name="radius" units="micrometre" initial_value="12"/>
      <variable name="length" units="micrometre" initial_value="74"/>
      <variable name="V_i" units="micrometre3" public_interface="out"/>
      <variable name="V_Cell" units="micrometre3"/>
      <variable name="V_i_ratio" units="dimensionless"/>
      <variable name="V_ds_ratio" units="dimensionless" initial_value="0.1"/>
      <variable name="V_rel_ratio" units="dimensionless" initial_value="0.1"/>
      <variable name="V_e_ratio" units="dimensionless" initial_value="0.4"/>
      <variable name="V_up_ratio" units="dimensionless" initial_value="0.01"/>
      <variable name="Kdecay" units="per_second" initial_value="10"/>
      <variable name="i_up" units="millimolar_per_second" public_interface="in"/>
      <variable name="i_trans" units="millimolar_per_second" public_interface="in"/>
      <variable name="i_rel" units="millimolar_per_second" public_interface="in"/>
      <variable name="i_NaCa_cyt" units="nanoA" public_interface="in"/>
      <variable name="i_NaCa_ds" units="nanoA" public_interface="in"/>
      <variable name="i_Ca_L_Ca_cyt" units="nanoA" public_interface="in"/>
      <variable name="i_Ca_L_Ca_ds" units="nanoA" public_interface="in"/>
      <variable name="i_b_Ca" units="nanoA" public_interface="in"/>
      <variable name="F" units="coulomb_per_mole" public_interface="in"/>
      <variable name="time" units="second" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>V_Cell</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <cn cellml:units="dimensionless">3.141592654</cn>
                  <apply>
                     <power/>
                     <apply>
                        <divide/>
                        <ci>radius</ci>
                        <cn cellml:units="dimensionless">1000</cn>
                     </apply>
                     <cn cellml:units="dimensionless">2</cn>
                  </apply>
                  <ci>length</ci>
               </apply>
               <cn cellml:units="dimensionless">1000</cn>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>V_i_ratio</ci>
            <apply>
               <minus/>
               <apply>
                  <minus/>
                  <apply>
                     <minus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <ci>V_e_ratio</ci>
                  </apply>
                  <ci>V_up_ratio</ci>
               </apply>
               <ci>V_rel_ratio</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>V_i</ci>
            <apply>
               <times/>
               <ci>V_Cell</ci>
               <ci>V_i_ratio</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>Ca_i</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <minus/>
                  <apply>
                     <minus/>
                     <apply>
                        <plus/>
                        <apply>
                           <times/>
                           <apply>
                              <divide/>
                              <apply>
                                 <minus/>
                                 <cn cellml:units="mA_nA">1</cn>
                              </apply>
                              <apply>
                                 <times/>
                                 <cn cellml:units="dimensionless">2</cn>
                                 <cn cellml:units="litre_micrometre3">1</cn>
                                 <ci>V_i</ci>
                                 <ci>F</ci>
                              </apply>
                           </apply>
                           <apply>
                              <minus/>
                              <apply>
                                 <minus/>
                                 <apply>
                                    <plus/>
                                    <ci>i_Ca_L_Ca_cyt</ci>
                                    <ci>i_b_Ca</ci>
                                 </apply>
                                 <apply>
                                    <times/>
                                    <cn cellml:units="dimensionless">2</cn>
                                    <ci>i_NaCa_cyt</ci>
                                 </apply>
                              </apply>
                              <apply>
                                 <times/>
                                 <cn cellml:units="dimensionless">2</cn>
                                 <ci>i_NaCa_ds</ci>
                              </apply>
                           </apply>
                        </apply>
                        <apply>
                           <times/>
                           <ci>Ca_ds</ci>
                           <ci>V_ds_ratio</ci>
                           <ci>Kdecay</ci>
                        </apply>
                        <apply>
                           <divide/>
                           <apply>
                              <times/>
                              <ci>i_rel</ci>
                              <ci>V_rel_ratio</ci>
                           </apply>
                           <ci>V_i_ratio</ci>
                        </apply>
                     </apply>
                     <apply>
                        <diff/>
                        <bvar>
                           <ci>time</ci>
                        </bvar>
                        <ci>Ca_Calmod</ci>
                     </apply>
                  </apply>
                  <apply>
                     <diff/>
                     <bvar>
                        <ci>time</ci>
                     </bvar>
                     <ci>Ca_Trop</ci>
                  </apply>
               </apply>
               <ci>i_up</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>Ca_ds</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <divide/>
                  <apply>
                     <times/>
                     <apply>
                        <minus/>
                        <cn cellml:units="mA_nA">1</cn>
                     </apply>
                     <ci>i_Ca_L_Ca_ds</ci>
                  </apply>
                  <apply>
                     <times/>
                     <cn cellml:units="dimensionless">2</cn>
                     <cn cellml:units="litre_micrometre3">1</cn>
                     <ci>V_ds_ratio</ci>
                     <ci>V_i</ci>
                     <ci>F</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>Ca_ds</ci>
                  <ci>Kdecay</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>Ca_up</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <apply>
                     <divide/>
                     <ci>V_i_ratio</ci>
                     <ci>V_up_ratio</ci>
                  </apply>
                  <ci>i_up</ci>
               </apply>
               <ci>i_trans</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>Ca_rel</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <apply>
                     <divide/>
                     <ci>V_up_ratio</ci>
                     <ci>V_rel_ratio</ci>
                  </apply>
                  <ci>i_trans</ci>
               </apply>
               <ci>i_rel</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>Ca_Calmod</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>alpha_Calmod</ci>
                  <ci>Ca_i</ci>
                  <apply>
                     <minus/>
                     <ci>Calmod</ci>
                     <ci>Ca_Calmod</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>beta_Calmod</ci>
                  <ci>Ca_Calmod</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>Ca_Trop</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>alpha_Trop</ci>
                  <ci>Ca_i</ci>
                  <apply>
                     <minus/>
                     <ci>Trop</ci>
                     <ci>Ca_Trop</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>beta_Trop</ci>
                  <ci>Ca_Trop</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <group>
      <relationship_ref relationship="containment"/>
      <component_ref component="membrane">
         <component_ref component="fast_sodium_current">
            <component_ref component="fast_sodium_current_m_gate"/>
            <component_ref component="fast_sodium_current_h_gate"/>
         </component_ref>
         <component_ref component="L_type_Ca_channel">
            <component_ref component="L_type_Ca_channel_d_gate"/>
            <component_ref component="L_type_Ca_channel_f_gate"/>
         </component_ref>
         <component_ref component="rapid_delayed_rectifier_potassium_current">
            <component_ref component="rapid_delayed_rectifier_potassium_current_xr1_gate"/>
            <component_ref component="rapid_delayed_rectifier_potassium_current_xr2_gate"/>
         </component_ref>
         <component_ref component="slow_delayed_rectifier_potassium_current">
            <component_ref component="slow_delayed_rectifier_potassium_current_xs_gate"/>
         </component_ref>
         <component_ref component="time_independent_potassium_current"/>
         <component_ref component="persistent_sodium_current"/>
         <component_ref component="sodium_background_current"/>
         <component_ref component="calcium_background_current"/>
         <component_ref component="transient_outward_current">
            <component_ref component="transient_outward_current_s_gate"/>
            <component_ref component="transient_outward_current_r_gate"/>
         </component_ref>
         <component_ref component="sodium_potassium_pump"/>
         <component_ref component="sodium_calcium_exchanger"/>
         <component_ref component="sarcoplasmic_reticulum_calcium_pump"/>
         <component_ref component="calcium_translocation"/>
         <component_ref component="calcium_release"/>
         <component_ref component="extracellular_sodium_concentration"/>
         <component_ref component="intracellular_sodium_concentration"/>
         <component_ref component="extracellular_potassium_concentration"/>
         <component_ref component="intracellular_potassium_concentration"/>
         <component_ref component="extracellular_calcium_concentration"/>
         <component_ref component="intracellular_calcium_concentration"/>
         <component_ref component="reversal_potentials"/>
      </component_ref>
   </group>
   <group>
      <relationship_ref relationship="encapsulation"/>
      <component_ref component="fast_sodium_current">
         <component_ref component="fast_sodium_current_m_gate"/>
         <component_ref component="fast_sodium_current_h_gate"/>
      </component_ref>
      <component_ref component="L_type_Ca_channel">
         <component_ref component="L_type_Ca_channel_d_gate"/>
         <component_ref component="L_type_Ca_channel_f_gate"/>
      </component_ref>
      <component_ref component="transient_outward_current">
         <component_ref component="transient_outward_current_s_gate"/>
         <component_ref component="transient_outward_current_r_gate"/>
      </component_ref>
      <component_ref component="rapid_delayed_rectifier_potassium_current">
         <component_ref component="rapid_delayed_rectifier_potassium_current_xr1_gate"/>
         <component_ref component="rapid_delayed_rectifier_potassium_current_xr2_gate"/>
      </component_ref>
      <component_ref component="slow_delayed_rectifier_potassium_current">
         <component_ref component="slow_delayed_rectifier_potassium_current_xs_gate"/>
      </component_ref>
   </group>
   <connection>
      <map_components component_1="membrane" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="reversal_potentials" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="fast_sodium_current" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="L_type_Ca_channel" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="rapid_delayed_rectifier_potassium_current" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="slow_delayed_rectifier_potassium_current" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="time_independent_potassium_current" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="persistent_sodium_current" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="sodium_background_current" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="calcium_background_current" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="transient_outward_current" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="sodium_potassium_pump" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="sodium_calcium_exchanger" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="intracellular_sodium_concentration" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="intracellular_potassium_concentration" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="intracellular_calcium_concentration" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="sarcoplasmic_reticulum_calcium_pump" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="calcium_translocation" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="calcium_release" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="calcium_release" component_2="membrane"/>
      <map_variables variable_1="V" variable_2="V"/>
   </connection>
   <connection>
      <map_components component_1="membrane" component_2="reversal_potentials"/>
      <map_variables variable_1="R" variable_2="R"/>
      <map_variables variable_1="T" variable_2="T"/>
      <map_variables variable_1="F" variable_2="F"/>
   </connection>
   <connection>
      <map_components component_1="membrane" component_2="fast_sodium_current"/>
      <map_variables variable_1="V" variable_2="V"/>
      <map_variables variable_1="i_Na" variable_2="i_Na"/>
   </connection>
   <connection>
      <map_components component_1="membrane" component_2="L_type_Ca_channel"/>
      <map_variables variable_1="V" variable_2="V"/>
      <map_variables variable_1="i_Ca_L_Ca_cyt" variable_2="i_Ca_L_Ca_cyt"/>
      <map_variables variable_1="i_Ca_L_Ca_ds" variable_2="i_Ca_L_Ca_ds"/>
      <map_variables variable_1="i_Ca_L_Na_cyt" variable_2="i_Ca_L_Na_cyt"/>
      <map_variables variable_1="i_Ca_L_Na_ds" variable_2="i_Ca_L_Na_ds"/>
      <map_variables variable_1="i_Ca_L_K_cyt" variable_2="i_Ca_L_K_cyt"/>
      <map_variables variable_1="i_Ca_L_K_ds" variable_2="i_Ca_L_K_ds"/>
      <map_variables variable_1="R" variable_2="R"/>
      <map_variables variable_1="T" variable_2="T"/>
      <map_variables variable_1="F" variable_2="F"/>
   </connection>
   <connection>
      <map_components component_1="membrane" component_2="rapid_delayed_rectifier_potassium_current"/>
      <map_variables variable_1="V" variable_2="V"/>
      <map_variables variable_1="i_Kr" variable_2="i_Kr"/>
   </connection>
   <connection>
      <map_components component_1="membrane" component_2="slow_delayed_rectifier_potassium_current"/>
      <map_variables variable_1="V" variable_2="V"/>
      <map_variables variable_1="i_Ks" variable_2="i_Ks"/>
   </connection>
   <connection>
      <map_components component_1="membrane" component_2="transient_outward_current"/>
      <map_variables variable_1="V" variable_2="V"/>
      <map_variables variable_1="i_to" variable_2="i_to"/>
   </connection>
   <connection>
      <map_components component_1="membrane" component_2="time_independent_potassium_current"/>
      <map_variables variable_1="V" variable_2="V"/>
      <map_variables variable_1="i_K1" variable_2="i_K1"/>
      <map_variables variable_1="R" variable_2="R"/>
      <map_variables variable_1="T" variable_2="T"/>
      <map_variables variable_1="F" variable_2="F"/>
   </connection>
   <connection>
      <map_components component_1="membrane" component_2="persistent_sodium_current"/>
      <map_variables variable_1="V" variable_2="V"/>
      <map_variables variable_1="i_p_Na" variable_2="i_p_Na"/>
   </connection>
   <connection>
      <map_components component_1="membrane" component_2="sodium_background_current"/>
      <map_variables variable_1="V" variable_2="V"/>
      <map_variables variable_1="i_b_Na" variable_2="i_b_Na"/>
   </connection>
   <connection>
      <map_components component_1="membrane" component_2="calcium_background_current"/>
      <map_variables variable_1="V" variable_2="V"/>
      <map_variables variable_1="i_b_Ca" variable_2="i_b_Ca"/>
   </connection>
   <connection>
      <map_components component_1="membrane" component_2="sodium_potassium_pump"/>
      <map_variables variable_1="i_NaK" variable_2="i_NaK"/>
   </connection>
   <connection>
      <map_components component_1="membrane" component_2="sodium_calcium_exchanger"/>
      <map_variables variable_1="V" variable_2="V"/>
      <map_variables variable_1="i_NaCa_cyt" variable_2="i_NaCa_cyt"/>
      <map_variables variable_1="i_NaCa_ds" variable_2="i_NaCa_ds"/>
      <map_variables variable_1="R" variable_2="R"/>
      <map_variables variable_1="T" variable_2="T"/>
      <map_variables variable_1="F" variable_2="F"/>
   </connection>
   <connection>
      <map_components component_1="membrane" component_2="intracellular_sodium_concentration"/>
      <map_variables variable_1="F" variable_2="F"/>
   </connection>
   <connection>
      <map_components component_1="membrane" component_2="intracellular_potassium_concentration"/>
      <map_variables variable_1="F" variable_2="F"/>
   </connection>
   <connection>
      <map_components component_1="membrane" component_2="intracellular_calcium_concentration"/>
      <map_variables variable_1="F" variable_2="F"/>
   </connection>
   <connection>
      <map_components component_1="reversal_potentials" component_2="sodium_background_current"/>
      <map_variables variable_1="E_Na" variable_2="E_Na"/>
   </connection>
   <connection>
      <map_components component_1="reversal_potentials" component_2="persistent_sodium_current"/>
      <map_variables variable_1="E_Na" variable_2="E_Na"/>
   </connection>
   <connection>
      <map_components component_1="reversal_potentials" component_2="time_independent_potassium_current"/>
      <map_variables variable_1="E_K" variable_2="E_K"/>
   </connection>
   <connection>
      <map_components component_1="reversal_potentials" component_2="rapid_delayed_rectifier_potassium_current"/>
      <map_variables variable_1="E_K" variable_2="E_K"/>
   </connection>
   <connection>
      <map_components component_1="reversal_potentials" component_2="slow_delayed_rectifier_potassium_current"/>
      <map_variables variable_1="E_Ks" variable_2="E_Ks"/>
   </connection>
   <connection>
      <map_components component_1="reversal_potentials" component_2="transient_outward_current"/>
      <map_variables variable_1="E_K" variable_2="E_K"/>
   </connection>
   <connection>
      <map_components component_1="reversal_potentials" component_2="calcium_background_current"/>
      <map_variables variable_1="E_Ca" variable_2="E_Ca"/>
   </connection>
   <connection>
      <map_components component_1="reversal_potentials" component_2="fast_sodium_current"/>
      <map_variables variable_1="E_mh" variable_2="E_mh"/>
   </connection>
   <connection>
      <map_components component_1="reversal_potentials" component_2="extracellular_potassium_concentration"/>
      <map_variables variable_1="K_o" variable_2="K_o"/>
   </connection>
   <connection>
      <map_components component_1="reversal_potentials" component_2="extracellular_sodium_concentration"/>
      <map_variables variable_1="Na_o" variable_2="Na_o"/>
   </connection>
   <connection>
      <map_components component_1="reversal_potentials" component_2="intracellular_potassium_concentration"/>
      <map_variables variable_1="K_i" variable_2="K_i"/>
   </connection>
   <connection>
      <map_components component_1="reversal_potentials" component_2="intracellular_sodium_concentration"/>
      <map_variables variable_1="Na_i" variable_2="Na_i"/>
   </connection>
   <connection>
      <map_components component_1="reversal_potentials" component_2="extracellular_calcium_concentration"/>
      <map_variables variable_1="Ca_o" variable_2="Ca_o"/>
   </connection>
   <connection>
      <map_components component_1="reversal_potentials" component_2="intracellular_calcium_concentration"/>
      <map_variables variable_1="Ca_i" variable_2="Ca_i"/>
   </connection>
   <connection>
      <map_components component_1="time_independent_potassium_current" component_2="extracellular_potassium_concentration"/>
      <map_variables variable_1="K_o" variable_2="K_o"/>
   </connection>
   <connection>
      <map_components component_1="time_independent_potassium_current" component_2="intracellular_potassium_concentration"/>
      <map_variables variable_1="i_K1" variable_2="i_K1"/>
   </connection>
   <connection>
      <map_components component_1="rapid_delayed_rectifier_potassium_current" component_2="intracellular_potassium_concentration"/>
      <map_variables variable_1="i_Kr" variable_2="i_Kr"/>
   </connection>
   <connection>
      <map_components component_1="slow_delayed_rectifier_potassium_current" component_2="intracellular_potassium_concentration"/>
      <map_variables variable_1="i_Ks" variable_2="i_Ks"/>
   </connection>
   <connection>
      <map_components component_1="fast_sodium_current" component_2="intracellular_sodium_concentration"/>
      <map_variables variable_1="i_Na" variable_2="i_Na"/>
   </connection>
   <connection>
      <map_components component_1="sodium_background_current" component_2="intracellular_sodium_concentration"/>
      <map_variables variable_1="i_b_Na" variable_2="i_b_Na"/>
      <map_variables variable_1="Na_i" variable_2="Na_i"/>
   </connection>
   <connection>
      <map_components component_1="persistent_sodium_current" component_2="intracellular_sodium_concentration"/>
      <map_variables variable_1="i_p_Na" variable_2="i_p_Na"/>
   </connection>
   <connection>
      <map_components component_1="L_type_Ca_channel" component_2="intracellular_sodium_concentration"/>
      <map_variables variable_1="i_Ca_L_Na_cyt" variable_2="i_Ca_L_Na_cyt"/>
      <map_variables variable_1="i_Ca_L_Na_ds" variable_2="i_Ca_L_Na_ds"/>
      <map_variables variable_1="Na_i" variable_2="Na_i"/>
   </connection>
   <connection>
      <map_components component_1="L_type_Ca_channel" component_2="intracellular_calcium_concentration"/>
      <map_variables variable_1="i_Ca_L_Ca_cyt" variable_2="i_Ca_L_Ca_cyt"/>
      <map_variables variable_1="i_Ca_L_Ca_ds" variable_2="i_Ca_L_Ca_ds"/>
      <map_variables variable_1="Ca_i" variable_2="Ca_i"/>
      <map_variables variable_1="Ca_ds" variable_2="Ca_ds"/>
   </connection>
   <connection>
      <map_components component_1="L_type_Ca_channel" component_2="intracellular_potassium_concentration"/>
      <map_variables variable_1="i_Ca_L_K_cyt" variable_2="i_Ca_L_K_cyt"/>
      <map_variables variable_1="i_Ca_L_K_ds" variable_2="i_Ca_L_K_ds"/>
      <map_variables variable_1="K_i" variable_2="K_i"/>
   </connection>
   <connection>
      <map_components component_1="L_type_Ca_channel" component_2="extracellular_potassium_concentration"/>
      <map_variables variable_1="K_o" variable_2="K_o"/>
   </connection>
   <connection>
      <map_components component_1="L_type_Ca_channel" component_2="extracellular_sodium_concentration"/>
      <map_variables variable_1="Na_o" variable_2="Na_o"/>
   </connection>
   <connection>
      <map_components component_1="L_type_Ca_channel" component_2="extracellular_calcium_concentration"/>
      <map_variables variable_1="Ca_o" variable_2="Ca_o"/>
   </connection>
   <connection>
      <map_components component_1="calcium_background_current" component_2="intracellular_calcium_concentration"/>
      <map_variables variable_1="i_b_Ca" variable_2="i_b_Ca"/>
   </connection>
   <connection>
      <map_components component_1="sodium_potassium_pump" component_2="intracellular_sodium_concentration"/>
      <map_variables variable_1="i_NaK" variable_2="i_NaK"/>
      <map_variables variable_1="Na_i" variable_2="Na_i"/>
   </connection>
   <connection>
      <map_components component_1="sodium_potassium_pump" component_2="extracellular_potassium_concentration"/>
      <map_variables variable_1="K_o" variable_2="K_o"/>
   </connection>
   <connection>
      <map_components component_1="sodium_potassium_pump" component_2="intracellular_potassium_concentration"/>
      <map_variables variable_1="i_NaK" variable_2="i_NaK"/>
   </connection>
   <connection>
      <map_components component_1="transient_outward_current" component_2="intracellular_potassium_concentration"/>
      <map_variables variable_1="i_to" variable_2="i_to"/>
   </connection>
   <connection>
      <map_components component_1="sodium_calcium_exchanger" component_2="intracellular_calcium_concentration"/>
      <map_variables variable_1="i_NaCa_cyt" variable_2="i_NaCa_cyt"/>
      <map_variables variable_1="i_NaCa_ds" variable_2="i_NaCa_ds"/>
      <map_variables variable_1="Ca_i" variable_2="Ca_i"/>
      <map_variables variable_1="Ca_ds" variable_2="Ca_ds"/>
   </connection>
   <connection>
      <map_components component_1="sodium_calcium_exchanger" component_2="intracellular_sodium_concentration"/>
      <map_variables variable_1="i_NaCa_cyt" variable_2="i_NaCa_cyt"/>
      <map_variables variable_1="Na_i" variable_2="Na_i"/>
   </connection>
   <connection>
      <map_components component_1="sodium_calcium_exchanger" component_2="extracellular_sodium_concentration"/>
      <map_variables variable_1="Na_o" variable_2="Na_o"/>
   </connection>
   <connection>
      <map_components component_1="sodium_calcium_exchanger" component_2="extracellular_calcium_concentration"/>
      <map_variables variable_1="Ca_o" variable_2="Ca_o"/>
   </connection>
   <connection>
      <map_components component_1="sarcoplasmic_reticulum_calcium_pump" component_2="intracellular_calcium_concentration"/>
      <map_variables variable_1="i_up" variable_2="i_up"/>
      <map_variables variable_1="Ca_i" variable_2="Ca_i"/>
      <map_variables variable_1="Ca_up" variable_2="Ca_up"/>
   </connection>
   <connection>
      <map_components component_1="calcium_translocation" component_2="intracellular_calcium_concentration"/>
      <map_variables variable_1="i_trans" variable_2="i_trans"/>
      <map_variables variable_1="Ca_rel" variable_2="Ca_rel"/>
      <map_variables variable_1="Ca_up" variable_2="Ca_up"/>
   </connection>
   <connection>
      <map_components component_1="calcium_release" component_2="intracellular_calcium_concentration"/>
      <map_variables variable_1="i_rel" variable_2="i_rel"/>
      <map_variables variable_1="Ca_rel" variable_2="Ca_rel"/>
      <map_variables variable_1="Ca_i" variable_2="Ca_i"/>
      <map_variables variable_1="Ca_ds" variable_2="Ca_ds"/>
   </connection>
   <connection>
      <map_components component_1="rapid_delayed_rectifier_potassium_current" component_2="rapid_delayed_rectifier_potassium_current_xr1_gate"/>
      <map_variables variable_1="xr1" variable_2="xr1"/>
      <map_variables variable_1="time" variable_2="time"/>
      <map_variables variable_1="V" variable_2="V"/>
   </connection>
   <connection>
      <map_components component_1="rapid_delayed_rectifier_potassium_current" component_2="rapid_delayed_rectifier_potassium_current_xr2_gate"/>
      <map_variables variable_1="xr2" variable_2="xr2"/>
      <map_variables variable_1="time" variable_2="time"/>
      <map_variables variable_1="V" variable_2="V"/>
   </connection>
   <connection>
      <map_components component_1="slow_delayed_rectifier_potassium_current" component_2="slow_delayed_rectifier_potassium_current_xs_gate"/>
      <map_variables variable_1="xs" variable_2="xs"/>
      <map_variables variable_1="time" variable_2="time"/>
      <map_variables variable_1="V" variable_2="V"/>
   </connection>
   <connection>
      <map_components component_1="fast_sodium_current" component_2="fast_sodium_current_m_gate"/>
      <map_variables variable_1="m" variable_2="m"/>
      <map_variables variable_1="time" variable_2="time"/>
      <map_variables variable_1="V" variable_2="V"/>
   </connection>
   <connection>
      <map_components component_1="fast_sodium_current" component_2="fast_sodium_current_h_gate"/>
      <map_variables variable_1="h" variable_2="h"/>
      <map_variables variable_1="time" variable_2="time"/>
      <map_variables variable_1="V" variable_2="V"/>
   </connection>
   <connection>
      <map_components component_1="L_type_Ca_channel" component_2="L_type_Ca_channel_d_gate"/>
      <map_variables variable_1="d" variable_2="d"/>
      <map_variables variable_1="time" variable_2="time"/>
      <map_variables variable_1="V" variable_2="V"/>
   </connection>
   <connection>
      <map_components component_1="L_type_Ca_channel" component_2="L_type_Ca_channel_f_gate"/>
      <map_variables variable_1="CaChoncyt" variable_2="CaChoncyt"/>
      <map_variables variable_1="CaChonds" variable_2="CaChonds"/>
      <map_variables variable_1="KCaChoff" variable_2="KCaChoff"/>
      <map_variables variable_1="Kmdsinact" variable_2="Kmdsinact"/>
      <map_variables variable_1="Ca_i" variable_2="Ca_i"/>
      <map_variables variable_1="Ca_ds" variable_2="Ca_ds"/>
      <map_variables variable_1="time" variable_2="time"/>
      <map_variables variable_1="V" variable_2="V"/>
   </connection>
   <connection>
      <map_components component_1="transient_outward_current" component_2="transient_outward_current_s_gate"/>
      <map_variables variable_1="s" variable_2="s"/>
      <map_variables variable_1="time" variable_2="time"/>
      <map_variables variable_1="V" variable_2="V"/>
   </connection>
   <connection>
      <map_components component_1="transient_outward_current" component_2="transient_outward_current_r_gate"/>
      <map_variables variable_1="r" variable_2="r"/>
      <map_variables variable_1="time" variable_2="time"/>
      <map_variables variable_1="V" variable_2="V"/>
   </connection>
   <connection>
      <map_components component_1="intracellular_sodium_concentration" component_2="intracellular_calcium_concentration"/>
      <map_variables variable_1="V_i" variable_2="V_i"/>
   </connection>
   <connection>
      <map_components component_1="intracellular_potassium_concentration" component_2="intracellular_calcium_concentration"/>
      <map_variables variable_1="V_i" variable_2="V_i"/>
   </connection>
</model>