Location: Wang, Sobie, 2008 @ 313ea109136e / wang_sobie_2008.cellml

Author:
Hanne <Hanne@hanne-nielsens-macbook.local>
Date:
2009-12-15 11:47:27+13:00
Desc:
Added images in ai and svg format, removed non pub med references
Permanent Source URI:
https://models.cellml.org/workspace/wang_sobie_2008/rawfile/313ea109136ea42097fc104189095e41e1707e32/wang_sobie_2008.cellml

<?xml version='1.0' encoding='utf-8'?>
<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#" xmlns:ns7="http://www.cellml.org/metadata/simulation/1.0#" cmeta:id="wang_model_2008" name="wang_model_2008">
<documentation xmlns="http://cellml.org/tmp-documentation">
    <article>
      <articleinfo>
        <title>Mathematical model of the neonatal mouse ventricular action potential</title>
        <author>
          <firstname>Catherine</firstname>
          <surname>Lloyd</surname>
          <affiliation>
            <shortaffil>ABI</shortaffil>
          </affiliation>
        </author>
      </articleinfo>
      <section id="sec_status">
        <title>Model Status</title>
        <para>
            This CellML model is known to run in PCEnv and COR to recreate the published results.  The units have been checked and are consistent.
          </para>
      </section>
      <sect1 id="sec_structure">
        <title>Model Structure</title>
        <para>
ABSTRACT: Therapies for heart disease are based largely on our understanding of the adult myocardium. The dramatic differences in action potential (AP) shape between neonatal and adult cardiac myocytes, however, indicate that a different set of molecular interactions in neonatal myocytes necessitates different treatment for newborns. Computational modeling is useful for synthesizing data to determine how interactions between components lead to systems-level behavior, but this technique has not been used extensively to study neonatal heart cell function. We created a mathematical model of the neonatal (day 1) mouse myocyte by modifying, on the basis of experimental data, the densities and/or formulations of ion transport mechanisms in an adult cell model. The new model reproduces the characteristic AP shape of neonatal cells, with a brief plateau phase and longer duration than the adult (action potential duration at 80% repolarization = 60.1 vs. 12.6 ms). The simulation results are consistent with experimental data, including 1) decreased density and altered inactivation of transient outward K+ currents, 2) increased delayed rectifier K+ currents, 3) Ca2+ entry through T-type as well as L-type Ca2+ channels, 4) increased Ca2+ influx through Na+/Ca2+ exchange, and 5) Ca2+ transients resulting from transmembrane Ca2+ entry rather than release from the sarcoplasmic reticulum (SR). Simulations performed with the model generated novel predictions, including increased SR Ca2+ leak and elevated intracellular Na+ concentration in neonatal compared with adult myocytes. This new model can therefore be used for testing hypotheses and obtaining a better quantitative understanding of differences between neonatal and adult physiology.
</para>
        
        <para>
