Rendering of the source text

<?xml version='1.0' encoding='utf-8'?>
<!--  

AUTHOR :  John Davidson
          Bioengineering Institute
          The University of Auckland
          

--><model xmlns="http://www.cellml.org/cellml/1.0#" xmlns:cmeta="http://www.cellml.org/metadata/1.0#" xmlns:dc="http://purl.org/dc/elements/1.1/" xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#" xmlns:bqs="http://www.cellml.org/bqs/1.0#" xmlns:cellml="http://www.cellml.org/cellml/1.0#" xmlns:dcterms="http://purl.org/dc/terms/" xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" cmeta:id="Slow_Twitch" name="Slow_Twitch">

<documentation xmlns="http://cellml.org/tmp-documentation">
<article>
  <articleinfo>
  <title>A mathematical model of fatigue in skeletal muscle force contraction</title>
  <author>
    <firstname>Catherine</firstname>
          <surname>Lloyd</surname>
    <affiliation>
      <shortaffil>Auckland Bioengineering Institute, The University of Auckland</shortaffil>
    </affiliation>
  </author>
</articleinfo>
  <section id="sec_status">
    <title>Model Status</title>
    <para>
        This cellML model is known to run in OpenCell, JSim and COR to recreate the published results.  The units have been checked in COR and they are consistent in as much as they are known to be equivalent.  This particular version of the CellML model represents slow twitch muscle.
          </para>
  </section>

<sect1 id="sec_structure">
  <title>Model Structure</title>

<para>
Abstract: The ability for muscle to repeatedly generate force is limited by fatigue. The cellular mechanisms behind
muscle fatigue are complex and potentially include breakdown at many points along the excitation–contraction pathway. In this paper we construct a mathematical model of the skeletal muscle excitation–contraction pathway based on the cellular biochemical events that link excitation to contraction. The model includes descriptions of membrane voltage, calcium cycling and crossbridge dynamics and was parameterised and validated using the response characteristics of mouse skeletal muscle to a range of electrical stimuli. This model was used to uncover the complexities of skeletal muscle fatigue. We also parameterised our model to describe force kinetics in fast and slow twitch fibre types, which have a number of biochemical and biophysical differences. How these differences interact to generate different force/fatigue responses in fast- and slow- twitch fibres is not well understood and we used our modelling approach to bring new insights to this relationship.
</para>

<informalfigure float="0">
<mediaobject>
  <imageobject>
    <objectinfo>
      <title>model diagram</title>
    </objectinfo>
    <imagedata fileref="shorten_2007a.png"/>
  </imageobject>
</mediaobject>
<caption>Schematic diagram of the muscle membrane electrophysiological model. There are descriptions of K+, Cl- and Na+ transport in the sarcolemma and t-tubules membranes. The model includes descriptions of an inward rectifier K+ current (IIR), a delayed rectifier K+ current (IDR), a Na+ current (INa), a Na+/K+ pump (INaK) and a Cl- current (ICl). Currents in the t-tubules are denoted by a subscript t.</caption>
</informalfigure>

<informalfigure float="0">
<mediaobject>
  <imageobject>
    <objectinfo>
      <title>model diagram</title>
    </objectinfo>
    <imagedata fileref="shorten_2007b.png"/>
  </imageobject>
</mediaobject>
<caption>Schematic diagram of the Ca2+ dynamics model. Ca2+ binds to parvalbumin and ATP myoplasmic buffers, is transported and buffered to calsequestrin in the SR and bind to tropinin to generate force.</caption>
</informalfigure>

<informalfigure float="0">
<mediaobject>
  <imageobject>
    <objectinfo>
      <title>model diagram</title>
    </objectinfo>
    <imagedata fileref="shorten_2007d.png"/>
  </imageobject>
</mediaobject>
<caption>Schematic diagram of the crossbridge dynamics model. Two Ca2+ ions must bind to the troponin-tropomyosin regulatory unit (RU) to permit crossbridge cycling (XB) kinetics. Open circles denote Ca2+ binding sites, the black circles represent the number of Ca2+ ions bound to troponin and the black bar represents the tropomyosin block.</caption>
</informalfigure>

<informalfigure float="0">
<mediaobject>
  <imageobject>
    <objectinfo>
      <title>model diagram</title>
    </objectinfo>
    <imagedata fileref="shorten_2007c.png"/>
  </imageobject>
</mediaobject>
<caption>Schematic diagram of the phosphate dynamics model. Phosphate is generated during crossbridge cycling and SR Ca2+ pumping, feeds back to slow crossbridge cycling and is transported into the SR to precipitate with Ca2+ and reduce Ca2+ release from the SR.</caption>
</informalfigure>

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

<para>
A mathematical model of fatigue in skeletal muscle force contraction, Paul R. Shorten, Paul O’Callaghan, John B. Davidson and Tanya K. Soboleva, 2007, <emphasis>Journal of Muscle Research and Cell Motility</emphasis>, 28, 293-313. <ulink url="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=pubmed&amp;cmd=Retrieve&amp;dopt=AbstractPlus&amp;list_uids=18080210&amp;query_hl=1&amp;itool=pubmed_docsum">PubMed ID: 18080210</ulink>
</para>

</sect1>
</article>
</documentation> 
    
<!-- units definition


-->    
    
<units name="millisecond">
  <unit units="second" prefix="milli"/>
</units>

<units name="per_millisecond">
  <unit units="second" prefix="milli" exponent="-1"/>
</units>

<units name="millivolt">
  <unit units="volt" prefix="milli"/>
</units>

<units name="milliS_per_cm2">
  <unit units="siemens" prefix="milli"/>
  <unit units="metre" prefix="centi" exponent="-2"/>
</units>

<units name="microF_per_cm2">
  <unit units="farad" prefix="micro"/>
  <unit units="metre" prefix="centi" exponent="-2"/>
</units>

<units name="microA_per_cm2">
  <unit units="ampere" prefix="micro"/>
  <unit units="metre" prefix="centi" exponent="-2"/>
</units>

<units name="ohm_cm2">
  <unit units="ohm"/>
  <unit units="metre" prefix="centi" exponent="2"/>
</units>

<units name="centi_metre">
  <unit units="metre" prefix="centi"/>
</units>

<units name="C_per_mol">
  <unit units="coulomb"/>
  <unit units="mole" exponent="-1"/>
</units>

<units name="per_millisecond_per_millivolt">
  <unit units="second" prefix="milli" exponent="-1"/>
  <unit units="volt" prefix="milli" exponent="-1"/>
</units>

<units name="milliJ_per_degreeK_per_mol">
  <unit units="joule" prefix="milli"/>
  <unit units="kelvin" exponent="-1"/>
  <unit units="mole" exponent="-1"/>
</units>

<units name="degreeK">
  <unit units="kelvin"/>
</units>

<units name="milliM">
  <unit units="mole" prefix="milli"/>
  <unit units="litre" exponent="-1"/>
</units>
<units name="M">
  <unit units="mole"/>
  <unit units="litre" exponent="-1"/>
</units>

<!--<units name="milliM2">
  <unit prefix="micro" units="mole" exponent="2"/>
  <unit units="litre" exponent="-2" />
</units>-->

<units name="micro_mol_per_cm2_per_second">
  <unit units="mole" prefix="micro"/>
  <unit units="metre" prefix="centi" exponent="-2"/>
  <unit units="second" exponent="-1"/>
</units>

<!--
Stern and Rios units
-->

<units name="micrometre3_per_millisecond">
  <unit units="metre" prefix="micro" exponent="3"/>
  <unit units="second" prefix="milli" exponent="-1"/>
</units>

<units name="per_micromolar2_per_millisecond">
  <unit units="mole" prefix="micro" exponent="-2"/>
  <unit units="litre" exponent="2"/>
  <unit units="second" prefix="milli" exponent="-1"/>
</units>

<units name="per_micromolar_per_millisecond">
  <unit units="mole" prefix="micro" exponent="-1"/>
  <unit units="litre"/>
  <unit units="second" prefix="milli" exponent="-1"/>
</units>

<units name="micromolar">
  <unit units="mole" prefix="micro"/>
  <unit units="litre" exponent="-1"/>
</units>

<!-- 
Razumova units
-->

<units name="micromolar_per_millisecond_micrometre3">
  <unit units="mole" prefix="micro"/>
  <unit units="litre" exponent="-1"/>
  <unit units="second" prefix="milli" exponent="-1"/>
  <unit units="metre" prefix="micro" exponent="3"/>
</units>

<units name="micrometre">
  <unit units="metre" prefix="micro"/>
</units>

<units name="micrometre3">
  <unit units="metre" prefix="micro" exponent="3"/>
</units>

<units name="milliV_milliM">
  <unit units="volt" prefix="milli"/>
  <unit units="mole" prefix="milli"/>
  <unit units="litre" exponent="-1"/>
</units>

<units name="milliM2">  
  <unit units="mole" prefix="milli" exponent="2"/>
  <unit units="litre" exponent="-2"/>
</units>

<units name="per_milliM3_per_millisecond">  
  <unit units="mole" prefix="milli" exponent="-3"/>
  <unit units="litre" exponent="3"/>
  <unit units="second" prefix="milli" exponent="-1"/>
</units>

<units name="per_milliM2_per_millisecond">  
  <unit units="mole" prefix="milli" exponent="-2"/>
  <unit units="litre" exponent="2"/>
  <unit units="second" prefix="milli" exponent="-1"/>
</units>

<!--  components
-->

<component name="cell">

  <variable units="millisecond" public_interface="out" name="time"/>
  
</component>



<component name="wal_environment">

  <variable units="microF_per_cm2" name="C_m" initial_value="0.58"/>
  <variable units="dimensionless" name="gam" initial_value="2.79"/>
  <variable units="ohm_cm2" name="R_a" initial_value="150"/>
  <variable units="centi_metre" name="tsi" initial_value="0.000001"/>
  <variable units="centi_metre" name="tsi2" initial_value="0.0025"/>
  <variable units="centi_metre" name="tsi3" initial_value="0.0005"/>
  <variable units="C_per_mol" private_interface="out" name="FF" initial_value="96485"/>
  <variable units="millisecond" name="tau_K" initial_value="559"/>
  <variable units="millisecond" name="tau_Na" initial_value="559"/>
  <variable units="dimensionless" name="f_T" initial_value="0.00174"/>
  <variable units="millisecond" name="tau_K2" initial_value="40229.885"/>
  <variable units="millisecond" name="tau_Na2" initial_value="40229.885"/>
  <variable units="microA_per_cm2" name="I_K_rest" initial_value="0.34"/>
  <variable units="microA_per_cm2" name="I_Na_rest" initial_value="-0.43"/>
  <variable units="per_millisecond" private_interface="out" name="alpha_h_bar" initial_value="0.0081"/>
  <variable units="per_millisecond_per_millivolt" private_interface="out" name="alpha_m_bar" initial_value="0.288"/>
  <variable units="per_millisecond_per_millivolt" private_interface="out" name="alpha_n_bar" initial_value="0.0131"/>
  <variable units="per_millisecond" private_interface="out" name="beta_h_bar" initial_value="4.38"/>
  <variable units="per_millisecond" private_interface="out" name="beta_m_bar" initial_value="1.38"/>
  <variable units="per_millisecond" private_interface="out" name="beta_n_bar" initial_value="0.067"/>
  <variable units="millivolt" private_interface="out" name="V_m" initial_value="-46"/>
  <variable units="millivolt" private_interface="out" name="V_n" initial_value="-40"/>
  <variable units="millivolt" private_interface="out" name="V_h" initial_value="-45"/>
  <variable units="millivolt" private_interface="out" name="V_a" initial_value="70"/>
  <variable units="millivolt" private_interface="out" name="V_S_inf" initial_value="-68"/>
  <variable units="millivolt" private_interface="out" name="V_h_K_inf" initial_value="-40"/>
  <variable units="millivolt" private_interface="out" name="A_a" initial_value="150"/>
  <variable units="millivolt" private_interface="out" name="A_S_inf" initial_value="7.1"/>
  <variable units="millivolt" private_interface="out" name="A_h_K_inf" initial_value="7.5"/>
  <variable units="millivolt" private_interface="out" name="K_alpha_h" initial_value="14.7"/>
  <variable units="millivolt" private_interface="out" name="K_beta_h" initial_value="9"/>
  <variable units="millivolt" private_interface="out" name="K_alpha_m" initial_value="10"/>
  <variable units="millivolt" private_interface="out" name="K_alpha_n" initial_value="7"/> 
  <variable units="millivolt" private_interface="out" name="K_beta_m" initial_value="18"/>
  <variable units="millivolt" private_interface="out" name="K_beta_n" initial_value="40"/>
  <variable units="milliJ_per_degreeK_per_mol" private_interface="out" name="RR" initial_value="8314.41"/>
  <variable units="degreeK" private_interface="out" name="TT" initial_value="293"/>
  <variable units="milliS_per_cm2" private_interface="out" name="g_Cl_bar" initial_value="3.275"/>
  <variable units="milliS_per_cm2" private_interface="out" name="g_K_bar" initial_value="10.8"/>
  <variable units="milliS_per_cm2" private_interface="out" name="g_Na_bar" initial_value="134"/>
  <variable units="milliS_per_cm2" private_interface="out" name="G_K" initial_value="1.85"/>
  <variable units="dimensionless" private_interface="out" name="del" initial_value="0.4"/>
  <variable units="milliM2" private_interface="out" name="K_K" initial_value="950"/>
  <variable units="milliM2" private_interface="out" name="K_S" initial_value="1"/>
  <variable units="milliM" private_interface="out" name="K_m_K" initial_value="1"/>
  <variable units="milliM" private_interface="out" name="K_m_Na" initial_value="13"/>
  <variable units="milliM" private_interface="out" name="S_i" initial_value="10"/>
  <variable units="micro_mol_per_cm2_per_second" private_interface="out" name="J_NaK_bar" initial_value="0.0001656"/>
  <variable units="millivolt" private_interface="out" name="V_tau" initial_value="70"/>
  
  
  
  <variable units="microA_per_cm2" private_interface="out" name="I_T"/>
  <variable units="millivolt" public_interface="out" private_interface="out" name="vS" initial_value="-79.974"/>
  <variable units="millivolt" public_interface="out" private_interface="out" name="vT" initial_value="-80.2"/>
  <variable units="microA_per_cm2" name="I_ionic_s"/>
  <variable units="microA_per_cm2" name="I_ionic_t"/>
  <variable units="milliM" private_interface="out" name="K_t" initial_value="5.9"/>
  <variable units="milliM" private_interface="out" name="K_i" initial_value="150.9"/>
  <variable units="milliM" private_interface="out" name="K_e" initial_value="5.9"/>
  <variable units="milliM" private_interface="out" name="Na_i" initial_value="12.7"/>
  <variable units="milliM" private_interface="out" name="Na_t" initial_value="133"/>
  <variable units="milliM" private_interface="out" name="Na_e" initial_value="133"/>
  <variable units="millivolt" private_interface="out" name="E_K"/>
  <variable units="millivolt" private_interface="out" name="E_K_t"/>
  <variable units="milliM" private_interface="out" name="Cl_i"/>
  <variable units="milliM" private_interface="out" name="Cl_o"/>
  <variable units="milliM" private_interface="out" name="Cl_i_t"/>
  <variable units="milliM" private_interface="out" name="Cl_o_t"/>
  <variable units="milliV_milliM" private_interface="out" name="J_K"/>
  <variable units="milliV_milliM" private_interface="out" name="J_K_t"/>
  <variable units="dimensionless" private_interface="out" name="eta_Cl" initial_value="0.1"/>
  <variable units="dimensionless" private_interface="out" name="eta_IR" initial_value="1.0"/>
  <variable units="dimensionless" private_interface="out" name="eta_DR" initial_value="0.45"/>
  <variable units="dimensionless" private_interface="out" name="eta_Na" initial_value="0.1"/>
  <variable units="dimensionless" private_interface="out" name="eta_NaK" initial_value="0.1"/>
  
  
  
  
   
  <variable units="microA_per_cm2" private_interface="in" name="I_Cl"/>
  <variable units="microA_per_cm2" private_interface="in" name="I_IR"/>
  <variable units="microA_per_cm2" private_interface="in" name="I_DR"/>
  <variable units="microA_per_cm2" private_interface="in" name="I_Na"/>
  <variable units="microA_per_cm2" private_interface="in" name="I_NaK"/>
  <variable units="microA_per_cm2" private_interface="in" name="I_Cl_t"/>
  <variable units="microA_per_cm2" private_interface="in" name="I_IR_t"/>
  <variable units="microA_per_cm2" private_interface="in" name="I_DR_t"/>
  <variable units="microA_per_cm2" private_interface="in" name="I_Na_t"/>
  <variable units="microA_per_cm2" private_interface="in" name="I_NaK_t"/>
  
  <variable units="microA_per_cm2" name="I_HH"/>
  <variable units="millisecond" public_interface="in" private_interface="out" name="time"/>
  
  
<math xmlns="http://www.w3.org/1998/Math/MathML"> 
<apply id="I_T_Eqn">
        <eq/>

  <ci> I_T </ci>
  
  <apply>
          <times/> 
  
    <apply>
            <divide/> <!-- unit fix 5_08 -->
      <cn cellml:units="millivolt"> 1000 </cn>
      <cn cellml:units="volt"> 1 </cn>
    </apply>
    <apply>
            <divide/>
  
      <apply>
              <minus/>
    
        <ci> vS </ci>
      
        <ci> vT </ci>
      
      </apply>
    
      <ci> R_a </ci>
      
    </apply>
     
  </apply>
 
</apply>
 
</math> 


<math xmlns="http://www.w3.org/1998/Math/MathML"> 
<apply id="vS_diff_calculation">
        <eq/>
  
  <apply>
          <diff/>
    
    <bvar>
            <ci> time </ci>
          </bvar>
 
    <ci> vS </ci>
    
  </apply> 
  
  <apply>
          <minus/> 
    
    <apply>
            <divide/> 
      
      <apply>
              <plus/> 
	
        <ci> I_ionic_s </ci>
	  
        <ci> I_T </ci>
	  
      </apply>
	
      <ci> C_m </ci>
	
    </apply>
      
  </apply>
    
</apply> 
  
</math> 


<math xmlns="http://www.w3.org/1998/Math/MathML"> 
<apply id="dvT_Eqn">
        <eq/>

  <apply>
          <diff/>
    
    <bvar>
            <ci> time </ci>
          </bvar>
 
    <ci> vT </ci>
    
  </apply> 
 
  <apply>
          <minus/>
   
    <apply>
            <divide/>
    
      <apply>
              <minus/>
    
        <ci> I_ionic_t </ci>

        <apply>
                <divide/>
  
          <ci> I_T </ci>
 
          <ci> gam </ci>

        </apply>
  
      </apply>

      <ci> C_m </ci>

    </apply>
  
  </apply>

</apply> 
 
</math> 


<math xmlns="http://www.w3.org/1998/Math/MathML"> 
<apply id="I_ionic_s_Eqn">
        <eq/> 

  <ci> I_ionic_s </ci> 
  
  <apply>
          <plus/> 
    
    <ci> I_Cl </ci>
 
    <ci> I_IR </ci>
  
    <ci> I_DR </ci>

    <ci> I_Na </ci>

    <ci> I_NaK </ci>
    
    <apply>
            <minus/>
    
    	<ci> I_HH </ci>
	
    </apply>
    
  </apply>
    
</apply> 
 
</math> 
 
 

<math xmlns="http://www.w3.org/1998/Math/MathML"> 
<apply id="I_ionic_t_Eqn">
        <eq/> 

  <ci> I_ionic_t </ci> 
 
  <apply>
          <plus/> 
         
    <ci> I_Cl_t </ci>
 
    <ci> I_IR_t </ci>

    <ci> I_DR_t </ci>

    <ci> I_Na_t </ci>

    <ci> I_NaK_t </ci>

  </apply>
  
</apply>
  
</math> 
 





      
<math xmlns="http://www.w3.org/1998/Math/MathML"> 
<apply id="dK_i_Eqn">
        <eq/>
 
  <apply>
          <diff/>
    
    <bvar>
            <ci> time </ci>
          </bvar>
 
    <ci> K_i </ci>
    
  </apply>
  
  <apply>
          <minus/>
  
  
    <apply>
            <times/>
    
      <apply>
              <minus/>
      
        <ci> f_T </ci>
	
      </apply>
  
      <apply>
              <divide/> 
 
        <apply>
                <plus/>
  
          <ci> I_IR_t </ci>
 
          <ci> I_DR_t </ci>
	  
	  <ci> I_K_rest </ci>

          <apply>
                  <times/> 

            <apply>
                    <minus/> 
 
              <cn cellml:units="dimensionless"> 2 </cn>
  
            </apply>

            <ci> I_NaK_t </ci>

          </apply>
  
        </apply>  
  
        <apply>
                <times/> 
          <apply>
                  <divide/>
            <cn cellml:units="milliM"> 1000 </cn>
            <cn cellml:units="M"> 1 </cn>
          </apply>
          <ci> FF </ci>

          <ci> tsi </ci>

        </apply>

      </apply>
    
    </apply> 
 
    <apply>
            <divide/> 
 
      <apply>
              <plus/>
  
        <ci> I_IR </ci>
 
        <ci> I_DR </ci>
	
	<ci> I_K_rest </ci>
	
	<apply>
                <times/>
	
	  <cn cellml:units="dimensionless"> -2 </cn>
	  
	  <ci> I_NaK </ci>
	  
	</apply>
 
      </apply>
  
      <apply>
              <times/>
        <apply>
                <divide/>
          <cn cellml:units="milliM"> 1000 </cn>
	  <cn cellml:units="M"> 1 </cn>
        </apply>
	<ci> FF </ci>
	
	<ci> tsi2 </ci>
	
      </apply>

    </apply>

  </apply>
  
</apply>
 
</math>  
   
  
<math xmlns="http://www.w3.org/1998/Math/MathML"> 
<apply id="dK_t_Eqn">
        <eq/>
 
  <apply>
          <diff/>
    
    <bvar>
            <ci> time </ci>
          </bvar>
 
    <ci> K_t </ci>
    
  </apply>
  
  <apply>
          <minus/>
   
    <apply>
            <divide/> 
 
      <apply>
              <plus/>
  
        <ci> I_IR_t </ci>
 
        <ci> I_DR_t </ci>
        
	<ci> I_K_rest </ci>
        
	<apply>
                <times/> 

          <apply>
                  <minus/> 
 
            <cn cellml:units="dimensionless"> 2 </cn>
  
          </apply>

          <ci> I_NaK_t </ci>

        </apply>
  
      </apply>  
  
      <apply>
              <times/> 
        <apply>
                <divide/>
          <cn cellml:units="milliM"> 1000 </cn>
          <cn cellml:units="M"> 1 </cn>
        </apply>
        <ci> FF </ci>

        <ci> tsi </ci>

      </apply>

 
    </apply> 
 
    <apply>
            <divide/> 
 
      <apply>
              <minus/>
  
        <ci> K_t </ci>
 
        <ci> K_e </ci>
 
      </apply>
  
      <ci> tau_K </ci>

    </apply>

  </apply>
  
</apply>
 
</math>

<math xmlns="http://www.w3.org/1998/Math/MathML"> 
<apply id="dK_e_Eqn">
        <eq/>
 
  <apply>
          <diff/>
    
    <bvar>
            <ci> time </ci>
          </bvar>
 
    <ci> K_e </ci>
    
  </apply>
  
  <apply>
          <plus/>
   
    <apply>
            <divide/> 
 
      <apply>
              <plus/>
  
        <ci> I_IR </ci>
 
        <ci> I_DR </ci>

	<ci> I_K_rest </ci>
        
        <apply>
                <times/> 

          <apply>
                  <minus/> 
 
            <cn cellml:units="dimensionless"> 2 </cn>
  
          </apply>

          <ci> I_NaK </ci>

        </apply>
  
      </apply>  
  
      <apply>
              <times/> 
        <apply>
                <divide/>
          <cn cellml:units="milliM"> 1000 </cn>
          <cn cellml:units="M"> 1 </cn>
        </apply>
        <ci> FF </ci>

        <ci> tsi3 </ci>

      </apply>

 
    </apply> 
 
    <apply>
            <divide/> 
 
      <apply>
              <minus/>
  
        <ci> K_t </ci>
 
        <ci> K_e </ci>
 
      </apply>
  
      <ci> tau_K2 </ci>

    </apply>

  </apply>
  
</apply>
 
</math>
 




  

      
<math xmlns="http://www.w3.org/1998/Math/MathML"> 
<apply id="dNa_i_Eqn">
        <eq/>
 
  <apply>
          <diff/>
    
    <bvar>
            <ci> time </ci>
          </bvar>
 
    <ci> Na_i </ci>
    
  </apply>
  
  <apply>
          <minus/>
  
  
    <apply>
            <times/>
    
      <apply>
              <minus/>
      
        <ci> f_T </ci>
	
      </apply>
  
      <apply>
              <divide/> 
 
        <apply>
                <plus/>
  
          <ci> I_Na_t </ci>
 
	  <ci> I_Na_rest </ci>
        
          <apply>
                  <times/> 

            <cn cellml:units="dimensionless"> 3 </cn>
  
            <ci> I_NaK_t </ci>

          </apply>
  
        </apply>  
  
        <apply>
                <times/> 
          <apply>
                  <divide/>
            <cn cellml:units="milliM"> 1000 </cn>
            <cn cellml:units="M"> 1 </cn>
          </apply>
          <ci> FF </ci>

          <ci> tsi </ci>

        </apply>

      </apply>
    
    </apply> 
 
    <apply>
            <divide/> 
 
      <apply>
              <plus/>
  
        <ci> I_Na </ci>
	
	<ci> I_Na_rest </ci>
        
	<apply>
                <times/>
	
	  <cn cellml:units="dimensionless"> 3 </cn>
	  
	  <ci> I_NaK </ci>
	  
	</apply>
 
      </apply>
  
      <apply>
              <times/>
        <apply>
                <divide/>
          <cn cellml:units="milliM"> 1000 </cn>
	  <cn cellml:units="M"> 1 </cn>
        </apply>
	<ci> FF </ci>
	
	<ci> tsi2 </ci>
	
      </apply>

    </apply>

  </apply>
  
</apply>
 
