Location: Li, Smith, 2009 @ 06e9c20c9cee / C57BL7_WT_model_version1.cellml

Author:
jlaw060 <jlaw060@bioeng203.bioeng.auckland.ac.nz>
Date:
2009-06-23 15:26:28+12:00
Desc:
renamed cellml file, changed cmeta ID of model
Permanent Source URI:
http://models.cellml.org/workspace/li_smith_2009/rawfile/06e9c20c9ceecf9538965391e9fd03f960d06c2e/C57BL7_WT_model_version1.cellml

<?xml version="1.0"?>
<!--
This CellML file was generated on 19/06/2009 at 09:10:40 using:

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

CellML 1.0 was used to generate this cellular model
http://www.CellML.org/
--><model xmlns="http://www.cellml.org/cellml/1.0#" xmlns:cmeta="http://www.cellml.org/metadata/1.0#" xml:base="file:///people/jlaw060/Desktop/C57BL7_WT_model_version1.cellml" cmeta:id="Li_Smith_2009" name="Li_Smith_2009">

<documentation xmlns="http://cellml.org/tmp-documentation">
<article>
  <articleinfo>
  <title>G-Protein Specificity In Synaptic Signalling</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>
</para>
  </section>
  <sect1 id="sec_structure">
<title>Model Structure</title>

<para>
A model of Ca<superscript>2+</superscript> dynamics in mouse ventricular myocytes from the C57BL/6 strain has been developed, based on experimental data at the same temperature and under consistent experimental conditions and protocols. 
</para>
<para>
Parameterisation of the NCX, SERCA and LCC formulations were done sequentially. The model of NCX, based on the Luo and Rudy 1994 formulation, was first fitted from the decay of the caffeine-induced [Ca<superscript>2+</superscript>]i transient, which allowed us to infer the flux through SERCA from the decay of the field-stimulated [Ca<superscript>2+</superscript>]i transient. Frequency-dependence of the flux through SERCA was accounted for using  Hund and Rudy 2002 formulation of Ca<superscript>2+</superscript>/calmodulin kinase II regulatory pathway, with parameters fitted to experimentally measured SERCA flux at 0.5, 1 and 3Hz pacing frequencies.  The L-type Ca<superscript>2+</superscript> channel model, with slight modifications to the Hinch 2004 formulation, was fitted to data from voltage-clamp experiments. All these component models were then incorporated into the whole-cell AP model of Bondarenko et al., with parameters of RyR adjusted to match the experimentally observed time-to-peak of the [Ca<superscript>2+</superscript>]i transient at 3Hz, and parameters
for the background Ca<superscript>2+</superscript> current and PMCA adjusted to match experimentally observed diastolic [Ca<superscript>2+</superscript>]i levels.
</para>

<informalfigure float="0" id="fig_reaction_diagram">
<mediaobject>
  <imageobject>
    <objectinfo>
      <title>reaction diagram</title>
    </objectinfo>
    <imagedata fileref="bondarenko_2004.png"/>
  </imageobject>
</mediaobject>
<caption>Schematic diagram of the Bondarenko 2004 model.</caption>
</informalfigure>

<sect2>
<title><emphasis>Associated data</emphasis></title>

<para><ulink url="Caffeine.zip">Caffeine</ulink> - Caffeine-induced Ca<superscript>2+</superscript> transients from 18 individual myocytes</para>
<para>Each file contains the time course of caffeine-induced Ca<superscript>2+</superscript> transients from a myocyte. The traces were exported from IonOptix for the duration of the caffeine application.</para>
<para>Related parameter: K<subscript>ncx</subscript></para>
<para>Fitting process: 1. The decay of individual transients were aligned and averaged. 
2. The rate of decay of the average transient, after taking into account of intracellular buffering, was used to calculate J<subscript>ncx</subscript>. J<subscript>ncx</subscript>  was then used to fit K<subscript>ncx</subscript> in the Luo and Rudy 1994 model.</para>


<para><ulink url="3Hz.zip">3Hz</ulink> - Field-stimulated Ca<superscript>2+</superscript> transients at 3Hz from 19 individual myocytes</para>
<para>Each file contains the time course of a Ca<superscript>2+</superscript> transient, field stimulated at 3Hz, from a myocyte. The time course was exported from IonOptix for a duration of 1s starting from the first stimulus.</para>
<para>Related parameter: K<subscript>m, serca</subscript>, V<subscript>serca</subscript> at 3Hz</para>
<para>Fitting process: 1. The individual transients were averaged. The decay of the average transient, from 70ms after-peak until the end of the period was analysed. 
2. J<subscript>ncx</subscript> was inferred from previously fitted NCX model, and the difference between the total Ca<superscript>2+</superscript> flux and J<subscript>ncx</subscript> was taken as the net SR Ca<superscript>2+</superscript> uptake J<subscript>up</subscript>. 
3 A small SR leak (0.008μM/ms) was included so the flux of Ca<superscript>2+</superscript> through SERCA J<subscript>serca</subscript> = J<subscript>up</subscript>+ J<subscript>leak</subscript>
4. J<subscript>serca</subscript> was used to fit the K<subscript>m,serca</subscript> and V<subscript>serca</subscript> at 3Hz.</para>

<para><ulink url="1Hz.zip">1Hz</ulink> - Field-stimulated Ca<superscript>2+</superscript> transients at 1Hz from 7 individual myocytes</para>
<para>Each file contains the time course of a Ca<superscript>2+</superscript> transient, field stimulated at 1Hz, from a myocyte.</para>
<para>Related parameter: Estimated V<subscript>serca</subscript> at 1Hz</para>
<para>1. The derivation of J<subscript>serca</subscript> at 1Hz is similar to above. 
2. J<subscript>serca</subscript> was used to constrain V<subscript>serca</subscript> at 1Hz
</para>

<para><ulink url="0.5Hz.zip">0.5Hz</ulink> - Field-stimulated Ca<superscript>2+</superscript> transients at 0.5Hz from 7 individual myocytes</para>
<para>Each file contains the time course of a Ca<superscript>2+</superscript> transient, field stimulated at 0.5Hz, from a myocyte</para>
<para>Related parameter: Estimated V<subscript>serca</subscript> at 0.5Hz</para>
<para>1. The derivation of J<subscript>serca</subscript> at 0.5Hz is similar to above. 
2. J<subscript>serca</subscript> was used to constrain V<subscript>serca</subscript> at 0.5Hz</para>

<para><ulink url="Ica.zip">Ica</ulink> - Time courses of transmembrane currents from 5 individual myocytes.</para>
<para>Each file contains the time course of experimentally measured L-type Ca<superscript>2+</superscript> current from one myocyte, measured in nA. </para>
<para>PCaL, ∆VL, VL, ΦL, tL, τL, KL, a, b, p, q</para>
<para>Fitting process: 1. The time course of transmembrane currents individual myocytes were averaged for each test potential between -20 and 20mV (-20, -10, 0, 10, 20mV). 
2. Averaged currents were converted into current densities by division by cell membrane capacitance (138.66pF).
3. Simulated sum of late K<superscript>+</superscript> currents, from the Bondarenko et al 2004 formulation, for test potentials -20, -10, 0, 10, 20mV were subtracted from the decay phase of each current.
4. The decay of the current at each test potential, after correcting for late K<superscript>+</superscript> currents, was fitted to a single exponential function, giving a time constant of decay tau value.
5. The IV relation and tau values at each test potential were used to fit the model of LCC.
</para>