The complete original paper reference is cited below:
</para>
        <para>
          Mathematical model of the neonatal mouse ventricular action potential, Linda J. Wang and Eric A. Sobie, 2008,<emphasis>American Journal of Physiology: Heart and Circulatory Physiology</emphasis>, 294, H2565-H2575. <ulink url="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?cmd=Retrieve&amp;amp;db=PubMed&amp;amp;list_uids=18408122&amp;amp;dopt=Abstract">PubMed ID: 18408122</ulink>
        </para>
		
		       <informalfigure float="0" id="fig_cell_diagram">
          <mediaobject>
            <imageobject>
              <objectinfo>
                <title>cell diagram</title>
              </objectinfo>
              <imagedata fileref="wang_2008.png"/>
            </imageobject>
          </mediaobject>
          <caption>Schematic diagram of the neonatal mouse model.</caption>
        </informalfigure> 		
      </sect1>
    </article>
  </documentation>
   
   
   
   
   
   <units name="millisecond">
      <unit units="second" prefix="milli"/>
   </units>
   <units name="per_millisecond">
      <unit units="second" prefix="milli" exponent="-1"/>
   </units>
   <units name="micromolar_per_millisecond">
      <unit units="micromolar"/>
      <unit units="millisecond" exponent="-1"/>
   </units>
   <units name="per_micromolar_millisecond">
      <unit units="micromolar" exponent="-1"/>
      <unit units="millisecond" exponent="-1"/>
   </units>
   <units name="micromolar3_per_millisecond">
      <unit units="micromolar" exponent="-3"/>
      <unit units="millisecond" exponent="-1"/>
   </units>
   <units name="micromolar4_per_millisecond">
      <unit units="micromolar" exponent="-4"/>
      <unit units="millisecond" exponent="-1"/>
   </units>
   <units name="millivolt">
      <unit units="volt" prefix="milli"/>
   </units>
   <units name="per_millivolt">
      <unit units="volt" prefix="milli" exponent="-1"/>
   </units>
   <units name="millivolt2">
      <unit units="volt" prefix="milli" exponent="2"/>
   </units>
   <units name="per_millivolt_millisecond">
      <unit units="millivolt" exponent="-1"/>
      <unit units="millisecond" exponent="-1"/>
   </units>
   <units name="milliS_per_microF">
      <unit units="siemens" prefix="milli"/>
      <unit units="farad" prefix="micro" exponent="-1"/>
   </units>
   <units name="picoF">
      <unit units="farad" prefix="pico"/>
   </units>
   <units name="per_picoF">
      <unit units="picoF" exponent="-1"/>
   </units>
   <units name="microF_per_cm2">
      <unit units="farad" prefix="micro"/>
      <unit units="metre" prefix="centi" exponent="-2"/>
   </units>
   <units name="picoA_per_picoF">
      <unit units="ampere" prefix="pico"/>
      <unit units="farad" prefix="pico" exponent="-1"/>
   </units>
   <units name="micromolar">
      <unit units="mole" prefix="micro"/>
      <unit units="litre" exponent="-1"/>
   </units>
   <units name="joule_per_mole_kelvin">
      <unit units="joule"/>
      <unit units="mole" exponent="-1"/>
      <unit units="kelvin" exponent="-1"/>
   </units>
   <units name="coulomb_per_millimole">
      <unit units="coulomb"/>
      <unit units="mole" prefix="milli" exponent="-1"/>
   </units>
   <units name="cm2">
      <unit units="metre" prefix="centi" exponent="2"/>
   </units>
   <units name="cm2_per_microF">
      <unit units="metre" prefix="centi" exponent="2"/>
	  <unit units="farad" prefix="micro" exponent="-1"/>
   </units>
   <units name="cm_per_second">
      <unit units="metre" prefix="centi"/>
	  <unit units="second" exponent="-1"/>
   </units>
   <units name="microlitre">
      <unit units="litre" prefix="micro"/>
   </units>
   
   <component name="environment">
      <variable units="millisecond" public_interface="out" name="time"/>
   </component>
   
   <component name="membrane">
      <variable units="millivolt" public_interface="out" name="V" initial_value="-80.6475"/>
      <variable units="microF_per_cm2" public_interface="out" name="Cm" initial_value="1"/>
      <variable units="microlitre" public_interface="out" name="Vmyo" initial_value="2.2826e-6"/>
      <variable units="microlitre" public_interface="out" name="VJSR" initial_value="0.12e-8"/>
      <variable units="microlitre" public_interface="out" name="VNSR" initial_value="2.098e-7"/>
      <variable units="microlitre" public_interface="out" name="Vss" initial_value="3.0734e-8"/>
      <variable units="cm2" public_interface="out" name="Acap" initial_value="1.5410e-5"/>
      <variable units="micromolar" public_interface="out" name="Ko" initial_value="5400"/>
      <variable units="micromolar" public_interface="out" name="Nao" initial_value="140000"/>
      <variable units="micromolar" public_interface="out" name="Cao" initial_value="1000"/>
	  <variable units="micromolar" public_interface="out" name="Clo" initial_value="132000"/>
      <variable units="micromolar" public_interface="out" name="Cli" initial_value="30000"/>
      <variable units="joule_per_mole_kelvin" public_interface="out" name="R" initial_value="8.314"/>
      <variable units="kelvin" public_interface="out" name="T" initial_value="298"/>
      <variable units="coulomb_per_millimole" public_interface="out" name="F" initial_value="96.5"/>
      
	  <variable units="picoA_per_picoF" name="i_stim"/>
	  <variable units="millisecond" name="stim_start" initial_value="20"/>
      <variable units="millisecond" name="stim_end" initial_value="500000"/>
      <variable units="millisecond" name="stim_period" initial_value="2000.0"/>
      <variable units="millisecond" name="stim_duration" initial_value="0.5"/>
      <variable units="picoA_per_picoF" name="stim_amplitude" initial_value="-80"/>
      
	  <variable units="millisecond" public_interface="in" name="time"/>
      <variable units="picoA_per_picoF" public_interface="in" name="i_CaL"/>
	  <variable units="picoA_per_picoF" public_interface="in" name="i_CaT"/>
      <variable units="picoA_per_picoF" public_interface="in" name="i_pCa"/>
      <variable units="picoA_per_picoF" public_interface="in" name="i_NaCa"/>
      <variable units="picoA_per_picoF" public_interface="in" name="i_Cab"/>
      <variable units="picoA_per_picoF" public_interface="in" name="i_Na"/>
      <variable units="picoA_per_picoF" public_interface="in" name="i_Nab"/>
      <variable units="picoA_per_picoF" public_interface="in" name="i_NaK"/>
      <variable units="picoA_per_picoF" public_interface="in" name="i_Kto_f"/>
      <variable units="picoA_per_picoF" public_interface="in" name="i_K1"/>
      <variable units="picoA_per_picoF" public_interface="in" name="i_Ks"/>
      <variable units="picoA_per_picoF" public_interface="in" name="i_Kur"/>
      <variable units="picoA_per_picoF" public_interface="in" name="i_Kss"/>
      <variable units="picoA_per_picoF" public_interface="in" name="i_ClCa"/>
      <variable units="picoA_per_picoF" public_interface="in" name="i_Kr"/>
      
      
	  <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_stim</ci>
            <piecewise>
               <piece>
                  <ci>stim_amplitude</ci>
                  <apply>
                     <and/>
                     <apply>
                        <geq/>
                        <ci>time</ci>
                        <ci>stim_start</ci>
                     </apply>
                     <apply>
                        <leq/>
                        <ci>time</ci>
                        <ci>stim_end</ci>
                     </apply>
                     <apply>
                        <leq/>
                        <apply>
                           <minus/>
                           <apply>
                              <minus/>
                              <ci>time</ci>
                              <ci>stim_start</ci>
                           </apply>
                           <apply>
                              <times/>
                              <apply>
                                 <floor/>
                                 <apply>
                                    <divide/>
                                    <apply>
                                       <minus/>
                                       <ci>time</ci>
                                       <ci>stim_start</ci>
                                    </apply>
                                    <ci>stim_period</ci>
                                 </apply>
                              </apply>
                              <ci>stim_period</ci>
                           </apply>
                        </apply>
                        <ci>stim_duration</ci>
                     </apply>
                  </apply>
               </piece>
               <otherwise>
                  <cn cellml:units="picoA_per_picoF">0</cn>
               </otherwise>
            </piecewise>
         </apply>
         
		 <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>V</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <plus/>
                  <ci>i_CaL</ci>
				  <ci>i_CaT</ci>
                  <ci>i_pCa</ci>
                  <ci>i_NaCa</ci>
                  <ci>i_Cab</ci>
                  <ci>i_Na</ci>
                  <ci>i_Nab</ci>
                  <ci>i_NaK</ci>
                  <ci>i_Kto_f</ci>
                  <ci>i_K1</ci>
                  <ci>i_Ks</ci>
                  <ci>i_Kur</ci>
                  <ci>i_Kss</ci>
                  <ci>i_Kr</ci>
                  <ci>i_ClCa</ci>
                  <ci>i_stim</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   
   <component name="calcium_concentration">
      <variable units="micromolar" public_interface="out" name="Cai" initial_value="0.2049"/>
      <variable units="micromolar" public_interface="out" name="Cass" initial_value="0.1867"/>
      <variable units="micromolar" public_interface="out" name="CaJSR" initial_value="705.5417"/>
      <variable units="micromolar" public_interface="out" name="CaNSR" initial_value="703.3063"/>
      
	  <variable units="dimensionless" name="Bi"/>
      <variable units="dimensionless" name="Bss"/>
      <variable units="dimensionless" name="BJSR"/>
      <variable units="micromolar" name="CMDN_tot" initial_value="25"/>
      <variable units="micromolar" name="CSQN_tot" initial_value="15000"/>
      <variable units="micromolar" name="Km_CMDN" initial_value="0.238"/>
      <variable units="micromolar" name="Km_CSQN" initial_value="800"/>
      
	  <variable units="cm2" public_interface="in" name="Acap"/>
      <variable units="microlitre" public_interface="in" name="Vmyo"/>
      <variable units="microlitre" public_interface="in" name="VJSR"/>
      <variable units="microlitre" public_interface="in" name="Vss"/>
      <variable units="microlitre" public_interface="in" name="VNSR"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="microF_per_cm2" public_interface="in" name="Cm"/>
      <variable units="millisecond" public_interface="in" name="time"/>
      <variable units="coulomb_per_millimole" public_interface="in" name="F"/>
      <variable units="picoA_per_picoF" public_interface="in" name="i_Cab"/>
      <variable units="picoA_per_picoF" public_interface="in" name="i_CaL"/>
	  <variable units="picoA_per_picoF" public_interface="in" name="i_CaT"/>
      <variable units="picoA_per_picoF" public_interface="in" name="i_NaCa"/>
      <variable units="picoA_per_picoF" public_interface="in" name="i_pCa"/>
      <variable units="micromolar_per_millisecond" public_interface="in" name="J_leak"/>
      <variable units="micromolar_per_millisecond" public_interface="in" name="J_rel"/>
      <variable units="micromolar_per_millisecond" public_interface="in" name="J_up"/>
      <variable units="micromolar_per_millisecond" public_interface="in" name="J_tr"/>
      <variable units="micromolar_per_millisecond" public_interface="in" name="J_trpn"/>
      <variable units="micromolar_per_millisecond" public_interface="in" name="J_xfer"/>
      
	  <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>Cai</ci>
            </apply>
            <apply>
               <times/>
               <ci>Bi</ci>
               <apply>
                  <minus/>
                  <apply>
                     <plus/>
                     <ci>J_leak</ci>
                     <ci>J_xfer</ci>
                  </apply>
                  <apply>
                     <plus/>
                     <ci>J_up</ci>
                     <ci>J_trpn</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         
		 <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>Cass</ci>
            </apply>
            <apply>
               <times/>
               <ci>Bss</ci>
               <apply>
                  <minus/>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <ci>J_rel</ci>
                        <ci>VJSR</ci>
                     </apply>
                     <ci>Vss</ci>
                  </apply>
                  <apply>
                     <plus/>
                     <apply>
                        <divide/>
                        <apply>
                           <times/>
                           <ci>J_xfer</ci>
                           <ci>Vmyo</ci>
                        </apply>
                        <ci>Vss</ci>
                     </apply>
                     <apply>
                        <times/>
						<apply>
                           <minus/>
						   <apply>
                              <plus/>
							  <ci>i_CaL</ci>
							  <ci>i_Cab</ci>
							  <ci>i_pCa</ci>
							  <ci>i_CaT</ci>
						   </apply>
						   <apply>
                              <times/>
						      <cn cellml:units="dimensionless">2</cn>
							  <ci>i_CaT</ci>
						   </apply>
						</apply>
						<apply>
                           <divide/>
                           <apply>
                              <times/>
                              <ci>Acap</ci>
                              <ci>Cm</ci>
                           </apply>
                           <apply>
                              <times/>
                              <cn cellml:units="dimensionless">2</cn>
                              <ci>Vss</ci>
                              <ci>F</ci>
						   </apply>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         
		 <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>CaJSR</ci>
            </apply>
            <apply>
               <times/>
               <ci>BJSR</ci>
               <apply>
                  <minus/>
                  <ci>J_xfer</ci>
                  <ci>J_rel</ci>
               </apply>
            </apply>
         </apply>
         
		 <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>CaNSR</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <divide/>
                  <apply>
                     <times/>
                     <apply>
                        <minus/>
                        <ci>J_up</ci>
                        <ci>J_leak</ci>
                     </apply>
                     <ci>Vmyo</ci>
                  </apply>
                  <ci>VNSR</ci>
               </apply>
               <apply>
                  <divide/>
                  <apply>
                     <times/>
                     <ci>J_tr</ci>
                     <ci>VJSR</ci>
                  </apply>
                  <ci>VNSR</ci>
               </apply>
            </apply>
         </apply>
         
		 <apply>
            <eq/>
            <ci>Bi</ci>
            <apply>
               <power/>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <ci>CMDN_tot</ci>
                        <ci>Km_CMDN</ci>
                     </apply>
                     <apply>
                        <power/>
                        <apply>
                           <plus/>
                           <ci>Km_CMDN</ci>
                           <ci>Cai</ci>
                        </apply>
                        <cn cellml:units="dimensionless">2</cn>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <minus/>
                  <cn cellml:units="dimensionless">1</cn>
               </apply>
            </apply>
         </apply>
         
		 <apply>
            <eq/>
            <ci>Bss</ci>
            <apply>
               <power/>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <ci>CMDN_tot</ci>
                        <ci>Km_CMDN</ci>
                     </apply>
                     <apply>
                        <power/>
                        <apply>
                           <plus/>
                           <ci>Km_CMDN</ci>
                           <ci>Cass</ci>
                        </apply>
                        <cn cellml:units="dimensionless">2</cn>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <minus/>
                  <cn cellml:units="dimensionless">1</cn>
               </apply>
            </apply>
         </apply>
         
		 <apply>
            <eq/>
            <ci>BJSR</ci>
            <apply>
               <power/>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <ci>CSQN_tot</ci>
                        <ci>Km_CSQN</ci>
                     </apply>
                     <apply>
                        <power/>
                        <apply>
                           <plus/>
                           <ci>Km_CSQN</ci>
                           <ci>CaJSR</ci>
                        </apply>
                        <cn cellml:units="dimensionless">2</cn>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <minus/>
                  <cn cellml:units="dimensionless">1</cn>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   
   <component name="calcium_fluxes">
      <variable units="micromolar_per_millisecond" public_interface="out" name="J_leak"/>
      <variable units="micromolar_per_millisecond" public_interface="out" name="J_rel"/>
      <variable units="micromolar_per_millisecond" public_interface="out" name="J_up"/>
      <variable units="micromolar_per_millisecond" public_interface="out" name="J_tr"/>
      <variable units="micromolar_per_millisecond" public_interface="out" name="J_trpn"/>
      <variable units="micromolar_per_millisecond" public_interface="out" name="J_xfer"/>
      <variable units="per_micromolar_millisecond" public_interface="out" name="k_plus_htrpn" initial_value="0.00237"/>
      <variable units="per_millisecond" public_interface="out" name="k_minus_htrpn" initial_value="3.2e-5"/>
      <variable units="per_micromolar_millisecond" public_interface="out" name="k_plus_ltrpn" initial_value="0.0327"/>
      <variable units="per_millisecond" public_interface="out" name="k_minus_ltrpn" initial_value="0.0196"/>
      
	  <variable units="dimensionless" name="P_RyR" initial_value="2.8603e-14"/>
      <variable units="per_millisecond" name="v1" initial_value="0.45"/>
      <variable units="millisecond" name="tau_tr" initial_value="20"/>
      <variable units="per_millisecond" name="v2" initial_value="2.088e-5"/>
      <variable units="millisecond" name="tau_xfer" initial_value="8"/>
      <variable units="micromolar_per_millisecond" name="v3" initial_value="0.09"/>
      <variable units="micromolar" name="Km_up" initial_value="0.5"/>
      
	  <variable units="micromolar" public_interface="in" name="LTRPN_tot"/>
      <variable units="micromolar" public_interface="in" name="HTRPN_tot"/>
      <variable units="micromolar" public_interface="in" name="LTRPN_Ca"/>
      <variable units="micromolar" public_interface="in" name="HTRPN_Ca"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="millisecond" public_interface="in" name="time"/>
      <variable units="picoA_per_picoF" public_interface="in" name="i_CaL_max"/>
      <variable units="picoA_per_picoF" public_interface="in" name="i_CaL"/>
      <variable units="micromolar" public_interface="in" name="Cai"/>
      <variable units="micromolar" public_interface="in" name="Cass"/>
      <variable units="micromolar" public_interface="in" name="CaJSR"/>
      <variable units="micromolar" public_interface="in" name="CaNSR"/>
      <variable units="dimensionless" public_interface="in" name="P_O1"/>
      <variable units="dimensionless" public_interface="in" name="P_O2"/>
      
	  <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>J_rel</ci>
            <apply>
               <times/>
               <ci>v1</ci>
               <apply>
                  <plus/>
                  <ci>P_O1</ci>
                  <ci>P_O2</ci>
               </apply>
               <apply>
                  <minus/>
                  <ci>CaJSR</ci>
                  <ci>Cass</ci>
               </apply>
               <ci>P_RyR</ci>
            </apply>
         </apply>
         
		 <apply>
            <eq/>
            <ci>J_tr</ci>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>CaNSR</ci>
                  <ci>CaJSR</ci>
               </apply>
               <ci>tau_tr</ci>
            </apply>
         </apply>
         
		 <apply>
            <eq/>
            <ci>J_xfer</ci>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>Cass</ci>
                  <ci>Cai</ci>
               </apply>
               <ci>tau_xfer</ci>
            </apply>
         </apply>
         
		 <apply>
            <eq/>
            <ci>J_leak</ci>
            <apply>
               <times/>
               <ci>v2</ci>
               <apply>
                  <minus/>
                  <ci>CaNSR</ci>
                  <ci>Cai</ci>
               </apply>
            </apply>
         </apply>
        
		 <apply>
            <eq/>
            <ci>J_up</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <ci>v3</ci>
                  <apply>
                     <power/>
                     <ci>Cai</ci>
                     <cn cellml:units="dimensionless">2</cn>
                  </apply>
               </apply>
               <apply>
                  <plus/>
                  <apply>
                     <power/>
                     <ci>Km_up</ci>
                     <cn cellml:units="dimensionless">2</cn>
                  </apply>
                  <apply>
                     <power/>
                     <ci>Cai</ci>
                     <cn cellml:units="dimensionless">2</cn>
                  </apply>
               </apply>
            </apply>
         </apply>
        
  		 <apply>
            <eq/>
            <ci>J_trpn</ci>
            <apply>
               <minus/>
               <apply>
                  <plus/>
                  <apply>
                     <times/>
                     <ci>k_plus_htrpn</ci>
                     <ci>Cai</ci>
                     <apply>
                        <minus/>
                        <ci>HTRPN_tot</ci>
                        <ci>HTRPN_Ca</ci>
                     </apply>
                  </apply>
                  <apply>
                     <times/>
                     <ci>k_plus_ltrpn</ci>
                     <ci>Cai</ci>
                     <apply>
                        <minus/>
                        <ci>LTRPN_tot</ci>
                        <ci>LTRPN_Ca</ci>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <plus/>
                  <apply>
                     <times/>
                     <ci>k_minus_htrpn</ci>
                     <ci>HTRPN_Ca</ci>
                  </apply>
                  <apply>
                     <times/>
                     <ci>k_minus_ltrpn</ci>
                     <ci>LTRPN_Ca</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         
		 <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>P_RyR</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <apply>
                     <minus/>
                     <cn cellml:units="per_millisecond">0.04</cn>
                  </apply>
                  <ci>P_RyR</ci>
               </apply>
               <apply>
                  <times/>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <cn cellml:units="per_millisecond">0.1</cn>
                        <ci>i_CaL</ci>
                     </apply>
                     <ci>i_CaL_max</ci>
                  </apply>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <minus/>
                           <apply>
                              <power/>
                              <apply>
                                 <minus/>
                                 <ci>V</ci>
                                 <cn cellml:units="millivolt">5</cn>
                              </apply>
                              <cn cellml:units="dimensionless">2</cn>
                           </apply>
                        </apply>
                        <cn cellml:units="millivolt2">648</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   
   <component name="calcium_buffering">
      <variable units="micromolar" public_interface="out" name="LTRPN_tot" initial_value="35"/>
      <variable units="micromolar" public_interface="out" name="HTRPN_tot" initial_value="70"/>
      <variable units="micromolar" public_interface="out" name="LTRPN_Ca" initial_value="8.9220"/>
      <variable units="micromolar" public_interface="out" name="HTRPN_Ca" initial_value="66.0407"/>
      
	  <variable units="millisecond" public_interface="in" name="time"/>
      <variable units="micromolar" public_interface="in" name="Cai"/>
      <variable units="per_micromolar_millisecond" public_interface="in" name="k_plus_htrpn"/>
      <variable units="per_millisecond" public_interface="in" name="k_minus_htrpn"/>
      <variable units="per_micromolar_millisecond" public_interface="in" name="k_plus_ltrpn"/>
      <variable units="per_millisecond" public_interface="in" name="k_minus_ltrpn"/>
      
	  <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>LTRPN_Ca</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>k_plus_ltrpn</ci>
                  <ci>Cai</ci>
                  <apply>
                     <minus/>
                     <ci>LTRPN_tot</ci>
                     <ci>LTRPN_Ca</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>k_minus_ltrpn</ci>
                  <ci>LTRPN_Ca</ci>
               </apply>
            </apply>
         </apply>
         
		 <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>HTRPN_Ca</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>k_plus_htrpn</ci>
                  <ci>Cai</ci>
                  <apply>
                     <minus/>
                     <ci>HTRPN_tot</ci>
                     <ci>HTRPN_Ca</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>k_minus_htrpn</ci>
                  <ci>HTRPN_Ca</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   
   <component name="ryanodine_receptors">
      <variable units="dimensionless" public_interface="out" name="P_O1" initial_value="1.1e-3"/>
      <variable units="dimensionless" public_interface="out" name="P_O2" initial_value="2.9799e-8"/>
      <variable units="dimensionless" public_interface="out" name="P_C1"/>
      <variable units="dimensionless" public_interface="out" name="P_C2" initial_value="0.1003"/>
      
	  <variable units="micromolar4_per_millisecond" name="k_plus_a" initial_value="0.00608"/>
      <variable units="per_millisecond" name="k_minus_a" initial_value="0.07125"/>
      <variable units="micromolar3_per_millisecond" name="k_plus_b" initial_value="0.00405"/>
      <variable units="per_millisecond" name="k_minus_b" initial_value="0.965"/>
      <variable units="per_millisecond" name="k_plus_c" initial_value="0.009"/>
      <variable units="per_millisecond" name="k_minus_c" initial_value="0.0008"/>
      <variable units="dimensionless" name="m" initial_value="3"/>
      <variable units="dimensionless" name="n" initial_value="4"/>
      
	  <variable units="millisecond" public_interface="in" name="time"/>
      <variable units="micromolar" public_interface="in" name="Cass"/>
      
	  <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>P_O1</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <plus/>
                  <apply>
                     <times/>
                     <ci>k_plus_a</ci>
                     <apply>
                        <power/>
                        <ci>Cass</ci>
                        <ci>n</ci>
                     </apply>
                     <ci>P_C1</ci>
                  </apply>
                  <apply>
                     <times/>
                     <ci>k_minus_b</ci>
                     <ci>P_O2</ci>
                  </apply>
                  <apply>
                     <times/>
                     <ci>k_minus_c</ci>
                     <ci>P_C2</ci>
                  </apply>
               </apply>
               <apply>
                  <plus/>
                  <apply>
                     <times/>
                     <ci>k_minus_a</ci>
                     <ci>P_O1</ci>
                  </apply>
                  <apply>
                     <times/>
                     <ci>k_plus_b</ci>
                     <apply>
                        <power/>
                        <ci>Cass</ci>
                        <ci>m</ci>
                     </apply>
                     <ci>P_O1</ci>
                  </apply>
                  <apply>
                     <times/>
                     <ci>k_plus_c</ci>
                     <ci>P_O1</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         
		 <apply>
            <eq/>
            <ci>P_C1</ci>
            <apply>
               <minus/>
               <cn cellml:units="dimensionless">1</cn>
               <apply>
                  <plus/>
                  <ci>P_C2</ci>
                  <ci>P_O1</ci>
                  <ci>P_O2</ci>
               </apply>
            </apply>
         </apply>
         
		 <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>P_O2</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>k_plus_b</ci>
                  <apply>
                     <power/>
                     <ci>Cass</ci>
                     <ci>m</ci>
                  </apply>
                  <ci>P_O1</ci>
               </apply>
               <apply>
                  <times/>
                  <ci>k_minus_b</ci>
                  <ci>P_O2</ci>
               </apply>
            </apply>
         </apply>
         
		 <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>P_C2</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>k_plus_c</ci>
                  <ci>P_O1</ci>
               </apply>
               <apply>
                  <times/>
                  <ci>k_minus_c</ci>
                  <ci>P_C2</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   
   <component name="L_type_calcium_current">
      <variable units="picoA_per_picoF" public_interface="out" name="i_CaL"/>
      <variable units="picoA_per_picoF" public_interface="out" name="i_CaL_max" initial_value="7"/>
      <variable units="millivolt" public_interface="out" name="E_CaL" initial_value="63"/>
      
	  <variable units="milliS_per_microF" name="g_CaL" initial_value="0.19019"/>
      <variable units="dimensionless" name="O" initial_value="4.4776e-12"/>
      <variable units="dimensionless" name="C1"/>
      <variable units="dimensionless" name="C2" initial_value="1.6996e-4"/>
      <variable units="dimensionless" name="C3" initial_value="1.1754e-8"/>
      <variable units="dimensionless" name="C4" initial_value="6.1239e-10"/>
      <variable units="dimensionless" name="I1" initial_value="1.5450e-9"/>
      <variable units="dimensionless" name="I2" initial_value="6.4226e-8"/>
      <variable units="dimensionless" name="I3" initial_value="6.5201e-7"/>
      <variable units="per_millisecond" name="alpha"/>
      <variable units="per_millisecond" name="beta"/>
      <variable units="per_millisecond" name="gamma"/>
      <variable units="per_millisecond" name="Kpcf"/>
      <variable units="per_millisecond" name="Kpcb" initial_value="0.0005"/>
      <variable units="per_millisecond" name="Kpc_max" initial_value="0.23324"/>
      <variable units="micromolar" name="Kpc_half" initial_value="20"/>
      
	  <variable units="millisecond" public_interface="in" name="time"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="micromolar" public_interface="in" name="Cass"/>
     
 	  <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_CaL</ci>
            <apply>
               <times/>
               <ci>g_CaL</ci>
               <ci>O</ci>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>E_CaL</ci>
               </apply>
            </apply>
         </apply>
         
		 <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>O</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <plus/>
                  <apply>
                     <times/>
                     <ci>alpha</ci>
                     <ci>C4</ci>
                  </apply>
                  <apply>
                     <times/>
                     <ci>Kpcb</ci>
                     <ci>I1</ci>
                  </apply>
                  <apply>
                     <times/>
                     <cn cellml:units="dimensionless">0.001</cn>
                     <apply>
                        <minus/>
                        <apply>
                           <times/>
                           <ci>alpha</ci>
                           <ci>I2</ci>
                        </apply>
                        <apply>
                           <times/>
                           <ci>Kpcf</ci>
                           <ci>O</ci>
                        </apply>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <plus/>
                  <apply>
                     <times/>
                     <cn cellml:units="dimensionless">4</cn>
                     <ci>beta</ci>
                     <ci>O</ci>
                  </apply>
                  <apply>
                     <times/>
                     <ci>gamma</ci>
                     <ci>O</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
        
		 <apply>
            <eq/>
            <ci>C1</ci>
            <apply>
               <minus/>
               <cn cellml:units="dimensionless">1</cn>
               <apply>
                  <plus/>
                  <ci>O</ci>
                  <ci>C2</ci>
                  <ci>C2</ci>
                  <ci>C3</ci>
                  <ci>C4</ci>
                  <ci>I1</ci>
                  <ci>I2</ci>
                  <ci>I3</ci>
               </apply>
            </apply>
         </apply>
         
		 <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>C2</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <plus/>
                  <apply>
                     <times/>
                     <cn cellml:units="dimensionless">4</cn>
                     <ci>alpha</ci>
                     <ci>C1</ci>
                  </apply>
                  <apply>
                     <times/>
                     <cn cellml:units="dimensionless">2</cn>
                     <ci>beta</ci>
                     <ci>C3</ci>
                  </apply>
               </apply>
               <apply>
                  <plus/>
                  <apply>
                     <times/>
                     <ci>beta</ci>
                     <ci>C2</ci>
                  </apply>
                  <apply>
                     <times/>
                     <cn cellml:units="dimensionless">3</cn>
                     <ci>alpha</ci>
                     <ci>C2</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         
		 <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>C3</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <plus/>
                  <apply>
                     <times/>
                     <cn cellml:units="dimensionless">3</cn>
                     <ci>alpha</ci>
                     <ci>C2</ci>
                  </apply>
                  <apply>
                     <times/>
                     <cn cellml:units="dimensionless">3</cn>
                     <ci>beta</ci>
                     <ci>C4</ci>
                  </apply>
               </apply>
               <apply>
                  <plus/>
                  <apply>
                     <times/>
                     <cn cellml:units="dimensionless">2</cn>
                     <ci>beta</ci>
                     <ci>C3</ci>
                  </apply>
                  <apply>
                     <times/>
                     <cn cellml:units="dimensionless">2</cn>
                     <ci>alpha</ci>
                     <ci>C3</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         
		 <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>C4</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <plus/>
                  <apply>
                     <times/>
                     <cn cellml:units="dimensionless">2</cn>
                     <ci>alpha</ci>
                     <ci>C3</ci>
                  </apply>
                  <apply>
                     <times/>
                     <cn cellml:units="dimensionless">4</cn>
                     <ci>beta</ci>
                     <ci>O</ci>
                  </apply>
                  <apply>
                     <times/>
                     <cn cellml:units="millisecond">0.01</cn>
                     <apply>
                        <minus/>
                        <apply>
                           <times/>
                           <cn cellml:units="dimensionless">4</cn>
                           <ci>Kpcb</ci>
                           <ci>beta</ci>
                           <ci>I1</ci>
                        </apply>
                        <apply>
                           <times/>
                           <ci>alpha</ci>
                           <ci>gamma</ci>
                           <ci>C4</ci>
                        </apply>
                     </apply>
                  </apply>
                  <apply>
                     <times/>
                     <cn cellml:units="dimensionless">0.002</cn>
                     <apply>
                        <minus/>
                        <apply>
                           <times/>
                           <cn cellml:units="dimensionless">4</cn>
                           <ci>beta</ci>
                           <ci>I2</ci>
                        </apply>
                        <apply>
                           <times/>
                           <ci>Kpcf</ci>
                           <ci>C4</ci>
                        </apply>
                     </apply>
                  </apply>
                  <apply>
                     <times/>
                     <cn cellml:units="millisecond">4</cn>
                     <ci>beta</ci>
                     <ci>Kpcb</ci>
                     <ci>I3</ci>
                  </apply>
               </apply>
               <apply>
                  <plus/>
                  <apply>
                     <times/>
                     <cn cellml:units="dimensionless">3</cn>
                     <ci>beta</ci>
                     <ci>C4</ci>
                  </apply>
                  <apply>
                     <times/>
                     <ci>alpha</ci>
                     <ci>C4</ci>
                  </apply>
                  <apply>
                     <times/>
                     <cn cellml:units="millisecond">1</cn>
                     <ci>gamma</ci>
                     <ci>Kpcf</ci>
                     <ci>C4</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         
		 <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>I1</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <plus/>
                  <apply>
                     <times/>
                     <ci>gamma</ci>
                     <ci>O</ci>
                  </apply>
                  <apply>
                     <times/>
                     <cn cellml:units="dimensionless">0.001</cn>
                     <apply>
                        <minus/>
                        <apply>
                           <times/>
                           <ci>alpha</ci>
                           <ci>I3</ci>
                        </apply>
                        <apply>
                           <times/>
                           <ci>Kpcf</ci>
                           <ci>I1</ci>
                        </apply>
                     </apply>
                  </apply>
                  <apply>
                     <times/>
                     <cn cellml:units="millisecond">0.01</cn>
                     <apply>
                        <minus/>
                        <apply>
                           <times/>
                           <ci>alpha</ci>
                           <ci>gamma</ci>
                           <ci>C4</ci>
                        </apply>
                        <apply>
                           <times/>
                           <cn cellml:units="dimensionless">4</cn>
                           <ci>beta</ci>
                           <ci>Kpcf</ci>
                           <ci>I1</ci>
                        </apply>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>Kpcb</ci>
                  <ci>I1</ci>
               </apply>
            </apply>
         </apply>
         
		 <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>I2</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <plus/>
                  <apply>
                     <times/>
                     <cn cellml:units="dimensionless">0.001</cn>
                     <apply>
                        <minus/>
                        <apply>
                           <times/>
                           <ci>Kpcf</ci>
                           <ci>O</ci>
                        </apply>
                        <apply>
                           <times/>
                           <ci>alpha</ci>
                           <ci>I2</ci>
                        </apply>
                     </apply>
                  </apply>
                  <apply>
                     <times/>
                     <ci>Kpcb</ci>
                     <ci>I3</ci>
                  </apply>
                  <apply>
                     <times/>
                     <cn cellml:units="dimensionless">0.002</cn>
                     <apply>
                        <minus/>
                        <apply>
                           <times/>
                           <ci>Kpcf</ci>
                           <ci>C4</ci>
                        </apply>
                        <apply>
                           <times/>
                           <cn cellml:units="dimensionless">4</cn>
                           <ci>beta</ci>
                           <ci>I2</ci>
                        </apply>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>gamma</ci>
                  <ci>I2</ci>
               </apply>
            </apply>
         </apply>
         
		 <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>I3</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <plus/>
                  <apply>
                     <times/>
                     <cn cellml:units="dimensionless">0.001</cn>
                     <apply>
                        <minus/>
                        <apply>
                           <times/>
                           <ci>Kpcf</ci>
                           <ci>I1</ci>
                        </apply>
                        <apply>
                           <times/>
                           <ci>alpha</ci>
                           <ci>I3</ci>
                        </apply>
                     </apply>
                  </apply>
                  <apply>
                     <times/>
                     <ci>gamma</ci>
                     <ci>I2</ci>
                  </apply>
                  <apply>
                     <times/>
                     <cn cellml:units="millisecond">1</cn>
                     <ci>gamma</ci>
                     <ci>Kpcf</ci>
                     <ci>C4</ci>
                  </apply>
               </apply>
               <apply>
                  <plus/>
                  <apply>
                     <times/>
                     <cn cellml:units="millisecond">4</cn>
                     <ci>beta</ci>
                     <ci>Kpcb</ci>
                     <ci>I3</ci>
                  </apply>
                  <apply>
                     <times/>
                     <ci>Kpcb</ci>
                     <ci>I3</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         
		 <apply>
            <eq/>
            <ci>alpha</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <cn cellml:units="per_millisecond">0.4</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <plus/>
                           <ci>V</ci>
                           <cn cellml:units="millivolt">12</cn>
                        </apply>
                        <cn cellml:units="millivolt">10</cn>
                     </apply>
                  </apply>
                  <apply>
                     <minus/>
                     <apply>
                        <plus/>
                        <cn cellml:units="dimensionless">1</cn>
                        <apply>
                           <times/>
                           <cn cellml:units="dimensionless">0.7</cn>
                           <apply>
                              <exp/>
                              <apply>
                                 <divide/>
                                 <apply>
                                    <minus/>
                                    <apply>
                                       <power/>
                                       <apply>
                                          <plus/>
                                          <ci>V</ci>
                                          <cn cellml:units="millivolt">40</cn>
                                       </apply>
                                       <cn cellml:units="dimensionless">2</cn>
                                    </apply>
                                 </apply>
                                 <cn cellml:units="millivolt2">10</cn>
                              </apply>
                           </apply>
                        </apply>
                     </apply>
                     <apply>
                        <times/>
                        <cn cellml:units="dimensionless">0.75</cn>
                        <apply>
                           <exp/>
                           <apply>
                              <divide/>
                              <apply>
                                 <minus/>
                                 <apply>
                                    <power/>
                                    <apply>
                                       <plus/>
                                       <ci>V</ci>
                                       <cn cellml:units="millivolt">20</cn>
                                    </apply>
                                    <cn cellml:units="dimensionless">2</cn>
                                 </apply>
                              </apply>
                              <cn cellml:units="millivolt2">400</cn>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <times/>
                     <cn cellml:units="dimensionless">0.12</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <plus/>
                              <ci>V</ci>
                              <cn cellml:units="millivolt">12</cn>
                           </apply>
                           <cn cellml:units="millivolt">10</cn>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         
		 <apply>
            <eq/>
            <ci>beta</ci>
            <apply>
               <times/>
               <cn cellml:units="per_millisecond">0.05</cn>
               <apply>
                  <exp/>
                  <apply>
                     <divide/>
                     <apply>
                        <minus/>
                        <apply>
                           <plus/>
                           <ci>V</ci>
                           <cn cellml:units="millivolt">12</cn>
                        </apply>
                     </apply>
                     <cn cellml:units="millivolt">13</cn>
                  </apply>
               </apply>
            </apply>
         </apply>
         
		 <apply>
            <eq/>
            <ci>gamma</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <ci>Kpc_max</ci>
                  <ci>Cass</ci>
               </apply>
               <apply>
                  <plus/>
                  <ci>Kpc_half</ci>
                  <ci>Cass</ci>
               </apply>
            </apply>
         </apply>
         
		 <apply>
            <eq/>
            <ci>Kpcf</ci>
            <apply>
               <times/>
               <cn cellml:units="per_millisecond">13</cn>
               <apply>
                  <minus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <minus/>
                           <apply>
                              <power/>
                              <apply>
                                 <plus/>
                                 <ci>V</ci>
                                 <cn cellml:units="millivolt">14.5</cn>
                              </apply>
                              <cn cellml:units="dimensionless">2</cn>
                           </apply>
                        </apply>
                        <cn cellml:units="millivolt2">100</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   
   <component name="T_type_calcium_current">
      <variable units="picoA_per_picoF" public_interface="out" name="i_CaT"/>
      
	  <variable units="millivolt" name="E_CaT" initial_value="50"/>
      <variable units="milliS_per_microF" name="g_CaT" initial_value="0.055"/>
      
	  <variable units="millisecond" public_interface="in" private_interface="out" name="time"/>
      <variable units="millivolt" public_interface="in" private_interface="out" name="V"/>
      <variable units="dimensionless" private_interface="in" name="b"/>
      <variable units="dimensionless" private_interface="in" name="g"/>
      
	  <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_CaT</ci>
            <apply>
               <times/>
               <ci>g_CaT</ci>
               <ci>b</ci>
               <ci>g</ci>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>E_CaT</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   
   <component name="T_type_calcium_current_b_gate">
      <variable units="dimensionless" public_interface="out" name="b" initial_value="0.0047"/>
      <variable units="dimensionless" name="b_infinity"/>
      <variable units="millisecond" name="tau_b"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="millisecond" public_interface="in" name="time"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>b_infinity</ci>
            <apply>
               <divide/>
               <cn cellml:units="dimensionless">1</cn>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
						<apply>
                           <minus/>
						   <apply>
                              <plus/>
					          <ci>V</ci>
                              <cn cellml:units="millivolt">48</cn>
                           </apply>
						</apply>
						<cn cellml:units="millivolt">6.1</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
		 
		 <apply>
            <eq/>
            <ci>tau_b</ci>
            <apply>
               <plus/>
			   <cn cellml:units="millisecond">0.1</cn>
               <apply>
                  <divide/>
                  <cn cellml:units="millisecond">5.4</cn>
                  <apply>
                     <plus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
						   <apply>
                              <plus/>
					          <ci>V</ci>
                              <cn cellml:units="millivolt">100</cn>
                           </apply>
						   <cn cellml:units="millivolt">6.6</cn>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
        
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>b</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>b_infinity</ci>
				  <ci>b</ci>
               </apply>
               <ci>tau_b</ci>
            </apply>
         </apply>
      </math>
   </component>
   
   <component name="T_type_calcium_current_g_gate">
      <variable units="dimensionless" public_interface="out" name="g" initial_value="0.9020"/>
      <variable units="dimensionless" name="g_infinity"/>
      <variable units="millisecond" name="tau_g"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="millisecond" public_interface="in" name="time"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>g_infinity</ci>
            <apply>
               <divide/>
               <cn cellml:units="dimensionless">1</cn>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
						<apply>
                           <plus/>
					       <ci>V</ci>
                           <cn cellml:units="millivolt">66</cn>
                        </apply>
						<cn cellml:units="millivolt">6.6</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
		 
		 <apply>
            <eq/>
            <ci>tau_g</ci>
            <apply>
               <plus/>
			   <cn cellml:units="millisecond">8</cn>
               <apply>
                  <divide/>
                  <cn cellml:units="millisecond">32</cn>
                  <apply>
                     <plus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
						   <apply>
                              <plus/>
					          <ci>V</ci>
                              <cn cellml:units="millivolt">65</cn>
                           </apply>
						   <cn cellml:units="millivolt">5</cn>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
        
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>g</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>g_infinity</ci>
				  <ci>g</ci>
               </apply>
               <ci>tau_g</ci>
            </apply>
         </apply>
      </math>
   </component>
   
   <component name="calcium_pump_current">
      <variable units="picoA_per_picoF" public_interface="out" name="i_pCa"/>
      
	  <variable units="picoA_per_picoF" name="i_pCa_max" initial_value="0.2"/>
      <variable units="micromolar" name="Km_pCa" initial_value="0.5"/>
      
	  <variable units="micromolar" public_interface="in" name="Cai"/>
      
	  <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_pCa</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <ci>i_pCa_max</ci>
                  <apply>
                     <power/>
                     <ci>Cai</ci>
                     <cn cellml:units="dimensionless">2</cn>
                  </apply>
               </apply>
               <apply>
                  <plus/>
                  <apply>
                     <power/>
                     <ci>Km_pCa</ci>
                     <cn cellml:units="dimensionless">2</cn>
                  </apply>
                  <apply>
                     <power/>
                     <ci>Cai</ci>
                     <cn cellml:units="dimensionless">2</cn>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
  
   <component name="sodium_calcium_exchange_current">
      <variable units="picoA_per_picoF" public_interface="out" name="i_NaCa"/>
     
   	  <variable units="picoA_per_picoF" name="k_NaCa" initial_value="907.68"/>
      <variable units="micromolar" name="K_mNa" initial_value="87500"/>
      <variable units="micromolar" name="K_mCa" initial_value="1380"/>
      <variable units="dimensionless" name="k_sat" initial_value="0.1"/>
      <variable units="dimensionless" name="eta" initial_value="0.35"/>
      
	  <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="joule_per_mole_kelvin" public_interface="in" name="R"/>
      <variable units="kelvin" public_interface="in" name="T"/>
      <variable units="coulomb_per_millimole" public_interface="in" name="F"/>
      <variable units="micromolar" public_interface="in" name="Nai"/>
      <variable units="micromolar" public_interface="in" name="Nao"/>
      <variable units="micromolar" public_interface="in" name="Cai"/>
      <variable units="micromolar" public_interface="in" name="Cao"/>
      
	  <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_NaCa</ci>
            <apply>
               <times/>
               <apply>
                  <divide/>
                  <apply>
                     <times/>
                     <apply>
                        <divide/>
                        <apply>
                           <times/>
                           <apply>
                              <divide/>
                              <apply>
                                 <times/>
                                 <ci>k_NaCa</ci>
                                 <cn cellml:units="dimensionless">1</cn>
                              </apply>
                              <apply>
                                 <plus/>
                                 <apply>
                                    <power/>
                                    <ci>K_mNa</ci>
                                    <cn cellml:units="dimensionless">3</cn>
                                 </apply>
                                 <apply>
                                    <power/>
                                    <ci>Nao</ci>
                                    <cn cellml:units="dimensionless">3</cn>
                                 </apply>
                              </apply>
                           </apply>
                           <cn cellml:units="dimensionless">1</cn>
                        </apply>
                        <apply>
                           <plus/>
                           <ci>K_mCa</ci>
                           <ci>Cao</ci>
                        </apply>
                     </apply>
                     <cn cellml:units="dimensionless">1</cn>
                  </apply>
                  <apply>
                     <plus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <apply>
                        <times/>
                        <ci>k_sat</ci>
                        <apply>
                           <exp/>
                           <apply>
                              <divide/>
                              <apply>
                                 <times/>
                                 <apply>
                                    <minus/>
                                    <ci>eta</ci>
                                    <cn cellml:units="dimensionless">1</cn>
                                 </apply>
                                 <ci>V</ci>
                                 <ci>F</ci>
                              </apply>
                              <apply>
                                 <times/>
                                 <ci>R</ci>
                                 <ci>T</ci>
                              </apply>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <minus/>
                  <apply>
                     <times/>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <times/>
                              <ci>eta</ci>
                              <ci>V</ci>
                              <ci>F</ci>
                           </apply>
                           <apply>
                              <times/>
                              <ci>R</ci>
                              <ci>T</ci>
                           </apply>
                        </apply>
                     </apply>
                     <apply>
                        <power/>
                        <ci>Nai</ci>
                        <cn cellml:units="dimensionless">3</cn>
                     </apply>
                     <ci>Cao</ci>
                  </apply>
                  <apply>
                     <times/>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <times/>
                              <apply>
                                 <minus/>
                                 <ci>eta</ci>
                                 <cn cellml:units="dimensionless">1</cn>
                              </apply>
                              <ci>V</ci>
                              <ci>F</ci>
                           </apply>
                           <apply>
                              <times/>
                              <ci>R</ci>
                              <ci>T</ci>
                           </apply>
                        </apply>
                     </apply>
                     <apply>
                        <power/>
                        <ci>Nao</ci>
                        <cn cellml:units="dimensionless">3</cn>
                     </apply>
                     <ci>Cai</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   
   <component name="calcium_background_current">
      <variable units="picoA_per_picoF" public_interface="out" name="i_Cab"/>
      
	  <variable units="milliS_per_microF" name="g_Cab" initial_value="0.00025"/>
      <variable units="millivolt" name="E_CaN"/>
      
	  <variable units="millivolt" public_interface="in" private_interface="out" name="V"/>
      <variable units="joule_per_mole_kelvin" public_interface="in" name="R"/>
      <variable units="kelvin" public_interface="in" name="T"/>
      <variable units="coulomb_per_millimole" public_interface="in" name="F"/>
      <variable units="micromolar" public_interface="in" name="Cai"/>
      <variable units="micromolar" public_interface="in" name="Cao"/>
	  
	  <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_Cab</ci>
            <apply>
               <times/>
               <ci>g_Cab</ci>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>E_CaN</ci>
               </apply>
            </apply>
         </apply>
		 
		 <apply>
            <eq/>
            <ci>E_CaN</ci>
            <apply>
               <times/>
               <apply>
                  <divide/>
                  <apply>
                     <times/>
                     <ci>R</ci>
                     <ci>T</ci>
                  </apply>
                  <apply>
                     <times/>
                     <cn cellml:units="dimensionless">2</cn>
                     <ci>F</ci>
                  </apply>
               </apply>
               <apply>
                  <ln/>
                  <apply>
                     <divide/>
                     <ci>Cao</ci>
                     <ci>Cai</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   
   <component name="sodium_concentration">
      <variable units="micromolar" public_interface="out" name="Nai" initial_value="2.1747e4"/>
      
	  <variable units="millisecond" public_interface="in" name="time"/>
      <variable units="coulomb_per_millimole" public_interface="in" name="F"/>
      <variable units="microlitre" public_interface="in" name="Vmyo"/>
      <variable units="microF_per_cm2" public_interface="in" name="Cm"/>
      <variable units="cm2" public_interface="in" name="Acap"/>
      <variable units="picoA_per_picoF" public_interface="in" name="i_Na"/>
      <variable units="picoA_per_picoF" public_interface="in" name="i_NaCa"/>
      <variable units="picoA_per_picoF" public_interface="in" name="i_NaK"/>
      <variable units="picoA_per_picoF" public_interface="in" name="i_Nab"/>
      
	  <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>Nai</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <apply>
                     <minus/>
                     <apply>
                        <plus/>
                        <ci>i_Na</ci>
                        <ci>i_Nab</ci>
                        <apply>
                           <times/>
                           <cn cellml:units="dimensionless">3</cn>
                           <ci>i_NaK</ci>
                        </apply>
                        <apply>
                           <times/>
                           <cn cellml:units="dimensionless">3</cn>
                           <ci>i_NaCa</ci>
                        </apply>
                     </apply>
                  </apply>
                  <ci>Acap</ci>
                  <ci>Cm</ci>
               </apply>
               <apply>
                  <times/>
                  <ci>Vmyo</ci>
                  <ci>F</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   
   <component name="fast_sodium_current">
      <variable units="picoA_per_picoF" public_interface="out" name="i_Na"/>
      <variable units="millivolt" public_interface="out" name="E_Na"/>
      
	  <variable units="milliS_per_microF" name="g_Na" initial_value="10"/>
      
	  <variable units="millisecond" public_interface="in" private_interface="out" name="time"/>
      <variable units="millivolt" public_interface="in" private_interface="out" name="V"/>
      <variable units="joule_per_mole_kelvin" public_interface="in" name="R"/>
      <variable units="kelvin" public_interface="in" name="T"/>
      <variable units="coulomb_per_millimole" public_interface="in" name="F"/>
      <variable units="micromolar" public_interface="in" name="Nao"/>
      <variable units="micromolar" public_interface="in" name="Nai"/>
	  
	  <variable units="dimensionless" private_interface="in" name="m"/>
      <variable units="dimensionless" private_interface="in" name="h"/>
      <variable units="dimensionless" private_interface="in" name="j"/>
      
	  <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_Na</ci>
            <apply>
               <times/>
               <ci>g_Na</ci>
			   <apply>
                  <power/>
                  <ci>m</ci>
				  <cn cellml:units="dimensionless">3</cn>
			   </apply>
			   <ci>h</ci>
			   <ci>j</ci>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>E_Na</ci>
               </apply>
            </apply>
         </apply>
         
		 <apply>
            <eq/>
            <ci>E_Na</ci>
            <apply>
               <times/>
               <apply>
                  <divide/>
                  <apply>
                     <times/>
                     <ci>R</ci>
                     <ci>T</ci>
                  </apply>
                  <ci>F</ci>
               </apply>
               <apply>
                  <ln/>
                  <apply>
                     <divide/>
                     <ci>Nao</ci>
                     <ci>Nai</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   
   <component name="fast_sodium_current_m_gate">
      <variable units="dimensionless" public_interface="out" name="m" initial_value="0.0032"/>
      <variable units="per_millisecond" name="alpha_m"/>
      <variable units="per_millisecond" name="beta_m"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="millisecond" public_interface="in" name="time"/>
      
	  <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>alpha_m</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <cn cellml:units="per_millivolt_millisecond">0.32</cn>
                  <apply>
                     <plus/>
                     <ci>V</ci>
					 <cn cellml:units="millivolt">47.13</cn>
                  </apply>
               </apply>
			   <apply>
                  <minus/>
                  <cn cellml:units="dimensionless">1</cn>
				  <apply>
                     <exp/>
					 <apply>
                        <times/>
                        <cn cellml:units="per_millivolt">-0.1</cn>
                        <apply>
                           <plus/>
                           <ci>V</ci>
					       <cn cellml:units="millivolt">47.13</cn>
                        </apply>
                     </apply>
				  </apply>
               </apply>
            </apply>
         </apply>
         
		 <apply>
            <eq/>
            <ci>beta_m</ci>
            <apply>
               <times/>
               <cn cellml:units="per_millisecond">0.08</cn>
               <apply>
                  <exp/>
                  <apply>
                     <divide/>
                     <apply>
                        <minus/>
                        <ci>V</ci>
                     </apply>
                     <cn cellml:units="millivolt">11</cn>
                  </apply>
               </apply>
            </apply>
         </apply>
         
		 <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>m</ci>
            </apply>
            <apply>
               <minus/>
			   <ci>alpha_m</ci>
               <apply>
                  <times/>
                  <apply>
                     <plus/>
                     <ci>alpha_m</ci>
                     <ci>beta_m</ci>
                  </apply>
				  <ci>m</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   
   <component name="fast_sodium_current_h_gate">
      <variable units="dimensionless" public_interface="out" name="h" initial_value="0.9606"/>
      <variable units="per_millisecond" name="alpha_h"/>
      <variable units="per_millisecond" name="beta_h"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="millisecond" public_interface="in" name="time"/>
      
	  <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>alpha_h</ci>
            <piecewise>
               <piece>
                  <apply>
                     <times/>
                     <cn cellml:units="per_millisecond">0.135</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <plus/>
                              <cn cellml:units="millivolt">80</cn>
                              <ci>V</ci>
                           </apply>
                           <apply>
                              <minus/>
                              <cn cellml:units="millivolt">6.8</cn>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
                  <apply>
                     <lt/>
                     <ci>V</ci>
                     <apply>
                        <minus/>
                        <cn cellml:units="millivolt">40</cn>
                     </apply>
                  </apply>
               </piece>
               <otherwise>
                  <cn cellml:units="per_millisecond">0</cn>
               </otherwise>
            </piecewise>
         </apply>
         
		 <apply>
            <eq/>
            <ci>beta_h</ci>
            <piecewise>
               <piece>
                  <apply>
                     <plus/>
                     <apply>
                        <times/>
                        <cn cellml:units="per_millisecond">3.56</cn>
                        <apply>
                           <exp/>
                           <apply>
                              <times/>
                              <cn cellml:units="per_millivolt">0.079</cn>
                              <ci>V</ci>
                           </apply>
                        </apply>
                     </apply>
                     <apply>
                        <times/>
                        <cn cellml:units="per_millisecond">310000</cn>
                        <apply>
                           <exp/>
                           <apply>
                              <times/>
                              <cn cellml:units="per_millivolt">0.35</cn>
                              <ci>V</ci>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
                  <apply>
                     <lt/>
                     <ci>V</ci>
                     <apply>
                        <minus/>
                        <cn cellml:units="millivolt">40</cn>
                     </apply>
                  </apply>
               </piece>
               <otherwise>
                  <apply>
                     <divide/>
                     <cn cellml:units="per_millisecond">1</cn>
                     <apply>
                        <times/>
                        <cn cellml:units="dimensionless">0.13</cn>
                        <apply>
                           <plus/>
                           <cn cellml:units="dimensionless">1</cn>
                           <apply>
                              <exp/>
                              <apply>
                                 <divide/>
                                 <apply>
                                    <plus/>
                                    <ci>V</ci>
                                    <cn cellml:units="millivolt">10.66</cn>
                                 </apply>
                                 <apply>
                                    <minus/>
                                    <cn cellml:units="millivolt">11.1</cn>
                                 </apply>
                              </apply>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </otherwise>
            </piecewise>
         </apply>
         
		 <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>h</ci>
            </apply>
            <apply>
               <minus/>
			   <ci>alpha_h</ci>
               <apply>
                  <times/>
                  <apply>
                     <plus/>
                     <ci>alpha_h</ci>
                     <ci>beta_h</ci>
                  </apply>
				  <ci>h</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   
   <component name="fast_sodium_current_j_gate">
      <variable units="dimensionless" public_interface="out" name="j" initial_value="0.9744"/>
      <variable units="per_millisecond" name="alpha_j"/>
      <variable units="per_millisecond" name="beta_j"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="millisecond" public_interface="in" name="time"/>
      
	  <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>alpha_j</ci>
            <piecewise>
               <piece>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <apply>
                           <minus/>
                           <apply>
                              <plus/>
                              <apply>
                                 <times/>
                                 <cn cellml:units="per_millivolt_millisecond">127140</cn>
                                 <apply>
                                    <exp/>
                                    <apply>
                                       <times/>
                                       <cn cellml:units="per_millivolt">0.2444</cn>
                                       <ci>V</ci>
                                    </apply>
                                 </apply>
                              </apply>
                              <apply>
                                 <times/>
                                 <cn cellml:units="per_millivolt_millisecond">3.474e-5</cn>
                                 <apply>
                                    <exp/>
                                    <apply>
                                       <times/>
                                       <apply>
                                          <minus/>
                                          <cn cellml:units="per_millivolt">0.04391</cn>
                                       </apply>
                                       <ci>V</ci>
                                    </apply>
                                 </apply>
                              </apply>
                           </apply>
                        </apply>
                        <apply>
                           <plus/>
                           <ci>V</ci>
                           <cn cellml:units="millivolt">37.78</cn>
                        </apply>
                     </apply>
                     <apply>
                        <plus/>
                        <cn cellml:units="dimensionless">1</cn>
                        <apply>
                           <exp/>
                           <apply>
                              <times/>
                              <cn cellml:units="per_millivolt">0.311</cn>
                              <apply>
                                 <plus/>
                                 <ci>V</ci>
                                 <cn cellml:units="millivolt">79.23</cn>
                              </apply>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
                  <apply>
                     <lt/>
                     <ci>V</ci>
                     <apply>
                        <minus/>
                        <cn cellml:units="millivolt">40</cn>
                     </apply>
                  </apply>
               </piece>
               <otherwise>
                  <cn cellml:units="per_millisecond">0</cn>
               </otherwise>
            </piecewise>
         </apply>
         
		 <apply>
            <eq/>
            <ci>beta_j</ci>
            <piecewise>
               <piece>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <cn cellml:units="per_millisecond">0.1212</cn>
                        <apply>
                           <exp/>
                           <apply>
                              <times/>
                              <apply>
                                 <minus/>
                                 <cn cellml:units="per_millivolt">0.01052</cn>
                              </apply>
                              <ci>V</ci>
                           </apply>
                        </apply>
                     </apply>
                     <apply>
                        <plus/>
                        <cn cellml:units="dimensionless">1</cn>
                        <apply>
                           <exp/>
                           <apply>
                              <times/>
                              <apply>
                                 <minus/>
                                 <cn cellml:units="per_millivolt">0.1378</cn>
                              </apply>
                              <apply>
                                 <plus/>
                                 <ci>V</ci>
                                 <cn cellml:units="millivolt">40.14</cn>
                              </apply>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
                  <apply>
                     <lt/>
                     <ci>V</ci>
                     <apply>
                        <minus/>
                        <cn cellml:units="millivolt">40</cn>
                     </apply>
                  </apply>
               </piece>
               <otherwise>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <cn cellml:units="per_millisecond">0.3</cn>
                        <apply>
                           <exp/>
                           <apply>
                              <times/>
                              <apply>
                                 <minus/>
                                 <cn cellml:units="per_millivolt">2.535e-7</cn>
                              </apply>
                              <ci>V</ci>
                           </apply>
                        </apply>
                     </apply>
                     <apply>
                        <plus/>
                        <cn cellml:units="dimensionless">1</cn>
                        <apply>
                           <exp/>
                           <apply>
                              <times/>
                              <apply>
                                 <minus/>
                                 <cn cellml:units="per_millivolt">0.1</cn>
                              </apply>
                              <apply>
                                 <plus/>
                                 <ci>V</ci>
                                 <cn cellml:units="millivolt">32</cn>
                              </apply>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </otherwise>
            </piecewise>
         </apply>
         
		 <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>j</ci>
            </apply>
            <apply>
               <minus/>
			   <ci>alpha_j</ci>
               <apply>
                  <times/>
                  <apply>
                     <plus/>
                     <ci>alpha_j</ci>
                     <ci>beta_j</ci>
                  </apply>
				  <ci>j</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   
   <component name="sodium_background_current">
      <variable units="picoA_per_picoF" public_interface="out" name="i_Nab"/>
      
	  <variable units="milliS_per_microF" name="g_Nab" initial_value="0.0026"/>
      
	  <variable units="millivolt" public_interface="in" name="E_Na"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      
	  <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_Nab</ci>
            <apply>
               <times/>
               <ci>g_Nab</ci>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>E_Na</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   
   <component name="potassium_concentration">
      <variable units="micromolar" public_interface="out" name="Ki" initial_value="1.3645e5"/>
      
	  <variable units="millisecond" public_interface="in" name="time"/>
      <variable units="coulomb_per_millimole" public_interface="in" name="F"/>
      <variable units="microlitre" public_interface="in" name="Vmyo"/>
      <variable units="microF_per_cm2" public_interface="in" name="Cm"/>
      <variable units="cm2" public_interface="in" name="Acap"/>
      <variable units="picoA_per_picoF" public_interface="in" name="i_Kto_f"/>
      <variable units="picoA_per_picoF" public_interface="in" name="i_K1"/>
      <variable units="picoA_per_picoF" public_interface="in" name="i_Ks"/>
      <variable units="picoA_per_picoF" public_interface="in" name="i_Kss"/>
      <variable units="picoA_per_picoF" public_interface="in" name="i_Kur"/>
      <variable units="picoA_per_picoF" public_interface="in" name="i_Kr"/>
      <variable units="picoA_per_picoF" public_interface="in" name="i_NaK"/>
      
	  <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>Ki</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <apply>
                     <minus/>
                     <apply>
                        <minus/>
                        <apply>
                           <plus/>
                           <ci>i_Kto_f</ci>
                           <ci>i_K1</ci>
                           <ci>i_Ks</ci>
                           <ci>i_Kss</ci>
                           <ci>i_Kur</ci>
                           <ci>i_Kr</ci>
                        </apply>
                        <apply>
                           <times/>
                           <cn cellml:units="dimensionless">2</cn>
                           <ci>i_NaK</ci>
                        </apply>
                     </apply>
                  </apply>
                  <ci>Acap</ci>
                  <ci>Cm</ci>
               </apply>
               <apply>
                  <times/>
                  <ci>Vmyo</ci>
                  <ci>F</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   
   <component name="fast_transient_outward_potassium_current">
      <variable units="picoA_per_picoF" public_interface="out" name="i_Kto_f"/>
      <variable units="millivolt" public_interface="out" name="E_K"/>
      
	  <variable units="milliS_per_microF" name="g_Kto_f" initial_value="0.1017"/>
      <variable units="dimensionless" name="ato_f" initial_value="0.0032"/>
      <variable units="dimensionless" name="ito_f" initial_value="0.7530"/>
	  <variable units="dimensionless" name="ito_f_infinity"/>
	  <variable units="millisecond" name="tau_ito_f"/>
      <variable units="per_millisecond" name="alpha_a"/>
      <variable units="per_millisecond" name="beta_a"/>
      <variable units="per_millisecond" name="alpha_i"/>
      <variable units="per_millisecond" name="beta_i"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="joule_per_mole_kelvin" public_interface="in" name="R"/>
      <variable units="kelvin" public_interface="in" name="T"/>
      <variable units="coulomb_per_millimole" public_interface="in" name="F"/>
      <variable units="micromolar" public_interface="in" name="Ki"/>
      <variable units="micromolar" public_interface="in" name="Ko"/>
      <variable units="millisecond" public_interface="in" name="time"/>
      
	  <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_Kto_f</ci>
            <apply>
               <times/>
               <ci>g_Kto_f</ci>
               <apply>
                  <power/>
                  <ci>ato_f</ci>
                  <cn cellml:units="dimensionless">6.5</cn>
               </apply>
               <ci>ito_f</ci>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>E_K</ci>
               </apply>
            </apply>
         </apply>
         
		 <apply>
            <eq/>
            <ci>E_K</ci>
            <apply>
               <times/>
               <apply>
                  <divide/>
                  <apply>
                     <times/>
                     <ci>R</ci>
                     <ci>T</ci>
                  </apply>
                  <ci>F</ci>
               </apply>
               <apply>
                  <ln/>
                  <apply>
                     <divide/>
                     <ci>Ko</ci>
                     <ci>Ki</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         
		 <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>ato_f</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>alpha_a</ci>
                  <apply>
                     <minus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <ci>ato_f</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>beta_a</ci>
                  <ci>ato_f</ci>
               </apply>
            </apply>
         </apply>
         
		 <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>ito_f</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
				  <ci>ito_f_infinity</ci>
				  <ci>ito_f</ci>
               </apply>
			   <ci>tau_ito_f</ci>
            </apply>
		 </apply>
		 
		 <apply>
            <eq/>
			<ci>ito_f_infinity</ci>
			<apply>
               <divide/>
			   <ci>alpha_i</ci>
			   <apply>
                  <plus/>
			      <ci>alpha_i</ci>
				  <ci>beta_i</ci>
			   </apply>
			</apply>
		 </apply>
			
	     <apply>
            <eq/>
            <ci>tau_ito_f</ci>
            <apply>
               <divide/>
			   <cn cellml:units="millisecond">1</cn>
               <apply>
                  <plus/>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <cn cellml:units="dimensionless">0.000152</cn>
                        <apply>
                           <exp/>
                           <apply>
                              <divide/>
                              <apply>
                                 <minus/>
                                 <apply>
                                    <plus/>
                                    <ci>V</ci>
                                    <cn cellml:units="millivolt">13.5</cn>
                                 </apply>
                              </apply>
                              <cn cellml:units="millivolt">7</cn>
                           </apply>
                        </apply>
                     </apply>
					 <apply>
                        <plus/>
						<apply>
                           <times/>
                           <cn cellml:units="dimensionless">0.067083</cn>
                           <apply>
                              <exp/>
                              <apply>
                                 <divide/>
                                 <apply>
                                    <minus/>
                                    <apply>
                                       <plus/>
                                       <ci>V</ci>
                                       <cn cellml:units="millivolt">33.5</cn>
                                    </apply>
                                 </apply>
                                 <cn cellml:units="millivolt">7</cn>
                              </apply>
                           </apply>
						</apply>
						<cn cellml:units="dimensionless">1</cn>
                     </apply>
				  </apply>
				  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <cn cellml:units="dimensionless">0.00095</cn>
                        <apply>
                           <exp/>
                           <apply>
                              <divide/>
                              <apply>
                                 <plus/>
                                 <ci>V</ci>
                                 <cn cellml:units="millivolt">33.5</cn>
                              </apply>
                              <cn cellml:units="millivolt">7</cn>
                           </apply>
                        </apply>
                     </apply>
					 <apply>
                        <plus/>
						<apply>
                           <times/>
                           <cn cellml:units="dimensionless">0.051335</cn>
                           <apply>
                              <exp/>
                              <apply>
                                 <divide/>
                                 <apply>
                                    <plus/>
                                    <ci>V</ci>
                                    <cn cellml:units="millivolt">33.5</cn>
                                 </apply>
                                 <cn cellml:units="millivolt">7</cn>
                              </apply>
                           </apply>
						</apply>
						<cn cellml:units="dimensionless">1</cn>
                     </apply>
				  </apply>
               </apply>
			</apply>
         </apply>
		 
		 <apply>
            <eq/>
            <ci>alpha_a</ci>
            <apply>
               <times/>
               <cn cellml:units="per_millisecond">0.18064</cn>
               <apply>
                  <exp/>
                  <apply>
                     <times/>
                     <cn cellml:units="per_millivolt">0.03577</cn>
                     <apply>
                        <plus/>
                        <ci>V</ci>
                        <cn cellml:units="millivolt">30</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         
		 <apply>
            <eq/>
            <ci>beta_a</ci>
            <apply>
               <times/>
               <cn cellml:units="per_millisecond">0.3956</cn>
               <apply>
                  <exp/>
                  <apply>
                     <times/>
                     <apply>
                        <minus/>
                        <cn cellml:units="per_millivolt">0.06237</cn>
                     </apply>
                     <apply>
                        <plus/>
                        <ci>V</ci>
                        <cn cellml:units="millivolt">30</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         
		 <apply>
            <eq/>
            <ci>alpha_i</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <cn cellml:units="per_millisecond">0.000152</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <minus/>
                           <apply>
                              <minus/>
                              <ci>V</ci>
                              <cn cellml:units="millivolt">3.81</cn>
                           </apply>
                        </apply>
                        <cn cellml:units="millivolt">15.75</cn>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <plus/>
                  <apply>
                     <times/>
                     <cn cellml:units="dimensionless">0.0067083</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <minus/>
                              <apply>
                                 <plus/>
                                 <ci>V</ci>
                                 <cn cellml:units="millivolt">132.05</cn>
                              </apply>
                           </apply>
                           <cn cellml:units="millivolt">15.75</cn>
                        </apply>
                     </apply>
                  </apply>
                  <cn cellml:units="dimensionless">1</cn>
               </apply>
            </apply>
         </apply>
         
		 <apply>
            <eq/>
            <ci>beta_i</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <cn cellml:units="per_millisecond">0.00095</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <plus/>
                           <ci>V</ci>
                           <cn cellml:units="millivolt">132.05</cn>
                        </apply>
                        <cn cellml:units="millivolt">15.75</cn>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <plus/>
                  <apply>
                     <times/>
                     <cn cellml:units="dimensionless">0.051335</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <plus/>
                              <ci>V</ci>
                              <cn cellml:units="millivolt">132.05</cn>
                           </apply>
                           <cn cellml:units="millivolt">15.75</cn>
                        </apply>
                     </apply>
                  </apply>
                  <cn cellml:units="dimensionless">1</cn>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   
   <component name="ultra_rapidly_activating_delayed_rectifier_potassium_current">
      <variable units="picoA_per_picoF" public_interface="out" name="i_Kur"/>
      <variable units="dimensionless" public_interface="out" name="ass"/>
	  
	  <variable units="milliS_per_microF" name="g_Kur" initial_value="0.0048"/>
      <variable units="dimensionless" name="aur" initial_value="5.2471e-4"/>
      <variable units="dimensionless" name="iur" initial_value="0.9866"/>
      <variable units="millisecond" name="tau_aur"/>
      <variable units="millisecond" name="tau_iur"/>
      <variable units="dimensionless" name="iss"/>
	  
	  <variable units="millivolt" public_interface="in" name="E_K"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="millisecond" public_interface="in" name="time"/>
      
	  <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_Kur</ci>
            <apply>
               <times/>
               <ci>g_Kur</ci>
               <ci>aur</ci>
               <ci>iur</ci>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>E_K</ci>
               </apply>
            </apply>
         </apply>
         
		 <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>aur</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>ass</ci>
                  <ci>aur</ci>
               </apply>
               <ci>tau_aur</ci>
            </apply>
         </apply>
         
		 <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>iur</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>iss</ci>
                  <ci>iur</ci>
               </apply>
               <ci>tau_iur</ci>
            </apply>
         </apply>
		 
		 <apply>
            <eq/>
            <ci>ass</ci>
            <apply>
               <divide/>
               <cn cellml:units="dimensionless">1</cn>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <minus/>
                           <apply>
                              <plus/>
                              <ci>V</ci>
                              <cn cellml:units="millivolt">22.5</cn>
                           </apply>
                        </apply>
                        <cn cellml:units="millivolt">7.7</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         
		 <apply>
            <eq/>
            <ci>tau_aur</ci>
            <apply>
               <plus/>
               <apply>
                  <times/>
                  <cn cellml:units="millisecond">0.493</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <times/>
                        <apply>
                           <minus/>
                           <cn cellml:units="per_millivolt">0.0629</cn>
                        </apply>
                        <ci>V</ci>
                     </apply>
                  </apply>
               </apply>
               <cn cellml:units="millisecond">2.058</cn>
            </apply>
         </apply>
		 
		 <apply>
            <eq/>
            <ci>iss</ci>
            <apply>
               <divide/>
               <cn cellml:units="dimensionless">1</cn>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <plus/>
                           <ci>V</ci>
                           <cn cellml:units="millivolt">45.2</cn>
                        </apply>
                        <cn cellml:units="millivolt">5.7</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         
		 <apply>
            <eq/>
            <ci>tau_iur</ci>
            <apply>
               <minus/>
               <cn cellml:units="millisecond">1200</cn>
               <apply>
                  <divide/>
                  <cn cellml:units="millisecond">170</cn>
                  <apply>
                     <plus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <plus/>
                              <ci>V</ci>
                              <cn cellml:units="millivolt">45.2</cn>
                           </apply>
                           <cn cellml:units="millivolt">5.7</cn>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   
   <component name="non_inactivating_steady_state_potassium_current">
      <variable units="picoA_per_picoF" public_interface="out" name="i_Kss"/>
      
	  <variable units="milliS_per_microF" name="g_Kss" initial_value="0.015"/>
      <variable units="dimensionless" name="aKss" initial_value="0.8500"/>
      <variable units="millisecond" name="tau_Kss"/>
      
	  <variable units="millivolt" public_interface="in" name="E_K"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="millisecond" public_interface="in" name="time"/>
      <variable units="dimensionless" public_interface="in" name="ass"/>
      
	  <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_Kss</ci>
            <apply>
               <times/>
               <ci>g_Kss</ci>
               <ci>aKss</ci>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>E_K</ci>
               </apply>
            </apply>
         </apply>
         
		 <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>aKss</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>ass</ci>
                  <ci>aKss</ci>
               </apply>
               <ci>tau_Kss</ci>
            </apply>
         </apply>
         
		 <apply>
            <eq/>
            <ci>tau_Kss</ci>
            <apply>
               <plus/>
               <apply>
                  <times/>
                  <cn cellml:units="millisecond">39.3</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <times/>
                        <apply>
                           <minus/>
                           <cn cellml:units="per_millivolt">0.0862</cn>
                        </apply>
                        <ci>V</ci>
                     </apply>
                  </apply>
               </apply>
               <cn cellml:units="millisecond">13.17</cn>
            </apply>
         </apply>
      </math>
   </component>
   
   <component name="time_independent_potassium_current">
      <variable units="picoA_per_picoF" public_interface="out" name="i_K1"/>
      
	  <variable units="milliS_per_microF" name="g_K1" initial_value="0.235"/>
	  
	  <variable units="micromolar" public_interface="in" name="Ko"/>
      <variable units="millivolt" public_interface="in" name="E_K"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="millisecond" public_interface="in" name="time"/>
      
	  <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_K1</ci>
            <apply>
               <times/>
			   <ci>g_K1</ci>
			   <apply>
                  <divide/>
                  <ci>Ko</ci>
                  <apply>
                     <plus/>
                     <ci>Ko</ci>
                     <cn cellml:units="micromolar">210</cn>
                  </apply>
               </apply>
			   <apply>
                  <divide/>
                  <apply>
                     <minus/>
                     <ci>V</ci>
                     <ci>E_K</ci>
                  </apply>
                  <apply>
                     <plus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <times/>
                           <cn cellml:units="per_millivolt">0.0896</cn>
                           <apply>
                              <minus/>
                              <ci>V</ci>
                              <ci>E_K</ci>
						   </apply>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   
   <component name="slow_delayed_rectifier_potassium_current">
      <variable units="picoA_per_picoF" public_interface="out" name="i_Ks"/>
      
	  <variable units="milliS_per_microF" name="g_Ks" initial_value="0.046"/>
      <variable units="dimensionless" name="nKs" initial_value="0.0026"/>
      <variable units="per_millisecond" name="alpha_n"/>
      <variable units="per_millisecond" name="beta_n"/>
      
	  <variable units="millivolt" public_interface="in" name="E_K"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="millisecond" public_interface="in" name="time"/>
      
	  <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_Ks</ci>
            <apply>
               <times/>
               <ci>g_Ks</ci>
               <apply>
                  <power/>
                  <ci>nKs</ci>
                  <cn cellml:units="dimensionless">2</cn>
               </apply>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>E_K</ci>
               </apply>
            </apply>
         </apply>
         
		 <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>nKs</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>alpha_n</ci>
                  <apply>
                     <minus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <ci>nKs</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>beta_n</ci>
                  <ci>nKs</ci>
               </apply>
            </apply>
         </apply>
         
		 <apply>
            <eq/>
            <ci>alpha_n</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <cn cellml:units="per_millivolt_millisecond">0.00000481333</cn>
                  <apply>
                     <plus/>
                     <ci>V</ci>
                     <cn cellml:units="millivolt">26.5</cn>
                  </apply>
               </apply>
               <apply>
                  <minus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <times/>
                        <apply>
                           <minus/>
                           <cn cellml:units="per_millivolt">0.128</cn>
                        </apply>
                        <apply>
                           <plus/>
                           <ci>V</ci>
                           <cn cellml:units="millivolt">26.5</cn>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
       
   	     <apply>
            <eq/>
            <ci>beta_n</ci>
            <apply>
               <times/>
               <cn cellml:units="per_millisecond">0.0000953333</cn>
               <apply>
                  <exp/>
                  <apply>
                     <times/>
                     <apply>
                        <minus/>
                        <cn cellml:units="per_millivolt">0.038</cn>
                     </apply>
                     <apply>
                        <plus/>
                        <ci>V</ci>
                        <cn cellml:units="millivolt">26.5</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
  
   <component name="rapid_delayed_rectifier_potassium_current">
      <variable units="picoA_per_picoF" public_interface="out" name="i_Kr"/>
      
	  <variable units="milliS_per_microF" name="g_Kr" initial_value="01.17"/>
      <variable units="dimensionless" name="O_K" initial_value="2.1e-3"/>
      <variable units="dimensionless" name="C_K1" initial_value="1.1e-3"/>
      <variable units="dimensionless" name="C_K2" initial_value="8.3641e-4"/>
      <variable units="dimensionless" name="C_K0"/>
      <variable units="dimensionless" name="I_K" initial_value="4.3522e-4"/>
      <variable units="per_millisecond" name="alpha_a0"/>
      <variable units="per_millisecond" name="beta_a0"/>
      <variable units="per_millisecond" name="kb" initial_value="0.036778"/>
      <variable units="per_millisecond" name="kf" initial_value="0.023761"/>
      <variable units="per_millisecond" name="alpha_a1"/>
      <variable units="per_millisecond" name="beta_a1"/>
      <variable units="per_millisecond" name="alpha_i"/>
      <variable units="per_millisecond" name="beta_i"/>
      
	  <variable units="millisecond" public_interface="in" name="time"/>
      <variable units="millivolt" public_interface="in" name="V"/>
	  <variable units="millivolt" public_interface="in" name="E_K"/>
      
	  <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_Kr</ci>
            <apply>
               <times/>
               <ci>g_Kr</ci>
               <ci>O_K</ci>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>E_K</ci>
               </apply>
            </apply>
         </apply>
         
		 <apply>
            <eq/>
            <ci>C_K0</ci>
            <apply>
               <minus/>
               <cn cellml:units="dimensionless">1</cn>
               <apply>
                  <plus/>
                  <ci>C_K1</ci>
                  <ci>C_K2</ci>
                  <ci>O_K</ci>
                  <ci>I_K</ci>
               </apply>
            </apply>
         </apply>
         
		 <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>C_K2</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <plus/>
                  <apply>
                     <times/>
                     <ci>kf</ci>
                     <ci>C_K1</ci>
                  </apply>
                  <apply>
                     <times/>
                     <ci>beta_a1</ci>
                     <ci>O_K</ci>
                  </apply>
               </apply>
               <apply>
                  <plus/>
                  <apply>
                     <times/>
                     <ci>kb</ci>
                     <ci>C_K2</ci>
                  </apply>
                  <apply>
                     <times/>
                     <ci>alpha_a1</ci>
                     <ci>C_K2</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         
		 <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>C_K1</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <plus/>
                  <apply>
                     <times/>
                     <ci>alpha_a0</ci>
                     <ci>C_K0</ci>
                  </apply>
                  <apply>
                     <times/>
                     <ci>kb</ci>
                     <ci>C_K2</ci>
                  </apply>
               </apply>
               <apply>
                  <plus/>
                  <apply>
                     <times/>
                     <ci>beta_a0</ci>
                     <ci>C_K1</ci>
                  </apply>
                  <apply>
                     <times/>
                     <ci>kf</ci>
                     <ci>C_K1</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         
		 <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>O_K</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <plus/>
                  <apply>
                     <times/>
                     <ci>alpha_a1</ci>
                     <ci>C_K2</ci>
                  </apply>
                  <apply>
                     <times/>
                     <ci>beta_i</ci>
                     <ci>I_K</ci>
                  </apply>
               </apply>
               <apply>
                  <plus/>
                  <apply>
                     <times/>
                     <ci>beta_a1</ci>
                     <ci>O_K</ci>
                  </apply>
                  <apply>
                     <times/>
                     <ci>alpha_i</ci>
                     <ci>O_K</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         
		 <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>I_K</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>alpha_i</ci>
                  <ci>O_K</ci>
               </apply>
               <apply>
                  <times/>
                  <ci>beta_i</ci>
                  <ci>I_K</ci>
               </apply>
            </apply>
         </apply>
         
		 <apply>
            <eq/>
            <ci>alpha_a0</ci>
            <apply>
               <times/>
               <cn cellml:units="per_millisecond">0.022348</cn>
               <apply>
                  <exp/>
                  <apply>
                     <times/>
                     <cn cellml:units="per_millivolt">0.01176</cn>
                     <ci>V</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         
		 <apply>
            <eq/>
            <ci>beta_a0</ci>
            <apply>
               <times/>
               <cn cellml:units="per_millisecond">0.047002</cn>
               <apply>
                  <exp/>
                  <apply>
                     <times/>
                     <apply>
                        <minus/>
                        <cn cellml:units="per_millivolt">0.0631</cn>
                     </apply>
                     <ci>V</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         
		 <apply>
            <eq/>
            <ci>alpha_a1</ci>
            <apply>
               <times/>
               <cn cellml:units="per_millisecond">0.013733</cn>
               <apply>
                  <exp/>
                  <apply>
                     <times/>
                     <cn cellml:units="per_millivolt">0.038198</cn>
                     <ci>V</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
        
 		 <apply>
            <eq/>
            <ci>beta_a1</ci>
            <apply>
               <times/>
               <cn cellml:units="per_millisecond">0.0000689</cn>
               <apply>
                  <exp/>
                  <apply>
                     <times/>
                     <apply>
                        <minus/>
                        <cn cellml:units="per_millivolt">0.04178</cn>
                     </apply>
                     <ci>V</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         
		 <apply>
            <eq/>
            <ci>alpha_i</ci>
            <apply>
               <times/>
               <cn cellml:units="per_millisecond">0.090821</cn>
               <apply>
                  <exp/>
                  <apply>
                     <times/>
                     <cn cellml:units="per_millivolt">0.023391</cn>
                     <apply>
                        <plus/>
                        <ci>V</ci>
                        <cn cellml:units="millivolt">5</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         
		 <apply>
            <eq/>
            <ci>beta_i</ci>
            <apply>
               <times/>
               <cn cellml:units="per_millisecond">0.006497</cn>
               <apply>
                  <exp/>
                  <apply>
                     <times/>
                     <apply>
                        <minus/>
                        <cn cellml:units="per_millivolt">0.03268</cn>
                     </apply>
                     <apply>
                        <plus/>
                        <ci>V</ci>
                        <cn cellml:units="millivolt">5</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   
   <component name="sodium_potassium_pump_current">
      <variable units="picoA_per_picoF" public_interface="out" name="i_NaK"/>
      
	  <variable units="picoA_per_picoF" name="i_NaK_max" initial_value="0.88"/>
      <variable units="micromolar" name="Km_Nai" initial_value="21000"/>
      <variable units="micromolar" name="Km_Ko" initial_value="1500"/>
      <variable units="dimensionless" name="f_NaK"/>
      <variable units="dimensionless" name="sigma"/>
      
	  <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="joule_per_mole_kelvin" public_interface="in" name="R"/>
      <variable units="kelvin" public_interface="in" name="T"/>
      <variable units="coulomb_per_millimole" public_interface="in" name="F"/>
      <variable units="micromolar" public_interface="in" name="Nai"/>
      <variable units="micromolar" public_interface="in" name="Nao"/>
      <variable units="micromolar" public_interface="in" name="Ko"/>
     
 	  <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_NaK</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <ci>i_NaK_max</ci>
                        <ci>f_NaK</ci>
                        <cn cellml:units="dimensionless">1</cn>
                     </apply>
                     <apply>
                        <plus/>
                        <cn cellml:units="dimensionless">1</cn>
                        <apply>
                           <power/>
                           <apply>
                              <divide/>
                              <ci>Km_Nai</ci>
                              <ci>Nai</ci>
                           </apply>
                           <cn cellml:units="dimensionless">1.5</cn>
                        </apply>
                     </apply>
                  </apply>
                  <ci>Ko</ci>
               </apply>
               <apply>
                  <plus/>
                  <ci>Ko</ci>
                  <ci>Km_Ko</ci>
               </apply>
            </apply>
         </apply>
         
		 <apply>
            <eq/>
            <ci>f_NaK</ci>
            <apply>
               <divide/>
               <cn cellml:units="dimensionless">1</cn>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <times/>
                     <cn cellml:units="dimensionless">0.1245</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <times/>
                              <apply>
                                 <minus/>
                                 <cn cellml:units="dimensionless">0.1</cn>
                              </apply>
                              <ci>V</ci>
                              <ci>F</ci>
                           </apply>
                           <apply>
                              <times/>
                              <ci>R</ci>
                              <ci>T</ci>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
                  <apply>
                     <times/>
                     <cn cellml:units="dimensionless">0.0365</cn>
                     <ci>sigma</ci>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <times/>
                              <apply>
                                 <minus/>
                                 <ci>V</ci>
                              </apply>
                              <ci>F</ci>
                           </apply>
                           <apply>
                              <times/>
                              <ci>R</ci>
                              <ci>T</ci>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         
		 <apply>
            <eq/>
            <ci>sigma</ci>
            <apply>
               <times/>
               <apply>
                  <divide/>
                  <cn cellml:units="dimensionless">1</cn>
                  <cn cellml:units="dimensionless">7</cn>
               </apply>
               <apply>
                  <minus/>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <ci>Nao</ci>
                        <cn cellml:units="micromolar">67300</cn>
                     </apply>
                  </apply>
                  <cn cellml:units="dimensionless">1</cn>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   
   <component name="calcium_activated_chloride_current">
      <variable units="picoA_per_picoF" public_interface="out" name="i_ClCa"/>
      
	  <variable units="cm_per_second" name="P_ClCa" initial_value="2.74e-7"/>
      <variable units="dimensionless" name="f_ClCa"/>
      <variable units="micromolar" name="Km_Cl" initial_value="4"/>
	  <variable units="dimensionless" name="n_ClCa" initial_value="3"/>
      
	  <variable units="millivolt" public_interface="in" name="V"/>
	  <variable units="joule_per_mole_kelvin" public_interface="in" name="R"/>
      <variable units="kelvin" public_interface="in" name="T"/>
      <variable units="coulomb_per_millimole" public_interface="in" name="F"/>
      <variable units="micromolar" public_interface="in" name="Cli"/>
	  <variable units="micromolar" public_interface="in" name="Clo"/>
	  <variable units="micromolar" public_interface="in" name="Cass"/>
      
	  <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>f_ClCa</ci>
            <apply>
               <power/>
			   <apply>
                  <divide/>
				  <ci>Cass</ci>
                  <apply>
                     <plus/>
                     <ci>Km_Cl</ci>
					 <ci>Cass</ci>
				  </apply>
			   </apply>
			   <ci>n_ClCa</ci>
			</apply>
         </apply>
					 
		 <apply>
            <eq/>
            <ci>i_ClCa</ci>
            <apply>
               <times/>
			   <cn cellml:units="cm2_per_microF">1</cn>
			   <ci>P_ClCa</ci>
               <ci>f_ClCa</ci>
               <apply>
                  <divide/>
                  <apply>
                     <times/>
                     <ci>V</ci>
					 <apply>
                        <power/>
                        <ci>F</ci>
						<cn cellml:units="dimensionless">2</cn>
                     </apply>
				  </apply>
				  <apply>
                     <times/>
                     <ci>R</ci>
                     <ci>T</ci>
                  </apply>
               </apply>
			   <apply>
                  <divide/>
				  <apply>
                     <minus/>
                     <apply>
                        <times/>
						<ci>Clo</ci>
						<apply>
                           <exp/>
                           <apply>
                              <divide/>
                              <apply>
                                 <times/>
								 <ci>V</ci>
                                 <ci>F</ci>
                              </apply>
							  <apply>
                                 <times/>
								 <ci>R</ci>
                                 <ci>T</ci>
                              </apply>
						   </apply>
						</apply>
				     </apply>
					 <ci>Cli</ci>
				  </apply>
                  <apply>
                     <minus/>
					 <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <times/>
					          <ci>V</ci>
                              <ci>F</ci>
                           </apply>
						   <apply>
                              <times/>
							  <ci>R</ci>
                              <ci>T</ci>
                           </apply>
						</apply>
				     </apply>
                     <cn cellml:units="dimensionless">1</cn>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   
   <group>
      <relationship_ref relationship="encapsulation"/>
      <component_ref component="fast_sodium_current">
         <component_ref component="fast_sodium_current_m_gate"/>
         <component_ref component="fast_sodium_current_h_gate"/>
         <component_ref component="fast_sodium_current_j_gate"/>
      </component_ref>
      <component_ref component="T_type_calcium_current">
         <component_ref component="T_type_calcium_current_b_gate"/>
         <component_ref component="T_type_calcium_current_g_gate"/>
      </component_ref>
   </group>
   
   <connection>
      <map_components component_2="environment" component_1="membrane"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="environment" component_1="calcium_concentration"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="environment" component_1="calcium_fluxes"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="environment" component_1="calcium_buffering"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="environment" component_1="ryanodine_receptors"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="environment" component_1="L_type_calcium_current"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="environment" component_1="T_type_calcium_current"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="environment" component_1="sodium_concentration"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="environment" component_1="fast_sodium_current"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="environment" component_1="potassium_concentration"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="environment" component_1="fast_transient_outward_potassium_current"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="environment" component_1="slow_delayed_rectifier_potassium_current"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="environment" component_1="time_independent_potassium_current"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="environment" component_1="ultra_rapidly_activating_delayed_rectifier_potassium_current"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="environment" component_1="non_inactivating_steady_state_potassium_current"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="environment" component_1="rapid_delayed_rectifier_potassium_current"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="L_type_calcium_current" component_1="membrane"/>
      <map_variables variable_2="i_CaL" variable_1="i_CaL"/>
      <map_variables variable_2="V" variable_1="V"/>
   </connection>
   <connection>
      <map_components component_2="T_type_calcium_current" component_1="membrane"/>
      <map_variables variable_2="i_CaT" variable_1="i_CaT"/>
      <map_variables variable_2="V" variable_1="V"/>
   </connection>
   <connection>
      <map_components component_2="calcium_pump_current" component_1="membrane"/>
      <map_variables variable_2="i_pCa" variable_1="i_pCa"/>
   </connection>
   <connection>
      <map_components component_2="sodium_calcium_exchange_current" component_1="membrane"/>
      <map_variables variable_2="i_NaCa" variable_1="i_NaCa"/>
      <map_variables variable_2="V" variable_1="V"/>
      <map_variables variable_2="R" variable_1="R"/>
      <map_variables variable_2="T" variable_1="T"/>
      <map_variables variable_2="F" variable_1="F"/>
      <map_variables variable_2="Cao" variable_1="Cao"/>
      <map_variables variable_2="Nao" variable_1="Nao"/>
   </connection>
   <connection>
      <map_components component_2="calcium_background_current" component_1="membrane"/>
      <map_variables variable_2="i_Cab" variable_1="i_Cab"/>
      <map_variables variable_2="V" variable_1="V"/>
      <map_variables variable_2="R" variable_1="R"/>
      <map_variables variable_2="T" variable_1="T"/>
      <map_variables variable_2="F" variable_1="F"/>
      <map_variables variable_2="Cao" variable_1="Cao"/>
   </connection>
   <connection>
      <map_components component_2="fast_sodium_current" component_1="membrane"/>
      <map_variables variable_2="i_Na" variable_1="i_Na"/>
      <map_variables variable_2="V" variable_1="V"/>
      <map_variables variable_2="R" variable_1="R"/>
      <map_variables variable_2="T" variable_1="T"/>
      <map_variables variable_2="F" variable_1="F"/>
      <map_variables variable_2="Nao" variable_1="Nao"/>
   </connection>
   <connection>
      <map_components component_2="sodium_background_current" component_1="membrane"/>
      <map_variables variable_2="i_Nab" variable_1="i_Nab"/>
      <map_variables variable_2="V" variable_1="V"/>
   </connection>
   <connection>
      <map_components component_2="sodium_potassium_pump_current" component_1="membrane"/>
      <map_variables variable_2="i_NaK" variable_1="i_NaK"/>
      <map_variables variable_2="V" variable_1="V"/>
      <map_variables variable_2="R" variable_1="R"/>
      <map_variables variable_2="T" variable_1="T"/>
      <map_variables variable_2="F" variable_1="F"/>
      <map_variables variable_2="Nao" variable_1="Nao"/>
	  <map_variables variable_2="Ko" variable_1="Ko"/>
   </connection>
   <connection>
      <map_components component_2="fast_transient_outward_potassium_current" component_1="membrane"/>
      <map_variables variable_2="i_Kto_f" variable_1="i_Kto_f"/>
      <map_variables variable_2="V" variable_1="V"/>
      <map_variables variable_2="R" variable_1="R"/>
      <map_variables variable_2="T" variable_1="T"/>
      <map_variables variable_2="F" variable_1="F"/>
      <map_variables variable_2="Ko" variable_1="Ko"/>
   </connection>
   <connection>
      <map_components component_2="time_independent_potassium_current" component_1="membrane"/>
      <map_variables variable_2="i_K1" variable_1="i_K1"/>
      <map_variables variable_2="V" variable_1="V"/>
      <map_variables variable_2="Ko" variable_1="Ko"/>
   </connection>
   <connection>
      <map_components component_2="slow_delayed_rectifier_potassium_current" component_1="membrane"/>
      <map_variables variable_2="i_Ks" variable_1="i_Ks"/>
      <map_variables variable_2="V" variable_1="V"/>
   </connection>
   <connection>
      <map_components component_2="ultra_rapidly_activating_delayed_rectifier_potassium_current" component_1="membrane"/>
      <map_variables variable_2="i_Kur" variable_1="i_Kur"/>
      <map_variables variable_2="V" variable_1="V"/>
   </connection>
   <connection>
      <map_components component_2="non_inactivating_steady_state_potassium_current" component_1="membrane"/>
      <map_variables variable_2="i_Kss" variable_1="i_Kss"/>
      <map_variables variable_2="V" variable_1="V"/>
   </connection>
   <connection>
      <map_components component_2="rapid_delayed_rectifier_potassium_current" component_1="membrane"/>
      <map_variables variable_2="i_Kr" variable_1="i_Kr"/>
      <map_variables variable_2="V" variable_1="V"/>
   </connection>
   <connection>
      <map_components component_2="calcium_activated_chloride_current" component_1="membrane"/>
      <map_variables variable_2="i_ClCa" variable_1="i_ClCa"/>
      <map_variables variable_2="V" variable_1="V"/>
	  <map_variables variable_2="R" variable_1="R"/>
      <map_variables variable_2="T" variable_1="T"/>
      <map_variables variable_2="F" variable_1="F"/>
	  <map_variables variable_2="Clo" variable_1="Clo"/>
      <map_variables variable_2="Cli" variable_1="Cli"/>
   </connection>
   <connection>
      <map_components component_2="calcium_concentration" component_1="membrane"/>
      <map_variables variable_2="Cm" variable_1="Cm"/>
      <map_variables variable_2="V" variable_1="V"/>
      <map_variables variable_2="Acap" variable_1="Acap"/>
      <map_variables variable_2="Vmyo" variable_1="Vmyo"/>
      <map_variables variable_2="VNSR" variable_1="VNSR"/>
      <map_variables variable_2="VJSR" variable_1="VJSR"/>
      <map_variables variable_2="Vss" variable_1="Vss"/>
      <map_variables variable_2="F" variable_1="F"/>
   </connection>
   <connection>
      <map_components component_2="calcium_fluxes" component_1="membrane"/>
      <map_variables variable_2="V" variable_1="V"/>
   </connection>
   <connection>
      <map_components component_2="sodium_concentration" component_1="membrane"/>
      <map_variables variable_2="Cm" variable_1="Cm"/>
      <map_variables variable_2="Acap" variable_1="Acap"/>
      <map_variables variable_2="Vmyo" variable_1="Vmyo"/>
      <map_variables variable_2="F" variable_1="F"/>
   </connection>
   <connection>
      <map_components component_2="potassium_concentration" component_1="membrane"/>
      <map_variables variable_2="Cm" variable_1="Cm"/>
      <map_variables variable_2="Acap" variable_1="Acap"/>
      <map_variables variable_2="Vmyo" variable_1="Vmyo"/>
      <map_variables variable_2="F" variable_1="F"/>
   </connection>
   <connection>
      <map_components component_2="calcium_fluxes" component_1="calcium_concentration"/>
      <map_variables variable_2="Cass" variable_1="Cass"/>
      <map_variables variable_2="CaJSR" variable_1="CaJSR"/>
      <map_variables variable_2="Cai" variable_1="Cai"/>
      <map_variables variable_2="CaNSR" variable_1="CaNSR"/>
      <map_variables variable_2="J_rel" variable_1="J_rel"/>
      <map_variables variable_2="J_tr" variable_1="J_tr"/>
      <map_variables variable_2="J_xfer" variable_1="J_xfer"/>
      <map_variables variable_2="J_leak" variable_1="J_leak"/>
      <map_variables variable_2="J_up" variable_1="J_up"/>
      <map_variables variable_2="J_trpn" variable_1="J_trpn"/>
   </connection>
   <connection>
      <map_components component_2="calcium_buffering" component_1="calcium_concentration"/>
      <map_variables variable_2="Cai" variable_1="Cai"/>
   </connection>
   <connection>
      <map_components component_2="ryanodine_receptors" component_1="calcium_concentration"/>
      <map_variables variable_2="Cass" variable_1="Cass"/>
   </connection>
   <connection>
      <map_components component_2="L_type_calcium_current" component_1="calcium_concentration"/>
      <map_variables variable_2="i_CaL" variable_1="i_CaL"/>
      <map_variables variable_2="Cass" variable_1="Cass"/>
   </connection>
   <connection>
      <map_components component_2="T_type_calcium_current" component_1="calcium_concentration"/>
      <map_variables variable_2="i_CaT" variable_1="i_CaT"/>
   </connection>
   <connection>
      <map_components component_2="calcium_pump_current" component_1="calcium_concentration"/>
      <map_variables variable_2="i_pCa" variable_1="i_pCa"/>
      <map_variables variable_2="Cai" variable_1="Cai"/>
   </connection>
   <connection>
      <map_components component_2="sodium_calcium_exchange_current" component_1="calcium_concentration"/>
      <map_variables variable_2="i_NaCa" variable_1="i_NaCa"/>
      <map_variables variable_2="Cai" variable_1="Cai"/>
   </connection>
   <connection>
      <map_components component_2="calcium_background_current" component_1="calcium_concentration"/>
      <map_variables variable_2="i_Cab" variable_1="i_Cab"/>
      <map_variables variable_2="Cai" variable_1="Cai"/>
   </connection>
   <connection>
      <map_components component_2="calcium_activated_chloride_current" component_1="calcium_concentration"/>
      <map_variables variable_2="Cass" variable_1="Cass"/>
   </connection>
   <connection>
      <map_components component_2="calcium_buffering" component_1="calcium_fluxes"/>
      <map_variables variable_2="LTRPN_tot" variable_1="LTRPN_tot"/>
      <map_variables variable_2="HTRPN_tot" variable_1="HTRPN_tot"/>
      <map_variables variable_2="LTRPN_Ca" variable_1="LTRPN_Ca"/>
      <map_variables variable_2="HTRPN_Ca" variable_1="HTRPN_Ca"/>
      <map_variables variable_2="k_plus_htrpn" variable_1="k_plus_htrpn"/>
      <map_variables variable_2="k_minus_htrpn" variable_1="k_minus_htrpn"/>
      <map_variables variable_2="k_plus_ltrpn" variable_1="k_plus_ltrpn"/>
      <map_variables variable_2="k_minus_ltrpn" variable_1="k_minus_ltrpn"/>
   </connection>
   <connection>
      <map_components component_2="L_type_calcium_current" component_1="calcium_fluxes"/>
      <map_variables variable_2="i_CaL_max" variable_1="i_CaL_max"/>
      <map_variables variable_2="i_CaL" variable_1="i_CaL"/>
   </connection>
   <connection>
      <map_components component_2="ryanodine_receptors" component_1="calcium_fluxes"/>
      <map_variables variable_2="P_O1" variable_1="P_O1"/>
      <map_variables variable_2="P_O2" variable_1="P_O2"/>
   </connection>
   <connection>
      <map_components component_2="fast_sodium_current" component_1="sodium_concentration"/>
      <map_variables variable_2="Nai" variable_1="Nai"/>
      <map_variables variable_2="i_Na" variable_1="i_Na"/>
   </connection>
   <connection>
      <map_components component_2="sodium_background_current" component_1="sodium_concentration"/>
      <map_variables variable_2="i_Nab" variable_1="i_Nab"/>
   </connection>
   <connection>
      <map_components component_2="sodium_calcium_exchange_current" component_1="sodium_concentration"/>
      <map_variables variable_2="Nai" variable_1="Nai"/>
      <map_variables variable_2="i_NaCa" variable_1="i_NaCa"/>
   </connection>
   <connection>
      <map_components component_2="sodium_potassium_pump_current" component_1="sodium_concentration"/>
      <map_variables variable_2="Nai" variable_1="Nai"/>
      <map_variables variable_2="i_NaK" variable_1="i_NaK"/>
   </connection>
   <connection>
      <map_components component_2="sodium_background_current" component_1="fast_sodium_current"/>
      <map_variables variable_2="E_Na" variable_1="E_Na"/>
   </connection>
   <connection>
      <map_components component_2="fast_transient_outward_potassium_current" component_1="potassium_concentration"/>
      <map_variables variable_2="Ki" variable_1="Ki"/>
      <map_variables variable_2="i_Kto_f" variable_1="i_Kto_f"/>
   </connection>
   <connection>
      <map_components component_2="time_independent_potassium_current" component_1="potassium_concentration"/>
      <map_variables variable_2="i_K1" variable_1="i_K1"/>
   </connection>
   <connection>
      <map_components component_2="slow_delayed_rectifier_potassium_current" component_1="potassium_concentration"/>
      <map_variables variable_2="i_Ks" variable_1="i_Ks"/>
   </connection>
   <connection>
      <map_components component_2="rapid_delayed_rectifier_potassium_current" component_1="potassium_concentration"/>
      <map_variables variable_2="i_Kr" variable_1="i_Kr"/>
   </connection>
   <connection>
      <map_components component_2="ultra_rapidly_activating_delayed_rectifier_potassium_current" component_1="potassium_concentration"/>
      <map_variables variable_2="i_Kur" variable_1="i_Kur"/>
   </connection>
   <connection>
      <map_components component_2="sodium_potassium_pump_current" component_1="potassium_concentration"/>
      <map_variables variable_2="i_NaK" variable_1="i_NaK"/>
   </connection>
   <connection>
      <map_components component_2="non_inactivating_steady_state_potassium_current" component_1="ultra_rapidly_activating_delayed_rectifier_potassium_current"/>
      <map_variables variable_2="ass" variable_1="ass"/>
   </connection>
   <connection>
      <map_components component_2="non_inactivating_steady_state_potassium_current" component_1="potassium_concentration"/>
      <map_variables variable_2="i_Kss" variable_1="i_Kss"/>
   </connection>
   <connection>
      <map_components component_2="time_independent_potassium_current" component_1="fast_transient_outward_potassium_current"/>
      <map_variables variable_2="E_K" variable_1="E_K"/>
   </connection>
   <connection>
      <map_components component_2="slow_delayed_rectifier_potassium_current" component_1="fast_transient_outward_potassium_current"/>
      <map_variables variable_2="E_K" variable_1="E_K"/>
   </connection>
   <connection>
      <map_components component_2="ultra_rapidly_activating_delayed_rectifier_potassium_current" component_1="fast_transient_outward_potassium_current"/>
      <map_variables variable_2="E_K" variable_1="E_K"/>
   </connection>
   <connection>
      <map_components component_2="non_inactivating_steady_state_potassium_current" component_1="fast_transient_outward_potassium_current"/>
      <map_variables variable_2="E_K" variable_1="E_K"/>
   </connection>
   <connection>
      <map_components component_2="rapid_delayed_rectifier_potassium_current" component_1="fast_transient_outward_potassium_current"/>
      <map_variables variable_2="E_K" variable_1="E_K"/>
   </connection>
   <connection>
      <map_components component_2="fast_sodium_current_m_gate" component_1="fast_sodium_current"/>
      <map_variables variable_2="m" variable_1="m"/>
      <map_variables variable_2="time" variable_1="time"/>
      <map_variables variable_2="V" variable_1="V"/>
   </connection>
   <connection>
      <map_components component_2="fast_sodium_current_h_gate" component_1="fast_sodium_current"/>
      <map_variables variable_2="h" variable_1="h"/>
      <map_variables variable_2="time" variable_1="time"/>
      <map_variables variable_2="V" variable_1="V"/>
   </connection>
   <connection>
      <map_components component_2="fast_sodium_current_j_gate" component_1="fast_sodium_current"/>
      <map_variables variable_2="j" variable_1="j"/>
      <map_variables variable_2="time" variable_1="time"/>
      <map_variables variable_2="V" variable_1="V"/>
   </connection>
   <connection>
      <map_components component_2="T_type_calcium_current_b_gate" component_1="T_type_calcium_current"/>
      <map_variables variable_2="b" variable_1="b"/>
      <map_variables variable_2="time" variable_1="time"/>
      <map_variables variable_2="V" variable_1="V"/>
   </connection>
   <connection>
      <map_components component_2="T_type_calcium_current_g_gate" component_1="T_type_calcium_current"/>
      <map_variables variable_2="g" variable_1="g"/>
      <map_variables variable_2="time" variable_1="time"/>
      <map_variables variable_2="V" variable_1="V"/>
   </connection>