</math>  
   
  
<math xmlns="http://www.w3.org/1998/Math/MathML"> 
<apply id="dNa_t_Eqn">
        <eq/>
 
  <apply>
          <diff/>
    
    <bvar>
            <ci> time </ci>
          </bvar>
 
    <ci> Na_t </ci>
    
  </apply>
  
  <apply>
          <minus/>
   
    <apply>
            <divide/> 
 
      <apply>
              <plus/>
  
        <ci> I_Na_t </ci>

	<ci> I_Na_rest </ci>
        
        <apply>
                <times/> 

          <cn cellml:units="dimensionless"> 3 </cn>
  
          <ci> I_NaK_t </ci>

        </apply>
  
      </apply>  
  
      <apply>
              <times/> 
        <apply>
                <divide/>
          <cn cellml:units="milliM"> 1000 </cn>
          <cn cellml:units="M"> 1 </cn>
        </apply>
        <ci> FF </ci>

        <ci> tsi </ci>

      </apply>

 
    </apply> 
 
    <apply>
            <divide/> 
 
      <apply>
              <minus/>
  
        <ci> Na_t </ci>
 
        <ci> Na_e </ci>
 
      </apply>
  
      <ci> tau_Na </ci>

    </apply>

  </apply>
  
</apply>
 
</math>

<math xmlns="http://www.w3.org/1998/Math/MathML"> 
<apply id="dNa_e_Eqn">
        <eq/>
 
  <apply>
          <diff/>
    
    <bvar>
            <ci> time </ci>
          </bvar>
 
    <ci> Na_e </ci>
    
  </apply>
  
  <apply>
          <plus/>
   
    <apply>
            <divide/> 
 
      <apply>
              <plus/>
  
        <ci> I_Na </ci>
 
	<ci> I_Na_rest </ci>
        
        <apply>
                <times/> 

          <cn cellml:units="dimensionless"> 3 </cn>
  
          <ci> I_NaK </ci>

        </apply>
  
      </apply>  
  
      <apply>
              <times/> 
        <apply>
                <divide/>
          <cn cellml:units="milliM"> 1000 </cn>
          <cn cellml:units="M"> 1 </cn>
        </apply>
        <ci> FF </ci>

        <ci> tsi3 </ci>

      </apply>

 
    </apply> 
 
    <apply>
            <divide/> 
 
      <apply>
              <minus/>
  
        <ci> Na_t </ci>
 
        <ci> Na_e </ci>
 
      </apply>
  
      <ci> tau_Na2 </ci>

    </apply>

  </apply>
  
</apply>
 
</math>
 

<math xmlns="http://www.w3.org/1998/Math/MathML"> 
<apply id="E_K_eqn">
        <eq/>
  
  <ci> E_K </ci>
  
  <apply>
          <times/>
  
    <apply>
            <divide/>
  
      <apply>
              <times/>
    
        <ci> RR </ci>
      
        <ci> TT </ci>
      
      </apply>
    
      <ci> FF </ci>
      
    </apply>
      
    <apply>
            <ln/>
      
      <apply>
              <divide/>
      
        <ci> K_e </ci>
	
	<ci> K_i </ci>
	
      </apply>
      
    </apply>
    
  </apply>
  
</apply>

</math>


<math xmlns="http://www.w3.org/1998/Math/MathML"> 
<apply id="E_K_t_eqn">
        <eq/>
  
  <ci> E_K_t </ci>
  
  <apply>
          <times/>
  
    <apply>
            <divide/>
  
      <apply>
              <times/>
    
        <ci> RR </ci>
      
        <ci> TT </ci>
      
      </apply>
              
      <ci> FF </ci>
     
    </apply>
      
    <apply>
            <ln/>
      
      <apply>
              <divide/>
	  
	<ci> K_t </ci>
	
	<ci> K_i </ci>
	
      </apply>
      
    </apply>
    
  </apply>
  
</apply>

</math>



<math xmlns="http://www.w3.org/1998/Math/MathML"> 
<apply id="Cl_i_eqn">
        <eq/>
  
  <ci> Cl_i </ci>
  
  <apply>
          <divide/>
  
    <cn cellml:units="milliM"> 156.5 </cn>
    
    <apply>
            <plus/>
      
      <cn cellml:units="dimensionless"> 5 </cn>
      
      <apply>
              <exp/>
      
        <apply>
                <divide/>
	
	  <apply>
                  <times/>
	  
	    <apply>
                    <minus/>
	    
	      <ci> FF </ci>
	      
	    </apply>
	    
	    <ci> E_K </ci>
	    
	  </apply>
	  
	  <apply>
                  <times/>
	  
	    <ci> RR </ci>
	    
	    <ci> TT </ci>
	    
	  </apply>
	  
	</apply>
	
      </apply>
    
    </apply>
  
  </apply>
  
</apply>

</math>


<math xmlns="http://www.w3.org/1998/Math/MathML"> 
<apply id="Cl_o_eqn">
        <eq/>

  <ci> Cl_o </ci>
  
  <apply>
          <minus/>
  
    <cn cellml:units="milliM"> 156.5 </cn>
    
    <apply>
            <times/>
    
      <cn cellml:units="dimensionless"> 5 </cn>
      
      <ci> Cl_i </ci>
    
    </apply>
    
  </apply>
  
</apply>

</math>

<math xmlns="http://www.w3.org/1998/Math/MathML"> 
<apply id="Cl_i_t_eqn">
        <eq/>
  
  <ci> Cl_i_t </ci>
  
  <apply>
          <divide/>
  
    <cn cellml:units="milliM"> 156.5 </cn>
    
    <apply>
            <plus/>
      
      <cn cellml:units="dimensionless"> 5 </cn>
      
      <apply>
              <exp/>
      
        <apply>
                <divide/>
	
	  <apply>
                  <times/>
	  
	    <apply>
                    <minus/>
	    
	      <ci> FF </ci>
	      
	    </apply>
	    
	    <ci> E_K_t </ci>
	    
	  </apply>
	  
	  <apply>
                  <times/>
	  
	    <ci> RR </ci>
	    
	    <ci> TT </ci>
	    
	  </apply>
	  
	</apply>
	
      </apply>
    
    </apply>
  
  </apply>
  
</apply>

</math>


<math xmlns="http://www.w3.org/1998/Math/MathML"> 
<apply id="Cl_o_t_eqn">
        <eq/>

  <ci> Cl_o_t </ci>
  
  <apply>
          <minus/>
  
    <cn cellml:units="milliM"> 156.5 </cn>
    
    <apply>
            <times/>
    
      <cn cellml:units="dimensionless"> 5 </cn>
      
      <ci> Cl_i_t </ci>
    
    </apply>
    
  </apply>
  
</apply>

</math>

<math xmlns="http://www.w3.org/1998/Math/MathML"> 
<apply id="J_K_eqn">
        <eq/>

  <ci> J_K </ci>
  
  <apply>
          <times/>
  
    <ci> vS </ci>
  
    <apply>
            <divide/>
  
      <apply>
              <minus/>
    
        <ci> K_i </ci>
      
        <apply>
                <times/>
      
          <ci> K_e </ci>
	
	  <apply>
                  <exp/>
	
	    <apply>
                    <divide/>
	  
	      <apply>
                      <times/>
	    
	        <cn cellml:units="dimensionless"> -1 </cn>
	      
	        <ci> FF </ci>
	      
	        <ci> vS </ci>
	      
	      </apply>
	    
	      <apply>
                      <times/>
	    
	        <ci> RR </ci>
	      
	        <ci> TT </ci>
	      
	      </apply>
	    
	    </apply>
	  
	  </apply>
	
        </apply>
      
      </apply>
    
      <apply>
              <minus/>
    
        <cn cellml:units="dimensionless"> 1 </cn>
      
        <apply>
                <exp/>
      
          <apply>
                  <divide/>
	  
	    <apply>
                    <times/>
	    
	      <cn cellml:units="dimensionless"> -1 </cn>
	      
	      <ci> FF </ci>
	      
	      <ci> vS </ci>
	      
	    </apply>
	    
	    <apply>
                    <times/>
	    
	      <ci> RR </ci>
	      
	      <ci> TT </ci>
	      
	    </apply>
	    
	  </apply>
	
        </apply>
      
      </apply>
    
    </apply>
  
  </apply>
  
</apply>

</math>
	    
	      

	    
<math xmlns="http://www.w3.org/1998/Math/MathML"> 
<apply id="J_K_t_eqn">
        <eq/>

  <ci> J_K_t </ci>
  
  <apply>
          <times/>
  
    <ci> vT </ci>
  
    <apply>
            <divide/>
  
      <apply>
              <minus/>
    
        <ci> K_i </ci>
      
        <apply>
                <times/>
      
          <ci> K_t </ci>
	
	  <apply>
                  <exp/>
	
	    <apply>
                    <divide/>
	  
	      <apply>
                      <times/>
	    
	        <cn cellml:units="dimensionless"> -1 </cn>
	      
	        <ci> FF </ci>
	      
	        <ci> vT </ci>
	      
	      </apply>
	    
	      <apply>
                      <times/>
	    
	        <ci> RR </ci>
	      
	        <ci> TT </ci>
	      
	      </apply>
	    
	    </apply>
	  
	  </apply>
	
        </apply>
      
      </apply>
    
      <apply>
              <minus/>
    
        <cn cellml:units="dimensionless"> 1 </cn>
      
        <apply>
                <exp/>
      
          <apply>
                  <divide/>
	  
	    <apply>
                    <times/>
	    
	      <cn cellml:units="dimensionless"> -1 </cn>
	      
	      <ci> FF </ci>
	      
	      <ci> vT </ci>
	      
	    </apply>
	    
	    <apply>
                    <times/>
	    
	      <ci> RR </ci>
	      
	      <ci> TT </ci>
	      
	    </apply>
	    
	  </apply>
	
        </apply>
      
      </apply>
    
    </apply>
  
  </apply>
  
</apply>

</math>
	      



<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply id="I_HH_Eqn">
        <eq/> 

  <ci> I_HH </ci>
  
  <piecewise>
  
    <piece>
    
      <cn cellml:units="microA_per_cm2"> 150 </cn>
	    
      <apply>
              <and/>
    
	<apply>
                <geq/>
	
	  <ci> time </ci>
	  
	  <cn cellml:units="millisecond"> 0 </cn>
	  
	</apply>
	
	<apply>
                <lt/>
	
	  <ci> time </ci>
	  
	  <cn cellml:units="millisecond"> 0.5 </cn>
	  
	</apply>
	
      </apply>
      
    </piece>
    
    <piece>
    
      <cn cellml:units="microA_per_cm2"> 150 </cn>
	    
      <apply>
              <and/>
    
	<apply>
                <geq/>
	
	  <ci> time </ci>
	  
	  <cn cellml:units="millisecond"> 50 </cn>
	  
	</apply>
	
	<apply>
                <lt/>
	
	  <ci> time </ci>
	  
	  <cn cellml:units="millisecond"> 50.5 </cn>
	  
	</apply>
	
      </apply>
      
    </piece>
    
    <piece>
    
      <cn cellml:units="microA_per_cm2"> 150 </cn>
    
      <apply>
              <and/>
    	
	<apply>
                <geq/>
	
	  <ci> time </ci>
	  
	  <cn cellml:units="millisecond"> 100 </cn>
	  
	</apply>
	
	<apply>
                <lt/>
	
	  <ci> time </ci>
	  
	  <cn cellml:units="millisecond"> 100.5 </cn>
	  
	</apply>
	
      </apply>
      
    </piece>
    <piece>
    
      <cn cellml:units="microA_per_cm2"> 150 </cn>
    
      <apply>
              <and/>
    	
	<apply>
                <geq/>
	
	  <ci> time </ci>
	  
	  <cn cellml:units="millisecond"> 150 </cn>
	  
	</apply>
	
	<apply>
                <lt/>
	
	  <ci> time </ci>
	  
	  <cn cellml:units="millisecond"> 150.5 </cn>
	  
	</apply>
	
      </apply>
      
    </piece>
    <piece>
    
      <cn cellml:units="microA_per_cm2"> 150 </cn>
    
      <apply>
              <and/>
    	
	<apply>
                <geq/>
	
	  <ci> time </ci>
	  
	  <cn cellml:units="millisecond"> 200 </cn>
	  
	</apply>
	
	<apply>
                <lt/>
	
	  <ci> time </ci>
	  
	  <cn cellml:units="millisecond"> 200.5 </cn>
	  
	</apply>
	
      </apply>
      
    </piece>
    <piece>
    
      <cn cellml:units="microA_per_cm2"> 150 </cn>
    
      <apply>
              <and/>
    	
	<apply>
                <geq/>
	
	  <ci> time </ci>
	  
	  <cn cellml:units="millisecond"> 250 </cn>
	  
	</apply>
	
	<apply>
                <lt/>
	
	  <ci> time </ci>
	  
	  <cn cellml:units="millisecond"> 250.5 </cn>
	  
	</apply>
	
      </apply>
      
    </piece>
    <piece>
    
      <cn cellml:units="microA_per_cm2"> 150 </cn>
    
      <apply>
              <and/>
    	
	<apply>
                <geq/>
	
	  <ci> time </ci>
	  
	  <cn cellml:units="millisecond"> 300 </cn>
	  
	</apply>
	
	<apply>
                <lt/>
	
	  <ci> time </ci>
	  
	  <cn cellml:units="millisecond"> 300.5 </cn>
	  
	</apply>
	
      </apply>
      
    </piece>
    <piece>
    
      <cn cellml:units="microA_per_cm2"> 150 </cn>
    
      <apply>
              <and/>
    	
	<apply>
                <geq/>
	
	  <ci> time </ci>
	  
	  <cn cellml:units="millisecond"> 350 </cn>
	  
	</apply>
	
	<apply>
                <lt/>
	
	  <ci> time </ci>
	  
	  <cn cellml:units="millisecond"> 350.5 </cn>
	  
	</apply>
	
      </apply>
      
    </piece>
    <piece>
    
      <cn cellml:units="microA_per_cm2"> 150 </cn>
    
      <apply>
              <and/>
    	
	<apply>
                <geq/>
	
	  <ci> time </ci>
	  
	  <cn cellml:units="millisecond"> 400 </cn>
	  
	</apply>
	
	<apply>
                <lt/>
	
	  <ci> time </ci>
	  
	  <cn cellml:units="millisecond"> 400.5 </cn>
	  
	</apply>
	
      </apply>
      
    </piece>
    
    
    <otherwise>
    
      <cn cellml:units="microA_per_cm2"> 0 </cn>
      
    </otherwise>
    
  </piecewise>
  
</apply>

</math>
         




 
</component>


<component xmlns:cmeta="cmeta" name="sarco_Cl_channel">

  <variable units="dimensionless" name="a"/>
  <variable units="millivolt" public_interface="in" name="vS"/>
  <variable units="millivolt" public_interface="in" name="V_a"/>
  <variable units="millivolt" public_interface="in" name="A_a"/>
  <variable units="milliV_milliM" name="J_Cl"/>
  <variable units="milliM" public_interface="in" name="Cl_i"/>
  <variable units="milliM" public_interface="in" name="Cl_o"/>
  <variable units="C_per_mol" public_interface="in" name="FF"/>
  <variable units="milliJ_per_degreeK_per_mol" public_interface="in" name="RR"/>
  <variable units="degreeK" public_interface="in" name="TT"/>
  <variable units="milliS_per_cm2" public_interface="in" name="g_Cl_bar"/>
  <variable units="milliS_per_cm2" name="g_Cl"/>
  <variable units="microA_per_cm2" public_interface="out" name="I_Cl"/>
      
  <variable units="millisecond" public_interface="in" name="time"/>


  
<math xmlns="http://www.w3.org/1998/Math/MathML"> 
<apply id="a_Eqn">
        <eq/> 

  <ci> a </ci> 
 
  <apply>
          <divide/> 
  
    <cn cellml:units="dimensionless"> 1 </cn> 
 
    <apply>
            <plus/> 

      <cn cellml:units="dimensionless"> 1 </cn> 
 
      <apply>
              <exp/> 

        <apply>
                <divide/>
	
	  <apply>
                  <minus/>
	  	
	    <ci> vS </ci>
	    
	    <ci> V_a </ci> 
 
          </apply>
 
          <ci> A_a </ci>

        </apply>
  
      </apply>
 
    </apply>
  
  </apply>
  
</apply>
  
</math>  


<math xmlns="http://www.w3.org/1998/Math/MathML"> 
<apply id="J_Cl_eqn">
        <eq/>

  <ci> J_Cl </ci>
  
  <apply>
          <times/>
  
    <ci> vS </ci>
  
    <apply>
            <divide/>
  
      <apply>
              <minus/>
    
        <ci> Cl_i </ci>
      
        <apply>
                <times/>
      
          <ci> Cl_o </ci>
	
	  <apply>
                  <exp/>
	
	    <apply>
                    <divide/>
	  
	      <apply>
                      <times/>
	    
	        <ci> FF </ci>
	      
	        <ci> vS </ci>
	      
	      </apply>
	    
	      <apply>
                      <times/>
	    
	        <ci> RR </ci>
	      
	        <ci> TT </ci>
	      
	      </apply>
	    
	    </apply>
	  
	  </apply>
	
        </apply>
      
      </apply>
    
      <apply>
              <minus/>
    
        <cn cellml:units="dimensionless"> 1 </cn>
      
        <apply>
                <exp/>
      
          <apply>
                  <divide/>
	  
	    <apply>
                    <times/>
	    	      
	      <ci> FF </ci>
	      
	      <ci> vS </ci>
	      
	    </apply>
	    
	    <apply>
                    <times/>
	    
	      <ci> RR </ci>
	      
	      <ci> TT </ci>
	      
	    </apply>
	    
	  </apply>
	
        </apply>
      
      </apply>
    
    </apply>
  
  </apply>
  
</apply>

</math>
	


<math xmlns="http://www.w3.org/1998/Math/MathML"> 
<apply id="g_Cl_Eqn">
        <eq/> 

  <ci> g_Cl </ci>
  
  <apply>
          <times/>
   
    <ci> g_Cl_bar </ci>
 
    <apply>
            <power/> 

      <ci> a </ci>
 
      <cn cellml:units="dimensionless"> 4 </cn> 

    </apply>
  
  </apply>
  
</apply>
  
</math> 


<math xmlns="http://www.w3.org/1998/Math/MathML"> 
<apply id="I_Cl_Eqn">
        <eq/> 

  <ci> I_Cl </ci> 
 
  <apply>
          <times/> 
  
    <ci> g_Cl </ci>
 
    <apply>
            <divide/> 

      <ci> J_Cl </ci>
  
      <cn cellml:units="milliM"> 45 </cn>

    </apply>
  
  </apply>
  
</apply> 
 
</math>   
  

  
</component>


<component xmlns:cmeta="cmeta" name="sarco_IR_channel">

  <variable units="milliM" name="K_R"/>
  <variable units="milliM" public_interface="in" name="K_e"/>
  <variable units="dimensionless" public_interface="in" name="del"/>
  <variable units="millivolt" public_interface="in" name="E_K"/>
  <variable units="millivolt" public_interface="in" name="vS"/>
  <variable units="C_per_mol" public_interface="in" name="FF"/>
  <variable units="milliJ_per_degreeK_per_mol" public_interface="in" name="RR"/>
  <variable units="degreeK" public_interface="in" name="TT"/>
  <variable units="milliS_per_cm2" name="g_IR_bar"/>
  <variable units="milliS_per_cm2" public_interface="in" name="G_K"/>
  <variable units="milliM2" public_interface="in" name="K_K"/>
  <variable units="dimensionless" name="y"/>
  <variable units="milliM2" public_interface="in" name="K_S"/>
  <variable units="milliM" public_interface="in" name="S_i"/>
  <variable units="milliS_per_cm2" name="g_IR"/>
  <variable units="microA_per_cm2" public_interface="out" name="I_IR"/>
  <variable units="milliV_milliM" public_interface="in" name="J_K"/>
   
  <variable units="millisecond" public_interface="in" private_interface="out" name="time"/>
 

<math xmlns="http://www.w3.org/1998/Math/MathML"> 
<apply id="K_R_Eqn">
        <eq/> 

  <ci> K_R </ci> 
 
  <apply>
          <times/>
   
    <ci> K_e </ci>
 
    <apply>
            <exp/> 

      <apply>
              <times/>
  
        <apply>
                <times/>
  
          <apply>
                  <minus/>
  
            <ci> del </ci>
 
          </apply>

          <ci> E_K </ci>

        </apply>
  
        <apply>
                <divide/>
 
          <ci> FF </ci>
 
          <apply>
                  <times/> 

            <ci> RR </ci>
 
            <ci> TT </ci>

          </apply>
  
        </apply>
  
      </apply>
  
    </apply>
 
  </apply>
  
</apply>
  
</math> 

<math xmlns="http://www.w3.org/1998/Math/MathML"> 
<apply id="g_IR_bar_Eqn">
        <eq/> 

  <ci> g_IR_bar </ci> 
 
  <apply>
          <times/>
   
    <ci> G_K </ci>
 
    <apply>
            <divide/>
 
      <apply>
              <power/>
  
        <ci> K_R </ci>
 
        <cn cellml:units="dimensionless"> 2 </cn> 

      </apply>
  
      <apply>
              <plus/>
 
        <ci> K_K </ci>
 
        <apply>
                <power/>
 
          <ci> K_R </ci>
 
          <cn cellml:units="dimensionless"> 2 </cn> 

        </apply>
  
      </apply>
  
    </apply>
   
  </apply>
  
</apply>
  
</math>  




<math xmlns="http://www.w3.org/1998/Math/MathML"> 
<apply id="y_Eqn">
        <eq/> 
 
  <ci> y </ci>
  
  <apply>
          <minus/>
  
    <cn cellml:units="dimensionless"> 1 </cn>
  
    <apply>
            <power/> 
 
      <apply>
              <plus/>
   
        <cn cellml:units="dimensionless"> 1 </cn>
    
         <apply>
                <divide/>
    
          <apply>
                  <times/>
	  
            <ci> K_S </ci>
	    
	    <apply>
                    <plus/>
	    
	      <cn cellml:units="dimensionless"> 1 </cn>
	      
	      <apply>
                      <divide/>
	      
	        <apply>
                        <power/>
		
	          <ci> K_R </ci>
		  
	          <cn cellml:units="dimensionless"> 2 </cn>
		  
	        </apply>
		
	        <ci> K_K </ci>
		
	      </apply>
	      
	    </apply>
	    
          </apply>
  
          
	  
            <apply>
                  <times/>
	    
              <apply>
                    <power/>
	      
	        <ci> S_i </ci>
		
	        <cn cellml:units="dimensionless"> 2 </cn>
		
	      </apply>
	      
              <apply>
                    <exp/>
	      
                <apply>
                      <divide/>
		
                  <apply>
                        <times/>
		  
	            <cn cellml:units="dimensionless"> 2 </cn>
		    
	            <apply>
                          <minus/>
		    
	              <cn cellml:units="dimensionless"> 1 </cn>
		      
	              <ci> del </ci>
		      
	            </apply>
		    
	            <ci> vS </ci>
		    
	            <ci> FF </ci>
		    
	          </apply>
		  
	          <apply>
                        <times/>
		  
	            <ci> RR </ci>
		    
	            <ci> TT </ci>
		    
	          </apply>
		  
                </apply>
		
	      </apply>
	      
	    </apply>

          
      
        </apply>
    
      </apply>
      
      <cn cellml:units="dimensionless"> -1 </cn>
    
    </apply>
  
  </apply>
 
</apply>
  
</math> 


<math xmlns="http://www.w3.org/1998/Math/MathML"> 
<apply id="g_IR_Eqn">
        <eq/> 

  <ci> g_IR </ci> 
 
  <apply>
          <times/> 
  
    <ci> g_IR_bar </ci>
 
    <ci> y </ci>

  </apply>
  
</apply>
  
</math> 


<math xmlns="http://www.w3.org/1998/Math/MathML"> 
<apply id="I_IR_Eqn">
        <eq/>
 
  <ci> I_IR </ci> 
 
  <apply>
          <times/> 
  
    <ci> g_IR </ci>
    
    <piecewise>
  
      <piece>
           
        <cn cellml:units="dimensionless"> 1 </cn>
            
        <apply>
                <gt/>
              
          <ci> J_K </ci>
              
	  <cn cellml:units="milliV_milliM"> 0 </cn>
            
        </apply>
          
      </piece>
      
      <otherwise>
            
        <cn cellml:units="dimensionless"> 0 </cn>
                
      </otherwise>
        
    </piecewise>
    
    <apply>
            <divide/> 

      <ci> J_K </ci>
 
      <cn cellml:units="milliM"> 50 </cn>

    </apply>
  
  </apply>
  
</apply>
  
</math>



</component>


<component xmlns:cmeta="cmeta" name="sarco_DR_channel">

  
  <variable units="per_millisecond" name="alpha_n"/>
  <variable units="per_millisecond_per_millivolt" public_interface="in" name="alpha_n_bar"/>
  <variable units="millivolt" public_interface="in" name="V_n"/>
  <variable units="millivolt" public_interface="in" name="K_alpha_n"/>
  <variable units="per_millisecond" name="beta_n"/>
  <variable units="per_millisecond" public_interface="in" name="beta_n_bar"/>
  <variable units="millivolt" public_interface="in" name="K_beta_n"/>
  <variable units="dimensionless" name="h_K_inf"/>
  <variable units="millivolt" public_interface="in" name="V_h_K_inf"/>
  <variable units="millivolt" public_interface="in" name="A_h_K_inf"/>
  <variable units="millisecond" name="tau_h_K"/>
  <variable units="dimensionless" name="n" initial_value="0.009466"/>
  <variable units="dimensionless" name="h_K" initial_value="0.9952"/>
  <variable units="milliS_per_cm2" name="g_DR"/>
  <variable units="milliS_per_cm2" public_interface="in" name="g_K_bar"/>
  <variable units="microA_per_cm2" public_interface="out" name="I_DR"/>
  <variable units="milliV_milliM" public_interface="in" name="J_K"/>
  
  <variable units="millisecond" public_interface="in" name="time"/>
  
  <variable units="millivolt" public_interface="in" name="vS"/>
  

<math xmlns="http://www.w3.org/1998/Math/MathML"> 
<apply id="alpha_n_Eqn">
        <eq/> 

  <ci> alpha_n </ci> 
 
  <apply>
          <times/>
   
    <ci> alpha_n_bar </ci>
 
    <apply>
            <divide/> 

      <apply>
              <minus/>
     
        <ci> vS </ci>
	  
	<ci> V_n </ci>
	
      </apply>
 
      <apply>
              <minus/>
 
        <cn cellml:units="dimensionless"> 1 </cn>
  
        <apply>
                <exp/> 
 
          <apply>
                  <minus/> 
 
            <apply>
                    <divide/> 
 
              <apply>
                      <minus/>
        
	        <ci> vS </ci>
	  
	        <ci> V_n </ci>
	  
              </apply>
 
              <ci> K_alpha_n </ci> 

            </apply>
  
          </apply>

        </apply>
 
      </apply>
  
    </apply>
  
  </apply>
  