<para><ulink url="Ica.zip">Ica/IV.xls</ulink> - Experimentally measured current-voltage relation for the L-type Ca<superscript>2+</superscript> current</para>
<para>Average current-voltage curve.</para>
<para>Related parameter: PCaL, ∆VL, VL, ΦL, tL, τL, KL, a, b, p, q</para>
<para>Fitting process: as for Ica
</para>
</sect2>
</sect1>
</article>
</documentation>

   <units name="millisecond">
      <unit prefix="milli" units="second"/>
   </units>
   <units name="per_millisecond">
      <unit exponent="-1" prefix="milli" units="second"/>
   </units>
   <units name="micromolar_per_millisecond">
      <unit units="micromolar"/>
      <unit exponent="-1" units="millisecond"/>
   </units>
   <units name="per_micromolar_millisecond">
      <unit exponent="-1" units="micromolar"/>
      <unit exponent="-1" units="millisecond"/>
   </units>
   <units name="micromolar3_per_millisecond">
      <unit exponent="-3" units="micromolar"/>
      <unit exponent="-1" units="millisecond"/>
   </units>
   <units name="micromolar4_per_millisecond">
      <unit exponent="-4" units="micromolar"/>
      <unit exponent="-1" units="millisecond"/>
   </units>
   <units name="millivolt">
      <unit prefix="milli" units="volt"/>
   </units>
   <units name="per_millivolt">
      <unit exponent="-1" prefix="milli" units="volt"/>
   </units>
   <units name="millivolt2">
      <unit exponent="2" prefix="milli" units="volt"/>
   </units>
   <units name="per_millivolt_millisecond">
      <unit exponent="-1" units="millivolt"/>
      <unit exponent="-1" units="millisecond"/>
   </units>
   <units name="milliS_per_microF">
      <unit prefix="milli" units="siemens"/>
      <unit exponent="-1" prefix="micro" units="farad"/>
   </units>
   <units name="picoF">
      <unit prefix="pico" units="farad"/>
   </units>
   <units name="per_picoF">
      <unit exponent="-1" units="picoF"/>
   </units>
   <units name="microF_per_cm2">
      <unit prefix="micro" units="farad"/>
      <unit exponent="-2" prefix="centi" units="metre"/>
   </units>
   <units name="picoA_per_picoF">
      <unit prefix="pico" units="ampere"/>
      <unit exponent="-1" prefix="pico" units="farad"/>
   </units>
   <units name="micromolar">
      <unit prefix="micro" units="mole"/>
      <unit exponent="-1" units="litre"/>
   </units>
   <units name="joule_per_mole_kelvin">
      <unit units="joule"/>
      <unit exponent="-1" units="mole"/>
      <unit exponent="-1" units="kelvin"/>
   </units>
   <units name="coulomb_per_millimole">
      <unit units="coulomb"/>
      <unit exponent="-1" prefix="milli" units="mole"/>
   </units>
   <units name="cm2">
      <unit exponent="2" prefix="centi" units="metre"/>
   </units>
   <units name="microlitre">
      <unit prefix="micro" units="litre"/>
   </units>
   <component name="environment">
      <variable name="time" public_interface="out" units="millisecond"/>
   </component>
   <component name="membrane">
      <variable initial_value="-85.64004" name="V" public_interface="out" units="millivolt"/>
      <variable initial_value="1" name="Cm" public_interface="out" units="microF_per_cm2"/>
      <variable initial_value="2.2e-5" name="Vmyo" public_interface="out" units="microlitre"/>
      <variable initial_value="1.022e-7" name="VJSR" public_interface="out" units="microlitre"/>
      <variable initial_value="1.786e-6" name="VNSR" public_interface="out" units="microlitre"/>
      <variable initial_value="1.264e-9" name="Vss" public_interface="out" units="microlitre"/>
      <variable initial_value="0.00013866" name="Acap" public_interface="out" units="cm2"/>
      <variable initial_value="5400" name="Ko" public_interface="out" units="micromolar"/>
      <variable initial_value="134000" name="Nao" public_interface="out" units="micromolar"/>
      <variable initial_value="1400" name="Cao" public_interface="out" units="micromolar"/>
      <variable initial_value="8.314" name="R" public_interface="out" units="joule_per_mole_kelvin"/>
      <variable initial_value="308" name="T" public_interface="out" units="kelvin"/>
      <variable initial_value="96.5" name="F" public_interface="out" units="coulomb_per_millimole"/>
      <variable name="time" public_interface="in" units="millisecond"/>
      <variable name="i_CaL" public_interface="in" units="picoA_per_picoF"/>
      <variable name="i_pCa" public_interface="in" units="picoA_per_picoF"/>
      <variable name="i_NaCa" public_interface="in" units="picoA_per_picoF"/>
      <variable name="i_Cab" public_interface="in" units="picoA_per_picoF"/>
      <variable name="i_Na" public_interface="in" units="picoA_per_picoF"/>
      <variable name="i_Nab" public_interface="in" units="picoA_per_picoF"/>
      <variable name="i_NaK" public_interface="in" units="picoA_per_picoF"/>
      <variable name="i_Kto_f" public_interface="in" units="picoA_per_picoF"/>
      <variable name="i_Kto_s" public_interface="in" units="picoA_per_picoF"/>
      <variable name="i_K1" public_interface="in" units="picoA_per_picoF"/>
      <variable name="i_Ks" public_interface="in" units="picoA_per_picoF"/>
      <variable name="i_Kur" public_interface="in" units="picoA_per_picoF"/>
      <variable name="i_Kss" public_interface="in" units="picoA_per_picoF"/>
      <variable name="i_ClCa" public_interface="in" units="picoA_per_picoF"/>
      <variable name="i_Kr" public_interface="in" units="picoA_per_picoF"/>
      <variable initial_value="100" name="stim_offset" units="millisecond"/>
      <variable initial_value="333.3333333" name="stim_period" units="millisecond"/>
      <variable initial_value="3" name="stim_duration" units="millisecond"/>
      <variable initial_value="-15" name="stim_amplitude" units="picoA_per_picoF"/>
      <variable name="i_Stim" public_interface="out" units="picoA_per_picoF"/>
      <variable name="past" units="millisecond"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>past</ci>
            <apply>
               <times/>
               <apply>
                  <floor/>
                  <apply>
                     <divide/>
                     <ci>time</ci>
                     <ci>stim_period</ci>
                  </apply>
               </apply>
               <ci>stim_period</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>i_Stim</ci>
            <piecewise>
               <piece>
                  <ci>stim_amplitude</ci>
                  <apply>
                     <and/>
                     <apply>
                        <geq/>
                        <apply>
                           <minus/>
                           <ci>time</ci>
                           <ci>past</ci>
                        </apply>
                        <ci>stim_offset</ci>
                     </apply>
                     <apply>
                        <leq/>
                        <apply>
                           <minus/>
                           <ci>time</ci>
                           <ci>past</ci>
                        </apply>
                        <apply>
                           <plus/>
                           <ci>stim_offset</ci>
                           <ci>stim_duration</ci>
                        </apply>
                     </apply>
                  </apply>
               </piece>
               <otherwise>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="picoA_per_picoF">0</cn>
               </otherwise>
            </piecewise>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>V</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <plus/>
                  <ci>i_CaL</ci>
                  <ci>i_pCa</ci>
                  <ci>i_NaCa</ci>
                  <ci>i_Cab</ci>
                  <ci>i_Na</ci>
                  <ci>i_Nab</ci>
                  <ci>i_NaK</ci>
                  <ci>i_Kto_f</ci>
                  <ci>i_Kto_s</ci>
                  <ci>i_K1</ci>
                  <ci>i_Ks</ci>
                  <ci>i_Kur</ci>
                  <ci>i_Kss</ci>
                  <ci>i_Kr</ci>
                  <ci>i_ClCa</ci>
                  <ci>i_Stim</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="calcium_concentration">
      <variable initial_value="0.1040595" name="Cai" public_interface="out" units="micromolar"/>
      <variable initial_value="0.1043777" name="Cass" public_interface="out" units="micromolar"/>
      <variable initial_value="730.0589" name="CaJSR" public_interface="out" units="micromolar"/>
      <variable initial_value="841.106" name="CaNSR" public_interface="out" units="micromolar"/>
      <variable name="Bi" units="dimensionless"/>
      <variable name="Bss" units="dimensionless"/>
      <variable name="BJSR" units="dimensionless"/>
      <variable initial_value="109" name="Bmax" units="micromolar"/>
      <variable initial_value="15000" name="CSQN_tot" units="micromolar"/>
      <variable initial_value="0.6" name="Kd" units="micromolar"/>
      <variable initial_value="800" name="Km_CSQN" units="micromolar"/>
      <variable name="Acap" public_interface="in" units="cm2"/>
      <variable name="Vmyo" public_interface="in" units="microlitre"/>
      <variable name="VJSR" public_interface="in" units="microlitre"/>
      <variable name="Vss" public_interface="in" units="microlitre"/>
      <variable name="VNSR" public_interface="in" units="microlitre"/>
      <variable name="V" public_interface="in" units="millivolt"/>
      <variable name="Cm" public_interface="in" units="microF_per_cm2"/>
      <variable name="time" public_interface="in" units="millisecond"/>
      <variable name="F" public_interface="in" units="coulomb_per_millimole"/>
      <variable name="i_Cab" public_interface="in" units="picoA_per_picoF"/>
      <variable name="i_CaL" public_interface="in" units="picoA_per_picoF"/>
      <variable name="i_NaCa" public_interface="in" units="picoA_per_picoF"/>
      <variable name="i_pCa" public_interface="in" units="picoA_per_picoF"/>
      <variable name="J_leak" public_interface="in" units="micromolar_per_millisecond"/>
      <variable name="J_rel" public_interface="in" units="micromolar_per_millisecond"/>
      <variable name="J_serca" public_interface="in" units="micromolar_per_millisecond"/>
      <variable name="J_tr" public_interface="in" units="micromolar_per_millisecond"/>
      <variable name="J_xfer" public_interface="in" units="micromolar_per_millisecond"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>Cai</ci>
            </apply>
            <apply>
               <times/>
               <ci>Bi</ci>
               <apply>
                  <minus/>
                  <apply>
                     <minus/>
                     <apply>
                        <plus/>
                        <ci>J_leak</ci>
                        <ci>J_xfer</ci>
                     </apply>
                     <apply>
                        <divide/>
                        <apply>
                           <times/>
                           <apply>
                              <minus/>
                              <apply>
                                 <plus/>
                                 <ci>i_Cab</ci>
                                 <ci>i_pCa</ci>
                              </apply>
                              <apply>
                                 <times/>
                                 <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">2</cn>
                                 <ci>i_NaCa</ci>
                              </apply>
                           </apply>
                           <ci>Acap</ci>
                           <ci>Cm</ci>
                        </apply>
                        <apply>
                           <times/>
                           <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">2</cn>
                           <ci>Vmyo</ci>
                           <ci>F</ci>
                        </apply>
                     </apply>
                  </apply>
                  <ci>J_serca</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>Cass</ci>
            </apply>
            <apply>
               <times/>
               <ci>Bss</ci>
               <apply>
                  <minus/>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <ci>J_rel</ci>
                        <ci>VJSR</ci>
                     </apply>
                     <ci>Vss</ci>
                  </apply>
                  <apply>
                     <plus/>
                     <apply>
                        <divide/>
                        <apply>
                           <times/>
                           <ci>J_xfer</ci>
                           <ci>Vmyo</ci>
                        </apply>
                        <ci>Vss</ci>
                     </apply>
                     <apply>
                        <divide/>
                        <apply>
                           <times/>
                           <ci>i_CaL</ci>
                           <ci>Acap</ci>
                           <ci>Cm</ci>
                        </apply>
                        <apply>
                           <times/>
                           <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">2</cn>
                           <ci>Vss</ci>
                           <ci>F</ci>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>CaJSR</ci>
            </apply>
            <apply>
               <times/>
               <ci>BJSR</ci>
               <apply>
                  <minus/>
                  <ci>J_tr</ci>
                  <ci>J_rel</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>CaNSR</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <divide/>
                  <apply>
                     <times/>
                     <apply>
                        <minus/>
                        <ci>J_serca</ci>
                        <ci>J_leak</ci>
                     </apply>
                     <ci>Vmyo</ci>
                  </apply>
                  <ci>VNSR</ci>
               </apply>
               <apply>
                  <divide/>
                  <apply>
                     <times/>
                     <ci>J_tr</ci>
                     <ci>VJSR</ci>
                  </apply>
                  <ci>VNSR</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>Bi</ci>
            <apply>
               <power/>
               <apply>
                  <plus/>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <ci>Bmax</ci>
                        <ci>Kd</ci>
                     </apply>
                     <apply>
                        <power/>
                        <apply>
                           <plus/>
                           <ci>Kd</ci>
                           <ci>Cai</ci>
                        </apply>
                        <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">2</cn>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <minus/>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>Bss</ci>
            <apply>
               <power/>
               <apply>
                  <plus/>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <ci>Bmax</ci>
                        <ci>Kd</ci>
                     </apply>
                     <apply>
                        <power/>
                        <apply>
                           <plus/>
                           <ci>Kd</ci>
                           <ci>Cass</ci>
                        </apply>
                        <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">2</cn>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <minus/>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>BJSR</ci>
            <apply>
               <power/>
               <apply>
                  <plus/>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <ci>CSQN_tot</ci>
                        <ci>Km_CSQN</ci>
                     </apply>
                     <apply>
                        <power/>
                        <apply>
                           <plus/>
                           <ci>Km_CSQN</ci>
                           <ci>CaJSR</ci>
                        </apply>
                        <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">2</cn>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <minus/>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="calcium_fluxes">
      <variable name="J_leak" public_interface="out" units="micromolar_per_millisecond"/>
      <variable name="J_rel" public_interface="out" units="micromolar_per_millisecond"/>
      <variable name="J_serca" public_interface="out" units="micromolar_per_millisecond"/>
      <variable name="J_tr" public_interface="out" units="micromolar_per_millisecond"/>
      <variable name="J_xfer" public_interface="out" units="micromolar_per_millisecond"/>
      <variable initial_value="2.290355e-9" name="P_RyR" units="dimensionless"/>
      <variable initial_value="4.5" name="v1" units="per_millisecond"/>
      <variable initial_value="20" name="tau_tr" units="millisecond"/>
      <variable initial_value="1e-5" name="v2" units="per_millisecond"/>
      <variable initial_value="8" name="tau_xfer" units="millisecond"/>
      <variable name="CaMKb" units="dimensionless"/>
      <variable initial_value="0.08989079" name="CaMKt" units="dimensionless"/>
      <variable name="CaMKa" units="dimensionless"/>
      <variable initial_value="0.05" name="on_rate" units="per_millisecond"/>
      <variable initial_value="0.00068" name="off_rate" units="per_millisecond"/>
      <variable name="vmup" units="micromolar_per_millisecond"/>
      <variable initial_value="0.412" name="Km_up" units="micromolar"/>
      <variable name="V" public_interface="in" units="millivolt"/>
      <variable name="time" public_interface="in" units="millisecond"/>
      <variable name="i_CaL_max" public_interface="in" units="picoA_per_picoF"/>
      <variable name="i_CaL" public_interface="in" units="picoA_per_picoF"/>
      <variable name="Cai" public_interface="in" units="micromolar"/>
      <variable name="Cass" public_interface="in" units="micromolar"/>
      <variable name="CaJSR" public_interface="in" units="micromolar"/>
      <variable name="CaNSR" public_interface="in" units="micromolar"/>
      <variable name="P_O1" public_interface="in" units="dimensionless"/>
      <variable name="P_O2" public_interface="in" units="dimensionless"/>
      <variable initial_value="0.3158" name="vmup_init" units="micromolar_per_millisecond"/>
      <variable initial_value="-0.09" name="P_ryr_const1" units="per_millisecond"/>
      <variable initial_value="-0.225" name="P_ryr_const2" units="per_millisecond"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>J_rel</ci>
            <apply>
               <times/>
               <ci>v1</ci>
               <apply>
                  <plus/>
                  <ci>P_O1</ci>
                  <ci>P_O2</ci>
               </apply>
               <apply>
                  <minus/>
                  <ci>CaJSR</ci>
                  <ci>Cass</ci>
               </apply>
               <ci>P_RyR</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>J_tr</ci>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>CaNSR</ci>
                  <ci>CaJSR</ci>
               </apply>
               <ci>tau_tr</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>J_xfer</ci>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>Cass</ci>
                  <ci>Cai</ci>
               </apply>
               <ci>tau_xfer</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>J_leak</ci>
            <apply>
               <times/>
               <ci>v2</ci>
               <apply>
                  <minus/>
                  <ci>CaNSR</ci>
                  <ci>Cai</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>CaMKb</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">0.05</cn>
                  <apply>
                     <minus/>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                     <ci>CaMKt</ci>
                  </apply>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
               </apply>
               <apply>
                  <plus/>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                  <apply>
                     <divide/>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="micromolar">0.7</cn>
                     <ci>Cass</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>CaMKt</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>on_rate</ci>
                  <ci>CaMKb</ci>
                  <apply>
                     <plus/>
                     <ci>CaMKb</ci>
                     <ci>CaMKt</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>off_rate</ci>
                  <ci>CaMKt</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>CaMKa</ci>
            <apply>
               <plus/>
               <ci>CaMKb</ci>
               <ci>CaMKt</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>vmup</ci>
            <apply>
               <times/>
               <apply>
                  <plus/>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">3.1512</cn>
                        <apply>
                           <power/>
                           <ci>CaMKa</ci>
                           <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">2.2062</cn>
                        </apply>
                     </apply>
                     <apply>
                        <plus/>
                        <apply>
                           <power/>
                           <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">0.1588</cn>
                           <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">2.2062</cn>
                        </apply>
                        <apply>
                           <power/>
                           <ci>CaMKa</ci>
                           <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">2.2062</cn>
                        </apply>
                     </apply>
                  </apply>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
               </apply>
               <ci>vmup_init</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>J_serca</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <ci>vmup</ci>
                  <apply>
                     <power/>
                     <ci>Cai</ci>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">2</cn>
                  </apply>
               </apply>
               <apply>
                  <plus/>
                  <apply>
                     <power/>
                     <ci>Km_up</ci>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">2</cn>
                  </apply>
                  <apply>
                     <power/>
                     <ci>Cai</ci>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">2</cn>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>P_RyR</ci>
            </apply>
            <apply>
               <plus/>
               <apply>
                  <times/>
                  <ci>P_ryr_const1</ci>
                  <ci>P_RyR</ci>
               </apply>
               <apply>
                  <times/>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <ci>P_ryr_const2</ci>
                        <ci>i_CaL</ci>
                     </apply>
                     <ci>i_CaL_max</ci>
                  </apply>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <minus/>
                           <apply>
                              <power/>
                              <apply>
                                 <minus/>
                                 <ci>V</ci>
                                 <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">5</cn>
                              </apply>
                              <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">2</cn>
                           </apply>
                        </apply>
                        <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt2">648</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="ryanodine_receptors">
      <variable initial_value="0.003825599" name="P_O1" public_interface="out" units="dimensionless"/>
      <variable initial_value="1.835831e-8" name="P_O2" public_interface="out" units="dimensionless"/>
      <variable name="P_C1" public_interface="out" units="dimensionless"/>
      <variable initial_value="0.3797679" name="P_C2" public_interface="out" units="dimensionless"/>
      <variable initial_value="0.006075" name="k_plus_a" units="micromolar4_per_millisecond"/>
      <variable initial_value="0.07125" name="k_minus_a" units="per_millisecond"/>
      <variable initial_value="0.00405" name="k_plus_b" units="micromolar3_per_millisecond"/>
      <variable initial_value="0.965" name="k_minus_b" units="per_millisecond"/>
      <variable initial_value="0.009" name="k_plus_c" units="per_millisecond"/>
      <variable initial_value="0.0008" name="k_minus_c" units="per_millisecond"/>
      <variable initial_value="3" name="m" units="dimensionless"/>
      <variable initial_value="4" name="n" units="dimensionless"/>
      <variable name="time" public_interface="in" units="millisecond"/>
      <variable name="Cass" public_interface="in" units="micromolar"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>P_O1</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <plus/>
                  <apply>
                     <times/>
                     <ci>k_plus_a</ci>
                     <apply>
                        <power/>
                        <ci>Cass</ci>
                        <ci>n</ci>
                     </apply>
                     <ci>P_C1</ci>
                  </apply>
                  <apply>
                     <times/>
                     <ci>k_minus_b</ci>
                     <ci>P_O2</ci>
                  </apply>
                  <apply>
                     <times/>
                     <ci>k_minus_c</ci>
                     <ci>P_C2</ci>
                  </apply>
               </apply>
               <apply>
                  <plus/>
                  <apply>
                     <times/>
                     <ci>k_minus_a</ci>
                     <ci>P_O1</ci>
                  </apply>
                  <apply>
                     <times/>
                     <ci>k_plus_b</ci>
                     <apply>
                        <power/>
                        <ci>Cass</ci>
                        <ci>m</ci>
                     </apply>
                     <ci>P_O1</ci>
                  </apply>
                  <apply>
                     <times/>
                     <ci>k_plus_c</ci>
                     <ci>P_O1</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>P_C1</ci>
            <apply>
               <minus/>
               <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
               <apply>
                  <plus/>
                  <ci>P_C2</ci>
                  <ci>P_O1</ci>
                  <ci>P_O2</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>P_O2</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>k_plus_b</ci>
                  <apply>
                     <power/>
                     <ci>Cass</ci>
                     <ci>m</ci>
                  </apply>
                  <ci>P_O1</ci>
               </apply>
               <apply>
                  <times/>
                  <ci>k_minus_b</ci>
                  <ci>P_O2</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>P_C2</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>k_plus_c</ci>
                  <ci>P_O1</ci>
               </apply>
               <apply>
                  <times/>
                  <ci>k_minus_c</ci>
                  <ci>P_C2</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="L_type_calcium_current">
      <variable name="i_CaL" public_interface="out" units="picoA_per_picoF"/>
      <variable initial_value="4" name="i_CaL_max" public_interface="out" units="picoA_per_picoF"/>
      <variable initial_value="19.1078" name="P_CaL" units="per_millisecond"/>
      <variable initial_value="4.373318e-6" name="O" units="dimensionless"/>
      <variable name="C" units="dimensionless"/>
      <variable initial_value="0.009171979" name="I" units="dimensionless"/>
      <variable initial_value="0.8876797" name="y_gate" units="dimensionless"/>
      <variable name="y_gate_inf" units="dimensionless"/>
      <variable name="y_gate_tau" units="millisecond"/>
      <variable name="Cass" public_interface="in" units="micromolar"/>
      <variable name="alpha_p" units="per_millisecond"/>
      <variable name="alpha_m" units="per_millisecond"/>
      <variable name="epsilon_p" units="per_micromolar_millisecond"/>
      <variable name="epsilon_m" units="per_millisecond"/>
      <variable initial_value="-2" name="V_L" units="millivolt"/>
      <variable initial_value="7.0671" name="delta_V_L" units="millivolt"/>
      <variable initial_value="1.1683" name="t_L" units="millisecond"/>
      <variable initial_value="1.6411" name="phi_L" units="dimensionless"/>
      <variable initial_value="0.07" name="a" units="dimensionless"/>
      <variable initial_value="14" name="b" units="dimensionless"/>
      <variable initial_value="972.9715" name="tau_L" units="millisecond"/>
      <variable initial_value="0.0964" name="K_L" units="micromolar"/>
      <variable name="expVL" units="dimensionless"/>
      <variable name="FVRT" units="dimensionless"/>
      <variable name="FVRT_Ca" units="dimensionless"/>
      <variable name="V" public_interface="in" units="millivolt"/>
      <variable name="time" public_interface="in" units="millisecond"/>
      <variable name="Cao" public_interface="in" units="micromolar"/>
      <variable name="T" public_interface="in" units="kelvin"/>
      <variable name="R" public_interface="in" units="joule_per_mole_kelvin"/>
      <variable name="Acap" public_interface="in" units="cm2"/>
      <variable name="Vss" public_interface="in" units="microlitre"/>
      <variable name="Cm" public_interface="in" units="microF_per_cm2"/>
      <variable name="F" public_interface="in" units="coulomb_per_millimole"/>
      <variable initial_value="6.6755" name="const5" units="millivolt"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>FVRT</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <ci>F</ci>
                  <ci>V</ci>
               </apply>
               <apply>
                  <times/>
                  <ci>R</ci>
                  <ci>T</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>FVRT_Ca</ci>
            <apply>
               <times/>
               <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">2</cn>
               <ci>FVRT</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>expVL</ci>
            <apply>
               <exp/>
               <apply>
                  <divide/>
                  <apply>
                     <minus/>
                     <ci>V</ci>
                     <ci>V_L</ci>
                  </apply>
                  <ci>delta_V_L</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>alpha_p</ci>
            <apply>
               <divide/>
               <ci>expVL</ci>
               <apply>
                  <times/>
                  <ci>t_L</ci>
                  <apply>
                     <plus/>
                     <ci>expVL</ci>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>alpha_m</ci>
            <apply>
               <divide/>
               <ci>phi_L</ci>
               <ci>t_L</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>epsilon_p</ci>
            <apply>
               <divide/>
               <apply>
                  <plus/>
                  <ci>expVL</ci>
                  <ci>a</ci>
               </apply>
               <apply>
                  <times/>
                  <ci>tau_L</ci>
                  <ci>K_L</ci>
                  <apply>
                     <plus/>
                     <ci>expVL</ci>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>epsilon_m</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <ci>b</ci>
                  <apply>
                     <plus/>
                     <ci>expVL</ci>
                     <ci>a</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>tau_L</ci>
                  <apply>
                     <plus/>
                     <apply>
                        <times/>
                        <ci>b</ci>
                        <ci>expVL</ci>
                     </apply>
                     <ci>a</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>y_gate_inf</ci>
            <apply>
               <plus/>
               <apply>
                  <divide/>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                  <apply>
                     <plus/>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <plus/>
                              <ci>V</ci>
                              <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">16.6577</cn>
                           </apply>
                           <ci>const5</ci>
                        </apply>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <divide/>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">0.1</cn>
                  <apply>
                     <plus/>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <plus/>
                              <apply>
                                 <minus/>
                                 <ci>V</ci>
                              </apply>
                              <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">40</cn>
                           </apply>
                           <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">6</cn>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>y_gate_tau</ci>
            <apply>
               <plus/>
               <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millisecond">20</cn>
               <apply>
                  <divide/>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millisecond">600</cn>
                  <apply>
                     <plus/>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <plus/>
                              <ci>V</ci>
                              <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">30</cn>
                           </apply>
                           <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">9.6</cn>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>C</ci>
            <apply>
               <minus/>
               <apply>
                  <minus/>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                  <ci>O</ci>
               </apply>
               <ci>I</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>O</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>alpha_p</ci>
                  <ci>C</ci>
               </apply>
               <apply>
                  <times/>
                  <ci>alpha_m</ci>
                  <ci>O</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>I</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>epsilon_p</ci>
                  <ci>Cass</ci>
                  <ci>C</ci>
               </apply>
               <apply>
                  <times/>
                  <ci>epsilon_m</ci>
                  <ci>I</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>y_gate</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>y_gate_inf</ci>
                  <ci>y_gate</ci>
               </apply>
               <ci>y_gate_tau</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>i_CaL</ci>
            <piecewise>
               <piece>
                  <apply>
                     <times/>
                     <apply>
                        <divide/>
                        <apply>
                           <times/>
                           <apply>
                              <divide/>
                              <apply>
                                 <times/>
                                 <apply>
                                    <minus/>
                                    <ci>P_CaL</ci>
                                 </apply>
                                 <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">2</cn>
                                 <ci>Vss</ci>
                                 <ci>F</ci>
                              </apply>
                              <apply>
                                 <times/>
                                 <ci>Acap</ci>
                                 <ci>Cm</ci>
                              </apply>
                           </apply>
                           <ci>O</ci>
                           <ci>y_gate</ci>
                           <ci>FVRT_Ca</ci>
                        </apply>
                        <apply>
                           <minus/>
                           <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                           <apply>
                              <exp/>
                              <apply>
                                 <minus/>
                                 <ci>FVRT_Ca</ci>
                              </apply>
                           </apply>
                        </apply>
                     </apply>
                     <apply>
                        <minus/>
                        <apply>
                           <times/>
                           <ci>Cao</ci>
                           <apply>
                              <exp/>
                              <apply>
                                 <minus/>
                                 <ci>FVRT_Ca</ci>
                              </apply>
                           </apply>
                        </apply>
                        <ci>Cass</ci>
                     </apply>
                  </apply>
                  <apply>
                     <gt/>
                     <apply>
                        <abs/>
                        <ci>FVRT_Ca</ci>
                     </apply>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">0.00001</cn>
                  </apply>
               </piece>
               <otherwise>
                  <apply>
                     <times/>
                     <apply>
                        <divide/>
                        <apply>
                           <times/>
                           <apply>
                              <divide/>
                              <apply>
                                 <times/>
                                 <apply>
                                    <minus/>
                                    <ci>P_CaL</ci>
                                 </apply>
                                 <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">2</cn>
                                 <ci>Vss</ci>
                                 <ci>F</ci>
                              </apply>
                              <apply>
                                 <times/>
                                 <ci>Acap</ci>
                                 <ci>Cm</ci>
                              </apply>
                           </apply>
                           <ci>O</ci>
                           <ci>y_gate</ci>
                           <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">0.00001</cn>
                        </apply>
                        <apply>
                           <minus/>
                           <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                           <apply>
                              <exp/>
                              <apply>
                                 <minus/>
                                 <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">0.00001</cn>
                              </apply>
                           </apply>
                        </apply>
                     </apply>
                     <apply>
                        <minus/>
                        <apply>
                           <times/>
                           <ci>Cao</ci>
                           <apply>
                              <exp/>
                              <apply>
                                 <minus/>
                                 <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">0.00001</cn>
                              </apply>
                           </apply>
                        </apply>
                        <ci>Cass</ci>
                     </apply>
                  </apply>
               </otherwise>
            </piecewise>
         </apply>
      </math>
   </component>
   <component name="calcium_pump_current">
      <variable name="i_pCa" public_interface="out" units="picoA_per_picoF"/>
      <variable initial_value="0.9" name="i_pCa_max" units="picoA_per_picoF"/>
      <variable initial_value="0.4" name="Km_pCa" units="micromolar"/>
      <variable name="Cai" public_interface="in" units="micromolar"/>
      <variable name="J_pCa" public_interface="out" units="micromolar_per_millisecond"/>
      <variable name="Vmyo" public_interface="in" units="microlitre"/>
      <variable name="Acap" public_interface="in" units="cm2"/>
      <variable name="Cm" public_interface="in" units="microF_per_cm2"/>
      <variable name="F" public_interface="in" units="coulomb_per_millimole"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_pCa</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <ci>i_pCa_max</ci>
                  <apply>
                     <power/>
                     <ci>Cai</ci>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">2</cn>
                  </apply>
               </apply>
               <apply>
                  <plus/>
                  <apply>
                     <power/>
                     <ci>Km_pCa</ci>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">2</cn>
                  </apply>
                  <apply>
                     <power/>
                     <ci>Cai</ci>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">2</cn>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>J_pCa</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <apply>
                     <minus/>
                     <ci>i_pCa</ci>
                  </apply>
                  <ci>Acap</ci>
                  <ci>Cm</ci>
               </apply>
               <apply>
                  <times/>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">2</cn>
                  <ci>Vmyo</ci>
                  <ci>F</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="sodium_calcium_exchange_current">
      <variable name="J_ncx" units="micromolar_per_millisecond"/>
      <variable name="i_NaCa" public_interface="out" units="picoA_per_picoF"/>
      <variable initial_value="772.8991" name="k_NaCa" units="picoA_per_picoF"/>
      <variable initial_value="86500" name="K_mNa" units="micromolar"/>
      <variable initial_value="1380" name="K_mCa" units="micromolar"/>
      <variable initial_value="0.1" name="k_sat" units="dimensionless"/>
      <variable initial_value="0.35" name="eta" units="dimensionless"/>
      <variable name="V" public_interface="in" units="millivolt"/>
      <variable name="R" public_interface="in" units="joule_per_mole_kelvin"/>
      <variable name="T" public_interface="in" units="kelvin"/>
      <variable name="F" public_interface="in" units="coulomb_per_millimole"/>
      <variable name="Nai" public_interface="in" units="micromolar"/>
      <variable name="Nao" public_interface="in" units="micromolar"/>
      <variable name="Cai" public_interface="in" units="micromolar"/>
      <variable name="Cao" public_interface="in" units="micromolar"/>
      <variable name="Vmyo" public_interface="in" units="microlitre"/>
      <variable name="Acap" public_interface="in" units="cm2"/>
      <variable name="Cm" public_interface="in" units="microF_per_cm2"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_NaCa</ci>
            <apply>
               <times/>
               <apply>
                  <divide/>
                  <apply>
                     <times/>
                     <apply>
                        <divide/>
                        <apply>
                           <times/>
                           <apply>
                              <divide/>
                              <apply>
                                 <times/>
                                 <ci>k_NaCa</ci>
                                 <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                              </apply>
                              <apply>
                                 <plus/>
                                 <apply>
                                    <power/>
                                    <ci>K_mNa</ci>
                                    <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">3</cn>
                                 </apply>
                                 <apply>
                                    <power/>
                                    <ci>Nao</ci>
                                    <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">3</cn>
                                 </apply>
                              </apply>
                           </apply>
                           <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                        </apply>
                        <apply>
                           <plus/>
                           <ci>K_mCa</ci>
                           <ci>Cao</ci>
                        </apply>
                     </apply>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                  </apply>
                  <apply>
                     <plus/>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                     <apply>
                        <times/>
                        <ci>k_sat</ci>
                        <apply>
                           <exp/>
                           <apply>
                              <divide/>
                              <apply>
                                 <times/>
                                 <apply>
                                    <minus/>
                                    <ci>eta</ci>
                                    <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                                 </apply>
                                 <ci>V</ci>
                                 <ci>F</ci>
                              </apply>
                              <apply>
                                 <times/>
                                 <ci>R</ci>
                                 <ci>T</ci>
                              </apply>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <minus/>
                  <apply>
                     <times/>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <times/>
                              <ci>eta</ci>
                              <ci>V</ci>
                              <ci>F</ci>
                           </apply>
                           <apply>
                              <times/>
                              <ci>R</ci>
                              <ci>T</ci>
                           </apply>
                        </apply>
                     </apply>
                     <apply>
                        <power/>
                        <ci>Nai</ci>
                        <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">3</cn>
                     </apply>
                     <ci>Cao</ci>
                  </apply>
                  <apply>
                     <times/>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <times/>
                              <apply>
                                 <minus/>
                                 <ci>eta</ci>
                                 <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                              </apply>
                              <ci>V</ci>
                              <ci>F</ci>
                           </apply>
                           <apply>
                              <times/>
                              <ci>R</ci>
                              <ci>T</ci>
                           </apply>
                        </apply>
                     </apply>
                     <apply>
                        <power/>
                        <ci>Nao</ci>
                        <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">3</cn>
                     </apply>
                     <ci>Cai</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>J_ncx</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <ci>i_NaCa</ci>
                  <ci>Acap</ci>
                  <ci>Cm</ci>
               </apply>
               <apply>
                  <times/>
                  <ci>Vmyo</ci>
                  <ci>F</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="calcium_background_current">
      <variable name="i_Cab" public_interface="out" units="picoA_per_picoF"/>
      <variable initial_value="0.00088" name="g_Cab" units="milliS_per_microF"/>
      <variable name="E_CaN" units="millivolt"/>
      <variable name="V" public_interface="in" units="millivolt"/>
      <variable name="R" public_interface="in" units="joule_per_mole_kelvin"/>
      <variable name="T" public_interface="in" units="kelvin"/>
      <variable name="F" public_interface="in" units="coulomb_per_millimole"/>
      <variable name="Cai" public_interface="in" units="micromolar"/>
      <variable name="Cao" public_interface="in" units="micromolar"/>
      <variable name="J_Cab" units="micromolar_per_millisecond"/>
      <variable name="Vmyo" public_interface="in" units="microlitre"/>
      <variable name="Acap" public_interface="in" units="cm2"/>
      <variable name="Cm" public_interface="in" units="microF_per_cm2"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_Cab</ci>
            <apply>
               <times/>
               <ci>g_Cab</ci>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>E_CaN</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>E_CaN</ci>
            <apply>
               <times/>
               <apply>
                  <divide/>
                  <apply>
                     <times/>
                     <ci>R</ci>
                     <ci>T</ci>
                  </apply>
                  <apply>
                     <times/>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">2</cn>
                     <ci>F</ci>
                  </apply>
               </apply>
               <apply>
                  <ln/>
                  <apply>
                     <divide/>
                     <ci>Cao</ci>
                     <ci>Cai</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>J_Cab</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <apply>
                     <minus/>
                     <ci>i_Cab</ci>
                  </apply>
                  <ci>Acap</ci>
                  <ci>Cm</ci>
               </apply>
               <apply>
                  <times/>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">2</cn>
                  <ci>Vmyo</ci>
                  <ci>F</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="sodium_concentration">
      <variable initial_value="16522.45" name="Nai" public_interface="out" units="micromolar"/>
      <variable name="time" public_interface="in" units="millisecond"/>
      <variable name="F" public_interface="in" units="coulomb_per_millimole"/>
      <variable name="Vmyo" public_interface="in" units="microlitre"/>
      <variable name="Cm" public_interface="in" units="microF_per_cm2"/>
      <variable name="Acap" public_interface="in" units="cm2"/>
      <variable name="i_Na" public_interface="in" units="picoA_per_picoF"/>
      <variable name="i_NaCa" public_interface="in" units="picoA_per_picoF"/>
      <variable name="i_NaK" public_interface="in" units="picoA_per_picoF"/>
      <variable name="i_Nab" public_interface="in" units="picoA_per_picoF"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>Nai</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <apply>
                     <minus/>
                     <apply>
                        <plus/>
                        <ci>i_Na</ci>
                        <ci>i_Nab</ci>
                        <apply>
                           <times/>
                           <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">3</cn>
                           <ci>i_NaK</ci>
                        </apply>
                        <apply>
                           <times/>
                           <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">3</cn>
                           <ci>i_NaCa</ci>
                        </apply>
                     </apply>
                  </apply>
                  <ci>Acap</ci>
                  <ci>Cm</ci>
               </apply>
               <apply>
                  <times/>
                  <ci>Vmyo</ci>
                  <ci>F</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="fast_sodium_current">
      <variable name="i_Na" public_interface="out" units="picoA_per_picoF"/>
      <variable name="E_Na" public_interface="out" units="millivolt"/>
      <variable initial_value="13" name="g_Na" units="milliS_per_microF"/>
      <variable initial_value="2.639399e-7" name="O_Na" units="dimensionless"/>
      <variable initial_value="0.0001581035" name="C_Na1" units="dimensionless"/>
      <variable initial_value="0.01702105" name="C_Na2" units="dimensionless"/>
      <variable name="C_Na3" units="dimensionless"/>
      <variable initial_value="0.00001799179" name="I1_Na" units="dimensionless"/>
      <variable initial_value="0.000005460299" name="I2_Na" units="dimensionless"/>
      <variable initial_value="0.0000556206" name="IF_Na" units="dimensionless"/>
      <variable initial_value="0.005985434" name="IC_Na2" units="dimensionless"/>
      <variable initial_value="0.2543133" name="IC_Na3" units="dimensionless"/>
      <variable name="alpha_Na11" units="per_millisecond"/>
      <variable name="beta_Na11" units="per_millisecond"/>
      <variable name="alpha_Na12" units="per_millisecond"/>
      <variable name="beta_Na12" units="per_millisecond"/>
      <variable name="alpha_Na13" units="per_millisecond"/>
      <variable name="beta_Na13" units="per_millisecond"/>
      <variable name="alpha_Na3" units="per_millisecond"/>
      <variable name="beta_Na3" units="per_millisecond"/>
      <variable name="alpha_Na2" units="per_millisecond"/>
      <variable name="beta_Na2" units="per_millisecond"/>
      <variable name="alpha_Na4" units="per_millisecond"/>
      <variable name="beta_Na4" units="per_millisecond"/>
      <variable name="alpha_Na5" units="per_millisecond"/>
      <variable name="beta_Na5" units="per_millisecond"/>
      <variable name="time" public_interface="in" units="millisecond"/>
      <variable name="V" public_interface="in" units="millivolt"/>
      <variable name="R" public_interface="in" units="joule_per_mole_kelvin"/>
      <variable name="T" public_interface="in" units="kelvin"/>
      <variable name="F" public_interface="in" units="coulomb_per_millimole"/>
      <variable name="Nao" public_interface="in" units="micromolar"/>
      <variable name="Nai" public_interface="in" units="micromolar"/>
      <variable name="Ko" public_interface="in" units="micromolar"/>
      <variable name="Ki" public_interface="in" units="micromolar"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_Na</ci>
            <apply>
               <times/>
               <ci>g_Na</ci>
               <ci>O_Na</ci>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>E_Na</ci>
               </apply>
            </apply>
         </apply>
         <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/>
                     <apply>
                        <plus/>
                        <apply>
                           <times/>
                           <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">0.9</cn>
                           <ci>Nao</ci>
                        </apply>
                        <apply>
                           <times/>
                           <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">0.1</cn>
                           <ci>Ko</ci>
                        </apply>
                     </apply>
                     <apply>
                        <plus/>
                        <apply>
                           <times/>
                           <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">0.9</cn>
                           <ci>Nai</ci>
                        </apply>
                        <apply>
                           <times/>
                           <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">0.1</cn>
                           <ci>Ki</ci>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>C_Na3</ci>
            <apply>
               <minus/>
               <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
               <apply>
                  <plus/>
                  <ci>O_Na</ci>
                  <ci>C_Na1</ci>
                  <ci>C_Na2</ci>
                  <ci>IF_Na</ci>
                  <ci>I1_Na</ci>
                  <ci>I2_Na</ci>
                  <ci>IC_Na2</ci>
                  <ci>IC_Na3</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>C_Na2</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <plus/>
                  <apply>
                     <times/>
                     <ci>alpha_Na11</ci>
                     <ci>C_Na3</ci>
                  </apply>
                  <apply>
                     <times/>
                     <ci>beta_Na12</ci>
                     <ci>C_Na1</ci>
                  </apply>
                  <apply>
                     <times/>
                     <ci>alpha_Na3</ci>
                     <ci>IC_Na2</ci>
                  </apply>
               </apply>
               <apply>
                  <plus/>
                  <apply>
                     <times/>
                     <ci>beta_Na11</ci>
                     <ci>C_Na2</ci>
                  </apply>
                  <apply>
                     <times/>
                     <ci>alpha_Na12</ci>
                     <ci>C_Na2</ci>
                  </apply>
                  <apply>
                     <times/>
                     <ci>beta_Na3</ci>
                     <ci>C_Na2</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>C_Na1</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <plus/>
                  <apply>
                     <times/>
                     <ci>alpha_Na12</ci>
                     <ci>C_Na2</ci>
                  </apply>
                  <apply>
                     <times/>
                     <ci>beta_Na13</ci>
                     <ci>O_Na</ci>
                  </apply>
                  <apply>
                     <times/>
                     <ci>alpha_Na3</ci>
                     <ci>IF_Na</ci>
                  </apply>
               </apply>
               <apply>
                  <plus/>
                  <apply>
                     <times/>
                     <ci>beta_Na12</ci>
                     <ci>C_Na1</ci>
                  </apply>
                  <apply>
                     <times/>
                     <ci>alpha_Na13</ci>
                     <ci>C_Na1</ci>
                  </apply>
                  <apply>
                     <times/>
                     <ci>beta_Na3</ci>
                     <ci>C_Na1</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>O_Na</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <plus/>
                  <apply>
                     <times/>
                     <ci>alpha_Na13</ci>
                     <ci>C_Na1</ci>
                  </apply>
                  <apply>
                     <times/>
                     <ci>beta_Na2</ci>
                     <ci>IF_Na</ci>
                  </apply>
               </apply>
               <apply>
                  <plus/>
                  <apply>
                     <times/>
                     <ci>beta_Na13</ci>
                     <ci>O_Na</ci>
                  </apply>
                  <apply>
                     <times/>
                     <ci>alpha_Na2</ci>
                     <ci>O_Na</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>IF_Na</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <plus/>
                  <apply>
                     <times/>
                     <ci>alpha_Na2</ci>
                     <ci>O_Na</ci>
                  </apply>
                  <apply>
                     <times/>
                     <ci>beta_Na3</ci>
                     <ci>C_Na1</ci>
                  </apply>
                  <apply>
                     <times/>
                     <ci>beta_Na4</ci>
                     <ci>I1_Na</ci>
                  </apply>
                  <apply>
                     <times/>
                     <ci>alpha_Na12</ci>
                     <ci>IC_Na2</ci>
                  </apply>
               </apply>
               <apply>
                  <plus/>
                  <apply>
                     <times/>
                     <ci>beta_Na2</ci>
                     <ci>IF_Na</ci>
                  </apply>
                  <apply>
                     <times/>
                     <ci>alpha_Na3</ci>
                     <ci>IF_Na</ci>
                  </apply>
                  <apply>
                     <times/>
                     <ci>alpha_Na4</ci>
                     <ci>IF_Na</ci>
                  </apply>
                  <apply>
                     <times/>
                     <ci>beta_Na12</ci>
                     <ci>IF_Na</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>I1_Na</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <plus/>
                  <apply>
                     <times/>
                     <ci>alpha_Na4</ci>
                     <ci>IF_Na</ci>
                  </apply>
                  <apply>
                     <times/>
                     <ci>beta_Na5</ci>
                     <ci>I2_Na</ci>
                  </apply>
               </apply>
               <apply>
                  <plus/>
                  <apply>
                     <times/>
                     <ci>beta_Na4</ci>
                     <ci>I1_Na</ci>
                  </apply>
                  <apply>
                     <times/>
                     <ci>alpha_Na5</ci>
                     <ci>I1_Na</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>I2_Na</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>alpha_Na5</ci>
                  <ci>I1_Na</ci>
               </apply>
               <apply>
                  <times/>
                  <ci>beta_Na5</ci>
                  <ci>I2_Na</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>IC_Na2</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <plus/>
                  <apply>
                     <times/>
                     <ci>alpha_Na11</ci>
                     <ci>IC_Na3</ci>
                  </apply>
                  <apply>
                     <times/>
                     <ci>beta_Na12</ci>
                     <ci>IF_Na</ci>
                  </apply>
                  <apply>
                     <times/>
                     <ci>beta_Na3</ci>
                     <ci>IC_Na2</ci>
                  </apply>
               </apply>
               <apply>
                  <plus/>
                  <apply>
                     <times/>
                     <ci>beta_Na11</ci>
                     <ci>IC_Na2</ci>
                  </apply>
                  <apply>
                     <times/>
                     <ci>alpha_Na12</ci>
                     <ci>IC_Na2</ci>
                  </apply>
                  <apply>
                     <times/>
                     <ci>alpha_Na3</ci>
                     <ci>IC_Na2</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>IC_Na3</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <plus/>
                  <apply>
                     <times/>
                     <ci>beta_Na11</ci>
                     <ci>IC_Na2</ci>
                  </apply>
                  <apply>
                     <times/>
                     <ci>beta_Na3</ci>
                     <ci>C_Na3</ci>
                  </apply>
               </apply>
               <apply>
                  <plus/>
                  <apply>
                     <times/>
                     <ci>alpha_Na11</ci>
                     <ci>IC_Na3</ci>
                  </apply>
                  <apply>
                     <times/>
                     <ci>alpha_Na3</ci>
                     <ci>IC_Na3</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>alpha_Na11</ci>
            <apply>
               <divide/>
               <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="per_millisecond">3.802</cn>
               <apply>
                  <plus/>
                  <apply>
                     <times/>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">0.1027</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <minus/>
                              <apply>
                                 <plus/>
                                 <ci>V</ci>
                                 <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">2.5</cn>
                              </apply>
                           </apply>
                           <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">17</cn>
                        </apply>
                     </apply>
                  </apply>
                  <apply>
                     <times/>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">0.2</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <minus/>
                              <apply>
                                 <plus/>
                                 <ci>V</ci>
                                 <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">2.5</cn>
                              </apply>
                           </apply>
                           <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">150</cn>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>alpha_Na12</ci>
            <apply>
               <divide/>
               <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="per_millisecond">3.802</cn>
               <apply>
                  <plus/>
                  <apply>
                     <times/>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">0.1027</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <minus/>
                              <apply>
                                 <plus/>
                                 <ci>V</ci>
                                 <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">2.5</cn>
                              </apply>
                           </apply>
                           <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">15</cn>
                        </apply>
                     </apply>
                  </apply>
                  <apply>
                     <times/>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">0.23</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <minus/>
                              <apply>
                                 <plus/>
                                 <ci>V</ci>
                                 <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">2.5</cn>
                              </apply>
                           </apply>
                           <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">150</cn>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>alpha_Na13</ci>
            <apply>
               <divide/>
               <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="per_millisecond">3.802</cn>
               <apply>
                  <plus/>
                  <apply>
                     <times/>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">0.1027</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <minus/>
                              <apply>
                                 <plus/>
                                 <ci>V</ci>
                                 <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">2.5</cn>
                              </apply>
                           </apply>
                           <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">12</cn>
                        </apply>
                     </apply>
                  </apply>
                  <apply>
                     <times/>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">0.25</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <minus/>
                              <apply>
                                 <plus/>
                                 <ci>V</ci>
                                 <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">2.5</cn>
                              </apply>
                           </apply>
                           <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">150</cn>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>beta_Na11</ci>
            <apply>
               <times/>
               <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="per_millisecond">0.1917</cn>
               <apply>
                  <exp/>
                  <apply>
                     <divide/>
                     <apply>
                        <minus/>
                        <apply>
                           <plus/>
                           <ci>V</ci>
                           <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">2.5</cn>
                        </apply>
                     </apply>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">20.3</cn>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>beta_Na12</ci>
            <apply>
               <times/>
               <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="per_millisecond">0.2</cn>
               <apply>
                  <exp/>
                  <apply>
                     <divide/>
                     <apply>
                        <minus/>
                        <apply>
                           <minus/>
                           <ci>V</ci>
                           <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">2.5</cn>
                        </apply>
                     </apply>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">20.3</cn>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>beta_Na13</ci>
            <apply>
               <times/>
               <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="per_millisecond">0.22</cn>
               <apply>
                  <exp/>
                  <apply>
                     <divide/>
                     <apply>
                        <minus/>
                        <apply>
                           <minus/>
                           <ci>V</ci>
                           <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">7.5</cn>
                        </apply>
                     </apply>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">20.3</cn>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>alpha_Na3</ci>
            <apply>
               <times/>
               <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" type="e-notation" cellml:units="per_millisecond">7<sep/>-7</cn>
               <apply>
                  <exp/>
                  <apply>
                     <divide/>
                     <apply>
                        <minus/>
                        <apply>
                           <plus/>
                           <ci>V</ci>
                           <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">7</cn>
                        </apply>
                     </apply>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">7.7</cn>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>beta_Na3</ci>
            <apply>
               <plus/>
               <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="per_millisecond">0.0084</cn>
               <apply>
                  <times/>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="per_millivolt_millisecond">0.00002</cn>
                  <apply>
                     <plus/>
                     <ci>V</ci>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">7</cn>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>alpha_Na2</ci>
            <apply>
               <divide/>
               <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="per_millisecond">1</cn>
               <apply>
                  <plus/>
                  <apply>
                     <times/>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">0.188495</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <minus/>
                              <apply>
                                 <plus/>
                                 <ci>V</ci>
                                 <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">7</cn>
                              </apply>
                           </apply>
                           <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">16.6</cn>
                        </apply>
                     </apply>
                  </apply>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">0.393956</cn>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>beta_Na2</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <ci>alpha_Na13</ci>
                  <ci>alpha_Na2</ci>
                  <ci>alpha_Na3</ci>
               </apply>
               <apply>
                  <times/>
                  <ci>beta_Na13</ci>
                  <ci>beta_Na3</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>alpha_Na4</ci>
            <apply>
               <divide/>
               <ci>alpha_Na2</ci>
               <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1000</cn>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>beta_Na4</ci>
            <ci>alpha_Na3</ci>
         </apply>
         <apply>
            <eq/>
            <ci>alpha_Na5</ci>
            <apply>
               <divide/>
               <ci>alpha_Na2</ci>
               <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">95000</cn>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>beta_Na5</ci>
            <apply>
               <divide/>
               <ci>alpha_Na3</ci>
               <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">50</cn>
            </apply>
         </apply>
      </math>
   </component>
   <component name="sodium_background_current">
      <variable name="i_Nab" public_interface="out" units="picoA_per_picoF"/>
      <variable initial_value="0.0026" name="g_Nab" units="milliS_per_microF"/>
      <variable name="E_Na" public_interface="in" units="millivolt"/>
      <variable name="V" public_interface="in" units="millivolt"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_Nab</ci>
            <apply>
               <times/>
               <ci>g_Nab</ci>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>E_Na</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="potassium_concentration">
      <variable initial_value="141474" name="Ki" public_interface="out" units="micromolar"/>
      <variable name="time" public_interface="in" units="millisecond"/>
      <variable name="F" public_interface="in" units="coulomb_per_millimole"/>
      <variable name="Vmyo" public_interface="in" units="microlitre"/>
      <variable name="Cm" public_interface="in" units="microF_per_cm2"/>
      <variable name="Acap" public_interface="in" units="cm2"/>
      <variable name="i_Kto_f" public_interface="in" units="picoA_per_picoF"/>
      <variable name="i_Kto_s" public_interface="in" units="picoA_per_picoF"/>
      <variable name="i_K1" public_interface="in" units="picoA_per_picoF"/>
      <variable name="i_Ks" public_interface="in" units="picoA_per_picoF"/>
      <variable name="i_Kss" public_interface="in" units="picoA_per_picoF"/>
      <variable name="i_Kur" public_interface="in" units="picoA_per_picoF"/>
      <variable name="i_Kr" public_interface="in" units="picoA_per_picoF"/>
      <variable name="i_NaK" public_interface="in" units="picoA_per_picoF"/>
      <variable name="i_Stim" public_interface="in" units="picoA_per_picoF"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>Ki</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <apply>
                     <minus/>
                     <apply>
                        <minus/>
                        <apply>
                           <plus/>
                           <ci>i_Stim</ci>
                           <ci>i_Kto_f</ci>
                           <ci>i_Kto_s</ci>
                           <ci>i_K1</ci>
                           <ci>i_Ks</ci>
                           <ci>i_Kss</ci>
                           <ci>i_Kur</ci>
                           <ci>i_Kr</ci>
                        </apply>
                        <apply>
                           <times/>
                           <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">2</cn>
                           <ci>i_NaK</ci>
                        </apply>
                     </apply>
                  </apply>
                  <ci>Acap</ci>
                  <ci>Cm</ci>
               </apply>
               <apply>
                  <times/>
                  <ci>Vmyo</ci>
                  <ci>F</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="fast_transient_outward_K_I">
      <variable name="i_Kto_f" public_interface="out" units="picoA_per_picoF"/>
      <variable name="E_K" public_interface="out" units="millivolt"/>
      <variable initial_value="0.4067" name="g_Kto_f" units="milliS_per_microF"/>
      <variable initial_value="0.001937245" name="ato_f" units="dimensionless"/>
      <variable initial_value="0.9999985" name="ito_f" units="dimensionless"/>
      <variable name="alpha_a" units="per_millisecond"/>
      <variable name="beta_a" units="per_millisecond"/>
      <variable name="alpha_i" units="per_millisecond"/>
      <variable name="beta_i" units="per_millisecond"/>
      <variable name="V" public_interface="in" units="millivolt"/>
      <variable name="R" public_interface="in" units="joule_per_mole_kelvin"/>
      <variable name="T" public_interface="in" units="kelvin"/>
      <variable name="F" public_interface="in" units="coulomb_per_millimole"/>
      <variable name="Ki" public_interface="in" units="micromolar"/>
      <variable name="Ko" public_interface="in" units="micromolar"/>
      <variable name="time" public_interface="in" units="millisecond"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_Kto_f</ci>
            <apply>
               <times/>
               <ci>g_Kto_f</ci>
               <apply>
                  <power/>
                  <ci>ato_f</ci>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">3</cn>
               </apply>
               <ci>ito_f</ci>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>E_K</ci>
               </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>Ko</ci>
                     <ci>Ki</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>ato_f</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>alpha_a</ci>
                  <apply>
                     <minus/>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                     <ci>ato_f</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>beta_a</ci>
                  <ci>ato_f</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>ito_f</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>alpha_i</ci>
                  <apply>
                     <minus/>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                     <ci>ito_f</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>beta_i</ci>
                  <ci>ito_f</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>alpha_a</ci>
            <apply>
               <times/>
               <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="per_millisecond">0.18064</cn>
               <apply>
                  <exp/>
                  <apply>
                     <times/>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="per_millivolt">0.03577</cn>
                     <apply>
                        <plus/>
                        <ci>V</ci>
                        <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">30</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>beta_a</ci>
            <apply>
               <times/>
               <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="per_millisecond">0.3956</cn>
               <apply>
                  <exp/>
                  <apply>
                     <times/>
                     <apply>
                        <minus/>
                        <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="per_millivolt">0.06237</cn>
                     </apply>
                     <apply>
                        <plus/>
                        <ci>V</ci>
                        <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">30</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>alpha_i</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="per_millisecond">0.000152</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <minus/>
                           <apply>
                              <plus/>
                              <ci>V</ci>
                              <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">13.5</cn>
                           </apply>
                        </apply>
                        <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">7</cn>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <plus/>
                  <apply>
                     <times/>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">0.0067083</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <minus/>
                              <apply>
                                 <plus/>
                                 <ci>V</ci>
                                 <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">33.5</cn>
                              </apply>
                           </apply>
                           <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">7</cn>
                        </apply>
                     </apply>
                  </apply>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>beta_i</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="per_millisecond">0.00095</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <plus/>
                           <ci>V</ci>
                           <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">33.5</cn>
                        </apply>
                        <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">7</cn>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <plus/>
                  <apply>
                     <times/>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">0.051335</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <plus/>
                              <ci>V</ci>
                              <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">33.5</cn>
                           </apply>
                           <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">7</cn>
                        </apply>
                     </apply>
                  </apply>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="slow_transient_outward_K_I">
      <variable name="i_Kto_s" public_interface="out" units="picoA_per_picoF"/>
      <variable name="ass" public_interface="out" units="dimensionless"/>
      <variable name="iss" public_interface="out" units="dimensionless"/>
      <variable initial_value="0" name="g_Kto_s" units="milliS_per_microF"/>
      <variable initial_value="0.02000568" name="ato_s" units="dimensionless"/>
      <variable initial_value="0.9308568" name="ito_s" units="dimensionless"/>
      <variable name="tau_ta_s" units="millisecond"/>
      <variable name="tau_ti_s" units="millisecond"/>
      <variable name="E_K" public_interface="in" units="millivolt"/>
      <variable name="V" public_interface="in" units="millivolt"/>
      <variable name="time" public_interface="in" units="millisecond"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_Kto_s</ci>
            <apply>
               <times/>
               <ci>g_Kto_s</ci>
               <ci>ato_s</ci>
               <ci>ito_s</ci>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>E_K</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>ato_s</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>ass</ci>
                  <ci>ato_s</ci>
               </apply>
               <ci>tau_ta_s</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>ito_s</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>iss</ci>
                  <ci>ito_s</ci>
               </apply>
               <ci>tau_ti_s</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>ass</ci>
            <apply>
               <divide/>
               <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
               <apply>
                  <plus/>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <minus/>
                           <apply>
                              <plus/>
                              <ci>V</ci>
                              <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">22.5</cn>
                           </apply>
                        </apply>
                        <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">7.7</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>iss</ci>
            <apply>
               <divide/>
               <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
               <apply>
                  <plus/>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <plus/>
                           <ci>V</ci>
                           <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">45.2</cn>
                        </apply>
                        <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">5.7</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>tau_ta_s</ci>
            <apply>
               <plus/>
               <apply>
                  <times/>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millisecond">0.493</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <times/>
                        <apply>
                           <minus/>
                           <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="per_millivolt">0.0629</cn>
                        </apply>
                        <ci>V</ci>
                     </apply>
                  </apply>
               </apply>
               <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millisecond">2.058</cn>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>tau_ti_s</ci>
            <apply>
               <plus/>
               <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millisecond">270</cn>
               <apply>
                  <divide/>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millisecond">1050</cn>
                  <apply>
                     <plus/>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <plus/>
                              <ci>V</ci>
                              <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">45.2</cn>
                           </apply>
                           <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">5.7</cn>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="time_independent_K_I">
      <variable name="i_K1" public_interface="out" units="picoA_per_picoF"/>
      <variable name="Ko" public_interface="in" units="micromolar"/>
      <variable name="E_K" public_interface="in" units="millivolt"/>
      <variable name="V" public_interface="in" units="millivolt"/>
      <variable name="time" public_interface="in" units="millisecond"/>
      <variable initial_value="0.2938" name="g_K1" units="milliS_per_microF"/>
      <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>Ko</ci>
                     </apply>
                     <apply>
                        <plus/>
                        <ci>Ko</ci>
                        <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="micromolar">210</cn>
                     </apply>
                  </apply>
                  <apply>
                     <minus/>
                     <ci>V</ci>
                     <ci>E_K</ci>
                  </apply>
               </apply>
               <apply>
                  <plus/>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <times/>
                        <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="per_millivolt">0.0896</cn>
                        <apply>
                           <minus/>
                           <ci>V</ci>
                           <ci>E_K</ci>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="slow_delayed_rectifier_K_I">
      <variable name="i_Ks" public_interface="out" units="picoA_per_picoF"/>
      <variable initial_value="0.00575" name="g_Ks" units="milliS_per_microF"/>
      <variable initial_value="0.002206261" name="nKs" units="dimensionless"/>
      <variable name="alpha_n" units="per_millisecond"/>
      <variable name="beta_n" units="per_millisecond"/>
      <variable name="E_K" public_interface="in" units="millivolt"/>
      <variable name="V" public_interface="in" units="millivolt"/>
      <variable name="time" public_interface="in" units="millisecond"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_Ks</ci>
            <apply>
               <times/>
               <ci>g_Ks</ci>
               <apply>
                  <power/>
                  <ci>nKs</ci>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">2</cn>
               </apply>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>E_K</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>nKs</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>alpha_n</ci>
                  <apply>
                     <minus/>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                     <ci>nKs</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>beta_n</ci>
                  <ci>nKs</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>alpha_n</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="per_millivolt_millisecond">0.00000481333</cn>
                  <apply>
                     <plus/>
                     <ci>V</ci>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">26.5</cn>
                  </apply>
               </apply>
               <apply>
                  <minus/>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <times/>
                        <apply>
                           <minus/>
                           <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="per_millivolt">0.128</cn>
                        </apply>
                        <apply>
                           <plus/>
                           <ci>V</ci>
                           <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">26.5</cn>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>beta_n</ci>
            <apply>
               <times/>
               <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="per_millisecond">0.0000953333</cn>
               <apply>
                  <exp/>
                  <apply>
                     <times/>
                     <apply>
                        <minus/>
                        <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="per_millivolt">0.038</cn>
                     </apply>
                     <apply>
                        <plus/>
                        <ci>V</ci>
                        <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">26.5</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="ultra_rapidly_activating_delayed_rectifier_K_I">
      <variable name="i_Kur" public_interface="out" units="picoA_per_picoF"/>
      <variable initial_value="0.16" name="g_Kur" units="milliS_per_microF"/>
      <variable initial_value="0.02000568" name="aur" units="dimensionless"/>
      <variable initial_value="0.9822006" name="iur" units="dimensionless"/>
      <variable name="tau_aur" units="millisecond"/>
      <variable name="tau_iur" units="millisecond"/>
      <variable name="E_K" public_interface="in" units="millivolt"/>
      <variable name="V" public_interface="in" units="millivolt"/>
      <variable name="time" public_interface="in" units="millisecond"/>
      <variable name="ass" public_interface="in" units="dimensionless"/>
      <variable name="iss" public_interface="in" units="dimensionless"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_Kur</ci>
            <apply>
               <times/>
               <ci>g_Kur</ci>
               <ci>aur</ci>
               <ci>iur</ci>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>E_K</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>aur</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>ass</ci>
                  <ci>aur</ci>
               </apply>
               <ci>tau_aur</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>iur</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>iss</ci>
                  <ci>iur</ci>
               </apply>
               <ci>tau_iur</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>tau_aur</ci>
            <apply>
               <plus/>
               <apply>
                  <times/>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millisecond">0.493</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <times/>
                        <apply>
                           <minus/>
                           <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="per_millivolt">0.0629</cn>
                        </apply>
                        <ci>V</ci>
                     </apply>
                  </apply>
               </apply>
               <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millisecond">2.058</cn>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>tau_iur</ci>
            <apply>
               <minus/>
               <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millisecond">1200</cn>
               <apply>
                  <divide/>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millisecond">170</cn>
                  <apply>
                     <plus/>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <plus/>
                              <ci>V</ci>
                              <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">45.2</cn>
                           </apply>
                           <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">5.7</cn>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="non_inactivating_steady_state_K_I">
      <variable name="i_Kss" public_interface="out" units="picoA_per_picoF"/>
      <variable initial_value="0.05" name="g_Kss" units="milliS_per_microF"/>
      <variable initial_value="0.8883113" name="aKss" units="dimensionless"/>
      <variable initial_value="1" name="iKss" units="dimensionless"/>
      <variable name="tau_Kss" units="millisecond"/>
      <variable name="E_K" public_interface="in" units="millivolt"/>
      <variable name="V" public_interface="in" units="millivolt"/>
      <variable name="time" public_interface="in" units="millisecond"/>
      <variable name="ass" public_interface="in" units="dimensionless"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_Kss</ci>
            <apply>
               <times/>
               <ci>g_Kss</ci>
               <ci>aKss</ci>
               <ci>iKss</ci>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>E_K</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>aKss</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>ass</ci>
                  <ci>aKss</ci>
               </apply>
               <ci>tau_Kss</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>iKss</ci>
            </apply>
            <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="per_millisecond">0</cn>
         </apply>
         <apply>
            <eq/>
            <ci>tau_Kss</ci>
            <apply>
               <plus/>
               <apply>
                  <times/>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millisecond">39.3</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <times/>
                        <apply>
                           <minus/>
                           <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="per_millivolt">0.0862</cn>
                        </apply>
                        <ci>V</ci>
                     </apply>
                  </apply>
               </apply>
               <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millisecond">13.17</cn>
            </apply>
         </apply>
      </math>
   </component>
   <component name="rapid_delayed_rectifier_K_I">
      <variable name="i_Kr" public_interface="out" units="picoA_per_picoF"/>
      <variable initial_value="0.078" name="g_Kr" units="milliS_per_microF"/>
      <variable initial_value="0.0004858865" name="O_K" units="dimensionless"/>
      <variable initial_value="0.0007799137" name="C_K1" units="dimensionless"/>
      <variable initial_value="0.0005301217" name="C_K2" units="dimensionless"/>
      <variable name="C_K0" units="dimensionless"/>
      <variable initial_value="0.00007519518" name="I_K" units="dimensionless"/>
      <variable name="alpha_a0" units="per_millisecond"/>
      <variable name="beta_a0" units="per_millisecond"/>
      <variable initial_value="0.036778" name="kb" units="per_millisecond"/>
      <variable initial_value="0.023761" name="kf" units="per_millisecond"/>
      <variable name="alpha_a1" units="per_millisecond"/>
      <variable name="beta_a1" units="per_millisecond"/>
      <variable name="alpha_i" units="per_millisecond"/>
      <variable name="beta_i" units="per_millisecond"/>
      <variable name="time" public_interface="in" units="millisecond"/>
      <variable name="V" public_interface="in" units="millivolt"/>
      <variable name="R" public_interface="in" units="joule_per_mole_kelvin"/>
      <variable name="T" public_interface="in" units="kelvin"/>
      <variable name="F" public_interface="in" units="coulomb_per_millimole"/>
      <variable name="Nao" public_interface="in" units="micromolar"/>
      <variable name="Nai" public_interface="in" units="micromolar"/>
      <variable name="Ko" public_interface="in" units="micromolar"/>
      <variable name="Ki" public_interface="in" units="micromolar"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_Kr</ci>
            <apply>
               <times/>
               <ci>g_Kr</ci>
               <ci>O_K</ci>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <apply>
                     <times/>
                     <apply>
                        <divide/>
                        <apply>
                           <times/>
                           <ci>R</ci>
                           <ci>T</ci>
                        </apply>
                        <ci>F</ci>
                     </apply>
                     <apply>
                        <ln/>
                        <apply>
                           <divide/>
                           <apply>
                              <plus/>
                              <apply>
                                 <times/>
                                 <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">0.98</cn>
                                 <ci>Ko</ci>
                              </apply>
                              <apply>
                                 <times/>
                                 <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">0.02</cn>
                                 <ci>Nao</ci>
                              </apply>
                           </apply>
                           <apply>
                              <plus/>
                              <apply>
                                 <times/>
                                 <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">0.98</cn>
                                 <ci>Ki</ci>
                              </apply>
                              <apply>
                                 <times/>
                                 <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">0.02</cn>
                                 <ci>Nai</ci>
                              </apply>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>C_K0</ci>
            <apply>
               <minus/>
               <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
               <apply>
                  <plus/>
                  <ci>C_K1</ci>
                  <ci>C_K2</ci>
                  <ci>O_K</ci>
                  <ci>I_K</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>C_K2</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <plus/>
                  <apply>
                     <times/>
                     <ci>kf</ci>
                     <ci>C_K1</ci>
                  </apply>
                  <apply>
                     <times/>
                     <ci>beta_a1</ci>
                     <ci>O_K</ci>
                  </apply>
               </apply>
               <apply>
                  <plus/>
                  <apply>
                     <times/>
                     <ci>kb</ci>
                     <ci>C_K2</ci>
                  </apply>
                  <apply>
                     <times/>
                     <ci>alpha_a1</ci>
                     <ci>C_K2</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>C_K1</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <plus/>
                  <apply>
                     <times/>
                     <ci>alpha_a0</ci>
                     <ci>C_K0</ci>
                  </apply>
                  <apply>
                     <times/>
                     <ci>kb</ci>
                     <ci>C_K2</ci>
                  </apply>
               </apply>
               <apply>
                  <plus/>
                  <apply>
                     <times/>
                     <ci>beta_a0</ci>
                     <ci>C_K1</ci>
                  </apply>
                  <apply>
                     <times/>
                     <ci>kf</ci>
                     <ci>C_K1</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>O_K</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <plus/>
                  <apply>
                     <times/>
                     <ci>alpha_a1</ci>
                     <ci>C_K2</ci>
                  </apply>
                  <apply>
                     <times/>
                     <ci>beta_i</ci>
                     <ci>I_K</ci>
                  </apply>
               </apply>
               <apply>
                  <plus/>
                  <apply>
                     <times/>
                     <ci>beta_a1</ci>
                     <ci>O_K</ci>
                  </apply>
                  <apply>
                     <times/>
                     <ci>alpha_i</ci>
                     <ci>O_K</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>I_K</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>alpha_i</ci>
                  <ci>O_K</ci>
               </apply>
               <apply>
                  <times/>
                  <ci>beta_i</ci>
                  <ci>I_K</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>alpha_a0</ci>
            <apply>
               <times/>
               <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="per_millisecond">0.022348</cn>
               <apply>
                  <exp/>
                  <apply>
                     <times/>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="per_millivolt">0.01176</cn>
                     <ci>V</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>beta_a0</ci>
            <apply>
               <times/>
               <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="per_millisecond">0.047002</cn>
               <apply>
                  <exp/>
                  <apply>
                     <times/>
                     <apply>
                        <minus/>
                        <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="per_millivolt">0.0631</cn>
                     </apply>
                     <ci>V</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>alpha_a1</ci>
            <apply>
               <times/>
               <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="per_millisecond">0.013733</cn>
               <apply>
                  <exp/>
                  <apply>
                     <times/>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="per_millivolt">0.038198</cn>
                     <ci>V</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>beta_a1</ci>
            <apply>
               <times/>
               <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="per_millisecond">0.0000689</cn>
               <apply>
                  <exp/>
                  <apply>
                     <times/>
                     <apply>
                        <minus/>
                        <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="per_millivolt">0.04178</cn>
                     </apply>
                     <ci>V</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>alpha_i</ci>
            <apply>
               <times/>
               <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="per_millisecond">0.090821</cn>
               <apply>
                  <exp/>
                  <apply>
                     <times/>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="per_millivolt">0.023391</cn>
                     <apply>
                        <plus/>
                        <ci>V</ci>
                        <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">5</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>beta_i</ci>
            <apply>
               <times/>
               <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="per_millisecond">0.006497</cn>
               <apply>
                  <exp/>
                  <apply>
                     <times/>
                     <apply>
                        <minus/>
                        <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="per_millivolt">0.03268</cn>
                     </apply>
                     <apply>
                        <plus/>
                        <ci>V</ci>
                        <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">5</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="sodium_potassium_pump_current">
      <variable name="i_NaK" public_interface="out" units="picoA_per_picoF"/>
      <variable initial_value="1.66" name="i_NaK_max" units="picoA_per_picoF"/>
      <variable initial_value="21000" name="Km_Nai" units="micromolar"/>
      <variable initial_value="1500" name="Km_Ko" units="micromolar"/>
      <variable name="f_NaK" units="dimensionless"/>
      <variable name="sigma" units="dimensionless"/>
      <variable name="V" public_interface="in" units="millivolt"/>
      <variable name="R" public_interface="in" units="joule_per_mole_kelvin"/>
      <variable name="T" public_interface="in" units="kelvin"/>
      <variable name="F" public_interface="in" units="coulomb_per_millimole"/>
      <variable name="Nai" public_interface="in" units="micromolar"/>
      <variable name="Nao" public_interface="in" units="micromolar"/>
      <variable name="Ko" public_interface="in" units="micromolar"/>
      <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>f_NaK</ci>
                        <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                     </apply>
                     <apply>
                        <plus/>
                        <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                        <apply>
                           <power/>
                           <apply>
                              <divide/>
                              <ci>Km_Nai</ci>
                              <ci>Nai</ci>
                           </apply>
                           <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1.5</cn>
                        </apply>
                     </apply>
                  </apply>
                  <ci>Ko</ci>
               </apply>
               <apply>
                  <plus/>
                  <ci>Ko</ci>
                  <ci>Km_Ko</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>f_NaK</ci>
            <apply>
               <divide/>
               <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
               <apply>
                  <plus/>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                  <apply>
                     <times/>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">0.1245</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <times/>
                              <apply>
                                 <minus/>
                                 <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">0.1</cn>
                              </apply>
                              <ci>V</ci>
                              <ci>F</ci>
                           </apply>
                           <apply>
                              <times/>
                              <ci>R</ci>
                              <ci>T</ci>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
                  <apply>
                     <times/>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">0.0365</cn>
                     <ci>sigma</ci>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <times/>
                              <apply>
                                 <minus/>
                                 <ci>V</ci>
                              </apply>
                              <ci>F</ci>
                           </apply>
                           <apply>
                              <times/>
                              <ci>R</ci>
                              <ci>T</ci>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>sigma</ci>
            <apply>
               <times/>
               <apply>
                  <divide/>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">7</cn>
               </apply>
               <apply>
                  <minus/>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <ci>Nao</ci>
                        <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="micromolar">67300</cn>
                     </apply>
                  </apply>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="calcium_activated_chloride_current">
      <variable name="i_ClCa" public_interface="out" units="picoA_per_picoF"/>
      <variable initial_value="10" name="g_ClCa" units="milliS_per_microF"/>
      <variable name="O_ClCa" units="dimensionless"/>
      <variable initial_value="-40" name="E_Cl" units="millivolt"/>
      <variable initial_value="10" name="Km_Cl" units="micromolar"/>
      <variable name="V" public_interface="in" units="millivolt"/>
      <variable name="Cai" public_interface="in" units="micromolar"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_ClCa</ci>
            <apply>
               <times/>
               <apply>
                  <divide/>
                  <apply>
                     <times/>
                     <ci>g_ClCa</ci>
                     <ci>O_ClCa</ci>
                     <ci>Cai</ci>
                  </apply>
                  <apply>
                     <plus/>
                     <ci>Cai</ci>
                     <ci>Km_Cl</ci>
                  </apply>
               </apply>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>E_Cl</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>O_ClCa</ci>
            <apply>
               <divide/>
               <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">0.2</cn>
               <apply>
                  <plus/>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <minus/>
                           <apply>
                              <minus/>
                              <ci>V</ci>
                              <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">46.7</cn>
                           </apply>
                        </apply>
                        <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">7.8</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="transmembrane_currents"/>
   <component name="intracellular_currents"/>
   <component name="Ions_n_reversal_potentials"/>
   <group>
      <relationship_ref relationship="containment"/>
      <component_ref component="membrane">
         <component_ref component="transmembrane_currents">
            <component_ref component="L_type_calcium_current"/>
            <component_ref component="calcium_pump_current"/>
            <component_ref component="sodium_calcium_exchange_current"/>
            <component_ref component="calcium_background_current"/>
            <component_ref component="fast_sodium_current"/>
            <component_ref component="sodium_background_current"/>
            <component_ref component="fast_transient_outward_K_I"/>
            <component_ref component="slow_transient_outward_K_I"/>
            <component_ref component="time_independent_K_I"/>
            <component_ref component="slow_delayed_rectifier_K_I"/>
            <component_ref component="ultra_rapidly_activating_delayed_rectifier_K_I"/>
            <component_ref component="non_inactivating_steady_state_K_I"/>
            <component_ref component="rapid_delayed_rectifier_K_I"/>
            <component_ref component="sodium_potassium_pump_current"/>
            <component_ref component="calcium_activated_chloride_current"/>
         </component_ref>
         <component_ref component="Ions_n_reversal_potentials">
            <component_ref component="calcium_concentration"/>
            <component_ref component="sodium_concentration"/>
            <component_ref component="potassium_concentration"/>
         </component_ref>
         <component_ref component="intracellular_currents">
            <component_ref component="ryanodine_receptors"/>
            <component_ref component="calcium_fluxes"/>
         </component_ref>
      </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="calcium_concentration" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="calcium_fluxes" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="ryanodine_receptors" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="L_type_calcium_current" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="sodium_concentration" 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="potassium_concentration" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="fast_transient_outward_K_I" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="slow_transient_outward_K_I" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="slow_delayed_rectifier_K_I" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="time_independent_K_I" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="ultra_rapidly_activating_delayed_rectifier_K_I" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="non_inactivating_steady_state_K_I" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="rapid_delayed_rectifier_K_I" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="membrane" component_2="L_type_calcium_current"/>
      <map_variables variable_1="i_CaL" variable_2="i_CaL"/>
      <map_variables variable_1="V" variable_2="V"/>
      <map_variables variable_1="F" variable_2="F"/>
      <map_variables variable_1="T" variable_2="T"/>
      <map_variables variable_1="Vss" variable_2="Vss"/>
      <map_variables variable_1="Acap" variable_2="Acap"/>
      <map_variables variable_1="R" variable_2="R"/>
      <map_variables variable_1="Cm" variable_2="Cm"/>
      <map_variables variable_1="Cao" variable_2="Cao"/>
   </connection>
   <connection>
      <map_components component_1="membrane" component_2="calcium_pump_current"/>
      <map_variables variable_1="i_pCa" variable_2="i_pCa"/>
      <map_variables variable_1="F" variable_2="F"/>
      <map_variables variable_1="Vmyo" variable_2="Vmyo"/>
      <map_variables variable_1="Acap" variable_2="Acap"/>
      <map_variables variable_1="Cm" variable_2="Cm"/>
   </connection>
   <connection>
      <map_components component_1="membrane" component_2="sodium_calcium_exchange_current"/>
      <map_variables variable_1="i_NaCa" variable_2="i_NaCa"/>
      <map_variables variable_1="V" variable_2="V"/>
      <map_variables variable_1="R" variable_2="R"/>
      <map_variables variable_1="T" variable_2="T"/>
      <map_variables variable_1="F" variable_2="F"/>
      <map_variables variable_1="Cao" variable_2="Cao"/>
      <map_variables variable_1="Nao" variable_2="Nao"/>
      <map_variables variable_1="Vmyo" variable_2="Vmyo"/>
      <map_variables variable_1="Acap" variable_2="Acap"/>
      <map_variables variable_1="Cm" variable_2="Cm"/>
   </connection>
   <connection>
      <map_components component_1="membrane" component_2="calcium_background_current"/>
      <map_variables variable_1="i_Cab" variable_2="i_Cab"/>
      <map_variables variable_1="V" variable_2="V"/>
      <map_variables variable_1="R" variable_2="R"/>
      <map_variables variable_1="T" variable_2="T"/>
      <map_variables variable_1="F" variable_2="F"/>
      <map_variables variable_1="Cao" variable_2="Cao"/>
      <map_variables variable_1="Vmyo" variable_2="Vmyo"/>
      <map_variables variable_1="Acap" variable_2="Acap"/>
      <map_variables variable_1="Cm" variable_2="Cm"/>
   </connection>
   <connection>
      <map_components component_1="membrane" component_2="fast_sodium_current"/>
      <map_variables variable_1="i_Na" variable_2="i_Na"/>
      <map_variables variable_1="V" variable_2="V"/>
      <map_variables variable_1="R" variable_2="R"/>
      <map_variables variable_1="T" variable_2="T"/>
      <map_variables variable_1="F" variable_2="F"/>
      <map_variables variable_1="Ko" variable_2="Ko"/>
      <map_variables variable_1="Nao" variable_2="Nao"/>
   </connection>
   <connection>
      <map_components component_1="membrane" component_2="sodium_background_current"/>
      <map_variables variable_1="i_Nab" variable_2="i_Nab"/>
      <map_variables variable_1="V" variable_2="V"/>
   </connection>
   <connection>
      <map_components component_1="membrane" component_2="sodium_potassium_pump_current"/>
      <map_variables variable_1="i_NaK" variable_2="i_NaK"/>
      <map_variables variable_1="V" variable_2="V"/>
      <map_variables variable_1="R" variable_2="R"/>
      <map_variables variable_1="T" variable_2="T"/>
      <map_variables variable_1="F" variable_2="F"/>
      <map_variables variable_1="Ko" variable_2="Ko"/>
      <map_variables variable_1="Nao" variable_2="Nao"/>
   </connection>
   <connection>
      <map_components component_1="membrane" component_2="fast_transient_outward_K_I"/>
      <map_variables variable_1="i_Kto_f" variable_2="i_Kto_f"/>
      <map_variables variable_1="V" variable_2="V"/>
      <map_variables variable_1="R" variable_2="R"/>
      <map_variables variable_1="T" variable_2="T"/>
      <map_variables variable_1="F" variable_2="F"/>
      <map_variables variable_1="Ko" variable_2="Ko"/>
   </connection>
   <connection>
      <map_components component_1="membrane" component_2="slow_transient_outward_K_I"/>
      <map_variables variable_1="i_Kto_s" variable_2="i_Kto_s"/>
      <map_variables variable_1="V" variable_2="V"/>
   </connection>
   <connection>
      <map_components component_1="membrane" component_2="time_independent_K_I"/>
      <map_variables variable_1="i_K1" variable_2="i_K1"/>
      <map_variables variable_1="V" variable_2="V"/>
      <map_variables variable_1="Ko" variable_2="Ko"/>
   </connection>
   <connection>
      <map_components component_1="membrane" component_2="slow_delayed_rectifier_K_I"/>
      <map_variables variable_1="i_Ks" variable_2="i_Ks"/>
      <map_variables variable_1="V" variable_2="V"/>
   </connection>
   <connection>
      <map_components component_1="membrane" component_2="ultra_rapidly_activating_delayed_rectifier_K_I"/>
      <map_variables variable_1="i_Kur" variable_2="i_Kur"/>
      <map_variables variable_1="V" variable_2="V"/>
   </connection>
   <connection>
      <map_components component_1="membrane" component_2="non_inactivating_steady_state_K_I"/>
      <map_variables variable_1="i_Kss" variable_2="i_Kss"/>
      <map_variables variable_1="V" variable_2="V"/>
   </connection>
   <connection>
      <map_components component_1="membrane" component_2="rapid_delayed_rectifier_K_I"/>
      <map_variables variable_1="i_Kr" variable_2="i_Kr"/>
      <map_variables variable_1="V" variable_2="V"/>
      <map_variables variable_1="R" variable_2="R"/>
      <map_variables variable_1="T" variable_2="T"/>
      <map_variables variable_1="F" variable_2="F"/>
      <map_variables variable_1="Ko" variable_2="Ko"/>
      <map_variables variable_1="Nao" variable_2="Nao"/>
   </connection>
   <connection>
      <map_components component_1="membrane" component_2="calcium_activated_chloride_current"/>
      <map_variables variable_1="i_ClCa" variable_2="i_ClCa"/>
      <map_variables variable_1="V" variable_2="V"/>
   </connection>
   <connection>
      <map_components component_1="membrane" component_2="calcium_concentration"/>
      <map_variables variable_1="Cm" variable_2="Cm"/>
      <map_variables variable_1="V" variable_2="V"/>
      <map_variables variable_1="Acap" variable_2="Acap"/>
      <map_variables variable_1="Vmyo" variable_2="Vmyo"/>
      <map_variables variable_1="VNSR" variable_2="VNSR"/>
      <map_variables variable_1="VJSR" variable_2="VJSR"/>
      <map_variables variable_1="Vss" variable_2="Vss"/>
      <map_variables variable_1="F" variable_2="F"/>
   </connection>
   <connection>
      <map_components component_1="membrane" component_2="calcium_fluxes"/>
      <map_variables variable_1="V" variable_2="V"/>
   </connection>
   <connection>
      <map_components component_1="membrane" component_2="sodium_concentration"/>
      <map_variables variable_1="Cm" variable_2="Cm"/>
      <map_variables variable_1="Acap" variable_2="Acap"/>
      <map_variables variable_1="Vmyo" variable_2="Vmyo"/>
      <map_variables variable_1="F" variable_2="F"/>
   </connection>
   <connection>
      <map_components component_1="membrane" component_2="potassium_concentration"/>
      <map_variables variable_1="Cm" variable_2="Cm"/>
      <map_variables variable_1="Acap" variable_2="Acap"/>
      <map_variables variable_1="Vmyo" variable_2="Vmyo"/>
      <map_variables variable_1="F" variable_2="F"/>
      <map_variables variable_1="i_Stim" variable_2="i_Stim"/>
   </connection>
   <connection>
      <map_components component_1="calcium_concentration" component_2="calcium_fluxes"/>
      <map_variables variable_1="Cass" variable_2="Cass"/>
      <map_variables variable_1="CaJSR" variable_2="CaJSR"/>
      <map_variables variable_1="Cai" variable_2="Cai"/>
      <map_variables variable_1="CaNSR" variable_2="CaNSR"/>
      <map_variables variable_1="J_rel" variable_2="J_rel"/>
      <map_variables variable_1="J_tr" variable_2="J_tr"/>
      <map_variables variable_1="J_xfer" variable_2="J_xfer"/>
      <map_variables variable_1="J_leak" variable_2="J_leak"/>
      <map_variables variable_1="J_serca" variable_2="J_serca"/>
   </connection>
   <connection>
      <map_components component_1="calcium_concentration" component_2="ryanodine_receptors"/>
      <map_variables variable_1="Cass" variable_2="Cass"/>
   </connection>
   <connection>
      <map_components component_1="calcium_concentration" component_2="L_type_calcium_current"/>
      <map_variables variable_1="i_CaL" variable_2="i_CaL"/>
      <map_variables variable_1="Cass" variable_2="Cass"/>
   </connection>
   <connection>
      <map_components component_1="calcium_concentration" component_2="calcium_pump_current"/>
      <map_variables variable_1="i_pCa" variable_2="i_pCa"/>
      <map_variables variable_1="Cai" variable_2="Cai"/>
   </connection>
   <connection>
      <map_components component_1="calcium_concentration" component_2="sodium_calcium_exchange_current"/>
      <map_variables variable_1="i_NaCa" variable_2="i_NaCa"/>
      <map_variables variable_1="Cai" variable_2="Cai"/>
   </connection>
   <connection>
      <map_components component_1="calcium_concentration" component_2="calcium_background_current"/>
      <map_variables variable_1="i_Cab" variable_2="i_Cab"/>
      <map_variables variable_1="Cai" variable_2="Cai"/>
   </connection>
   <connection>
      <map_components component_1="calcium_concentration" component_2="calcium_activated_chloride_current"/>
      <map_variables variable_1="Cai" variable_2="Cai"/>
   </connection>
   <connection>
      <map_components component_1="calcium_fluxes" component_2="L_type_calcium_current"/>
      <map_variables variable_1="i_CaL_max" variable_2="i_CaL_max"/>
      <map_variables variable_1="i_CaL" variable_2="i_CaL"/>
   </connection>
   <connection>
      <map_components component_1="calcium_fluxes" component_2="ryanodine_receptors"/>
      <map_variables variable_1="P_O1" variable_2="P_O1"/>
      <map_variables variable_1="P_O2" variable_2="P_O2"/>
   </connection>
   <connection>
      <map_components component_1="sodium_concentration" component_2="fast_sodium_current"/>
      <map_variables variable_1="Nai" variable_2="Nai"/>
      <map_variables variable_1="i_Na" variable_2="i_Na"/>
   </connection>
   <connection>
      <map_components component_1="sodium_concentration" component_2="sodium_background_current"/>
      <map_variables variable_1="i_Nab" variable_2="i_Nab"/>
   </connection>
   <connection>
      <map_components component_1="sodium_concentration" component_2="sodium_calcium_exchange_current"/>
      <map_variables variable_1="Nai" variable_2="Nai"/>
      <map_variables variable_1="i_NaCa" variable_2="i_NaCa"/>
   </connection>
   <connection>
      <map_components component_1="sodium_concentration" component_2="sodium_potassium_pump_current"/>
      <map_variables variable_1="Nai" variable_2="Nai"/>
      <map_variables variable_1="i_NaK" variable_2="i_NaK"/>
   </connection>
   <connection>
      <map_components component_1="fast_sodium_current" component_2="sodium_background_current"/>
      <map_variables variable_1="E_Na" variable_2="E_Na"/>
   </connection>
   <connection>
      <map_components component_1="fast_sodium_current" component_2="potassium_concentration"/>
      <map_variables variable_1="Ki" variable_2="Ki"/>
   </connection>
   <connection>
      <map_components component_1="potassium_concentration" component_2="fast_transient_outward_K_I"/>
      <map_variables variable_1="Ki" variable_2="Ki"/>
      <map_variables variable_1="i_Kto_f" variable_2="i_Kto_f"/>
   </connection>
   <connection>
      <map_components component_1="potassium_concentration" component_2="slow_transient_outward_K_I"/>
      <map_variables variable_1="i_Kto_s" variable_2="i_Kto_s"/>
   </connection>
   <connection>
      <map_components component_1="potassium_concentration" component_2="time_independent_K_I"/>
      <map_variables variable_1="i_K1" variable_2="i_K1"/>
   </connection>
   <connection>
      <map_components component_1="potassium_concentration" component_2="slow_delayed_rectifier_K_I"/>
      <map_variables variable_1="i_Ks" variable_2="i_Ks"/>
   </connection>
   <connection>
      <map_components component_1="potassium_concentration" component_2="rapid_delayed_rectifier_K_I"/>
      <map_variables variable_1="i_Kr" variable_2="i_Kr"/>
      <map_variables variable_1="Ki" variable_2="Ki"/>
   </connection>
   <connection>
      <map_components component_1="sodium_concentration" component_2="rapid_delayed_rectifier_K_I"/>
      <map_variables variable_1="Nai" variable_2="Nai"/>
   </connection>
   <connection>
      <map_components component_1="potassium_concentration" component_2="ultra_rapidly_activating_delayed_rectifier_K_I"/>
      <map_variables variable_1="i_Kur" variable_2="i_Kur"/>
   </connection>
   <connection>
      <map_components component_1="potassium_concentration" component_2="sodium_potassium_pump_current"/>
      <map_variables variable_1="i_NaK" variable_2="i_NaK"/>
   </connection>
   <connection>
      <map_components component_1="potassium_concentration" component_2="non_inactivating_steady_state_K_I"/>
      <map_variables variable_1="i_Kss" variable_2="i_Kss"/>
   </connection>
   <connection>
      <map_components component_1="fast_transient_outward_K_I" component_2="slow_transient_outward_K_I"/>
      <map_variables variable_1="E_K" variable_2="E_K"/>
   </connection>
   <connection>
      <map_components component_1="fast_transient_outward_K_I" component_2="time_independent_K_I"/>
      <map_variables variable_1="E_K" variable_2="E_K"/>
   </connection>
   <connection>
      <map_components component_1="fast_transient_outward_K_I" component_2="slow_delayed_rectifier_K_I"/>
      <map_variables variable_1="E_K" variable_2="E_K"/>
   </connection>
   <connection>
      <map_components component_1="fast_transient_outward_K_I" component_2="ultra_rapidly_activating_delayed_rectifier_K_I"/>
      <map_variables variable_1="E_K" variable_2="E_K"/>
   </connection>
   <connection>
      <map_components component_1="fast_transient_outward_K_I" component_2="non_inactivating_steady_state_K_I"/>
      <map_variables variable_1="E_K" variable_2="E_K"/>
   </connection>
   <connection>
      <map_components component_1="slow_transient_outward_K_I" component_2="ultra_rapidly_activating_delayed_rectifier_K_I"/>
      <map_variables variable_1="ass" variable_2="ass"/>
      <map_variables variable_1="iss" variable_2="iss"/>
   </connection>
   <connection>
      <map_components component_1="slow_transient_outward_K_I" component_2="non_inactivating_steady_state_K_I"/>
      <map_variables variable_1="ass" variable_2="ass"/>
   </connection>
<rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#"><rdf:Description rdf:about="#MyModel_2006"><simulation xmlns="http://www.cellml.org/metadata/simulation/1.0#"><rdf:Description rdf:nodeID="n1"/></simulation></rdf:Description><rdf:Description rdf:nodeID="n1"><boundIntervals xmlns="http://www.cellml.org/metadata/simulation/1.0#"><rdf:Description rdf:nodeID="n2"/></boundIntervals></rdf:Description><rdf:Description rdf:nodeID="n2"><rdf:first><rdf:Description rdf:nodeID="n3"/></rdf:first><rdf:rest><rdf:Description rdf:about="http://www.w3.org/1999/02/22-rdf-syntax-ns#nil"/></rdf:rest></rdf:Description><rdf:Description rdf:nodeID="n3"><pointDensity xmlns="http://www.cellml.org/metadata/simulation/1.0#nonstandard-" xml:lang="en">10000</pointDensity><maximumStepSize xmlns="http://www.cellml.org/metadata/simulation/1.0#" xml:lang="en">0.01</maximumStepSize><endingValue xmlns="http://www.cellml.org/metadata/simulation/1.0#" xml:lang="en">1000</endingValue></rdf:Description></rdf:RDF></model>