<rdf:RDF>
  <rdf:Seq rdf:about="rdf:#0b6e6fb8-9e67-4163-8f2a-125b2973d05e">
    <rdf:li rdf:resource="rdf:#6685a93f-c980-4fa3-b78e-cd112c11390e"/>
    <rdf:li rdf:resource="rdf:#19e7066e-d983-4408-b59a-bd31a1ef5189"/>
  </rdf:Seq>
  <rdf:Description rdf:about="rdf:#4d0294d5-dd3e-4677-88c1-b641f2d1d3a8">
    <dcterms:W3CDTF>2009-03-02T00:00:00+00:00</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#19e7066e-d983-4408-b59a-bd31a1ef5189">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#7d52f800-28d5-464b-9904-191c0468798c"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#d77fcb00-54a6-4792-b702-dde28a037fd5">
    <vCard:FN/>
  </rdf:Description>
  <rdf:Description rdf:about="">
    <dc:publisher/>
    <cmeta:comment rdf:resource="rdf:#50781470-62a2-46f6-a423-b72a35b34587"/>
    <dcterms:created rdf:resource="rdf:#4d0294d5-dd3e-4677-88c1-b641f2d1d3a8"/>
    <dc:creator rdf:resource="rdf:#caf953e9-03ae-44e2-a1f3-5edbd2d3c6c7"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#50781470-62a2-46f6-a423-b72a35b34587">
    <dc:creator rdf:resource="rdf:#6c03d419-e71e-403d-9074-48fa3afe93de"/>
    <rdf:value>This CellML model is known to run in PCEnv and COR to recreate the published results. The units have been checked and are consistent.</rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#6685a93f-c980-4fa3-b78e-cd112c11390e">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#8ea3fb69-ae8b-4e35-80e4-372ce3abed38"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#632aec73-3b6c-4e1c-a245-49a7b52af31d">
    <vCard:Orgname>The University of Auckland</vCard:Orgname>
    <vCard:Orgunit>Auckland Bioengineering Institute</vCard:Orgunit>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#6c03d419-e71e-403d-9074-48fa3afe93de">
    <vCard:FN>Catherine Lloyd</vCard:FN>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#09c152f3-5f44-41ce-bad8-175e6f27bf86">
    <dc:creator rdf:resource="rdf:#d77fcb00-54a6-4792-b702-dde28a037fd5"/>
    <rdf:value/>
  </rdf:Description>
  <rdf:Description rdf:about="#wang_model_2008">
    <bqs:reference rdf:resource="rdf:#d0769982-0d4d-4d65-a54e-462d65917001"/>