</apply>
  
</math> 

<math xmlns="http://www.w3.org/1998/Math/MathML"> 
<apply id="beta_n_Eqn">
        <eq/>
 
  <ci> beta_n </ci>
  
  <apply>
          <times/> 
  
    <ci> beta_n_bar </ci>
 
    <apply>
            <exp/>
 
      <apply>
              <minus/> 
 
        <apply>
                <divide/>
  
          <apply>
                  <minus/>
  	
	    <ci> vS </ci>
	  
	    <ci> V_n </ci>
	
          </apply>
 
          <ci> K_beta_n </ci> 

        </apply>
  
      </apply>

    </apply>
 
  </apply>
  
</apply>
  
</math> 

<math xmlns="http://www.w3.org/1998/Math/MathML"> 
<apply id="h_K_inf_Eqn">
        <eq/>
 
  <ci> h_K_inf </ci> 
 
  <apply>
          <divide/> 
  
    <cn cellml:units="dimensionless"> 1 </cn>
  
    <apply>
            <plus/> 

      <cn cellml:units="dimensionless"> 1 </cn>
  
      <apply>
              <exp/> 

        <apply>
                <divide/>
  
          <apply>
                  <minus/>
	  
	    <ci> vS </ci>
	    
	    <ci> V_h_K_inf </ci>
	    
	  </apply>
 
          <ci> A_h_K_inf </ci>

        </apply>
  
      </apply>
 
    </apply>
  
  </apply>
   
</apply>
  
</math>
 
    
<math xmlns="http://www.w3.org/1998/Math/MathML"> 
<apply id="tau_h_K_Eqn">
        <eq/> 

  <ci> tau_h_K </ci>
  
  <apply>
          <times/>
    
    <cn cellml:units="millisecond"> 1000 </cn>
  
    <apply>
            <exp/>
   
      <apply>
              <minus/>
  
        <apply>
                <divide/> 
 
          <apply>
                  <plus/>
	 
	    <ci> vS </ci>
	  
	    <cn cellml:units="millivolt"> 40 </cn>
	  
	  </apply>
 
          <cn cellml:units="millivolt"> 25.75 </cn>
	  
	</apply> 

      </apply>
  
    </apply>

  </apply>
 
</apply>
  
</math> 


<math xmlns="http://www.w3.org/1998/Math/MathML"> 
<apply id="dn_Eqn">
        <eq/> 
  <apply>
          <diff/>
    
    <bvar>
            <ci> time </ci>
          </bvar>
 
    <ci> n </ci>
    
  </apply> 

  <apply>
          <minus/>
   
    <apply>
            <times/>
  
      <ci> alpha_n </ci>
 
      <apply>
              <minus/> 

        <cn cellml:units="dimensionless"> 1 </cn>
  
        <ci> n </ci>
 
     </apply>
  
  </apply>
  
  <apply>
            <times/>
 
    <ci> beta_n </ci>
 
    <ci> n </ci>

  </apply>
  
 </apply>
  
</apply>
  
</math> 
 

<math xmlns="http://www.w3.org/1998/Math/MathML"> 
<apply id="dh_K_Eqn">
        <eq/> 

  <apply>
          <diff/>
    
    <bvar>
            <ci> time </ci>
          </bvar>
 
    <ci> h_K </ci>
    
  </apply> 
  
  <apply>
          <divide/>
   
    <apply>
            <minus/>
    
      <ci> h_K_inf </ci>
      
      <ci> h_K </ci>
      
    </apply>
 
    <ci> tau_h_K </ci>

  </apply>
  
</apply>
  
</math> 


 
<math xmlns="http://www.w3.org/1998/Math/MathML"> 
<apply id="g_DR_Eqn">
        <eq/> 

  <ci> g_DR </ci> 
 
  <apply>
          <times/>
   
    <apply>
            <times/>
  
      <ci> g_K_bar </ci>
 
      <apply>
              <power/>
        
	<ci> n </ci>
	
	<cn cellml:units="dimensionless"> 4 </cn>
	
      </apply>

    </apply>
  
    <ci> h_K </ci>

  </apply>
  
</apply>
  
</math>   
  

<math xmlns="http://www.w3.org/1998/Math/MathML"> 
<apply id="I_DR_Eqn">
        <eq/> 

  <ci> I_DR </ci> 
 
  <apply>
          <times/>
  
    <ci> g_DR </ci>
 
    <apply>
            <divide/>
 
      <ci> J_K </ci>
 
      <cn cellml:units="milliM"> 50 </cn>
 
    </apply>
  
  </apply>
  
</apply>
  
</math>   
  
    
</component>


<component xmlns:cmeta="cmeta" name="sarco_Na_channel">

  <variable units="millivolt" public_interface="in" name="vS"/>
  <variable units="milliM" public_interface="in" name="Na_i"/>
  <variable units="milliM" public_interface="in" name="Na_e"/>
  <variable units="per_millisecond" name="alpha_h"/>
  <variable units="per_millisecond" public_interface="in" name="alpha_h_bar"/>
  <variable units="millivolt" public_interface="in" name="V_h"/>
  <variable units="millivolt" public_interface="in" name="K_alpha_h"/>
  <variable units="per_millisecond" name="beta_h"/>
  <variable units="per_millisecond" public_interface="in" name="beta_h_bar"/>
  <variable units="millivolt" public_interface="in" name="K_beta_h"/>
  <variable units="per_millisecond" name="alpha_m"/>
  <variable units="per_millisecond_per_millivolt" public_interface="in" name="alpha_m_bar"/>
  <variable units="millivolt" public_interface="in" name="V_m"/>
  <variable units="millivolt" public_interface="in" name="K_alpha_m"/>
  <variable units="per_millisecond" name="beta_m"/>
  <variable units="per_millisecond" public_interface="in" name="beta_m_bar"/>
  <variable units="millivolt" public_interface="in" name="K_beta_m"/>
  <variable units="dimensionless" name="S_inf"/>
  <variable units="millivolt" public_interface="in" name="V_S_inf"/>
  <variable units="millivolt" public_interface="in" name="A_S_inf"/>
  <variable units="C_per_mol" public_interface="in" name="FF"/>
  <variable units="milliJ_per_degreeK_per_mol" public_interface="in" name="RR"/>
  <variable units="degreeK" public_interface="in" name="TT"/>
  <variable units="millisecond" name="tau_S"/>
  <variable units="millivolt" public_interface="in" name="V_tau"/>
  <variable units="dimensionless" name="m" initial_value="0.0358"/>
  <variable units="dimensionless" name="h" initial_value="0.4981"/>
  <variable units="dimensionless" name="S" initial_value="0.581"/>
  <variable units="milliS_per_cm2" name="g_Na"/>
  <variable units="milliS_per_cm2" public_interface="in" name="g_Na_bar"/>
  <variable units="milliV_milliM" name="J_Na"/>
  <variable units="microA_per_cm2" public_interface="out" name="I_Na"/>
  
  <variable units="millisecond" public_interface="in" name="time"/>

<math xmlns="http://www.w3.org/1998/Math/MathML"> 
<apply id="alpha_h_Eqn">
        <eq/> 

  <ci> alpha_h </ci> 
 
  <apply>
          <times/>
   
    <ci> alpha_h_bar </ci>
 
    <apply>
            <exp/> 

      <apply>
              <minus/>
  
        <apply>
                <divide/> 
 
          <apply>
                  <minus/>
	  
	    <ci> vS </ci>
	    
	    <ci> V_h </ci>
	    
	  </apply>
 
          <ci> K_alpha_h </ci> 

        </apply>
  
      </apply>

    </apply>
 
  </apply>
  
</apply> 
 
</math> 

<math xmlns="http://www.w3.org/1998/Math/MathML"> 
<apply id="beta_h_Eqn">
        <eq/> 

  <ci> beta_h </ci> 
 
  <apply>
          <divide/>
   
    <ci> beta_h_bar </ci>
 
    <apply>
            <plus/> 

      <cn cellml:units="dimensionless"> 1 </cn> 
 
      <apply>
              <exp/>
 
        <apply>
                <minus/>
  
          <apply>
                  <divide/> 
 
            <apply>
                    <minus/>
	  
	      <ci> vS </ci>
	    
	      <ci> V_h </ci>
	    
	    </apply>
 
            <ci> K_beta_h </ci> 

          </apply>
  
        </apply>

      </apply>
 
    </apply>
  
  </apply>
  
</apply>
  
</math>  


<math xmlns="http://www.w3.org/1998/Math/MathML"> 
<apply id="alpha_m_Eqn">
        <eq/> 

  <ci> alpha_m </ci>
  
  <apply>
          <times/> 
  
    <ci> alpha_m_bar </ci>
 
    <apply>
            <divide/>
 
      <apply>
              <minus/>
      
        <ci> vS </ci>
	
	<ci> V_m </ci>
	
      </apply>
 
      <apply>
              <minus/> 

        <cn cellml:units="dimensionless"> 1 </cn>
 
        <apply>
                <exp/>
  
          <apply>
                  <minus/> 
 
            <apply>
                    <divide/> 
 
              <apply>
                      <minus/>
      
                <ci> vS </ci>
	
	        <ci> V_m </ci>
	
              </apply>
 
              <ci> K_alpha_m </ci>
 
            </apply>

          </apply>
 
        </apply>

      </apply>
  
    </apply>
  
  </apply>
  
</apply> 
 
</math> 

<math xmlns="http://www.w3.org/1998/Math/MathML"> 
<apply id="beta_m_Eqn">
        <eq/> 

  <ci> beta_m </ci> 
 
  <apply>
          <times/>
   
    <ci> beta_m_bar </ci>
 
    <apply>
            <exp/> 

      <apply>
              <minus/>
  
        <apply>
                <divide/>
  
          <apply>
                  <minus/>
      
             <ci> vS </ci>
	
	     <ci> V_m </ci>
	
          </apply>
 
          <ci> K_beta_m </ci>
 
        </apply>
  
      </apply>

    </apply>
 
  </apply>
  
</apply>
  
</math> 


<math xmlns="http://www.w3.org/1998/Math/MathML"> 
<apply id="S_inf_Eqn">
        <eq/>
 
  <ci> S_inf </ci> 
 
  <apply>
          <divide/> 
  
    <cn cellml:units="dimensionless"> 1 </cn>
  
    <apply>
            <plus/> 

      <cn cellml:units="dimensionless"> 1 </cn>
  
      <apply>
              <exp/> 

        <apply>
                <divide/>
  
          <apply>
                  <minus/>
	  
	    <ci> vS </ci>
	    
	    <ci> V_S_inf </ci>
	    
	  </apply>
 
          <ci> A_S_inf </ci>

        </apply>
  
      </apply>
 
    </apply>
  
  </apply>
   
</apply>
  
</math>
 
    
<math xmlns="http://www.w3.org/1998/Math/MathML"> 
<apply id="tau_S_Eqn">
        <eq/> 

  <ci> tau_S </ci>
  
  <apply>
          <divide/>
  
    <!-- factor of 1000 added on 25/6/07  
    
    -->
    
    
    <cn cellml:units="millisecond"> 8571 </cn>
    
    
    <apply>
            <plus/>
    
      <cn cellml:units="dimensionless"> 0.2 </cn>
      
      <apply>
              <times/>
      
        <cn cellml:units="dimensionless"> 5.65 </cn>
	
	<apply>
                <power/>
	
	  <apply>
                  <divide/>
	
	    <apply>
                    <plus/>
	  
	      <ci> vS </ci>
	    
	      <ci> V_tau </ci>
	    
	    </apply>
	    
	    <cn cellml:units="millivolt"> 100 </cn>
	  
	  </apply>
	  
	  <cn cellml:units="dimensionless"> 2 </cn>
	  
	</apply>
	
      </apply>
      
    </apply>
    
  </apply>
 
</apply>
  
</math> 


<math xmlns="http://www.w3.org/1998/Math/MathML"> 
<apply id="J_Na_eqn">
        <eq/>

  <ci> J_Na </ci>
  
  <apply>
          <times/>
  
    <ci> vS </ci>
  
    <apply>
            <divide/>
  
      <apply>
              <minus/>
    
        <ci> Na_i </ci>
      
        <apply>
                <times/>
      
          <ci> Na_e </ci>
	
	  <apply>
                  <exp/>
	
	    <apply>
                    <divide/>
	  
	      <apply>
                      <times/>
	    
	        <cn cellml:units="dimensionless"> -1 </cn>
	    	
	        <ci> FF </ci>
	      
	        <ci> vS </ci>
	      
	      </apply>
	    
	      <apply>
                      <times/>
	    
	        <ci> RR </ci>
	      
	        <ci> TT </ci>
	      
	      </apply>
	    
	    </apply>
	  
	  </apply>
	
        </apply>
      
      </apply>
    
      <apply>
              <minus/>
    
        <cn cellml:units="dimensionless"> 1 </cn>
      
        <apply>
                <exp/>
      
          <apply>
                  <divide/>
	  
	    <apply>
                    <times/>
	    
	      <cn cellml:units="dimensionless"> -1 </cn>
	    	      
	      <ci> FF </ci>
	      
	      <ci> vS </ci>
	      
	    </apply>
	    
	    <apply>
                    <times/>
	    
	      <ci> RR </ci>
	      
	      <ci> TT </ci>
	      
	    </apply>
	    
	  </apply>
	
        </apply>
      
      </apply>
    
    </apply>
  
  </apply>
  
</apply>

</math>
	
<math xmlns="http://www.w3.org/1998/Math/MathML"> 
<apply id="dm_Eqn">
        <eq/> 

  <apply>
          <diff/>
    
    <bvar>
            <ci> time </ci>
          </bvar>
 
    <ci> m </ci>
    
  </apply>
 
  <apply>
          <minus/>
   
    <apply>
            <times/> 
 
      <ci> alpha_m </ci>
 
      <apply>
              <minus/>
 
        <cn cellml:units="dimensionless"> 1 </cn> 
      
          <ci> m </ci>
 
      </apply>
  
   </apply>
  
   <apply>
            <times/>
 
       <ci> beta_m </ci>
 
     <ci> m </ci>

   </apply>
  
 </apply>
  
</apply>
  
</math>   


<math xmlns="http://www.w3.org/1998/Math/MathML"> 
<apply id="dh_Eqn">
        <eq/> 

  <apply>
          <diff/>
    
    <bvar>
            <ci> time </ci>
          </bvar>
 
    <ci> h </ci>
    
  </apply>
  
  <apply>
          <minus/> 
  
    <apply>
            <times/> 
 
      <ci> alpha_h </ci>
 
      <apply>
              <minus/>

        <cn cellml:units="dimensionless"> 1 </cn>
  
        <ci> h </ci>
 
      </apply>
 
    </apply>
  
    <apply>
            <times/> 

      <ci> beta_h </ci>

      <ci> h </ci>

    </apply>
  
  </apply>
  
</apply>
  
</math>   
  
  
<math xmlns="http://www.w3.org/1998/Math/MathML"> 
<apply id="dS_Eqn">
        <eq/> 

  <apply>
          <diff/>
    
    <bvar>
            <ci> time </ci>
          </bvar>
 
    <ci> S </ci>
    
  </apply> 
  
    <apply>
          <divide/>
   
    <apply>
            <minus/>
    
      <ci> S_inf </ci>
      
      <ci> S </ci>
      
    </apply>
 
    <ci> tau_S </ci>

  </apply>
  
</apply>
  
</math>     
  
  
<math xmlns="http://www.w3.org/1998/Math/MathML"> 
<apply id="g_Na_Eqn">
        <eq/>
 
  <ci> g_Na </ci>
  
  <apply>
          <times/>
   
    <apply>
            <times/>
  
      <apply>
              <times/>
  
        <ci> g_Na_bar </ci>
 
        <apply>
                <power/>
	
	  <ci> m </ci>
	  
	  <cn cellml:units="dimensionless"> 3 </cn>
	  
	</apply>

      </apply>
  
      <ci> h </ci>

    </apply>
  
    <ci> S </ci>

  </apply>
  
</apply>
  
</math> 



<math xmlns="http://www.w3.org/1998/Math/MathML"> 
<apply id="I_Na_Eqn">
        <eq/> 

  <ci> I_Na </ci> 
 
  <apply>
          <times/> 
  
    <ci> g_Na </ci>
 
    <apply>
            <divide/> 

      <ci> J_Na </ci>
 
      <cn cellml:units="milliM"> 75 </cn>
 
    </apply>
  
  </apply>
  
</apply> 
 
</math>   
 
   
</component> 


<component xmlns:cmeta="cmeta" name="sarco_NaK_channel">


  <variable units="dimensionless" name="sig"/>
  <variable units="dimensionless" name="f1"/>
  <variable units="C_per_mol" public_interface="in" name="FF"/>
  <variable units="milliJ_per_degreeK_per_mol" public_interface="in" name="RR"/>
  <variable units="degreeK" public_interface="in" name="TT"/>
  <variable units="microA_per_cm2" name="I_NaK_bar"/>
  <variable units="milliM" public_interface="in" name="K_m_K"/>
  <variable units="milliM" public_interface="in" name="K_m_Na"/>
  <variable units="microA_per_cm2" public_interface="out" name="I_NaK"/>
  <variable units="micro_mol_per_cm2_per_second" public_interface="in" name="J_NaK_bar"/>    
  <variable units="milliM" public_interface="in" name="K_e"/>
  <variable units="milliM" public_interface="in" name="Na_i"/>
  <variable units="milliM" public_interface="in" name="Na_e"/>
  
  <variable units="millisecond" public_interface="in" name="time"/>
  
  <variable units="millivolt" public_interface="in" name="vS"/>
  

<math xmlns="http://www.w3.org/1998/Math/MathML"> 
<apply id="sig_Eqn">
        <eq/>

  <ci> sig </ci>
  
  <apply>
          <times/>
  
    <apply>
            <divide/>
  
      <cn cellml:units="dimensionless"> 1 </cn>
      
      <cn cellml:units="dimensionless"> 7 </cn>
      
    </apply>
    
    <apply>
            <minus/>
    
      <apply>
              <exp/>
      
        <apply>
                <divide/>
	
	  <ci> Na_e </ci>
	  
	  <cn cellml:units="milliM"> 67.3 </cn>
	  
	</apply>
	
      </apply>
      
      <cn cellml:units="dimensionless"> 1 </cn>
      
    </apply>
    
  </apply>
  
</apply>

</math>


  
<math xmlns="http://www.w3.org/1998/Math/MathML"> 
<apply id="f1_Eqn">
        <eq/> 

  <ci> f1 </ci>
  
  <apply>
          <power/> 
  
    <apply>
            <plus/>
 
      <cn cellml:units="dimensionless"> 1 </cn>
      
      <apply>
              <times/>
      
        <cn cellml:units="dimensionless"> 0.12 </cn>
	
	<apply>
                <exp/>
	
	  <apply>
                  <times/>
	  
	    <cn cellml:units="dimensionless"> -0.1 </cn>
	    
	    <ci> vS </ci>
	    
	    <apply>
                    <divide/>
	      
	      <ci> FF </ci>
	      
	      <apply>
                      <times/>
	      
	        <ci> RR </ci>
		
		<ci> TT </ci>
		
              </apply>
	      
	    </apply>
	    
	  </apply>
	  
	</apply>

      </apply>
	    
      <apply>
              <times/>
      
        <cn cellml:units="dimensionless"> 0.04 </cn>
	
	<ci> sig </ci>
	
	<apply>
                <exp/>
	
	  <apply>
                  <minus/>
	
	    <apply>
                    <times/>
	  
	      <ci> vS </ci>
	    
	      <apply>
                      <divide/>
	    
	        <ci> FF </ci>
	      
	        <apply>
                        <times/>
	      
	          <ci> RR </ci>
	      
	          <ci> TT </ci>
		  
		</apply>
	      
	      </apply>
	      
	    </apply>
	    
	  </apply>
      
        </apply>
      
      </apply>
  
    </apply>
    
    <cn cellml:units="dimensionless"> -1 </cn>

  </apply>
  
</apply> 

</math> 


<math xmlns="http://www.w3.org/1998/Math/MathML"> 
<apply id="I_NaK_bar_Eqn">
        <eq/>
 
  <ci> I_NaK_bar </ci>
  
  <apply>
          <times/>
   
    <ci> FF </ci>
 
    <apply>
            <divide/> 

      <ci> J_NaK_bar </ci>
 
      <apply>
              <times/> 

        <apply>
                <power/> 
 
          <apply>
                  <plus/>
	  
	    <cn cellml:units="dimensionless"> 1 </cn>
	  
	    <apply>
                    <divide/>
	    
	      <ci> K_m_K </ci>
	      
	      <ci> K_e </ci>
	      
	    </apply>
	    
	  </apply>
 
          <cn cellml:units="dimensionless"> 2 </cn> 

        </apply>
  
        <apply>
                <power/> 

          <apply>
                  <plus/>
	  
	    <cn cellml:units="dimensionless"> 1 </cn>
	  
	    <apply>
                    <divide/>
	    
	      <ci> K_m_Na </ci>
	      
	      <ci> Na_i </ci>
	      
	    </apply>
	    
	  </apply>
 
          <cn cellml:units="dimensionless"> 3 </cn> 

        </apply>
  
      </apply>
  
    </apply>
  
  </apply>
  
</apply>
 
</math>   
  


<math xmlns="http://www.w3.org/1998/Math/MathML"> 
<apply id="I_NaK_Eqn">
        <eq/> 

  <ci> I_NaK </ci>
  
  <apply>
          <times/>
   
    <ci> I_NaK_bar </ci>
 
    <ci> f1 </ci>

  </apply>
  
</apply> 

</math> 



</component>





<component xmlns:cmeta="cmeta" name="t_Cl_channel">

  <variable units="dimensionless" name="a_t"/>
  <variable units="millivolt" public_interface="in" name="vT"/>
  <variable units="millivolt" public_interface="in" name="V_a"/>
  <variable units="millivolt" public_interface="in" name="A_a"/>
  <variable units="milliV_milliM" name="J_Cl_t"/>
  <variable units="milliM" public_interface="in" name="Cl_i_t"/>
  <variable units="milliM" public_interface="in" name="Cl_o_t"/>
  <variable units="C_per_mol" public_interface="in" name="FF"/>
  <variable units="milliJ_per_degreeK_per_mol" public_interface="in" name="RR"/>
  <variable units="degreeK" public_interface="in" name="TT"/>
  <variable units="milliS_per_cm2" public_interface="in" name="g_Cl_bar"/>
  <variable units="microA_per_cm2" public_interface="out" name="I_Cl_t"/>
  <variable units="milliS_per_cm2" name="g_Cl_t"/>
  <variable units="dimensionless" public_interface="in" name="eta_Cl"/>    
  <variable units="millisecond" public_interface="in" private_interface="out" name="time"/>


  
<math xmlns="http://www.w3.org/1998/Math/MathML"> 
<apply id="a_t_Eqn">
        <eq/> 

  <ci> a_t </ci> 
 
  <apply>
          <divide/> 
  
    <cn cellml:units="dimensionless"> 1 </cn> 
 
    <apply>
            <plus/> 

      <cn cellml:units="dimensionless"> 1 </cn> 
 
      <apply>
              <exp/> 

        <apply>
                <divide/>
	
	  <apply>
                  <minus/>
	  	
	    <ci> vT </ci>
	    
	    <ci> V_a </ci> 
 
          </apply>
 
          <ci> A_a </ci>

        </apply>
  
      </apply>
 
    </apply>
  
  </apply>
  
</apply>
  
</math>  


<math xmlns="http://www.w3.org/1998/Math/MathML"> 
<apply id="J_Cl_t_eqn">
        <eq/>

  <ci> J_Cl_t </ci>
  
  <apply>
          <times/>
  
    <ci> vT </ci>
  
    <apply>
            <divide/>
  
      <apply>
              <minus/>
    
        <ci> Cl_i_t </ci>
      
        <apply>
                <times/>
      
          <ci> Cl_o_t </ci>
	
	  <apply>
                  <exp/>
	
	    <apply>
                    <divide/>
	  
	      <apply>
                      <times/>
	    
	        <ci> FF </ci>
	      
	        <ci> vT </ci>
	      
	      </apply>
	    
	      <apply>
                      <times/>
	    
	        <ci> RR </ci>
	      
	        <ci> TT </ci>
	      
	      </apply>
	    
	    </apply>
	  
	  </apply>
	
        </apply>
      
      </apply>
    
      <apply>
              <minus/>
    
        <cn cellml:units="dimensionless"> 1 </cn>
      
        <apply>
                <exp/>
      
          <apply>
                  <divide/>
	  
	    <apply>
                    <times/>
	    	      
	      <ci> FF </ci>
	      
	      <ci> vT </ci>
	      
	    </apply>
	    
	    <apply>
                    <times/>
	    
	      <ci> RR </ci>
	      
	      <ci> TT </ci>
	      
	    </apply>
	    
	  </apply>
	
        </apply>
      
      </apply>
    
    </apply>
  
  </apply>
  
</apply>

</math>
	


<math xmlns="http://www.w3.org/1998/Math/MathML"> 
<apply id="g_Cl_t_Eqn">
        <eq/> 

  <ci> g_Cl_t </ci>
  
  <apply>
          <times/>
   
    <ci> g_Cl_bar </ci>
 
    <apply>
            <power/> 

      <ci> a_t </ci>
 
      <cn cellml:units="dimensionless"> 4 </cn> 

    </apply>
  
  </apply>
  
</apply>
  
</math> 


<math xmlns="http://www.w3.org/1998/Math/MathML"> 
<apply id="I_Cl_t_Eqn">
        <eq/> 

  <ci> I_Cl_t </ci> 
 
  <apply>
          <times/> 
  
    <ci> eta_Cl </ci>
  
    <ci> g_Cl_t </ci>
 
    <apply>
            <divide/> 

      <ci> J_Cl_t </ci>
  
      <cn cellml:units="milliM"> 45 </cn>

    </apply>
  
  </apply>
  
</apply> 
 
</math>   
  

  
</component>


<component xmlns:cmeta="cmeta" name="t_IR_channel">

  <variable units="milliM" name="K_R_t"/>
  <variable units="milliM" public_interface="in" name="K_t"/>
  <variable units="dimensionless" public_interface="in" name="del"/>
  <variable units="millivolt" public_interface="in" name="E_K_t"/>
  <variable units="millivolt" public_interface="in" name="vT"/>
  <variable units="C_per_mol" public_interface="in" name="FF"/>
  <variable units="milliJ_per_degreeK_per_mol" public_interface="in" name="RR"/>
  <variable units="degreeK" public_interface="in" name="TT"/>
  <variable units="milliS_per_cm2" name="g_IR_bar_t"/>
  <variable units="milliS_per_cm2" public_interface="in" name="G_K"/>
  <variable units="milliM2" public_interface="in" name="K_K"/>
  <variable units="dimensionless" name="y_t"/>
  <variable units="milliM2" public_interface="in" name="K_S"/>
  <variable units="milliM" public_interface="in" name="S_i"/>
  <variable units="milliS_per_cm2" name="g_IR_t"/>
  <variable units="microA_per_cm2" public_interface="out" name="I_IR_t"/>
  <variable units="milliV_milliM" public_interface="in" name="J_K_t"/>
  <variable units="dimensionless" public_interface="in" name="eta_IR"/> 
  <variable units="millisecond" public_interface="in" private_interface="out" name="time"/>
 

<math xmlns="http://www.w3.org/1998/Math/MathML"> 
<apply id="K_R_t_Eqn">
        <eq/> 

  <ci> K_R_t </ci> 
 
  <apply>
          <times/>
   
    <ci> K_t </ci>
 
    <apply>
            <exp/> 

      <apply>
              <times/>
  
        <apply>
                <times/>
  
          <apply>
                  <minus/>
  
            <ci> del </ci>
 
          </apply>

          <ci> E_K_t </ci>

        </apply>
  
        <apply>
                <divide/>
 
          <ci> FF </ci>
 
          <apply>
                  <times/> 

            <ci> RR </ci>
 
            <ci> TT </ci>

          </apply>
  
        </apply>
  
      </apply>
  
    </apply>
 
  </apply>
  
</apply>
  
</math> 

<math xmlns="http://www.w3.org/1998/Math/MathML"> 
<apply id="g_IR_bar_t_Eqn">
        <eq/> 

  <ci> g_IR_bar_t </ci> 
 
  <apply>
          <times/>
   
    <ci> G_K </ci>
 
    <apply>
            <divide/>
 
      <apply>
              <power/>
  
        <ci> K_R_t </ci>
 
        <cn cellml:units="dimensionless"> 2 </cn> 

      </apply>
  
      <apply>
              <plus/>
 
        <ci> K_K </ci>
 
        <apply>
                <power/>
 
          <ci> K_R_t </ci>
 
          <cn cellml:units="dimensionless"> 2 </cn> 

        </apply>
  
      </apply>
  
    </apply>
   
  </apply>
  
</apply>
  
</math>  




<math xmlns="http://www.w3.org/1998/Math/MathML"> 
<apply id="y_t_Eqn">
        <eq/> 
 
  <ci> y_t </ci>
  
  <apply>
          <minus/>
  
    <cn cellml:units="dimensionless"> 1 </cn>
  
    <apply>
            <power/> 
 
      <apply>
              <plus/>
   
        <cn cellml:units="dimensionless"> 1 </cn>
    
         <apply>
                <divide/>
    
          <apply>
                  <times/>
	  
            <ci> K_S </ci>
	    
	    <apply>
                    <plus/>
	    
	      <cn cellml:units="dimensionless"> 1 </cn>
	      
	      <apply>
                      <divide/>
	      
	        <apply>
                        <power/>
		
	          <ci> K_R_t </ci>
		  
	          <cn cellml:units="dimensionless"> 2 </cn>
		  
	        </apply>
		
	        <ci> K_K </ci>
		
	      </apply>
	      
	    </apply>
	    
          </apply>
  
          
	  
            <apply>
                  <times/>
	    
              <apply>
                    <power/>
	      
	        <ci> S_i </ci>
		
	        <cn cellml:units="dimensionless"> 2 </cn>
		
	      </apply>
	      
              <apply>
                    <exp/>
	      
                <apply>
                      <divide/>
		
                  <apply>
                        <times/>
		  
	            <cn cellml:units="dimensionless"> 2 </cn>
		    
	            <apply>
                          <minus/>
		    
	              <cn cellml:units="dimensionless"> 1 </cn>
		      
	              <ci> del </ci>
		      
	            </apply>
		    
	            <ci> vT </ci>
		    
	            <ci> FF </ci>
		    
	          </apply>
		  
	          <apply>
                        <times/>
		  
	            <ci> RR </ci>
		    
	            <ci> TT </ci>
		    
	          </apply>
		  
                </apply>
		
	      </apply>
	      
	    </apply>

          
      
        </apply>
    
      </apply>
      
      <cn cellml:units="dimensionless"> -1 </cn>
    
    </apply>
  
  </apply>
 
</apply>
  
</math> 


<math xmlns="http://www.w3.org/1998/Math/MathML"> 
<apply id="g_IR_t_Eqn">
        <eq/> 

  <ci> g_IR_t </ci> 
 
  <apply>
          <times/> 
  
    <ci> g_IR_bar_t </ci>
 
    <ci> y_t </ci>

  </apply>
  
</apply>
  
</math> 

<math xmlns="http://www.w3.org/1998/Math/MathML"> 
<apply id="I_IR_t_Eqn">
        <eq/>
 
  <ci> I_IR_t </ci> 
 
  <apply>
          <times/> 
  
    <ci> eta_IR </ci>
  
    <ci> g_IR_t </ci>
    
    <!--<ci> H </ci>-->
 
    <apply>
            <divide/> 

      <ci> J_K_t </ci>
 
      <cn cellml:units="milliM"> 50 </cn>

    </apply>
  
  </apply>
  
</apply>
  
</math>


</component>


<component xmlns:cmeta="cmeta" name="t_DR_channel">

  
  <variable units="per_millisecond" name="alpha_n_t"/>
  <variable units="per_millisecond_per_millivolt" public_interface="in" name="alpha_n_bar"/>
  <variable units="millivolt" public_interface="in" name="V_n"/>
  <variable units="millivolt" public_interface="in" name="K_alpha_n"/>
  <variable units="per_millisecond" name="beta_n_t"/>
  <variable units="per_millisecond" public_interface="in" name="beta_n_bar"/>
  <variable units="millivolt" public_interface="in" name="K_beta_n"/>
  <variable units="dimensionless" name="h_K_inf_t"/>
  <variable units="millivolt" public_interface="in" name="V_h_K_inf"/>
  <variable units="millivolt" public_interface="in" name="A_h_K_inf"/>
  <variable units="millisecond" name="tau_h_K_t"/>
  <variable units="dimensionless" name="n_t" initial_value="0.009466"/>
  <variable units="dimensionless" name="h_K_t" initial_value="0.9952"/>
  <variable units="milliS_per_cm2" name="g_DR_t"/>
  <variable units="milliS_per_cm2" public_interface="in" name="g_K_bar"/>
  <variable units="microA_per_cm2" public_interface="out" name="I_DR_t"/>
  <variable units="milliV_milliM" public_interface="in" name="J_K_t"/>
  <variable units="dimensionless" public_interface="in" name="eta_DR"/>
  <variable units="millisecond" public_interface="in" name="time"/>
  
  <variable units="millivolt" public_interface="in" name="vT"/>
  

<math xmlns="http://www.w3.org/1998/Math/MathML"> 
<apply id="alpha_n_t_Eqn">
        <eq/> 

  <ci> alpha_n_t </ci> 
 
  <apply>
          <times/>
   
    <ci> alpha_n_bar </ci>
 
    <apply>
            <divide/> 

      <apply>
              <minus/>
     
        <ci> vT </ci>
	  
	<ci> V_n </ci>
	
      </apply>
 
      <apply>
              <minus/>
 
        <cn cellml:units="dimensionless"> 1 </cn>
  
        <apply>
                <exp/> 
 
          <apply>
                  <minus/> 
 
            <apply>
                    <divide/> 
 
              <apply>
                      <minus/>
        
	        <ci> vT </ci>
	  
	        <ci> V_n </ci>
	  
              </apply>
 
              <ci> K_alpha_n </ci> 

            </apply>
  
          </apply>

        </apply>
 
      </apply>
  
    </apply>
  
  </apply>
  
</apply>
  
</math> 

<math xmlns="http://www.w3.org/1998/Math/MathML"> 
<apply id="beta_n_t_Eqn">
        <eq/>
 
  <ci> beta_n_t </ci>
  
  <apply>
          <times/> 
  
    <ci> beta_n_bar </ci>
 
    <apply>
            <exp/>
 
      <apply>
              <minus/> 
 
        <apply>
                <divide/>
  
          <apply>
                  <minus/>
  	
	    <ci> vT </ci>
	  
	    <ci> V_n </ci>
	
          </apply>
 
          <ci> K_beta_n </ci> 

        </apply>
  
      </apply>

    </apply>
 
  </apply>
  
</apply>
  
</math> 

<math xmlns="http://www.w3.org/1998/Math/MathML"> 
<apply id="h_K_inf_t_Eqn">
        <eq/>
 
  <ci> h_K_inf_t </ci> 
 
  <apply>
          <divide/> 
  
    <cn cellml:units="dimensionless"> 1 </cn>
  
    <apply>
            <plus/> 

      <cn cellml:units="dimensionless"> 1 </cn>
  
      <apply>
              <exp/> 

        <apply>
                <divide/>
  
          <apply>
                  <minus/>
	  
	    <ci> vT </ci>
	    
	    <ci> V_h_K_inf </ci>
	    
	  </apply>
 
          <ci> A_h_K_inf </ci>

        </apply>
  
      </apply>
 
    </apply>
  
  </apply>
   
</apply>
  
</math>
 
    
<math xmlns="http://www.w3.org/1998/Math/MathML"> 
<apply id="tau_h_K_t_Eqn">
        <eq/> 

  <ci> tau_h_K_t </ci>
  
  <apply>
          <times/>
  
    <cn cellml:units="millisecond"> 1 </cn>
  
    <apply>
            <exp/>
   
      <apply>
              <minus/>
  
        <apply>
                <divide/> 
 
          <apply>
                  <plus/>
	 
	    <ci> vT </ci>
	  
	    <cn cellml:units="millivolt"> 40 </cn>
	  
	  </apply>
 
          <cn cellml:units="millivolt"> 25.75 </cn>
	  
	</apply> 

      </apply>
  
    </apply>

  </apply>
 
</apply>
  
</math> 


<math xmlns="http://www.w3.org/1998/Math/MathML"> 
<apply id="dn_t_Eqn">
        <eq/> 
  <apply>
          <diff/>
    
    <bvar>
            <ci> time </ci>
          </bvar>
 
    <ci> n_t </ci>
    
  </apply> 

  <apply>
          <minus/>
   
    <apply>
            <times/>
  
      <ci> alpha_n_t </ci>
 
      <apply>
              <minus/> 

        <cn cellml:units="dimensionless"> 1 </cn>
  
        <ci> n_t </ci>
 
     </apply>
  
  </apply>
  
  <apply>
            <times/>
 
    <ci> beta_n_t </ci>
 
    <ci> n_t </ci>

  </apply>
  
 </apply>
  
</apply>
  
</math> 
 

<math xmlns="http://www.w3.org/1998/Math/MathML"> 
<apply id="dh_K_t_Eqn">
        <eq/> 

  <apply>
          <diff/>
    
    <bvar>
            <ci> time </ci>
          </bvar>
 
    <ci> h_K_t </ci>
    
  </apply> 
  
  <apply>
          <divide/>
   
    <apply>
            <minus/>
    
      <ci> h_K_inf_t </ci>
      
      <ci> h_K_t </ci>
      
    </apply>
 
    <ci> tau_h_K_t </ci>

  </apply>
  
</apply>
  
</math> 


 
<math xmlns="http://www.w3.org/1998/Math/MathML"> 
<apply id="g_DR_t_Eqn">
        <eq/> 

  <ci> g_DR_t </ci> 
 
  <apply>
          <times/>
   
    <apply>
            <times/>
  
      <ci> g_K_bar </ci>
 
      <apply>
              <power/>
        
	<ci> n_t </ci>
	
	<cn cellml:units="dimensionless"> 4 </cn>
	
      </apply>

    </apply>
  
    <ci> h_K_t </ci>

  </apply>
  
</apply>
  
</math>   
  

<math xmlns="http://www.w3.org/1998/Math/MathML"> 
<apply id="I_DR_t_Eqn">
        <eq/> 

  <ci> I_DR_t </ci> 
 
  <apply>
          <times/>
  
    <ci> eta_DR </ci>
  
    <ci> g_DR_t </ci>
 
    <apply>
            <divide/>
 
      <ci> J_K_t </ci>
 
      <cn cellml:units="milliM"> 50 </cn>
 
    </apply>
  
  </apply>
  
</apply>
  
</math>   
  
    
</component>


<component xmlns:cmeta="cmeta" name="t_Na_channel">
  
  
  <variable units="millivolt" public_interface="in" name="vT"/>
  <variable units="milliM" public_interface="in" name="Na_i"/>
  <variable units="milliM" public_interface="in" name="Na_t"/>
  <variable units="per_millisecond" name="alpha_h_t"/>
  <variable units="per_millisecond" public_interface="in" name="alpha_h_bar"/>
  <variable units="millivolt" public_interface="in" name="V_h"/>
  <variable units="millivolt" public_interface="in" name="K_alpha_h"/>
  <variable units="per_millisecond" name="beta_h_t"/>
  <variable units="per_millisecond" public_interface="in" name="beta_h_bar"/>
  <variable units="millivolt" public_interface="in" name="K_beta_h"/>
  <variable units="per_millisecond" name="alpha_m_t"/>
  <variable units="per_millisecond_per_millivolt" public_interface="in" name="alpha_m_bar"/>
  <variable units="millivolt" public_interface="in" name="V_m"/>
  <variable units="millivolt" public_interface="in" name="K_alpha_m"/>
  <variable units="per_millisecond" name="beta_m_t"/>
  <variable units="per_millisecond" public_interface="in" name="beta_m_bar"/>
  <variable units="millivolt" public_interface="in" name="K_beta_m"/>
  <variable units="dimensionless" name="S_inf_t"/>
  <variable units="millivolt" public_interface="in" name="V_S_inf"/>
  <variable units="millivolt" public_interface="in" name="A_S_inf"/>
  <variable units="C_per_mol" public_interface="in" name="FF"/>
  <variable units="milliJ_per_degreeK_per_mol" public_interface="in" name="RR"/>
  <variable units="degreeK" public_interface="in" name="TT"/>
  <variable units="millisecond" name="tau_S_t"/>
  <variable units="millivolt" public_interface="in" name="V_tau"/>
  <variable units="dimensionless" name="m_t" initial_value="0.0358"/>
  <variable units="dimensionless" name="h_t" initial_value="0.4981"/>
  <variable units="dimensionless" name="S_t" initial_value="0.581"/>
  <variable units="milliS_per_cm2" name="g_Na_t"/>
  <variable units="milliS_per_cm2" public_interface="in" name="g_Na_bar"/>
  <variable units="milliV_milliM" name="J_Na_t"/>
  <variable units="microA_per_cm2" public_interface="out" name="I_Na_t"/>
  <variable units="dimensionless" public_interface="in" name="eta_Na"/>
  <variable units="millisecond" public_interface="in" name="time"/>

  

<math xmlns="http://www.w3.org/1998/Math/MathML"> 
<apply id="alpha_h_t_Eqn">
        <eq/> 

  <ci> alpha_h_t </ci> 
 
  <apply>
          <times/>
   
    <ci> alpha_h_bar </ci>
 
    <apply>
            <exp/> 

      <apply>
              <minus/>
  
        <apply>
                <divide/> 
 
          <apply>
                  <minus/>
	  
	    <ci> vT </ci>
	    
	    <ci> V_h </ci>
	    
	  </apply>
 
          <ci> K_alpha_h </ci> 

        </apply>
  
      </apply>

    </apply>
 
  </apply>
  
</apply> 
 
</math> 

<math xmlns="http://www.w3.org/1998/Math/MathML"> 
<apply id="beta_h_t_Eqn">
        <eq/> 

  <ci> beta_h_t </ci> 
 
  <apply>
          <divide/>
   
    <ci> beta_h_bar </ci>
 
    <apply>
            <plus/> 

      <cn cellml:units="dimensionless"> 1 </cn> 
 
      <apply>
              <exp/>
 
        <apply>
                <minus/>
  
          <apply>
                  <divide/> 
 
            <apply>
                    <minus/>
	  
	      <ci> vT </ci>
	    
	      <ci> V_h </ci>
	    
	    </apply>
 
            <ci> K_beta_h </ci> 

          </apply>
  
        </apply>

      </apply>
 
    </apply>
  
  </apply>
  
</apply>
  
</math>  


<math xmlns="http://www.w3.org/1998/Math/MathML"> 
<apply id="alpha_m_t_Eqn">
        <eq/> 

  <ci> alpha_m_t </ci>
  
  <apply>
          <times/> 
  
    <ci> alpha_m_bar </ci>
 
    <apply>
            <divide/>
 
      <apply>
              <minus/>
      
        <ci> vT </ci>
	
	<ci> V_m </ci>
	
      </apply>
 
      <apply>
              <minus/> 

        <cn cellml:units="dimensionless"> 1 </cn>
 
        <apply>
                <exp/>
  
          <apply>
                  <minus/> 
 
            <apply>
                    <divide/> 
 
              <apply>
                      <minus/>
      
                <ci> vT </ci>
	
	        <ci> V_m </ci>
	
              </apply>
 
              <ci> K_alpha_m </ci>
 
            </apply>

          </apply>
 
        </apply>

      </apply>
  
    </apply>
  
  </apply>
  
</apply> 
 
</math> 

<math xmlns="http://www.w3.org/1998/Math/MathML"> 
<apply id="beta_m_t_Eqn">
        <eq/> 

  <ci> beta_m_t </ci> 
 
  <apply>
          <times/>
   
    <ci> beta_m_bar </ci>
 
    <apply>
            <exp/> 

      <apply>
              <minus/>
  
        <apply>
                <divide/>
  
          <apply>
                  <minus/>
      
             <ci> vT </ci>
	
	     <ci> V_m </ci>
	
          </apply>
 
          <ci> K_beta_m </ci>
 
        </apply>
  
      </apply>

    </apply>
 
  </apply>
  
</apply>
  
</math> 


<math xmlns="http://www.w3.org/1998/Math/MathML"> 
<apply id="S_inf_t_Eqn">
        <eq/>
 
  <ci> S_inf_t </ci> 
 
  <apply>
          <divide/> 
  
    <cn cellml:units="dimensionless"> 1 </cn>
  
    <apply>
            <plus/> 

      <cn cellml:units="dimensionless"> 1 </cn>
  
      <apply>
              <exp/> 

        <apply>
                <divide/>
  
          <apply>
                  <minus/>
	  
	    <ci> vT </ci>
	    
	    <ci> V_S_inf </ci>
	    
	  </apply>
 
          <ci> A_S_inf </ci>

        </apply>
  
      </apply>
 
    </apply>
  
  </apply>
   
</apply>
  
</math>
 
    
<math xmlns="http://www.w3.org/1998/Math/MathML"> 
<apply id="tau_S_t_Eqn">
        <eq/> 

  <ci> tau_S_t </ci>
  
  <apply>
          <divide/>
  
    <!-- factor of 1000 added on 25/6/07  
    
    -->
    <cn cellml:units="millisecond"> 8571 </cn>
    
    <apply>
            <plus/>
    
      <cn cellml:units="dimensionless"> 0.2 </cn>
      
      <apply>
              <times/>
      
        <cn cellml:units="dimensionless"> 5.65 </cn>
	
	<apply>
                <power/>
	
	  <apply>
                  <divide/>
	
	    <apply>
                    <plus/>
	  
	      <ci> vT </ci>
	    
	      <ci> V_tau </ci>
	    
	    </apply>
	    
	    <cn cellml:units="millivolt"> 100 </cn>
	  
	  </apply>
	  
	  <cn cellml:units="dimensionless"> 2 </cn>
	  
	</apply>
	
      </apply>
      
    </apply>
    
  </apply>
 
</apply>
  
</math> 


<math xmlns="http://www.w3.org/1998/Math/MathML"> 
<apply id="J_Na_t_eqn">
        <eq/>

  <ci> J_Na_t </ci>
  
  <apply>
          <times/>
  
    <ci> vT </ci>
  
    <apply>
            <divide/>
  
      <apply>
              <minus/>
    
        <ci> Na_i </ci>
      
        <apply>
                <times/>
      
          <ci> Na_t </ci>
	
	  <apply>
                  <exp/>
	
	    <apply>
                    <divide/>
	  
	      <apply>
                      <times/>
	    
	        <cn cellml:units="dimensionless"> -1 </cn>
	    	
	        <ci> FF </ci>
	      
	        <ci> vT </ci>
	      
	      </apply>
	    
	      <apply>
                      <times/>
	    
	        <ci> RR </ci>
	      
	        <ci> TT </ci>
	      
	      </apply>
	    
	    </apply>
	  
	  </apply>
	
        </apply>
      
      </apply>
    
      <apply>
              <minus/>
    
        <cn cellml:units="dimensionless"> 1 </cn>
      
        <apply>
                <exp/>
      
          <apply>
                  <divide/>
	  
	    <apply>
                    <times/>
	    
	      <cn cellml:units="dimensionless"> -1 </cn>
	    	      
	      <ci> FF </ci>
	      
	      <ci> vT </ci>
	      
	    </apply>
	    
	    <apply>
                    <times/>
	    
	      <ci> RR </ci>
	      
	      <ci> TT </ci>
	      
	    </apply>
	    
	  </apply>
	
        </apply>
      
      </apply>
    
    </apply>
  
  </apply>
  
</apply>

</math>
	
<math xmlns="http://www.w3.org/1998/Math/MathML"> 
<apply id="dm_t_Eqn">
        <eq/> 

  <apply>
          <diff/>
    
    <bvar>
            <ci> time </ci>
          </bvar>
 
    <ci> m_t </ci>
    
  </apply>
 
  <apply>
          <minus/>
   
    <apply>
            <times/> 
 
      <ci> alpha_m_t </ci>
 
      <apply>
              <minus/>
 
        <cn cellml:units="dimensionless"> 1 </cn> 
      
          <ci> m_t </ci>
 
      </apply>
  
   </apply>
  
   <apply>
            <times/>
 
       <ci> beta_m_t </ci>
 
     <ci> m_t </ci>

   </apply>
  
 </apply>
  
</apply>
  
</math>   


<math xmlns="http://www.w3.org/1998/Math/MathML"> 
<apply id="dh_t_Eqn">
        <eq/> 

  <apply>
          <diff/>
    
    <bvar>
            <ci> time </ci>
          </bvar>
 
    <ci> h_t </ci>
    
  </apply>
  
  <apply>
          <minus/> 
  
    <apply>
            <times/> 
 
      <ci> alpha_h_t </ci>
 
      <apply>
              <minus/>

        <cn cellml:units="dimensionless"> 1 </cn>
  
        <ci> h_t </ci>
 
      </apply>
 
    </apply>
  
    <apply>
            <times/> 

      <ci> beta_h_t </ci>

      <ci> h_t </ci>

    </apply>
  
  </apply>
  
</apply>
  
</math>   
  
  
<math xmlns="http://www.w3.org/1998/Math/MathML"> 
<apply id="dS_t_Eqn">
        <eq/> 

  <apply>
          <diff/>
    
    <bvar>
            <ci> time </ci>
          </bvar>
 
    <ci> S_t </ci>
    
  </apply> 
  
    <apply>
          <divide/>
   
    <apply>
            <minus/>
    
      <ci> S_inf_t </ci>
      
      <ci> S_t </ci>
      
    </apply>
 
    <ci> tau_S_t </ci>

  </apply>
  
</apply>
  
</math>     
  
  
<math xmlns="http://www.w3.org/1998/Math/MathML"> 
<apply id="g_Na_t_Eqn">
        <eq/>
 
  <ci> g_Na_t </ci>
  
  <apply>
          <times/>
   
    <apply>
            <times/>
  
      <apply>
              <times/>
  
        <ci> g_Na_bar </ci>
 
        <apply>
                <power/>
	
	  <ci> m_t </ci>
	  
	  <cn cellml:units="dimensionless"> 3 </cn>
	  
	</apply>

      </apply>
  
      <ci> h_t </ci>

    </apply>
  
    <ci> S_t </ci>

  </apply>
  
</apply>
  
</math> 



<math xmlns="http://www.w3.org/1998/Math/MathML"> 
<apply id="I_Na_t_Eqn">
        <eq/> 

  <ci> I_Na_t </ci> 
 
  <apply>
          <times/> 
  
    <ci> eta_Na </ci>
  
    <ci> g_Na_t </ci>
 
    <apply>
            <divide/> 

      <ci> J_Na_t </ci>
 
      <cn cellml:units="milliM"> 75 </cn>
 
    </apply>
  
  </apply>
  
</apply> 
 
</math>   
 
   
</component> 


<component xmlns:cmeta="cmeta" name="t_NaK_channel">


  <variable units="dimensionless" name="sig_t"/>
  <variable units="dimensionless" name="f1_t"/>
  <variable units="C_per_mol" public_interface="in" name="FF"/>
  <variable units="milliJ_per_degreeK_per_mol" public_interface="in" name="RR"/>
  <variable units="degreeK" public_interface="in" name="TT"/>
  <variable units="microA_per_cm2" name="I_NaK_bar_t"/>
  <variable units="milliM" public_interface="in" name="K_m_K"/>
  <variable units="milliM" public_interface="in" name="K_m_Na"/>
  <variable units="microA_per_cm2" public_interface="out" name="I_NaK_t"/>
  <variable units="micro_mol_per_cm2_per_second" public_interface="in" name="J_NaK_bar"/>    
  <variable units="milliM" public_interface="in" name="K_t"/>
  <variable units="milliM" public_interface="in" name="Na_i"/>
  <variable units="milliM" public_interface="in" name="Na_t"/>
  <variable units="dimensionless" public_interface="in" name="eta_NaK"/>
  <variable units="millisecond" public_interface="in" name="time"/>
  
  <variable units="millivolt" public_interface="in" name="vT"/>
  

<math xmlns="http://www.w3.org/1998/Math/MathML"> 
<apply id="sig_t_Eqn">
        <eq/>

  <ci> sig_t </ci>
  
  <apply>
          <times/>
  
    <apply>
            <divide/>
  
      <cn cellml:units="dimensionless"> 1 </cn>
      
      <cn cellml:units="dimensionless"> 7 </cn>
      
    </apply>
    
    <apply>
            <minus/>
    
      <apply>
              <exp/>
      
        <apply>
                <divide/>
	
	  <ci> Na_t </ci>
	  
	  <cn cellml:units="milliM"> 67.3 </cn>
	  
	</apply>
	
      </apply>
      
      <cn cellml:units="dimensionless"> 1 </cn>
      
    </apply>
    
  </apply>
  
</apply>