<bqs:reference rdf:parseType="Resource">
  <dc:subject rdf:parseType="Resource">
    <bqs:subject_type>keyword</bqs:subject_type>
    <rdf:value>
      <rdf:Bag>
        <rdf:li>electrophysiology</rdf:li>
        <rdf:li>calcium dynamics</rdf:li>
        <rdf:li>myocyte</rdf:li>
  	<rdf:li>mouse</rdf:li>
 	<rdf:li>myocardium</rdf:li>
      </rdf:Bag>
    </rdf:value>
  </dc:subject>
</bqs:reference>
    <cmeta:comment rdf:resource="rdf:#09c152f3-5f44-41ce-bad8-175e6f27bf86"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#af27a0fa-7d17-43df-a49b-7e01d8a1c114">
    <dcterms:W3CDTF>2008-04-11 00:00</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#7d52f800-28d5-464b-9904-191c0468798c">
    <vCard:Given>Eric</vCard:Given>
    <vCard:Family>Sobie</vCard:Family>
    <vCard:Other>A</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#1c62e8db-6d60-4a8b-828d-016d4b882233">
    <dc:title>American Journal of Physiology Heart and Circulatory Physiology</dc:title>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#caf953e9-03ae-44e2-a1f3-5edbd2d3c6c7">
    <vCard:ORG rdf:resource="rdf:#632aec73-3b6c-4e1c-a245-49a7b52af31d"/>
    <vCard:EMAIL rdf:resource="rdf:#d25fefba-6be1-4fc9-80ef-ade4d3798d48"/>
    <vCard:N rdf:resource="rdf:#7cd164da-489d-43b3-ac26-575850200ade"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#8ea3fb69-ae8b-4e35-80e4-372ce3abed38">
    <vCard:Given>Linda</vCard:Given>
    <vCard:Family>Wang</vCard:Family>
    <vCard:Other>J</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#d25fefba-6be1-4fc9-80ef-ade4d3798d48">
    <rdf:type rdf:resource="http://imc.org/vCard/3.0#internet"/>
    <rdf:value>c.lloyd@auckland.ac.nz</rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#d0769982-0d4d-4d65-a54e-462d65917001">
    <bqs:Pubmed_id>18408122</bqs:Pubmed_id>
    <bqs:JournalArticle rdf:resource="rdf:#9f4e73fd-4d59-471e-a9d5-4efb22972e41"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#9f4e73fd-4d59-471e-a9d5-4efb22972e41">
    <dc:creator rdf:resource="rdf:#0b6e6fb8-9e67-4163-8f2a-125b2973d05e"/>
    <dc:title>Mathematical model of the neonatal mouse ventricular action potential</dc:title>
    <bqs:volume>294</bqs:volume>
    <bqs:first_page>H2565</bqs:first_page>
    <bqs:Journal rdf:resource="rdf:#1c62e8db-6d60-4a8b-828d-016d4b882233"/>
    <dcterms:issued rdf:resource="rdf:#af27a0fa-7d17-43df-a49b-7e01d8a1c114"/>
    <bqs:last_page>H2575</bqs:last_page>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#7cd164da-489d-43b3-ac26-575850200ade">
    <vCard:Given>Catherine</vCard:Given>
    <vCard:Family>Lloyd</vCard:Family>
    <vCard:Other>May</vCard:Other>
  </rdf:Description>
</rdf:RDF>
</model>