</math>


  
<math xmlns="http://www.w3.org/1998/Math/MathML"> 
<apply id="f1_t_Eqn">
        <eq/> 

  <ci> f1_t </ci>
  
  <apply>
          <power/> 
  
    <apply>
            <plus/>
 
      <cn cellml:units="dimensionless"> 1 </cn>
      
      <apply>
              <times/>
      
        <cn cellml:units="dimensionless"> 0.12 </cn>
	
	<apply>
                <exp/>
	
	  <apply>
                  <times/>
	  
	    <cn cellml:units="dimensionless"> -0.1 </cn>
	    
	    <ci> vT </ci>
	    
	    <apply>
                    <divide/>
	      
	      <ci> FF </ci>
	      
	      <apply>
                      <times/>
	      
	        <ci> RR </ci>
		
		<ci> TT </ci>
		
              </apply>
	      
	    </apply>
	    
	  </apply>
	  
	</apply>

      </apply>
	    
      <apply>
              <times/>
      
        <cn cellml:units="dimensionless"> 0.04 </cn>
	
	<ci> sig_t </ci>
	
	<apply>
                <exp/>
	
	  <apply>
                  <minus/>
	
	    <apply>
                    <times/>
	  
	      <ci> vT </ci>
	    
	      <apply>
                      <divide/>
	    
	        <ci> FF </ci>
	      
	        <apply>
                        <times/>
	      
	          <ci> RR </ci>
	      
	          <ci> TT </ci>
		  
		</apply>
	      
	      </apply>
	      
	    </apply>
	    
	  </apply>
      
        </apply>
      
      </apply>
  
    </apply>
    
    <cn cellml:units="dimensionless"> -1 </cn>

  </apply>
  
</apply> 

</math> 


<math xmlns="http://www.w3.org/1998/Math/MathML"> 
<apply id="I_NaK_bar_t_Eqn">
        <eq/>
 
  <ci> I_NaK_bar_t </ci>
  
  <apply>
          <times/>
   
    <ci> FF </ci>
 
    <apply>
            <divide/> 

      <ci> J_NaK_bar </ci>
 
      <apply>
              <times/> 

        <apply>
                <power/> 
 
          <apply>
                  <plus/>
	  
	    <cn cellml:units="dimensionless"> 1 </cn>
	  
	    <apply>
                    <divide/>
	    
	      <ci> K_m_K </ci>
	      
	      <ci> K_t </ci>
	      
	    </apply>
	    
	  </apply>
 
          <cn cellml:units="dimensionless"> 2 </cn> 

        </apply>
  
        <apply>
                <power/> 

          <apply>
                  <plus/>
	  
	    <cn cellml:units="dimensionless"> 1 </cn>
	  
	    <apply>
                    <divide/>
	    
	      <ci> K_m_Na </ci>
	      
	      <ci> Na_i </ci>
	      
	    </apply>
	    
	  </apply>
 
          <cn cellml:units="dimensionless"> 3 </cn> 

        </apply>
  
      </apply>
  
    </apply>
  
  </apply>
  
</apply>
 
</math>   
  


<math xmlns="http://www.w3.org/1998/Math/MathML"> 
<apply id="I_NaK_t_Eqn">
        <eq/> 

  <ci> I_NaK_t </ci>
  
  <apply>
          <times/>
  
    <ci> eta_NaK </ci>
   
    <ci> I_NaK_bar_t </ci>
 
    <ci> f1_t </ci>

  </apply>
  
</apply> 

</math> 



</component>
<!--

Stern and Rios Component


-->
<component name="sternrios">

  <!-- these variables are defined here and used in other components -->
  
  <variable units="dimensionless" public_interface="out" name="O_0" initial_value="0"/>
  <variable units="dimensionless" public_interface="out" name="O_1" initial_value="0"/>
  <variable units="dimensionless" public_interface="out" name="O_2" initial_value="0"/>
  <variable units="dimensionless" public_interface="out" name="O_3" initial_value="0"/>
  <variable units="dimensionless" public_interface="out" name="O_4" initial_value="0"/>
  
  <!-- these variables are defined here and only used internally -->
  
  <variable units="per_millisecond" name="k_L" initial_value="0.002"/>
  <variable units="per_millisecond" name="k_Lm" initial_value="1000"/>
  <variable units="dimensionless" name="f" initial_value="0.2"/>
  <variable units="per_millisecond" name="alpha1" initial_value="0.2"/>
  <variable units="millivolt" name="K" initial_value="4.5"/>
  <variable units="millivolt" name="Vbar" initial_value="-20"/>
 
  <variable units="dimensionless" name="C_0" initial_value="1"/>
  <variable units="dimensionless" name="C_1" initial_value="0"/>
  <variable units="dimensionless" name="C_2" initial_value="0"/>
  <variable units="dimensionless" name="C_3" initial_value="0"/>
  <variable units="dimensionless" name="C_4" initial_value="0"/>
  <variable units="per_millisecond" name="k_C"/>
  <variable units="per_millisecond" name="k_Cm"/>
  
  <!-- these variables are imported from other components -->
  
  <variable units="millisecond" public_interface="in" name="time"/>
  <variable units="millivolt" public_interface="in" name="vT"/>
     
<math xmlns="http://www.w3.org/1998/Math/MathML"> 
<apply id="k_C_Eqn">
        <eq/>

  <ci> k_C </ci>
  
  <apply>
          <times/>
  
    <cn cellml:units="dimensionless"> 0.5 </cn>
    
    <ci> alpha1 </ci>
    
    <apply>
            <exp/>
    
      <apply>
              <divide/>
      
        <apply>
                <minus/>
	
	  <ci> vT </ci>
	  
	  <ci> Vbar </ci>
	  
	</apply>
	
	<apply>
                <times/>
	
	  <cn cellml:units="dimensionless"> 8 </cn>
	  
	  <ci> K </ci>
	  
	</apply>
	
      </apply>
      
    </apply>
    
  </apply>
  
</apply>

</math>


<math xmlns="http://www.w3.org/1998/Math/MathML"> 
<apply id="k_Cm_Eqn">
        <eq/>

  <ci> k_Cm </ci>
  
  <apply>
          <times/>
  
    <cn cellml:units="dimensionless"> 0.5 </cn>
    
    <ci> alpha1 </ci>
    
    <apply>
            <exp/>
    
      <apply>
              <divide/>
      
        <apply>
                <minus/>
	
	  <ci> Vbar </ci>
	  
	  <ci> vT  </ci>
	  
	</apply>
	
	<apply>
                <times/>
	
	  <cn cellml:units="dimensionless"> 8 </cn>
	  
	  <ci> K </ci>
	  
	</apply>
	
      </apply>
      
    </apply>
    
  </apply>
  
</apply>

</math>

<math xmlns="http://www.w3.org/1998/Math/MathML"> 
<apply id="dC_0_Eqn">
        <eq/> 

  <apply>
          <diff/>
    
    <bvar>
            <ci> time </ci>
          </bvar>
 
    <ci> C_0 </ci>
    
  </apply>
  
  <apply>
          <plus/>
  
    <apply>
            <times/>
    
      <apply>
              <minus/>
     
        <ci> k_L </ci>
       
      </apply>
     
      <ci> C_0 </ci>
     
    </apply>
    
    <apply>
            <times/>
    
      <ci> k_Lm </ci>
      
      <ci> O_0 </ci>
     
    </apply> 
  
    <apply>
            <times/>
    
      <cn cellml:units="dimensionless"> -4 </cn>
    
      <ci> k_C </ci>
    
      <ci> C_0 </ci>
     
    </apply>
    
    <apply>
            <times/>
    
      <ci> k_Cm </ci>
    
      <ci> C_1 </ci>
     
    </apply>
    
  </apply>
  
</apply>

</math>



<math xmlns="http://www.w3.org/1998/Math/MathML"> 
<apply id="dO_0_Eqn">
        <eq/> 

  <apply>
          <diff/>
    
    <bvar>
            <ci> time </ci>
          </bvar>
 
    <ci> O_0 </ci>
    
  </apply>
  
  <apply>
          <plus/>
  
    <apply>
            <times/>
    
      <ci> k_L </ci>
    
      <ci> C_0 </ci>
     
    </apply>
    
    <apply>
            <times/>
    
      <apply>
              <minus/>
    
        <ci> k_Lm </ci>
	
      </apply>
      
      <ci> O_0 </ci>
     
    </apply>
    
    <apply>
            <divide/> 
  
      <apply>
              <times/>
    
        <cn cellml:units="dimensionless"> -4 </cn>
    
        <ci> k_C </ci>
	
        <ci> O_0 </ci>
	
      </apply>
     
      <ci> f </ci>

    </apply>
    
    <apply>
            <times/>
    
      <ci> f </ci>
    
      <ci> k_Cm </ci>
    
      <ci> O_1 </ci>
     
    </apply>
    
  </apply>
  
</apply>

</math>



<math xmlns="http://www.w3.org/1998/Math/MathML"> 
<apply id="dC_1_Eqn">
        <eq/> 

  <apply>
          <diff/>
    
    <bvar>
            <ci> time </ci>
          </bvar>
 
    <ci> C_1 </ci>
    
  </apply>
  
  <apply>
          <plus/>
  
    <apply>
            <times/>
    
      <cn cellml:units="dimensionless"> 4 </cn>
    
      <ci> k_C </ci>
    
      <ci> C_0 </ci>
     
    </apply>
    
    <apply>
            <times/>
    
      <apply>
              <minus/>
    
        <ci> k_Cm </ci>
	
      </apply>
      
      <ci> C_1 </ci>
     
    </apply>
    
    <apply>
            <divide/>
    
      <apply>
              <times/> 
  
        <apply>
                <minus/>
    
          <ci> k_L </ci>
    
        </apply>
	
	<ci> C_1 </ci>
   	
      </apply>
      
      <ci> f </ci>
	
    </apply>
    
    <apply>
            <times/>
    
      <ci> f </ci>
    
      <ci> k_Lm </ci>
    
      <ci> O_1 </ci>
     
    </apply>
    
    <apply>
            <times/>
    
      <cn cellml:units="dimensionless"> -3 </cn>
      
      <ci> k_C </ci>
      
      <ci> C_1 </ci>
      
    </apply>
    
    <apply>
            <times/>
    
      <cn cellml:units="dimensionless"> 2 </cn>
      
      <ci> k_Cm </ci>
      
      <ci> C_2 </ci>
      
    </apply>
    
  </apply>
  
</apply>

</math>



<math xmlns="http://www.w3.org/1998/Math/MathML"> 
<apply id="dO_1_Eqn">
        <eq/> 

  <apply>
          <diff/>
    
    <bvar>
            <ci> time </ci>
          </bvar>
 
    <ci> O_1 </ci>
    
  </apply>
  
  <apply>
          <plus/>
  
    <apply>
            <divide/>
    
      <apply>
              <times/>
    
        <ci> k_L </ci>
	  
        <ci> C_1 </ci>

      </apply>
     
      <ci> f </ci>
   
    </apply>
    
    <apply>
            <times/>
    
      <apply>
              <minus/>
    
        <ci> k_Lm </ci>
	
      </apply>
      
      <ci> f </ci>
      
      <ci> O_1 </ci>
     
    </apply>
    
    <apply>
            <divide/> 
  
      <apply>
              <times/>
      
        <cn cellml:units="dimensionless"> 4 </cn>
    
        <ci> k_C </ci>
	
	<ci> O_0 </ci>
  
      </apply>
     
      <ci> f </ci>
   
    </apply>
    
    <apply>
            <times/>
    
      <apply>
              <minus/>
    
        <ci> f </ci>
	
      </apply>
    
      <ci> k_Cm </ci>
    
      <ci> O_1 </ci>
     
    </apply>
    
    <apply>
            <divide/>
    
      <apply>
              <times/>
    
        <cn cellml:units="dimensionless"> -3 </cn>
      
        <ci> k_C </ci>
	
	<ci> O_1 </ci>
	
      </apply>
     
      <ci> f </ci>
   
    </apply>
    
    <apply>
            <times/>
    
      <cn cellml:units="dimensionless"> 2 </cn>
      
      <ci> f </ci>
      
      <ci> k_Cm </ci>
      
      <ci> O_2 </ci>
      
    </apply>
    
  </apply>
  
</apply>

</math>




<math xmlns="http://www.w3.org/1998/Math/MathML"> 
<apply id="dC_2_Eqn">
        <eq/> 

  <apply>
          <diff/>
    
    <bvar>
            <ci> time </ci>
          </bvar>
 
    <ci> C_2 </ci>
    
  </apply>
  
  <apply>
          <plus/>
  
    <apply>
            <times/>
    
      <cn cellml:units="dimensionless"> 3 </cn>
      
      <ci> k_C </ci>
    
      <ci> C_1 </ci>
	
    </apply>
    
    <apply>
            <times/>
    
      <cn cellml:units="dimensionless"> -2 </cn>
    
      <ci> k_Cm </ci>
   
      <ci> C_2 </ci>
     
    </apply>
    
    <apply>
            <divide/>
    
      <apply>
              <times/>
      
        <apply>
                <minus/>
      
          <ci> k_L </ci>
    
        </apply>
        
        <ci> C_2 </ci>
	
      </apply>
    
      <apply>
              <power/>
      
        <ci> f </ci>
	  
	<cn cellml:units="dimensionless"> 2 </cn>
	
      </apply>
      
    </apply>
    
    <apply>
            <times/>
    
      <apply>
              <power/>
    
        <ci> f </ci>
	
	<cn cellml:units="dimensionless"> 2 </cn>
	
      </apply>

      <ci> k_Lm </ci>
    
      <ci> O_2 </ci>
     
    </apply>
    
    <apply>
            <times/>
    
      <cn cellml:units="dimensionless"> -2 </cn>
      
      <ci> k_C </ci>
      
      <ci> C_2 </ci>
	 
    </apply>
    
    <apply>
            <times/>
    
      <cn cellml:units="dimensionless"> 3 </cn>
       
      <ci> k_Cm </ci>
      
      <ci> C_3 </ci>
      
    </apply>
    
  </apply>
  
</apply>

</math>




<math xmlns="http://www.w3.org/1998/Math/MathML"> 
<apply id="dO_2_Eqn">
        <eq/> 

  <apply>
          <diff/>
    
    <bvar>
            <ci> time </ci>
          </bvar>
 
    <ci> O_2 </ci>
    
  </apply>
  
  <apply>
          <plus/>
  
    <apply>
            <divide/>
  
      <apply>
              <times/>
    
        <cn cellml:units="dimensionless"> 3 </cn>
      
        <ci> k_C </ci>
	
        <ci> O_1 </ci>
	
      </apply>
      
      <ci> f </ci>
    
  </apply>
    
    <apply>
            <times/>
    
      <cn cellml:units="dimensionless"> -2 </cn>
      
      <ci> f </ci>
    
      <ci> k_Cm </ci>
   
      <ci> O_2 </ci>
     
    </apply>
    
    <apply>
            <divide/> 
    
      <apply>
              <times/>
        
        <ci> k_L </ci>
  
	<ci> C_2 </ci>
	
      </apply>

      <apply>
              <power/>
      
        <ci> f </ci>
	  
	<cn cellml:units="dimensionless"> 2 </cn>
	  
      </apply>
    
    </apply>
    
    <apply>
            <times/>
    
      <apply>
              <minus/>
         
        <ci> k_Lm </ci>
	  
      </apply>
    
      <apply>
              <power/>
      
        <ci> f </ci>
	
	<cn cellml:units="dimensionless"> 2 </cn>
	
      </apply>

      <ci> O_2 </ci>
     
    </apply>
    
    <apply>
            <divide/>
    
      <apply>
              <times/>
    
        <cn cellml:units="dimensionless"> -2 </cn>
      
        <ci> k_C </ci>
	
        <ci> O_2 </ci>
	
      </apply>
      
      <ci> f </ci>
      
    </apply>
    
    <apply>
            <times/>
    
      <cn cellml:units="dimensionless"> 3 </cn>
      
      <ci> f </ci>
       
      <ci> k_Cm </ci>
      
      <ci> O_3 </ci>
      
    </apply>
    
  </apply>
  
</apply>

</math>




<math xmlns="http://www.w3.org/1998/Math/MathML"> 
<apply id="dC_3_Eqn">
        <eq/> 

  <apply>
          <diff/>
    
    <bvar>
            <ci> time </ci>
          </bvar>
 
    <ci> C_3 </ci>
    
  </apply>
  
  <apply>
          <plus/>
  
    <apply>
            <times/>
  
      <cn cellml:units="dimensionless"> 2 </cn>
      
      <ci> k_C </ci>
 
      <ci> C_2 </ci>
	
    </apply>
    
    <apply>
            <times/>
    
      <cn cellml:units="dimensionless"> -3 </cn>
    
      <ci> k_Cm </ci>
   
      <ci> C_3 </ci>
     
    </apply>
    
    <apply>
            <divide/>
    
      <apply>
              <times/> 
  
        <apply>
                <minus/>
      
          <ci> k_L </ci>
    
        </apply>
    
	<ci> C_3 </ci>
	
      </apply>
  
      <apply>
              <power/>
      
        <ci> f </ci>
	  
	<cn cellml:units="dimensionless"> 3 </cn>
	
      </apply>
      
    </apply>
    
    <apply>
            <times/>
       
      <ci> k_Lm </ci>
	
      <apply>
              <power/>
      
        <ci> f </ci>
	
	<cn cellml:units="dimensionless"> 3 </cn>
	
      </apply>

      <ci> O_3 </ci>
     
    </apply>
    
    <apply>
            <times/>
    
      <apply>
              <minus/>
    
        <ci> k_C </ci>
	
      </apply>
      
      <ci> C_3 </ci>
	 
    </apply>
    
    <apply>
            <times/>
    
      <cn cellml:units="dimensionless"> 4 </cn>
       
      <ci> k_Cm </ci>
      
      <ci> C_4 </ci>
      
    </apply>
    
  </apply>
  
</apply>

</math>




<math xmlns="http://www.w3.org/1998/Math/MathML"> 
<apply id="dO_3_Eqn">
        <eq/> 

  <apply>
          <diff/>
    
    <bvar>
            <ci> time </ci>
          </bvar>
 
    <ci> O_3 </ci>
    
  </apply>
  
  <apply>
          <plus/>
  
    <apply>
            <divide/>
   
      <apply>
              <times/>
      
        <ci> k_L </ci>
     
        <ci> C_3 </ci>
	
      </apply>
	
      <apply>
              <power/>
	  
	<ci> f </ci>
	  
	<cn cellml:units="dimensionless"> 3 </cn>

      </apply>
	
    </apply>
    
    <apply>
            <times/>
    
      <apply>
              <minus/>
    
        <ci> k_Lm </ci>
	
      </apply>
      
      <apply>
              <power/>
	  
	<ci> f </ci>
	  
	<cn cellml:units="dimensionless"> 3 </cn>
	  
      </apply>
   
      <ci> O_3 </ci>
     
    </apply>
    
    <apply>
            <divide/> 
  
      <apply>
              <times/>
      
        <cn cellml:units="dimensionless"> 2 </cn>
      
        <ci> k_C </ci>
    
	<ci> O_2 </ci>

      </apply>
     
      <ci> f </ci>
  
    </apply>
    
    <apply>
            <times/>
       
      <cn cellml:units="dimensionless"> -3 </cn>
	
      <ci> k_Cm </ci>
   
      <ci> f </ci>
      
      <ci> O_3 </ci>
  
    </apply>
    
    <apply>
            <divide/>
    
      <apply>
              <times/>
    
        <apply>
                <minus/>
    
          <ci> k_C </ci>
	  
	</apply>
	 
      <ci> O_3 </ci>
	
      </apply>
      
      <ci> f </ci>
 	 
    </apply>
    
    <apply>
            <times/>
    
      <cn cellml:units="dimensionless"> 4 </cn>
      
      <ci> f </ci>
       
      <ci> k_Cm </ci>
      
      <ci> O_4 </ci>
      
    </apply>
    
  </apply>
  
</apply>

</math>




<math xmlns="http://www.w3.org/1998/Math/MathML"> 
<apply id="dC_4_Eqn">
        <eq/> 

  <apply>
          <diff/>
    
    <bvar>
            <ci> time </ci>
          </bvar>
 
    <ci> C_4 </ci>
    
  </apply>
  
  <apply>
          <plus/>
  
    <apply>
            <times/>
  
      <ci> k_C </ci>
      
      <ci> C_3 </ci>
	
    </apply>
    
    <apply>
            <times/>
    
      <cn cellml:units="dimensionless"> -4 </cn>
    
      <ci> k_Cm </ci>
      
      <ci> C_4 </ci>
     
    </apply>
    
    <apply>
            <divide/>
    
      <apply>
              <times/> 
  
        <apply>
                <minus/>
      
          <ci> k_L </ci>
    
        </apply>
      
	<ci> C_4 </ci>
	
      </apply>
   
      <apply>
              <power/>
      
        <ci> f </ci>
	  
	<cn cellml:units="dimensionless"> 4 </cn>
	
      </apply>
      
    </apply>
    
    <apply>
            <times/>
     
      <ci> k_Lm </ci>
      
      <apply>
              <power/>
      
        <ci> f </ci>
	  
        <cn cellml:units="dimensionless"> 4 </cn>
	  
      </apply>
      
      <ci> O_4 </ci>
  
    </apply>
   
  </apply>
  
</apply>

</math>




<math xmlns="http://www.w3.org/1998/Math/MathML"> 
<apply id="dO_4_Eqn">
        <eq/> 

  <apply>
          <diff/>
    
    <bvar>
            <ci> time </ci>
          </bvar>
 
    <ci> O_4 </ci>
    
  </apply>
  
  <apply>
          <plus/>
  
    <apply>
            <divide/>
  
      <apply>
              <times/>
      
        <ci> k_C </ci>
      
        <ci> O_3 </ci>
	
      </apply>
	
      <ci> f </ci>
      
    </apply>
    
    <apply>
            <times/>
    
      <cn cellml:units="dimensionless"> -4 </cn>
      
      <ci> f </ci>
    
      <ci> k_Cm </ci>
      
      <ci> O_4 </ci>
     
    </apply>
    
    <apply>
            <divide/> 
  
      <apply>
              <times/>
       
        <ci> k_L </ci>
      
	<ci> C_4 </ci>
	
      </apply>

      <apply>
              <power/>
      
        <ci> f </ci>
	  
	<cn cellml:units="dimensionless"> 4 </cn>
	
      </apply>
      
    </apply>
    
    <apply>
            <times/>
    
      <apply>
              <minus/>
     
        <ci> k_Lm </ci>
	
      </apply>
      
      <apply>
              <power/>
      
        <ci> f </ci>
	  
        <cn cellml:units="dimensionless"> 4 </cn>
	  
      </apply>
      
      <ci> O_4 </ci>
  
    </apply>
   
  </apply>
  
</apply>

</math>


</component>



<component name="razumova">
 
  
  
  <!-- these variables are defined here and only used internally -->
  <variable units="micromolar_per_millisecond_micrometre3" name="nu_SR" initial_value="2.4375"/>
  <variable units="micromolar" name="K_SR" initial_value="1"/>
  <variable units="micrometre3_per_millisecond" name="L_e" initial_value="0.00004"/>
  <variable units="micrometre3_per_millisecond" name="tau_R" initial_value="0.75"/>
  <variable units="micrometre3_per_millisecond" name="tau_SR_R" initial_value="0.75"/>
  <variable units="micrometre" name="L_x" initial_value="1.1"/>
  <variable units="micrometre" name="R_R" initial_value="0.5"/>
  
  <variable units="micrometre3" name="V_o"/>
  <variable units="micrometre3" name="V_1"/>
  <variable units="micrometre3" name="V_2"/>
  <variable units="micrometre3" name="V_SR"/>
  <variable units="micrometre3" name="V_SR1"/>
  <variable units="micrometre3" name="V_SR2"/>
  
  <variable units="per_micromolar_per_millisecond" name="k_T_on" initial_value="0.0885"/>
  <variable units="per_millisecond" name="k_T_off" initial_value="0.115"/>
  <variable units="micromolar" name="T_tot" initial_value="140"/>
  <variable units="per_micromolar_per_millisecond" name="k_P_on" initial_value="0"/>
  <variable units="per_millisecond" name="k_P_off" initial_value="0"/>
  <variable units="micromolar" name="P_tot" initial_value="1500"/>
  <variable units="per_micromolar_per_millisecond" name="k_Mg_on" initial_value="0"/>
  <variable units="per_millisecond" name="k_Mg_off" initial_value="0"/>
  <variable units="per_micromolar_per_millisecond" name="k_Cs_on" initial_value="0.000004"/>
  <variable units="per_millisecond" name="k_Cs_off" initial_value="0.005"/>
  <variable units="micromolar" name="Cs_tot" initial_value="31000"/>
  <variable units="per_micromolar_per_millisecond" name="k_CATP_on" initial_value="0.15"/>
  <variable units="per_millisecond" name="k_CATP_off" initial_value="30"/>
  <variable units="per_micromolar_per_millisecond" name="k_MATP_on" initial_value="0.0015"/>
  <variable units="per_millisecond" name="k_MATP_off" initial_value="0.15"/>
  <variable units="micrometre3_per_millisecond" name="tau_ATP" initial_value="0.375"/>
  <variable units="micrometre3_per_millisecond" name="tau_Mg" initial_value="1.5"/>
  <variable units="per_millisecond" name="k_0_on" initial_value="0"/>
  <variable units="per_millisecond" name="k_0_off" initial_value="0.15"/>
  <variable units="per_millisecond" name="k_Ca_on" initial_value="0.15"/>
  <variable units="per_millisecond" name="k_Ca_off" initial_value="0.05"/>
  <variable units="per_millisecond" name="f_o" initial_value="0.5"/>
  <variable units="per_millisecond" name="f_p" initial_value="5"/>
  <variable units="per_millisecond" name="h_o" initial_value="0.08"/>
  <variable units="per_millisecond" name="h_p" initial_value="0.06"/>
  <variable units="per_millisecond" name="g_o" initial_value="0.04"/>
  <variable units="per_millisecond" name="b_p" initial_value="0.00000394"/>
  <variable units="micrometre3_per_millisecond" name="k_p" initial_value="0.00000362"/>
  <variable units="per_milliM3_per_millisecond" name="A_p" initial_value="1"/>
  <variable units="per_milliM2_per_millisecond" name="B_p" initial_value="0.0001"/>
  <variable units="milliM2" name="PP" initial_value="6"/>
  
  <variable units="micromolar" name="T_0"/>
  <variable units="micromolar" public_interface="out" name="Ca_1" initial_value="0.1"/>
  <variable units="micromolar" name="Ca_SR1" initial_value="1500"/>
  <variable units="micromolar" name="Ca_2" initial_value="0.1"/>
  <variable units="micromolar" name="Ca_SR2" initial_value="1500"/>
  <variable units="micromolar" name="Ca_T_2" initial_value="25"/>
  <variable units="micromolar" name="Ca_P1" initial_value="615"/>
  <variable units="micromolar" name="Ca_P2" initial_value="615"/>
  <variable units="micromolar" name="Mg_P1" initial_value="811"/>
  <variable units="micromolar" name="Mg_P2" initial_value="811"/>
  <variable units="micromolar" name="Ca_Cs1" initial_value="16900"/>
  <variable units="micromolar" name="Ca_Cs2" initial_value="16900"/>
  <variable units="micromolar" name="Ca_ATP1" initial_value="0.4"/>
  <variable units="micromolar" name="Ca_ATP2" initial_value="0.4"/>
  <variable units="micromolar" name="Mg_ATP1" initial_value="7200"/>
  <variable units="micromolar" name="Mg_ATP2" initial_value="7200"/>
  <variable units="micromolar" name="ATP1" initial_value="799.6"/>
  <variable units="micromolar" name="ATP2" initial_value="799.6"/>
  <variable units="micromolar" name="Mg1" initial_value="1000"/>
  <variable units="micromolar" name="Mg2" initial_value="1000"/>
  <variable units="micromolar" name="Ca_CaT2" initial_value="3"/>
  <variable units="micromolar" name="D_0" initial_value="0.8"/>
  <variable units="micromolar" name="D_1" initial_value="1.2"/>
  <variable units="micromolar" name="D_2" initial_value="3"/>
  <variable units="micromolar" name="A_1" initial_value="0.3"/>
  <variable units="micromolar" name="A_2" initial_value="0.23"/>
  <variable units="milliM" name="P" initial_value="0.23"/>
  <variable units="milliM" name="P_SR" initial_value="0.23"/>
  <variable units="milliM" name="P_C_SR" initial_value="0.23"/>
  <variable units="micrometre3_per_millisecond" name="i2" initial_value="60"/>
 
  <!-- these variables are imported from other components -->
  
  <variable units="millisecond" public_interface="in" name="time"/> 
  <variable units="dimensionless" public_interface="in" name="O_0"/>
  <variable units="dimensionless" public_interface="in" name="O_1"/>
  <variable units="dimensionless" public_interface="in" name="O_2"/>
  <variable units="dimensionless" public_interface="in" name="O_3"/>
  <variable units="dimensionless" public_interface="in" name="O_4"/>
  
  


<!-- 

Sub (The Cannell Part)



-->


 
<math xmlns="http://www.w3.org/1998/Math/MathML"> 
<apply id="V_o_Eqn">
        <eq/>

  <ci> V_o </ci>   
  
  <apply>
          <times/>
  
    <cn cellml:units="dimensionless"> 0.95 </cn>
    
    <ci> L_x </ci>
    
    <pi/>
    
    <!--<apply><pi/>
    
    </apply> -->
    
    <apply>
            <power/>
    
      <ci> R_R </ci>
      
      <cn cellml:units="dimensionless"> 2 </cn> 
      
    </apply> 
  
  </apply>
  
</apply>

</math> 

<math xmlns="http://www.w3.org/1998/Math/MathML"> 
<apply id="V_1_Eqn">
        <eq/> 

  <ci> V_1 </ci>
    
  <apply>
          <times/>
  
    <cn cellml:units="dimensionless"> 0.01 </cn>
    
    <ci> V_o </ci>
    
  </apply>
  
</apply>

</math>
       
<math xmlns="http://www.w3.org/1998/Math/MathML"> 
<apply id="V_2_Eqn">
        <eq/> 

  <ci> V_2 </ci>
    
  <apply>
          <times/>
  
    <cn cellml:units="dimensionless"> 0.99 </cn>
    
    <ci> V_o </ci>
    
  </apply>
  
</apply>

</math>
        
<math xmlns="http://www.w3.org/1998/Math/MathML"> 
<apply id="V_SR_Eqn">
        <eq/>

  <ci> V_SR </ci>   
  
  <apply>
          <times/>
  
    <cn cellml:units="dimensionless"> 0.05 </cn>
    
    <ci> L_x </ci>
    
    <pi/>
    
    <!--<apply><pi/>
    
    </apply>--> 
    
    <apply>
            <power/>
    
      <ci> R_R </ci>
      
      <cn cellml:units="dimensionless"> 2 </cn> 
      
    </apply> 
  
  </apply>
  
</apply>

</math> 
      
<math xmlns="http://www.w3.org/1998/Math/MathML"> 
<apply id="V_SR1_Eqn">
        <eq/>

  <ci> V_SR1 </ci>
  
  <apply>
          <times/>
  
    <cn cellml:units="dimensionless"> 0.01 </cn> 
    
    <ci> V_SR </ci>
    
  </apply>
  
</apply>

</math>


<math xmlns="http://www.w3.org/1998/Math/MathML"> 
<apply id="V_SR2_Eqn">
        <eq/>

  <ci> V_SR2 </ci>
  
  <apply>
          <times/>
  
    <cn cellml:units="dimensionless"> 0.99 </cn> 
    
    <ci> V_SR </ci>
    
  </apply>
  
</apply>

</math>


<math xmlns="http://www.w3.org/1998/Math/MathML"> 
<apply id="T_0_Eqn">
        <eq/>

  <ci> T_0 </ci>
  
  <apply>
          <plus/>
  
    <ci> T_tot </ci>
    
    <apply>
            <minus/>
    
      <ci> Ca_T_2 </ci>
    
    </apply>
    
    <apply>
            <minus/>
    
      <ci> Ca_CaT2 </ci>
    
    </apply>
     
    <apply>
            <minus/>
    
      <ci> D_0 </ci>
    
    </apply> 
    
    <apply>
            <minus/>
    
      <ci> D_1 </ci>
    
    </apply> 
      
    <apply>
            <minus/> 
    
      <ci> D_2 </ci>
    
    </apply> 
      
    <apply>
            <minus/>
    
      <ci> A_1 </ci>
    
    </apply> 
    
    <apply>
            <minus/>
    
      <ci> A_2 </ci>
    
    </apply>
    
  </apply>
  
</apply>

</math>




<!--

Sub (The Razumova Part)


-->  
  


<math xmlns="http://www.w3.org/1998/Math/MathML"> 
<apply id="dCa1_Eqn">
        <eq/>

  <apply>
          <diff/>
    <bvar>
            <ci> time </ci>
          </bvar> 
    <ci> Ca_1 </ci>
  </apply> 
  <apply>
          <plus/> 
    <apply>
            <plus/> 
      <apply>
              <plus/> 
        <apply>
                <plus/> 
          <apply>
                  <minus/> 
            <apply>
                    <times/> 
               
                <apply>
                      <times/> 
                  <ci> i2 </ci>
                  <apply>
                        <plus/> 
                    <ci> O_0 </ci>
                    <ci> O_1 </ci>
                    <ci> O_2 </ci>
                    <ci> O_3 </ci>
                    <ci> O_4 </ci>
                  </apply>
                </apply>
                
              
              <apply>
                      <divide/> 
                <apply>
                        <minus/> 
                  <ci> Ca_SR1 </ci>
                  <ci> Ca_1 </ci>
                </apply>
                <ci> V_1 </ci>
              </apply>
            </apply>
            <apply>
                    <times/> 
              <ci> nu_SR </ci>
              <apply>
                      <divide/> 
                <apply>
                        <divide/> 
                  <ci> Ca_1 </ci>
                  <apply>
                          <plus/> 
                    <ci> Ca_1 </ci>
                    <ci> K_SR </ci>
                  </apply>
                </apply>
                <ci> V_1 </ci>
              </apply>
            </apply>
          </apply>
          <apply>
                  <times/> 
            <ci> L_e </ci>
            <apply>
                    <divide/> 
              <apply>
                      <minus/> 
                <ci> Ca_SR1 </ci>
                <ci> Ca_1 </ci>
             </apply>
              <ci> V_1 </ci>
            </apply>
          </apply>
        </apply>
        <apply>
                <times/> 
          <apply>
                  <minus/> 
            <ci> tau_R </ci>
          </apply>
          <apply>
                  <divide/> 
            <apply>
                    <minus/> 
              <ci> Ca_1 </ci>
              <ci> Ca_2 </ci>
            </apply>
            <ci> V_1 </ci>
          </apply>
        </apply>
      </apply>
      <apply>
              <minus/> 
        <apply>
                <plus/> 
          <apply>
                  <times/> 
            <apply>
                    <times/> 
              <ci> k_P_on </ci>
              <ci> Ca_1 </ci>
            </apply>
            <apply>
                    <plus/> 
              <apply>
                      <plus/> 
                <ci> P_tot </ci>
                <apply>
                        <minus/> 
                  <ci> Ca_P1 </ci>
                </apply>
              </apply>
              <apply>
                      <minus/> 
                <ci> Mg_P1 </ci>
              </apply>
            </apply>
          </apply>
          <apply>
                  <times/> 
            <apply>
                    <minus/> 
              <ci> k_P_off </ci>
            </apply>
            <ci> Ca_P1 </ci>
          </apply>
        </apply>
      </apply>
    </apply>
    <apply>
            <minus/> 
      <apply>
              <plus/> 
        <apply>
                <times/> 
          <apply>
                  <times/> 
            <ci> k_CATP_on </ci>
            <ci> Ca_1 </ci>
          </apply>
          <ci> ATP1 </ci> 
        </apply>
        <apply>
                <times/> 
          <apply>
                  <minus/> 
            <ci> k_CATP_off </ci>
          </apply>
          <ci> Ca_ATP1 </ci>
        </apply>
      </apply>
    </apply>
  </apply>
</apply> 
</math> 


<math xmlns="http://www.w3.org/1998/Math/MathML"> 
<apply id="dCaSR1_Eqn">
        <eq/>
  <apply>
          <diff/>
    <bvar>
            <ci> time </ci>
          </bvar>
    <ci> Ca_SR1 </ci> 
  </apply>
  <apply>
          <plus/> 
    <apply>
            <plus/> 
      <apply>
              <plus/> 
        <apply>
                <plus/> 
          <apply>
                  <times/> 
            <apply>
                    <minus/> 
               
                <apply>
                      <times/> 
                  <ci> i2 </ci>
                  <apply>
                        <plus/> 
                    <ci> O_0 </ci>
                    <ci> O_1 </ci>
                    <ci> O_2 </ci>
                    <ci> O_3 </ci>
                    <ci> O_4 </ci>
                  </apply>
                </apply>
                
              
            </apply>
            <apply>
                    <divide/> 
              <apply>
                      <minus/> 
                <ci> Ca_SR1 </ci>
                <ci> Ca_1 </ci>
              </apply>
              <ci> V_SR1 </ci>
            </apply>
          </apply>
          <apply>
                  <times/> 
            <ci> nu_SR </ci>
            <apply>
                    <divide/> 
              <apply>
                      <divide/> 
                <ci> Ca_1 </ci>
                <apply>
                        <plus/> 
                  <ci> Ca_1 </ci>
                  <ci> K_SR </ci>
                </apply>
              </apply>
              <ci> V_SR1 </ci>
            </apply>
          </apply>
        </apply>
        <apply>
                <times/> 
          <apply>
                  <minus/> 
            <ci> L_e </ci>
          </apply>
          <apply>
                  <divide/> 
            <apply>
                    <minus/> 
              <ci> Ca_SR1 </ci>
              <ci> Ca_1 </ci>
            </apply>
            <ci> V_SR1 </ci>
          </apply>
        </apply>
      </apply>
      <apply>
              <times/> 
        <apply>
                <minus/> 
          <ci> tau_SR_R </ci>
        </apply>
        <apply>
                <divide/> 
          <apply>
                  <minus/> 
            <ci> Ca_SR1 </ci>
            <ci> Ca_SR2 </ci>
          </apply>
          <ci> V_SR1 </ci>
        </apply>
      </apply>
    </apply>
    <apply>
            <minus/> 
      <apply>
              <plus/> 
        <apply>
                <times/> 
          <apply>
                  <times/> 
            <ci> k_Cs_on </ci>
            <ci> Ca_SR1 </ci>
          </apply>
          <apply>
                  <minus/> 
            <ci> Cs_tot </ci>
            <ci> Ca_Cs1 </ci>
          </apply>
        </apply>
        <apply>
                <times/> 
          <apply>
                  <minus/> 
            <ci> k_Cs_off </ci>
          </apply>
          <ci> Ca_Cs1 </ci>
        </apply>
      </apply>
    </apply>
  </apply>
</apply> 
</math> 

<math xmlns="http://www.w3.org/1998/Math/MathML"> 
<apply id="dCa2_Eqn">
        <eq/> 
  <apply>
          <diff/>
    <bvar>
            <ci> time </ci>
          </bvar>
    <ci> Ca_2 </ci>
  </apply> 
  <apply>
          <plus/> 
    <apply>
            <plus/> 
      <apply>
              <plus/> 
        <apply>
                <plus/> 
          <apply>
                  <plus/> 
            <apply>
                    <times/> 
              <apply>
                      <minus/> 
                <ci> nu_SR </ci>
              </apply>
              <apply>
                      <divide/> 
                <apply>
                        <divide/> 
                  <ci> Ca_2 </ci>
                  <apply>
                          <plus/> 
                    <ci> Ca_2 </ci>
                    <ci> K_SR </ci>
                  </apply>
                </apply>
                <ci> V_2 </ci>
              </apply>
            </apply>
            <apply>
                    <times/> 
              <ci> L_e </ci>
              <apply>
                      <divide/> 
                <apply>
                        <plus/> 
                  <ci> Ca_SR2 </ci>
                  <apply>
                          <minus/> 
                    <ci> Ca_2 </ci>
                  </apply>
                </apply>
                <ci> V_2 </ci>
              </apply>
            </apply>
          </apply>
          <apply>
                  <times/> 
            <ci> tau_R </ci>
            <apply>
                    <divide/> 
              <apply>
                      <minus/> 
                <ci> Ca_1 </ci>
                <ci> Ca_2 </ci>
              </apply>
              <ci> V_2 </ci>
            </apply>
          </apply>
        </apply>
        <apply>
                <minus/> 
          <apply>
                  <plus/> 
            <apply>
                    <plus/> 
              <apply>
                      <plus/> 
                <apply>
                        <plus/> 
                  <apply>
                          <plus/> 
                    <apply>
                            <plus/> 
                      <apply>
                              <plus/> 
                        <apply>
                                <times/> 
                          <ci> k_T_on </ci>
                          <ci> Ca_2 </ci>
                          <ci> T_0 </ci> 
                        </apply>
                        <apply>
                                <times/> 
                          <apply>
                                  <minus/> 
                            <ci> k_T_off </ci>
                          </apply>
                          <ci> Ca_T_2 </ci>
                        </apply>
                      </apply>
                      <apply>
                              <times/> 
                        <ci> k_T_on </ci>
                        <ci> Ca_2 </ci>
                        <ci> Ca_T_2 </ci>
                      </apply>
                    </apply>
                    <apply>
                            <times/> 
                      <apply>
                              <minus/> 
                        <ci> k_T_off </ci>
                      </apply>
                      <ci> Ca_CaT2 </ci>
                    </apply>
                  </apply>
                  <apply>
                          <times/> 
                    <ci> k_T_on </ci>
                    <ci> Ca_2 </ci>
                    <ci> D_0 </ci>
                  </apply>
                </apply>
                <apply>
                        <times/> 
                  <apply>
                          <minus/> 
                    <ci> k_T_off </ci>
                  </apply>
                  <ci> D_1 </ci>
                </apply>
              </apply>
              <apply>
                      <times/> 
                <ci> k_T_on </ci>
                <ci> Ca_2 </ci>
                <ci> D_1 </ci>
              </apply>
            </apply>
            <apply>
                    <times/> 
              <apply>
                      <minus/> 
                <ci> k_T_off </ci>
              </apply>
              <ci> D_2 </ci>
            </apply>
          </apply>
        </apply>
      </apply>
      <apply>
              <minus/> 
        <apply>
                <plus/> 
          <apply>
                  <times/> 
            <apply>
                    <times/> 
              <ci> k_P_on </ci>
              <ci> Ca_2 </ci>
            </apply>
          <apply>
                    <plus/> 
            <ci> P_tot </ci>
            <apply>
                      <minus/> 
              <ci> Ca_P2 </ci>
            </apply>
            <apply>
                      <minus/> 
              <ci> Mg_P2 </ci>
            </apply>
          </apply>
        </apply>
        <apply>
                  <times/> 
          <apply>
                    <minus/> 
            <ci> k_P_off </ci>
          </apply>
          <ci> Ca_P2 </ci>
        </apply>
      </apply>
    </apply>
  </apply>
  <apply>
            <minus/> 
    <apply>
              <plus/> 
      <apply>
                <times/> 
        <apply>
                  <times/> 
          <ci> k_CATP_on </ci>
          <ci> Ca_2 </ci>
        </apply>
        <ci> ATP2 </ci> 
      </apply>
      <apply>
                <times/> 
        <apply>
                  <minus/> 
          <ci> k_CATP_off </ci>
        </apply>
        <ci> Ca_ATP2 </ci>
      </apply>
    </apply>
  </apply>
</apply>
</apply> 
</math> 

<math xmlns="http://www.w3.org/1998/Math/MathML"> 
<apply id="dCa_SR2_Eqn">
        <eq/>
  <apply>
          <diff/>
    <bvar>
            <ci> time </ci>
          </bvar> 
    <ci> Ca_SR2 </ci> 
  </apply>
  
  <apply>
          <minus/>
  
  <apply>
            <plus/> 
    <apply>
              <plus/> 
      <apply>
                <plus/> 
        <apply>
                  <times/> 
          <ci> nu_SR </ci>
          <apply>
                    <divide/> 
            <apply>
                      <divide/> 
              <ci> Ca_2 </ci>
              <apply>
                        <plus/> 
                <ci> Ca_2 </ci>
                <ci> K_SR </ci>
              </apply>
            </apply>
            <ci> V_SR2 </ci>
          </apply>
        </apply>
        <apply>
                  <times/> 
          <apply>
                    <minus/> 
            <ci> L_e </ci>
          </apply>
          <apply>
                    <divide/> 
            <apply>
                      <plus/> 
              <ci> Ca_SR2 </ci>
              <apply>
                        <minus/> 
                <ci> Ca_2 </ci>
              </apply>
            </apply>
            <ci> V_SR2 </ci>
          </apply>
        </apply>
      </apply>
      <apply>
                <times/> 
        <ci> tau_SR_R </ci>
        <apply>
                  <divide/> 
          <apply>
                    <plus/> 
            <ci> Ca_SR1 </ci>
            <apply>
                      <minus/> 
              <ci> Ca_SR2 </ci>
            </apply>
          </apply>
          <ci> V_SR2 </ci>
        </apply>
      </apply>
    </apply>
    <apply>
              <minus/> 
      <apply>
                <plus/> 
        <apply>
                  <times/> 
          <apply>
                    <times/> 
            <ci> k_Cs_on </ci>
            <ci> Ca_SR2 </ci>
          </apply>
          <apply>
                    <plus/> 
            <ci> Cs_tot </ci>
            <apply>
                      <minus/> 
              <ci> Ca_Cs2 </ci>
            </apply>
          </apply>
        </apply>
        <apply>
                  <times/> 
          <apply>
                    <minus/> 
            <ci> k_Cs_off </ci>
          </apply>
          <ci> Ca_Cs2 </ci>
        </apply>
      </apply>
    </apply>
  </apply>
  
  
  <apply>
            <times/>
  
    <!-- with factor to get units correct-->
    <apply>
              <divide/>
      <cn cellml:units="milliM"> 1000 </cn>
	<cn cellml:units="M"> 1 </cn>
    </apply>

    
    <apply>
              <minus/>
      <apply>
                <times/>
        <ci> A_p </ci>
	<apply>
                  <minus/>
	  <apply>
                    <times/>
	    <ci> P_SR </ci>
	    <apply>
                      <divide/>
	      <cn cellml:units="M"> 0.001 </cn>
		<cn cellml:units="milliM"> 1 </cn>
	    </apply>
	    <ci> Ca_SR2 </ci>
	  </apply>
	  <ci> PP </ci>
	</apply>
	
	<piecewise>
  
          <piece>
	  
	    <cn cellml:units="dimensionless"> 1 </cn>
            
            <apply>
                      <gt/>
	    
              <apply>
                        <minus/>
	        <apply>
                          <times/>
	          <ci> P_SR </ci>
                <apply>
                            <divide/>
	            <cn cellml:units="M"> 0.001 </cn>
			<cn cellml:units="milliM"> 1 </cn>
		    </apply>
	          <ci> Ca_SR2 </ci>
	        </apply>
	        <ci> PP </ci>
              </apply>
	      
	      <cn cellml:units="milliM2"> 0 </cn>
	      
	    </apply>
	 
	  </piece>
      
          <otherwise>
            
            <cn cellml:units="dimensionless"> 0 </cn>
                
          </otherwise>
        
        </piecewise>
	<apply>
                  <divide/>	
	  <cn cellml:units="M"> 0.001 </cn>
	  <cn cellml:units="milliM"> 1 </cn>
      </apply>
	<ci> P_SR </ci>
	<ci> Ca_SR2 </ci>
      </apply>
      <apply>
                <times/>
        
      	
        <ci> B_p </ci>
	<ci> P_C_SR </ci>
	<apply>
                  <minus/>
	  <ci> PP </ci>
	  <apply>
                    <times/>
	    <ci> P_SR </ci>
          <apply>
                      <divide/>
	      <cn cellml:units="M"> 0.001 </cn>
            <cn cellml:units="milliM"> 1 </cn>
          </apply>
	    <ci> Ca_SR2 </ci>
	  </apply>
	</apply>
	
	
	<piecewise>
  
          <piece>
	  
	    <cn cellml:units="dimensionless"> 1 </cn>
            
            <apply>
                      <gt/>
	    
              <apply>
                        <minus/>
	        <ci> PP </ci>
	        <apply>
                          <times/>
	          <ci> P_SR </ci>
		    <apply>
                            <divide/>
	            <cn cellml:units="M"> 0.001 </cn>
			<cn cellml:units="milliM"> 1 </cn>
		    </apply>
	          <ci> Ca_SR2 </ci>
	        </apply>
	      </apply>
	      
	      <cn cellml:units="milliM2"> 0 </cn>
	      
	    </apply>
	 
	  </piece>
      
          <otherwise>
            
            <cn cellml:units="dimensionless"> 0 </cn>
                
          </otherwise>
        
        </piecewise>
	
	
	
	
	
	
      </apply>
    </apply>
  </apply>
  
  
  </apply>
  
</apply>
 
</math>

<math xmlns="http://www.w3.org/1998/Math/MathML"> 
<apply id="dCa_T_2_Eqn">
        <eq/>

  <apply>
          <diff/>
    <bvar>
            <ci> time </ci>
          </bvar>
    <ci> Ca_T_2 </ci>
  </apply> 
  <apply>
          <plus/> 
 <apply>
            <plus/> 
 <apply>
              <plus/> 
 <apply>
                <plus/> 
 <apply>
                  <plus/> 
 <apply>
                    <times/> 
 <apply>
                      <times/> 
 <ci> k_T_on </ci>
<ci> Ca_2 </ci>
  </apply>
<ci> T_0 </ci> 
  </apply>
<apply>
                    <times/> 
 <apply>
                      <minus/> 
 <ci> k_T_off </ci>
</apply>
<ci> Ca_T_2 </ci>
  </apply>
  </apply>
<apply>
                  <times/> 
 <apply>
                    <times/> 
 <apply>
                      <minus/> 
 <ci> k_T_on </ci>
</apply>
<ci> Ca_2 </ci>
  </apply>
<ci> Ca_T_2 </ci>
  </apply>
  </apply>
<apply>
                <times/> 
 <ci> k_T_off </ci>
<ci> Ca_CaT2 </ci>
  </apply>
  </apply>
<apply>
              <times/> 
 <apply>
                <minus/> 
 <ci> k_0_on </ci>
</apply>
<ci> Ca_T_2 </ci>
  </apply>
  </apply>
<apply>
            <times/> 
 <ci> k_0_off </ci>
<ci> D_1 </ci>
  </apply>
  </apply>
 </apply> 
 </math> 
 
 
 
<math xmlns="http://www.w3.org/1998/Math/MathML"> 
<apply id="dCa_P1_Eqn">
        <eq/>
  <apply>
          <diff/>
    <bvar>
            <ci> time </ci>
          </bvar> 
  <ci> Ca_P1 </ci>
  </apply> 
  <apply>
          <plus/> 
 <apply>
            <times/> 
 <apply>
              <times/> 
 <ci> k_P_on </ci>
<ci> Ca_1 </ci>
  </apply>
<apply>
              <plus/> 
 <apply>
                <plus/> 
 <ci> P_tot </ci>
<apply>
                  <minus/> 
 <ci> Ca_P1 </ci>
</apply>
  </apply>
<apply>
                <minus/> 
 <ci> Mg_P1 </ci>
</apply>
  </apply>
  </apply>
<apply>
            <times/> 
 <apply>
              <minus/> 
 <ci> k_P_off </ci>
</apply>
<ci> Ca_P1 </ci>
  </apply>
  </apply>
 </apply> 
 </math> 

<math xmlns="http://www.w3.org/1998/Math/MathML"> 
<apply id="dCa_P2_Eqn">
        <eq/>
<apply>
          <diff/>
<bvar>
            <ci> time </ci>
          </bvar> 
 <ci> Ca_P2 </ci> 
 </apply>
  <apply>
          <plus/> 
 <apply>
            <times/> 
 <apply>
              <times/> 
 <ci> k_P_on </ci>
<ci> Ca_2 </ci>
  </apply>
<apply>
              <plus/> 
 <apply>
                <plus/> 
 <ci> P_tot </ci>
<apply>
                  <minus/> 
 <ci> Ca_P2 </ci>
</apply>
  </apply>
<apply>
                <minus/> 
 <ci> Mg_P2 </ci>
</apply>
  </apply>
  </apply>
<apply>
            <times/> 
 <apply>
              <minus/> 
 <ci> k_P_off </ci>
</apply>
<ci> Ca_P2 </ci>
  </apply>
  </apply>
 </apply> 
 </math> 



<math xmlns="http://www.w3.org/1998/Math/MathML"> 
<apply id="dMg_P1_Eqn">
        <eq/>
  <apply>
          <diff/>
    <bvar>
            <ci> time </ci>
          </bvar>
    <ci> Mg_P1 </ci>
  </apply> 
  <apply>
          <plus/> 
    <apply>
            <times/> 
      <apply>
              <times/> 
        <ci> k_Mg_on </ci>
        <apply>
                <plus/>
          <ci> P_tot </ci>
          <apply>
                  <minus/>
            <ci> Ca_P1 </ci>
          </apply>
          <apply>
                  <minus/>
            <ci> Mg_P1 </ci>
          </apply>
        </apply>
      </apply>
      <ci> Mg1 </ci>
    </apply>
    <apply>
            <times/> 
      <apply>
              <minus/> 
        <ci> k_Mg_off </ci>
      </apply>
      <ci> Mg_P1 </ci>
    </apply>
  </apply>
</apply> 
</math> 

<math xmlns="http://www.w3.org/1998/Math/MathML"> 
<apply id="dMP2_Eqn">
        <eq/> 
  <apply>
          <diff/>
    <bvar>
            <ci> time </ci>
          </bvar>
 <ci> Mg_P2 </ci>
 </apply> 
  <apply>
          <plus/> 
 <apply>
            <times/> 
 <apply>
              <times/> 
 <ci> k_Mg_on </ci>
 <apply>
                <plus/>
          <ci> P_tot </ci>
          <apply>
                  <minus/>
            <ci> Ca_P2 </ci>
          </apply>
          <apply>
                  <minus/>
            <ci> Mg_P2 </ci>
          </apply>
        </apply>
  </apply>
<ci> Mg2 </ci>
  </apply>
<apply>
            <times/> 
 <apply>
              <minus/> 
 <ci> k_Mg_off </ci>
</apply>
<ci> Mg_P2 </ci>
  </apply>
  </apply>
 </apply> 
 </math> 



<math xmlns="http://www.w3.org/1998/Math/MathML"> 
<apply id="dCa_Cs1_Eqn">
        <eq/> 
 <apply>
          <diff/>
    <bvar>
            <ci> time </ci>
          </bvar>
 <ci> Ca_Cs1 </ci>
 </apply> 
  <apply>
          <plus/> 
 <apply>
            <times/> 
 <apply>
              <times/> 
 <ci> k_Cs_on </ci>
<ci> Ca_SR1 </ci>
  </apply>
<apply>
              <plus/> 
 <ci> Cs_tot </ci>
<apply>
                <minus/> 
 <ci> Ca_Cs1 </ci>
</apply>
  </apply>
  </apply>
<apply>
            <times/> 
 <apply>
              <minus/> 
 <ci> k_Cs_off </ci>
</apply>
<ci> Ca_Cs1 </ci>
  </apply>
  </apply>
 </apply> 
 </math> 

<math xmlns="http://www.w3.org/1998/Math/MathML"> 
<apply id="dCs_Cs2_Eqn">
        <eq/> 
 <apply>
          <diff/>
    <bvar>
            <ci> time </ci>
          </bvar>
 <ci> Ca_Cs2 </ci>
 </apply> 
  <apply>
          <plus/> 
 <apply>
            <times/> 
 <apply>
              <times/> 
 <ci> k_Cs_on </ci>
<ci> Ca_SR2 </ci>
  </apply>
<apply>
              <plus/> 
 <ci> Cs_tot </ci>
<apply>
                <minus/> 
 <ci> Ca_Cs2 </ci>
</apply>
  </apply>
  </apply>
<apply>
            <times/> 
 <apply>
              <minus/> 
 <ci> k_Cs_off </ci>
</apply>
<ci> Ca_Cs2 </ci>
  </apply>
  </apply>
 </apply> 
 </math> 



<math xmlns="http://www.w3.org/1998/Math/MathML"> 
<apply id="dCa_ATP1_Eqn">
        <eq/> 
 <apply>
          <diff/>
    <bvar>
            <ci> time </ci>
          </bvar>
 <ci> Ca_ATP1 </ci> 
 </apply>
  <apply>
          <plus/> 
 <apply>
            <plus/> 
 <apply>
              <times/> 
 <apply>
                <times/> 
 <ci> k_CATP_on </ci>
<ci> Ca_1 </ci>
  </apply>
<ci> ATP1 </ci> 
  </apply>
<apply>
              <times/> 
 <apply>
                <minus/> 
 <ci> k_CATP_off </ci>
</apply>
<ci> Ca_ATP1 </ci>
  </apply>
  </apply>
<apply>
            <times/> 
 <apply>
              <minus/> 
 <ci> tau_ATP </ci>
</apply>
<apply>
              <divide/> 
 <apply>
                <plus/> 
 <ci> Ca_ATP1 </ci>
<apply>
                  <minus/> 
 <ci> Ca_ATP2 </ci>
</apply>
  </apply>
<ci> V_1 </ci>
  </apply>
  </apply>
  </apply>
 </apply> 
 </math> 

<math xmlns="http://www.w3.org/1998/Math/MathML"> 
<apply id="dCa_ATP2_Eqn">
        <eq/>
<apply>
          <diff/>
    <bvar>
            <ci> time </ci>
          </bvar> 
 <ci> Ca_ATP2 </ci> 
 </apply>
  <apply>
          <plus/> 
 <apply>
            <plus/> 
 <apply>
              <times/> 
 <apply>
                <times/> 
 <ci> k_CATP_on </ci>
<ci> Ca_2 </ci>
  </apply>
<ci> ATP2 </ci> 
  </apply>
<apply>
              <times/> 
 <apply>
                <minus/> 
 <ci> k_CATP_off </ci>
</apply>
<ci> Ca_ATP2 </ci>
  </apply>
  </apply>
<apply>
            <times/> 
 <ci> tau_ATP </ci>
<apply>
              <divide/> 
 <apply>
                <plus/> 
 <ci> Ca_ATP1 </ci>
<apply>
                  <minus/> 
 <ci> Ca_ATP2 </ci>
</apply>
  </apply>
<ci> V_2 </ci>
  </apply>
  </apply>
  </apply>
 </apply> 
 </math> 


<math xmlns="http://www.w3.org/1998/Math/MathML"> 
<apply id="dMg_ATP1_Eqn">
        <eq/> 
 <apply>
          <diff/>
    <bvar>
            <ci> time </ci>
          </bvar> 
 <ci> Mg_ATP1 </ci> 
  </apply>
  <apply>
          <plus/> 
 <apply>
            <plus/> 
 <apply>
              <times/> 
 <apply>
                <times/> 
 <ci> k_MATP_on </ci>
<ci> Mg1 </ci>
  </apply>
<ci> ATP1 </ci> 
  </apply>
<apply>
              <times/> 
 <apply>
                <minus/> 
 <ci> k_MATP_off </ci>
</apply>
<ci> Mg_ATP1 </ci>
  </apply>
  </apply>
<apply>
            <times/> 
 <apply>
              <minus/> 
 <ci> tau_ATP </ci>
</apply>
<apply>
              <divide/> 
 <apply>
                <plus/> 
 <ci> Mg_ATP1 </ci>
<apply>
                  <minus/> 
 <ci> Mg_ATP2 </ci>
</apply>
  </apply>
<ci> V_1 </ci>
  </apply>
  </apply>
  </apply>
 </apply> 
 </math> 

<math xmlns="http://www.w3.org/1998/Math/MathML"> 
<apply id="dMg_ATP2_Eqn">
        <eq/> 
 <apply>
          <diff/>
    <bvar>
            <ci> time </ci>
          </bvar> 
 <ci> Mg_ATP2 </ci> 
  </apply>
  <apply>
          <plus/> 
 <apply>
            <plus/> 
 <apply>
              <times/> 
 <apply>
                <times/> 
 <ci> k_MATP_on </ci>
<ci> Mg2 </ci>
  </apply>
<ci> ATP2 </ci> 
  </apply>
<apply>
              <times/> 
 <apply>
                <minus/> 
 <ci> k_MATP_off </ci>
</apply>
<ci> Mg_ATP2 </ci>
  </apply>
  </apply>
<apply>
            <times/> 
 <ci> tau_ATP </ci>
<apply>
              <divide/> 
 <apply>
                <plus/> 
 <ci> Mg_ATP1 </ci>
<apply>
                  <minus/> 
 <ci> Mg_ATP2 </ci>
</apply>
  </apply>
<ci> V_2 </ci>
  </apply>
  </apply>
  </apply>
 </apply> 
 </math> 


<math xmlns="http://www.w3.org/1998/Math/MathML"> 
<apply id="dATP1_Eqn">
        <eq/> 
 <apply>
          <diff/>
    <bvar>
            <ci> time </ci>
          </bvar> 
 <ci> ATP1 </ci> 
  </apply>
  <apply>
          <plus/> 
 <apply>
            <plus/> 
 <apply>
              <minus/> 
 <apply>
                <plus/> 
 <apply>
                  <times/> 
 <apply>
                    <times/> 
 <ci> k_CATP_on </ci>
<ci> Ca_1 </ci>
  </apply>
<ci> ATP1 </ci> 
  </apply>
<apply>
                  <times/> 
 <apply>
                    <minus/> 
 <ci> k_CATP_off </ci>
</apply>
<ci> Ca_ATP1 </ci>
  </apply>
  </apply>
</apply>
<apply>
              <minus/> 
 <apply>
                <plus/> 
 <apply>
                  <times/> 
 <apply>
                    <times/> 
 <ci> k_MATP_on </ci>
<ci> Mg1 </ci>
  </apply>
<ci> ATP1 </ci> 
  </apply>
<apply>
                  <times/> 
 <apply>
                    <minus/> 
 <ci> k_MATP_off </ci>
</apply>
<ci> Mg_ATP1 </ci>
  </apply>
  </apply>
</apply>
  </apply>
<apply>
            <times/> 
 <apply>
              <minus/> 
 <ci> tau_ATP </ci>
</apply>
<apply>
              <divide/> 
 <apply>
                <plus/> 
 <ci> ATP1 </ci>
<apply>
                  <minus/> 
 <ci> ATP2 </ci>
</apply>
  </apply>
<ci> V_1 </ci>
  </apply>
  </apply>
  </apply>
 </apply> 
 </math> 

<math xmlns="http://www.w3.org/1998/Math/MathML"> 
<apply id="dATP2_Eqn">
        <eq/> 
  <apply>
          <diff/>
    <bvar>
            <ci> time </ci>
          </bvar> 
 <ci> ATP2 </ci> 
  </apply>
  <apply>
          <plus/> 
 <apply>
            <plus/> 
 <apply>
              <minus/> 
 <apply>
                <plus/> 
 <apply>
                  <times/> 
 <apply>
                    <times/> 
 <ci> k_CATP_on </ci>
<ci> Ca_2 </ci>
  </apply>
<ci> ATP2 </ci> 
  </apply>
<apply>
                  <times/> 
 <apply>
                    <minus/> 
 <ci> k_CATP_off </ci>
</apply>
<ci> Ca_ATP2 </ci>
  </apply>
  </apply>
</apply>
<apply>
              <minus/> 
 <apply>
                <plus/> 
 <apply>
                  <times/> 
 <apply>
                    <times/> 
 <ci> k_MATP_on </ci>
<ci> Mg2 </ci>
  </apply>
<ci> ATP2 </ci> 
  </apply>
<apply>
                  <times/> 
 <apply>
                    <minus/> 
 <ci> k_MATP_off </ci>
</apply>
<ci> Mg_ATP2 </ci>
  </apply>
  </apply>
</apply>
  </apply>
<apply>
            <times/> 
 <ci> tau_ATP </ci>
<apply>
              <divide/> 
 <apply>
                <plus/> 
 <ci> ATP1 </ci>
<apply>
                  <minus/> 
 <ci> ATP2 </ci>
</apply>
  </apply>
<ci> V_2 </ci>
  </apply>
  </apply>
  </apply>
 </apply> 
 </math> 


<math xmlns="http://www.w3.org/1998/Math/MathML"> 
<apply id="dMg1_Eqn">
        <eq/> 
 <apply>
          <diff/>
    <bvar>
            <ci> time </ci>
          </bvar> 
 <ci> Mg1 </ci> 
  </apply>
  <apply>
          <plus/> 
 <apply>
            <plus/> 
 <apply>
              <minus/> 
 <apply>
                <plus/> 
 <apply>
                  <times/> 
 <apply>
                    <times/> 
 <ci> k_Mg_on </ci>
 <apply>
                      <plus/>
          <ci> P_tot </ci>
          <apply>
                        <minus/>
            <ci> Ca_P1 </ci>
          </apply>
          <apply>
                        <minus/>
            <ci> Mg_P1 </ci>
          </apply>
        </apply>


  </apply>
<ci> Mg1 </ci>
  </apply>
<apply>
                  <times/> 
 <apply>
                    <minus/> 
 <ci> k_Mg_off </ci>
</apply>
<ci> Mg_P1 </ci>
  </apply>
  </apply>
</apply>
<apply>
              <minus/> 
 <apply>
                <plus/> 
 <apply>
                  <times/> 
 <apply>
                    <times/> 
 <ci> k_MATP_on </ci>
<ci> Mg1 </ci>
  </apply>
<ci> ATP1 </ci> 
  </apply>
<apply>
                  <times/> 
 <apply>
                    <minus/> 
 <ci> k_MATP_off </ci>
</apply>
<ci> Mg_ATP1 </ci>
  </apply>
  </apply>
</apply>
  </apply>
<apply>
            <times/> 
 <apply>
              <minus/> 
 <ci> tau_Mg </ci>
</apply>
<apply>
              <divide/> 
 <apply>
                <plus/> 
 <ci> Mg1 </ci>
<apply>
                  <minus/> 
 <ci> Mg2 </ci>
</apply>
  </apply>
<ci> V_1 </ci>
  </apply>
  </apply>
  </apply>
 </apply> 
 </math> 

<math xmlns="http://www.w3.org/1998/Math/MathML"> 
<apply id="dMg2_Eqn">
        <eq/> 
   <apply>
          <diff/>
    <bvar>
            <ci> time </ci>
          </bvar> 
  <ci> Mg2 </ci> 
  </apply>
  <apply>
          <plus/> 
 <apply>
            <plus/> 
 <apply>
              <minus/> 
 <apply>
                <plus/> 
 <apply>
                  <times/> 
 <apply>
                    <times/> 
 <ci> k_Mg_on </ci>
 <apply>
                      <plus/>
          <ci> P_tot </ci>
          <apply>
                        <minus/>
            <ci> Ca_P2 </ci>
          </apply>
          <apply>
                        <minus/>
            <ci> Mg_P2 </ci>
          </apply>
        </apply>

  </apply>
<ci> Mg2 </ci>
  </apply>
<apply>
                  <times/> 
 <apply>
                    <minus/> 
 <ci> k_Mg_off </ci>
</apply>
<ci> Mg_P2 </ci>
  </apply>
  </apply>
</apply>
<apply>
              <minus/> 
 <apply>
                <plus/> 
 <apply>
                  <times/> 
 <apply>
                    <times/> 
 <ci> k_MATP_on </ci>
<ci> Mg2 </ci>
  </apply>
<ci> ATP2 </ci> 
  </apply>
<apply>
                  <times/> 
 <apply>
                    <minus/> 
 <ci> k_MATP_off </ci>
</apply>
<ci> Mg_ATP2 </ci>
  </apply>
  </apply>
</apply>
  </apply>
<apply>
            <times/> 
 <ci> tau_Mg </ci>
<apply>
              <divide/> 
 <apply>
                <plus/> 
 <ci> Mg1 </ci>
<apply>
                  <minus/> 
 <ci> Mg2 </ci>
</apply>
  </apply>
<ci> V_2 </ci>
  </apply>
  </apply>
  </apply>
 </apply> 
 </math> 


<math xmlns="http://www.w3.org/1998/Math/MathML"> 
<apply id="dCa_CaT2_Eqn">
        <eq/> 
 <apply>
          <diff/>
    <bvar>
            <ci> time </ci>
          </bvar> 
 <ci> Ca_CaT2 </ci> 
  </apply>
  <apply>
          <plus/> 
 <apply>
            <plus/> 
 <apply>
              <plus/> 
 <apply>
                <times/> 
 <apply>
                  <times/> 
 <ci> k_T_on </ci>
<ci> Ca_2 </ci>
  </apply>
<ci> Ca_T_2 </ci>
  </apply>
<apply>
                <times/> 
 <apply>
                  <minus/> 
 <ci> k_T_off </ci>
</apply>
<ci> Ca_CaT2 </ci>
  </apply>
  </apply>
<apply>
              <times/> 
 <apply>
                <minus/> 
 <ci> k_Ca_on </ci>
</apply>
<ci> Ca_CaT2 </ci>
  </apply>
  </apply>
<apply>
            <times/> 
 <ci> k_Ca_off </ci>
<ci> D_2 </ci>
  </apply>
  </apply>
 </apply> 
 </math> 


<math xmlns="http://www.w3.org/1998/Math/MathML"> 
<apply id="dD_0_Eqn">
        <eq/> 
 <apply>
          <diff/>
    <bvar>
            <ci> time </ci>
          </bvar> 
 <ci> D_0 </ci> 
  </apply>
  <apply>
          <plus/> 
 <apply>
            <plus/> 
 <apply>
              <plus/> 
 <apply>
                <times/> 
 <apply>
                  <times/> 
 <apply>
                    <minus/> 
 <ci> k_T_on </ci>
</apply>
<ci> Ca_2 </ci>
  </apply>
<ci> D_0 </ci>
  </apply>
<apply>
                <times/> 
 <ci> k_T_off </ci>
<ci> D_1 </ci>
  </apply>
  </apply>
<apply>
              <times/> 
 <ci> k_0_on </ci>
<ci> T_0 </ci>
  </apply>
  </apply>
<apply>
            <times/> 
 <apply>
              <minus/> 
 <ci> k_0_off </ci>
</apply>
<ci> D_0 </ci>
  </apply>
  </apply>
 </apply> 
 </math> 

<math xmlns="http://www.w3.org/1998/Math/MathML"> 
<apply id="dD_1_Eqn">
        <eq/> 
 
  <apply>
          <diff/>
    <bvar>
            <ci> time </ci>
          </bvar>
    <ci> D_1 </ci> 
  </apply> 
  <apply>
          <plus/> 
    <apply>
            <plus/> 
      <apply>
              <plus/> 
        <apply>
                <plus/> 
          <apply>
                  <plus/> 
            <apply>
                    <times/> 
              <ci> k_T_on </ci>
              <ci> Ca_2 </ci>
              <ci> D_0 </ci>
            </apply>
          <apply>
                    <times/> 
             <apply>
                      <minus/> 
               <ci> k_T_off </ci>
             </apply>
             <ci> D_1 </ci>
           </apply>
         </apply>
         <apply>
                  <times/> 
           <ci> k_0_on </ci>
           <ci> Ca_T_2 </ci>
         </apply>
       </apply>
       <apply>
                <times/> 
         <apply>
                  <minus/> 
           <ci> k_0_off </ci>
         </apply>
         <ci> D_1 </ci>
       </apply>
     </apply>
     <apply>
              <times/> 
 <apply>
                <times/> 
 <apply>
                  <minus/> 
 <ci> k_T_on </ci>
</apply>
<ci> Ca_2 </ci>
  </apply>
<ci> D_1 </ci>
  </apply>
  </apply>
<apply>
            <times/> 
 <ci> k_T_off </ci>
<ci> D_2 </ci>
  </apply>
  </apply>
 </apply> 
 </math> 



<math xmlns="http://www.w3.org/1998/Math/MathML"> 
<apply id="dD_2_Eqn">
        <eq/> 
  <apply>
          <diff/>
    <bvar>
            <ci> time </ci>
          </bvar>
    <ci> D_2 </ci> 
  </apply>
  <apply>
          <plus/> 
    <apply>
            <plus/> 
      <apply>
              <plus/> 
        <apply>
                <plus/> 
          <apply>
                  <plus/> 
            <apply>
                    <plus/> 
              <apply>
                      <times/> 
                <ci> k_T_on </ci>
                <ci> Ca_2 </ci>
                <ci> D_1 </ci>
              </apply>
<apply>
                      <times/> 
 <apply>
                        <minus/> 
 <ci> k_T_off </ci>
</apply>
<ci> D_2 </ci>
  </apply>
  </apply>
<apply>
                    <times/> 
 <ci> k_Ca_on </ci>
<ci> Ca_CaT2 </ci>
  </apply>
  </apply>
<apply>
                  <times/> 
 <apply>
                    <minus/> 
 <ci> k_Ca_off </ci>
</apply>
<ci> D_2 </ci>
  </apply>
  </apply>
<apply>
                <times/> 
 <apply>
                  <minus/> 
 <ci> f_o </ci>
</apply>
<ci> D_2 </ci>
  </apply>
  </apply>
<apply>
              <times/> 
 <ci> f_p </ci>
<ci> A_1 </ci>
  </apply>
  </apply>
<apply>
            <times/> 
 <ci> g_o </ci>
<ci> A_2 </ci>
  </apply>
  </apply>
 </apply> 
 </math> 

 
 <math xmlns="http://www.w3.org/1998/Math/MathML"> 
<apply id="dA_1_Eqn">
        <eq/> 
  <apply>
          <diff/>
    <bvar>
            <ci> time </ci>
          </bvar>
    <ci> A_1 </ci> 
  </apply>
  <apply>
          <plus/> 
    <apply>
            <plus/> 
      <apply>
              <plus/> 
         
          <apply>
                <times/> 
            <ci> f_o </ci>
            <ci> D_2 </ci>
          
        </apply>
<apply>
                <times/> 
 <apply>
                  <minus/> 
 <ci> f_p </ci>
</apply>
<ci> A_1 </ci>
  </apply>
  </apply>
<apply>
              <times/> 
 <ci> h_p </ci>
<ci> A_2 </ci>
  </apply>
  </apply>
<apply>
            <times/> 
 <apply>
              <minus/> 
 <ci> h_o </ci>
</apply>
<ci> A_1 </ci>
  </apply>
  </apply>
 </apply> 
 </math> 

<math xmlns="http://www.w3.org/1998/Math/MathML"> 
<apply id="dA_2_Eqn">
        <eq/> 
 <apply>
          <diff/>
    <bvar>
            <ci> time </ci>
          </bvar>
 <ci> A_2 </ci> 
  </apply>
  <apply>
          <plus/> 
 <apply>
            <plus/> 
 <apply>
              <times/> 
 <apply>
                <minus/> 
 <ci> h_p </ci>
</apply>
<ci> A_2 </ci>
  </apply>
<apply>
              <times/> 
 <ci> h_o </ci>
<ci> A_1 </ci>
  </apply>
  </apply>
<apply>
            <times/> 
 <apply>
              <minus/> 
 <ci> g_o </ci>
</apply>
<ci> A_2 </ci>
  </apply>
  </apply>
 </apply> 
 </math> 
 



<math xmlns="http://www.w3.org/1998/Math/MathML"> 
<apply id="dP_Eqn">
        <eq/> 
  <apply>
          <diff/>
    <bvar>
            <ci> time </ci>
          </bvar>
    <ci> P </ci> 
  </apply>
  <apply>
          <plus/>
    <apply>
            <times/>
	<apply>
              <divide/>
        <cn cellml:units="M"> 0.001 </cn>
	  <cn cellml:units="milliM"> 1 </cn>
	</apply>
      <apply>
              <minus/>
        <apply>
                <times/>
	  <ci> h_o </ci>
	  <ci> A_1 </ci>
	</apply>
	<apply>
                <times/>
	  <ci> h_p </ci>
	  <ci> A_2 </ci>
	</apply>
      </apply>
    </apply>
    <apply>
            <times/>
      <cn cellml:units="dimensionless"> -1 </cn>
      <ci> b_p </ci>
      <ci> P </ci>
    </apply>
    <apply>
            <times/>
      <cn cellml:units="dimensionless"> -1 </cn>
      <ci> k_p </ci>
      <apply>
              <divide/>
        <apply>
                <minus/>
	  <ci> P </ci>
	  <ci> P_SR </ci>
	</apply>
	<ci> V_2 </ci>
      </apply>
    </apply>
  </apply>
</apply>
</math>
    

<math xmlns="http://www.w3.org/1998/Math/MathML"> 
<apply id="dP_SR_Eqn">
        <eq/> 
  <apply>
          <diff/>
    <bvar>
            <ci> time </ci>
          </bvar>
    <ci> P_SR </ci> 
  </apply> 
  <apply>
          <minus/>
    <apply>
            <times/>
      <ci> k_p </ci>
      <apply>
              <divide/>
        <apply>
                <minus/>
	  <ci> P </ci>
	  <ci> P_SR </ci>
	</apply>
	<ci> V_SR2 </ci>
      </apply>
    </apply>
    
    <!-- to get units correct -->
    <apply>
            <times/>
      <cn cellml:units="dimensionless"> 1 </cn>
    
    <apply>
              <minus/>
      <apply>
                <times/>
        <ci> A_p </ci>
	<apply>
                  <minus/>
	  <apply>
                    <times/>
	    <ci> P_SR </ci>
	    <apply>
                      <divide/>
	    	<cn cellml:units="M"> 0.001 </cn>
		<cn cellml:units="milliM"> 1 </cn>
	    </apply>
	    <ci> Ca_SR2 </ci>
	  </apply>
	  <ci> PP </ci>
	</apply>
	
	<piecewise>
  
          <piece>
           
            <cn cellml:units="dimensionless"> 1 </cn>
            
            <apply>
                      <gt/>
              
              <apply>
                        <minus/>
	        <apply>
                          <times/>
	          <ci> P_SR </ci>
		    <apply>
                            <divide/>
	          	<cn cellml:units="M"> 0.001 </cn>
			<cn cellml:units="milliM"> 1 </cn>
		    </apply>
	          <ci> Ca_SR2 </ci>
	        </apply>
	        <ci> PP </ci>
              </apply>
              
	      <cn cellml:units="milliM2"> 0 </cn>
            
            </apply>
          
          </piece>
      
          <otherwise>
            
            <cn cellml:units="dimensionless"> 0 </cn>
                
          </otherwise>
        
        </piecewise>
        
	
	<apply>
                  <divide/>
	  <cn cellml:units="M"> 0.001 </cn>
	  <cn cellml:units="milliM"> 1 </cn>
	</apply>
	<ci> P_SR </ci>
	<ci> Ca_SR2 </ci>
      </apply>
      <apply>
                <times/>
      
       
        <ci> B_p </ci>
	<ci> P_C_SR </ci>
	<apply>
                  <minus/>
	  <ci> PP </ci>
	  <apply>
                    <times/>
	    <ci> P_SR </ci>
	    <apply>
                      <divide/>
	    	<cn cellml:units="M"> 0.001 </cn>
		<cn cellml:units="milliM"> 1 </cn>
	    </apply>
	    <ci> Ca_SR2 </ci>
	  </apply>
	</apply>
	
	<piecewise>
  
          <piece>
           
            <cn cellml:units="dimensionless"> 1 </cn>
            
            <apply>
                      <gt/>
              
              <apply>
                        <minus/>
	        <ci> PP </ci>
	        <apply>
                          <times/>
	          <ci> P_SR </ci>
		    <apply>
                            <divide/>
	          	<cn cellml:units="M"> 0.001 </cn>
			<cn cellml:units="milliM"> 1 </cn>
		    </apply>
	          <ci> Ca_SR2 </ci>
	        </apply>
	      </apply>
              
	      <cn cellml:units="milliM2"> 0 </cn>
            
            </apply>
          
          </piece>
      
          <otherwise>
            
            <cn cellml:units="dimensionless"> 0 </cn>
                
          </otherwise>
        
        </piecewise>
	
	
	
	
      </apply>
    </apply>
    
    </apply>
    
  </apply>
</apply>
</math>

<math xmlns="http://www.w3.org/1998/Math/MathML"> 
<apply id="dP_C_SR_Eqn">
        <eq/> 
  <apply>
          <diff/>
    <bvar>
            <ci> time </ci>
          </bvar>
    <ci> P_C_SR </ci> 
  </apply>
  
  
  <!-- to get units correct -->
  <apply>
          <times/>
    <cn cellml:units="dimensionless"> 1 </cn>
    
   
  <apply>
            <minus/>
      <apply>
              <times/>
        <ci> A_p </ci>
	<apply>
                <minus/>
	  <apply>
                  <times/>
	    <ci> P_SR </ci>
	    <apply>
                    <divide/>
	    	<cn cellml:units="M"> 0.001 </cn>
		<cn cellml:units="milliM"> 1 </cn>
	    </apply>
	    <ci> Ca_SR2 </ci>
	  </apply>
	  <ci> PP </ci>
	</apply>
	
	<piecewise>
  
          <piece>
           
            <cn cellml:units="dimensionless"> 1 </cn>
            
            <apply>
                    <gt/>
              
              <apply>
                      <minus/>
	        <apply>
                        <times/>
	          <ci> P_SR </ci>
		    <apply>
                          <divide/>
	          	<cn cellml:units="M"> 0.001 </cn>
			<cn cellml:units="milliM"> 1 </cn>
		    </apply>
	          <ci> Ca_SR2 </ci>
	        </apply>
	        <ci> PP </ci>
              </apply>
              
	      <cn cellml:units="milliM2"> 0 </cn>
            
            </apply>
          
          </piece>
      
          <otherwise>
            
            <cn cellml:units="dimensionless"> 0 </cn>
                
          </otherwise>
        
        </piecewise>
	
        
	
	<apply>
                <divide/>
	  <cn cellml:units="M"> 0.001 </cn>
	  <cn cellml:units="milliM"> 1 </cn>
	</apply>
	<ci> P_SR </ci>
	<ci> Ca_SR2 </ci>
      </apply>
      <apply>
              <times/>
      
        <ci> B_p </ci>
	<ci> P_C_SR </ci>
	<apply>
                <minus/>
	  <ci> PP </ci>
	  <apply>
                  <times/>
	    <ci> P_SR </ci>
	    <apply>
                    <divide/>
	    	<cn cellml:units="M"> 0.001 </cn>
		<cn cellml:units="milliM"> 1 </cn>
	    </apply>
	    <ci> Ca_SR2 </ci>
	  </apply>
	</apply>
	
	<piecewise>
  
          <piece>
           
            <cn cellml:units="dimensionless"> 1 </cn>
            
            <apply>
                    <gt/>
              
              <apply>
                      <minus/>
	        <ci> PP </ci>
	        <apply>
                        <times/>
	          <ci> P_SR </ci>
		    <apply>
                          <divide/>
	          	<cn cellml:units="M"> 0.001 </cn>
			<cn cellml:units="milliM"> 1 </cn>
		    </apply>
	          <ci> Ca_SR2 </ci>
	        </apply>
	      </apply>
              
	      <cn cellml:units="milliM2"> 0 </cn>
            
            </apply>
          
          </piece>
      
          <otherwise>
            
            <cn cellml:units="dimensionless"> 0 </cn>
                
          </otherwise>
        
        </piecewise>
	
	
	
      </apply>
    </apply>
    
    </apply>
</apply>
</math>
    
  
        
</component>



<group>
  
  <relationship_ref relationship="containment"/>
  <component_ref component="cell">
    <component_ref component="wal_environment"/>
    <component_ref component="sternrios"/>
    <component_ref component="razumova"/>
  </component_ref>

</group>

<group>

  <relationship_ref relationship="encapsulation"/>
  
  <component_ref component="wal_environment">
  
    <component_ref component="sarco_Cl_channel"/>
    <component_ref component="sarco_IR_channel"/>
    <component_ref component="sarco_DR_channel"/>
    <component_ref component="sarco_Na_channel"/>
    <component_ref component="sarco_NaK_channel"/>
    <component_ref component="t_Cl_channel"/>
    <component_ref component="t_IR_channel"/>
    <component_ref component="t_DR_channel"/>
    <component_ref component="t_Na_channel"/>
    <component_ref component="t_NaK_channel"/>
  </component_ref>
  
</group>

<connection>

  <map_components component_2="wal_environment" component_1="cell"/>
  <map_variables variable_2="time" variable_1="time"/>
   
</connection>

<connection>

  <map_components component_2="sarco_Cl_channel" component_1="wal_environment"/>
  <map_variables variable_2="vS" variable_1="vS"/>
  <map_variables variable_2="Cl_i" variable_1="Cl_i"/>
  <map_variables variable_2="Cl_o" variable_1="Cl_o"/>
  <map_variables variable_2="V_a" variable_1="V_a"/>
  <map_variables variable_2="A_a" variable_1="A_a"/>
  <map_variables variable_2="FF" variable_1="FF"/>
  <map_variables variable_2="RR" variable_1="RR"/>
  <map_variables variable_2="TT" variable_1="TT"/>
  <map_variables variable_2="time" variable_1="time"/>
  <map_variables variable_2="g_Cl_bar" variable_1="g_Cl_bar"/>
  <map_variables variable_2="I_Cl" variable_1="I_Cl"/>
  
  
</connection>


<connection>

  <map_components component_2="sarco_IR_channel" component_1="wal_environment"/>
  <map_variables variable_2="vS" variable_1="vS"/>
  <map_variables variable_2="E_K" variable_1="E_K"/>
  <map_variables variable_2="K_e" variable_1="K_e"/>
  <map_variables variable_2="J_K" variable_1="J_K"/>
  <map_variables variable_2="del" variable_1="del"/>
  <map_variables variable_2="FF" variable_1="FF"/>
  <map_variables variable_2="RR" variable_1="RR"/>
  <map_variables variable_2="TT" variable_1="TT"/>
  <map_variables variable_2="time" variable_1="time"/>
  <map_variables variable_2="G_K" variable_1="G_K"/>
  <map_variables variable_2="K_K" variable_1="K_K"/>
  <map_variables variable_2="K_S" variable_1="K_S"/>
  <map_variables variable_2="S_i" variable_1="S_i"/>
  <map_variables variable_2="I_IR" variable_1="I_IR"/>
  
</connection>

<connection>

  <map_components component_2="sarco_DR_channel" component_1="wal_environment"/>
  <map_variables variable_2="vS" variable_1="vS"/>
  <map_variables variable_2="J_K" variable_1="J_K"/>
  <map_variables variable_2="alpha_n_bar" variable_1="alpha_n_bar"/>
  <map_variables variable_2="V_n" variable_1="V_n"/>
  <map_variables variable_2="K_alpha_n" variable_1="K_alpha_n"/>
  <map_variables variable_2="beta_n_bar" variable_1="beta_n_bar"/>
  <map_variables variable_2="K_beta_n" variable_1="K_beta_n"/>
  <map_variables variable_2="V_h_K_inf" variable_1="V_h_K_inf"/>
  <map_variables variable_2="A_h_K_inf" variable_1="A_h_K_inf"/>
  <map_variables variable_2="time" variable_1="time"/>
  <map_variables variable_2="g_K_bar" variable_1="g_K_bar"/>
  <map_variables variable_2="I_DR" variable_1="I_DR"/>
  
</connection>


<connection>

  <map_components component_2="sarco_Na_channel" component_1="wal_environment"/>
  <map_variables variable_2="vS" variable_1="vS"/>
  <map_variables variable_2="Na_i" variable_1="Na_i"/>
  <map_variables variable_2="Na_e" variable_1="Na_e"/>
  <map_variables variable_2="alpha_h_bar" variable_1="alpha_h_bar"/>
  <map_variables variable_2="V_h" variable_1="V_h"/>
  <map_variables variable_2="K_alpha_h" variable_1="K_alpha_h"/>
  <map_variables variable_2="beta_h_bar" variable_1="beta_h_bar"/>
  <map_variables variable_2="K_beta_h" variable_1="K_beta_h"/>
  <map_variables variable_2="alpha_m_bar" variable_1="alpha_m_bar"/>
  <map_variables variable_2="V_m" variable_1="V_m"/>
  <map_variables variable_2="K_alpha_m" variable_1="K_alpha_m"/>
  <map_variables variable_2="beta_m_bar" variable_1="beta_m_bar"/>
  <map_variables variable_2="K_beta_m" variable_1="K_beta_m"/>
  <map_variables variable_2="V_S_inf" variable_1="V_S_inf"/>
  <map_variables variable_2="A_S_inf" variable_1="A_S_inf"/>
  <map_variables variable_2="V_tau" variable_1="V_tau"/>
  <map_variables variable_2="FF" variable_1="FF"/>
  <map_variables variable_2="RR" variable_1="RR"/>
  <map_variables variable_2="TT" variable_1="TT"/>
  <map_variables variable_2="time" variable_1="time"/>
  <map_variables variable_2="g_Na_bar" variable_1="g_Na_bar"/>
  <map_variables variable_2="I_Na" variable_1="I_Na"/>
  
</connection>


<connection>

  <map_components component_2="sarco_NaK_channel" component_1="wal_environment"/>
  <map_variables variable_2="vS" variable_1="vS"/>
  <map_variables variable_2="Na_i" variable_1="Na_i"/>
  <map_variables variable_2="Na_e" variable_1="Na_e"/>
  <map_variables variable_2="K_e" variable_1="K_e"/>
  <map_variables variable_2="K_m_K" variable_1="K_m_K"/>
  <map_variables variable_2="K_m_Na" variable_1="K_m_Na"/>
  <map_variables variable_2="FF" variable_1="FF"/>
  <map_variables variable_2="RR" variable_1="RR"/>
  <map_variables variable_2="TT" variable_1="TT"/>
  <map_variables variable_2="time" variable_1="time"/>
  <map_variables variable_2="J_NaK_bar" variable_1="J_NaK_bar"/>
  <map_variables variable_2="I_NaK" variable_1="I_NaK"/>
  
</connection>

<connection>

  <map_components component_2="t_Cl_channel" component_1="wal_environment"/>
  <map_variables variable_2="vT" variable_1="vT"/>
  <map_variables variable_2="Cl_i_t" variable_1="Cl_i_t"/>
  <map_variables variable_2="Cl_o_t" variable_1="Cl_o_t"/>
  <map_variables variable_2="V_a" variable_1="V_a"/>
  <map_variables variable_2="A_a" variable_1="A_a"/>
  <map_variables variable_2="FF" variable_1="FF"/>
  <map_variables variable_2="RR" variable_1="RR"/>
  <map_variables variable_2="TT" variable_1="TT"/>
  <map_variables variable_2="time" variable_1="time"/>
  <map_variables variable_2="g_Cl_bar" variable_1="g_Cl_bar"/>
  <map_variables variable_2="I_Cl_t" variable_1="I_Cl_t"/>
  <map_variables variable_2="eta_Cl" variable_1="eta_Cl"/>
  
</connection>


<connection>

  <map_components component_2="t_IR_channel" component_1="wal_environment"/>
  <map_variables variable_2="vT" variable_1="vT"/>
  <map_variables variable_2="E_K_t" variable_1="E_K_t"/>
  <map_variables variable_2="K_t" variable_1="K_t"/>
  <map_variables variable_2="J_K_t" variable_1="J_K_t"/>
  <map_variables variable_2="del" variable_1="del"/>
  <map_variables variable_2="FF" variable_1="FF"/>
  <map_variables variable_2="RR" variable_1="RR"/>
  <map_variables variable_2="TT" variable_1="TT"/>
  <map_variables variable_2="time" variable_1="time"/>
  <map_variables variable_2="G_K" variable_1="G_K"/>
  <map_variables variable_2="K_K" variable_1="K_K"/>
  <map_variables variable_2="K_S" variable_1="K_S"/>
  <map_variables variable_2="S_i" variable_1="S_i"/>
  <map_variables variable_2="I_IR_t" variable_1="I_IR_t"/>
  <map_variables variable_2="eta_IR" variable_1="eta_IR"/>
  
</connection>

<connection>

  <map_components component_2="t_DR_channel" component_1="wal_environment"/>
  <map_variables variable_2="vT" variable_1="vT"/>
  <map_variables variable_2="J_K_t" variable_1="J_K_t"/>
  <map_variables variable_2="alpha_n_bar" variable_1="alpha_n_bar"/>
  <map_variables variable_2="V_n" variable_1="V_n"/>
  <map_variables variable_2="K_alpha_n" variable_1="K_alpha_n"/>
  <map_variables variable_2="beta_n_bar" variable_1="beta_n_bar"/>
  <map_variables variable_2="K_beta_n" variable_1="K_beta_n"/>
  <map_variables variable_2="V_h_K_inf" variable_1="V_h_K_inf"/>
  <map_variables variable_2="A_h_K_inf" variable_1="A_h_K_inf"/>
  <map_variables variable_2="time" variable_1="time"/>
  <map_variables variable_2="g_K_bar" variable_1="g_K_bar"/>
  <map_variables variable_2="I_DR_t" variable_1="I_DR_t"/>
  <map_variables variable_2="eta_DR" variable_1="eta_DR"/>
  
</connection>


<connection>

  <map_components component_2="t_Na_channel" component_1="wal_environment"/>
  <map_variables variable_2="vT" variable_1="vT"/>
  <map_variables variable_2="Na_i" variable_1="Na_i"/>
  <map_variables variable_2="Na_t" variable_1="Na_t"/>
  <map_variables variable_2="alpha_h_bar" variable_1="alpha_h_bar"/>
  <map_variables variable_2="V_h" variable_1="V_h"/>
  <map_variables variable_2="K_alpha_h" variable_1="K_alpha_h"/>
  <map_variables variable_2="beta_h_bar" variable_1="beta_h_bar"/>
  <map_variables variable_2="K_beta_h" variable_1="K_beta_h"/>
  <map_variables variable_2="alpha_m_bar" variable_1="alpha_m_bar"/>
  <map_variables variable_2="V_m" variable_1="V_m"/>
  <map_variables variable_2="K_alpha_m" variable_1="K_alpha_m"/>
  <map_variables variable_2="beta_m_bar" variable_1="beta_m_bar"/>
  <map_variables variable_2="K_beta_m" variable_1="K_beta_m"/>
  <map_variables variable_2="V_S_inf" variable_1="V_S_inf"/>
  <map_variables variable_2="A_S_inf" variable_1="A_S_inf"/>
  <map_variables variable_2="V_tau" variable_1="V_tau"/>
  <map_variables variable_2="FF" variable_1="FF"/>
  <map_variables variable_2="RR" variable_1="RR"/>
  <map_variables variable_2="TT" variable_1="TT"/>
  <map_variables variable_2="time" variable_1="time"/>
  <map_variables variable_2="g_Na_bar" variable_1="g_Na_bar"/>
  <map_variables variable_2="I_Na_t" variable_1="I_Na_t"/>
  <map_variables variable_2="eta_Na" variable_1="eta_Na"/>
  
</connection>


<connection>

  <map_components component_2="t_NaK_channel" component_1="wal_environment"/>
  <map_variables variable_2="vT" variable_1="vT"/>
  <map_variables variable_2="Na_i" variable_1="Na_i"/>
  <map_variables variable_2="Na_t" variable_1="Na_t"/>
  <map_variables variable_2="K_t" variable_1="K_t"/>
  <map_variables variable_2="K_m_K" variable_1="K_m_K"/>
  <map_variables variable_2="K_m_Na" variable_1="K_m_Na"/>
  <map_variables variable_2="FF" variable_1="FF"/>
  <map_variables variable_2="RR" variable_1="RR"/>
  <map_variables variable_2="TT" variable_1="TT"/>
  <map_variables variable_2="time" variable_1="time"/>
  <map_variables variable_2="J_NaK_bar" variable_1="J_NaK_bar"/>
  <map_variables variable_2="I_NaK_t" variable_1="I_NaK_t"/>
  <map_variables variable_2="eta_NaK" variable_1="eta_NaK"/>
  
</connection>

<connection>
  <map_components component_2="sternrios" component_1="cell"/>
  <map_variables variable_2="time" variable_1="time"/>
</connection>

<connection>
  <map_components component_2="sternrios" component_1="wal_environment"/>
  <map_variables variable_2="vT" variable_1="vT"/>
</connection>

<connection>
  <map_components component_2="razumova" component_1="cell"/>
  <map_variables variable_2="time" variable_1="time"/>
</connection>

<connection>
  <map_components component_2="razumova" component_1="sternrios"/>
  <map_variables variable_2="O_0" variable_1="O_0"/>
  <map_variables variable_2="O_1" variable_1="O_1"/>
  <map_variables variable_2="O_2" variable_1="O_2"/>
  <map_variables variable_2="O_3" variable_1="O_3"/>
  <map_variables variable_2="O_4" variable_1="O_4"/>
  
</connection>


<rdf:RDF>
  <rdf:Seq rdf:about="rdf:#8edcae8f-29c4-4bb0-952b-9e986673097e">
    <rdf:li rdf:resource="rdf:#c54bc728-a493-48e9-aec8-4bde6e8ceafe"/>
    <rdf:li rdf:resource="rdf:#aca58b7c-3db3-4845-bc58-868e073d2a3b"/>
    <rdf:li rdf:resource="rdf:#63d96787-7bc6-4dc6-8f68-1f481dd6873f"/>
    <rdf:li rdf:resource="rdf:#c1882186-e4b4-4198-98b3-992db03eb285"/>
  </rdf:Seq>
  <rdf:Description rdf:about="rdf:#4d626140-c5b5-4901-930e-c8928572a480">
    <dcterms:W3CDTF>2008-05-26T13:00:58+12:00</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#6cc6ac24-6032-4652-92bc-ee43edec87cd">
    <vCard:ORG rdf:resource="rdf:#9f7d3b71-3ae1-468d-82fe-b9e16afdaf16"/>
    <vCard:EMAIL rdf:resource="rdf:#5e507b7c-8add-40a2-a84a-4030125a8bd7"/>
    <vCard:N rdf:resource="rdf:#e338712e-acc1-4b1b-a7c6-7b7ddea351eb"/>
  </rdf:Description>
  <rdf:Description rdf:about="">
    <dc:title>A mathematical model of fatigue in skeletal muscle force contraction: slow-twitch muscle</dc:title>
	<dc:publisher>The University of Auckland, Auckland Bioengineering Institute</dc:publisher>
    <cmeta:comment rdf:resource="rdf:#571d3981-cf12-4435-9a1f-48a26d0fc5da"/>
    <dcterms:created rdf:resource="rdf:#9f27f24f-c56a-4fbf-bf78-f8ef1fea7e19"/>
    <dc:creator rdf:resource="rdf:#6cc6ac24-6032-4652-92bc-ee43edec87cd"/>
    <cmeta:modification rdf:resource="rdf:#6b3ac2a0-7bc5-45d3-82db-2c219d0a8b52"/>
    <cmeta:modification rdf:resource="rdf:#47f7404b-ec61-425c-8a6e-daea5eb6bfcb"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#865ab0fa-a280-4eff-9d0d-8fc518a43c61">
    <dcterms:W3CDTF>2008-06-16T09:33:31+12:00</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#9f7d3b71-3ae1-468d-82fe-b9e16afdaf16">
    <vCard:Orgname>The University of Auckland</vCard:Orgname>
    <vCard:Orgunit>Auckland Bioengineering Institute</vCard:Orgunit>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#d60d740c-86ba-4c05-8aa5-5e49d609432b">
    <vCard:Given>Paul</vCard:Given>
    <vCard:Family>Shorten</vCard:Family>
    <vCard:Other>R</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#fecc3d95-47e2-4e02-bd13-1395bc335bd3">
    <vCard:Given>Paul</vCard:Given>
    <vCard:Family>O'Callaghan</vCard:Family>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#6a395a10-250f-40c1-8413-3c336ccccc9b">
    <vCard:N rdf:resource="rdf:#b443f7bf-386c-47c9-9b04-e6e17a20a37b"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#9f27f24f-c56a-4fbf-bf78-f8ef1fea7e19">
    <dcterms:W3CDTF>2008-05-26T00:00:00+00:00</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#e589a0a5-59ae-445e-a771-307e70e42365">
    <vCard:FN>Catherine Lloyd</vCard:FN>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#b443f7bf-386c-47c9-9b04-e6e17a20a37b">
    <vCard:Given>John</vCard:Given>
    <vCard:Family>Davidson </vCard:Family>
    <vCard:Other>B</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#0cb47d69-6153-4d58-945f-ba3253f36797">
    <dc:title>Journal of Muscle Research and Cell Motility</dc:title>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#5e507b7c-8add-40a2-a84a-4030125a8bd7">
    <rdf:type rdf:resource="http://imc.org/vCard/3.0#internet"/>
    <rdf:value>j.davidson@auckland.ac.nz</rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#571d3981-cf12-4435-9a1f-48a26d0fc5da">
    <dc:creator rdf:resource="rdf:#e589a0a5-59ae-445e-a771-307e70e42365"/>
    <rdf:value>This cellML model is known to run in PCEnv (0.4), COR and JSim to recreate the published results.  The units have been checked and they are consistent in as much as they are known to be equivalent.  This particular version of the CellML model represents slow twitch muscle.</rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#9b86b023-0efd-4680-9ef3-ac349bcbe0a3">
    <vCard:Given>John</vCard:Given>
    <vCard:Family>Davidson</vCard:Family>
    <vCard:Other>B</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#c1882186-e4b4-4198-98b3-992db03eb285">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#a9476cd0-9aa4-4dcb-b34b-0d98dee82c03"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#47f7404b-ec61-425c-8a6e-daea5eb6bfcb">
    <dcterms:modified rdf:resource="rdf:#4d626140-c5b5-4901-930e-c8928572a480"/>
    <rdf:value>Removed some variables which were defined as constants but which were not used anywhere else in the model.</rdf:value>
    <cmeta:modifier rdf:resource="rdf:#e34db315-bb4c-4eda-b4dc-959698aa7396"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#316589ff-d141-41c6-b5b9-e7edf5359c82">
    <bqs:Pubmed_id>18080210</bqs:Pubmed_id>
    <bqs:JournalArticle rdf:resource="rdf:#4315812c-7ff0-44d4-ae0d-e6f112530b69"/>
  </rdf:Description>
  <rdf:Description rdf:about="#Slow_Twitch">
    <bqs:reference rdf:resource="rdf:#316589ff-d141-41c6-b5b9-e7edf5359c82"/>
<bqs:reference rdf:parseType="Resource">
  <dc:subject rdf:parseType="Resource">
    <bqs:subject_type>keyword</bqs:subject_type>
    <rdf:value>
      <rdf:Bag>
        <rdf:li>excitation-contraction coupling</rdf:li>
      </rdf:Bag>
    </rdf:value>
  </dc:subject>
</bqs:reference>
    <cmeta:comment rdf:resource="rdf:#10574c1e-e2d2-4c0d-8f81-5a74453cf83f"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#63d96787-7bc6-4dc6-8f68-1f481dd6873f">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#9b86b023-0efd-4680-9ef3-ac349bcbe0a3"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#10574c1e-e2d2-4c0d-8f81-5a74453cf83f">
    <dc:creator rdf:resource="rdf:#62797b8c-2602-4a74-b0f0-1938580d3a44"/>
    <rdf:value/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#62797b8c-2602-4a74-b0f0-1938580d3a44">
    <vCard:FN/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#4315812c-7ff0-44d4-ae0d-e6f112530b69">
    <dc:creator rdf:resource="rdf:#8edcae8f-29c4-4bb0-952b-9e986673097e"/>
    <dc:title>A mathematical model of fatigue in skeletal muscle force contraction</dc:title>
    <bqs:volume>28</bqs:volume>
    <bqs:first_page>293</bqs:first_page>
    <bqs:Journal rdf:resource="rdf:#0cb47d69-6153-4d58-945f-ba3253f36797"/>
    <dcterms:issued rdf:resource="rdf:#de31a6a3-0afc-420a-980f-9ddc941bb6be"/>
    <bqs:last_page>313</bqs:last_page>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#e34db315-bb4c-4eda-b4dc-959698aa7396">
    <vCard:N rdf:resource="rdf:#82f93635-0168-48e6-aed0-f2a91b1a8432"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#de31a6a3-0afc-420a-980f-9ddc941bb6be">
    <dcterms:W3CDTF>2007-00-00 00:00</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#a9476cd0-9aa4-4dcb-b34b-0d98dee82c03">
    <vCard:Given>Tanya</vCard:Given>
    <vCard:Family>Soboleva</vCard:Family>
    <vCard:Other>K</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#e338712e-acc1-4b1b-a7c6-7b7ddea351eb">
    <vCard:Given>John</vCard:Given>
    <vCard:Family>Davidson</vCard:Family>
    <vCard:Other>B</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#c54bc728-a493-48e9-aec8-4bde6e8ceafe">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#d60d740c-86ba-4c05-8aa5-5e49d609432b"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#6b3ac2a0-7bc5-45d3-82db-2c219d0a8b52">
    <dcterms:modified rdf:resource="rdf:#865ab0fa-a280-4eff-9d0d-8fc518a43c61"/>
    <rdf:value>Fixed units such that the model now runs in JSim without any need for automated unit correction.</rdf:value>
    <cmeta:modifier rdf:resource="rdf:#6a395a10-250f-40c1-8413-3c336ccccc9b"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#82f93635-0168-48e6-aed0-f2a91b1a8432">
    <vCard:Given>Catherine</vCard:Given>
    <vCard:Family>Lloyd</vCard:Family>
    <vCard:Other>May</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#aca58b7c-3db3-4845-bc58-868e073d2a3b">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#fecc3d95-47e2-4e02-bd13-1395bc335bd3"/>
  </rdf:Description>
</rdf:RDF>
</model>