Rendering of the source text

<?xml version="1.0" encoding="utf-8"?>
<!--
This CellML file was generated on 02/11/2009 at 14:48:14 using:

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

CellML 1.0 was used to generate this model
http://www.cellml.org/
-->
<model 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:dcterms="http://purl.org/dc/terms/" xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" name="inada_NH_model_2009" cmeta:id="inada_NH_model_2009" xmlns="http://www.cellml.org/cellml/1.0#" xmlns:cellml="http://www.cellml.org/cellml/1.0#" xmlns:cmeta="http://www.cellml.org/metadata/1.0#">
	<documentation xmlns="http://cellml.org/tmp-documentation">
		<article>
			<articleinfo>
				<title>One-dimensional mathematical model of the atrioventricular node including atrio-nodal, nodal, and nodal-his cells</title>
				<author>
					<firstname>Catherine</firstname>
					<surname>Lloyd</surname>
					<affiliation>
						<shortaffil>Bioengineering Institute, University of Auckland</shortaffil>
					</affiliation>
				</author>
			</articleinfo>
			<section id="sec_status">
				<title>Model Status</title>
				<para>
				This CellML model runs in both COR and PCEnv to reproduce the published results.  The units have been checked and they balance.  This version of the CellML model represents the nodal-His (NH) cell.  The CellML model has been based on both the published paper and the original C++ code.
				</para>
			</section>
			<sect1 id="sec_structure">
				<title>Model Structure</title>
				<para>
ABSTRACT: Mathematical models are a repository of knowledge as well as research and teaching tools. Although action potential models have been developed for most regions of the heart, there is no model for the atrioventricular node (AVN). We have developed action potential models for single atrio-nodal, nodal, and nodal-His cells. The models have the same action potential shapes and refractoriness as observed in experiments. Using these models, together with models for the sinoatrial node (SAN) and atrial muscle, we have developed a one-dimensional (1D) multicellular model including the SAN and AVN. The multicellular model has slow and fast pathways into the AVN and using it we have analyzed the rich behavior of the AVN. Under normal conditions, action potentials were initiated in the SAN center and then propagated through the atrium and AVN. The relationship between the AVN conduction time and the timing of a premature stimulus (conduction curve) is consistent with experimental data. After premature stimulation, atrioventricular nodal reentry could occur. After slow pathway ablation or block of the L-type Ca(2+) current, atrioventricular nodal reentry was abolished. During atrial fibrillation, the AVN limited the number of action potentials transmitted to the ventricle. In the absence of SAN pacemaking, the inferior nodal extension acted as the pacemaker. In conclusion, we have developed what we believe is the first detailed mathematical model of the AVN and it shows the typical physiological and pathophysiological characteristics of the tissue. The model can be used as a tool to analyze the complex structure and behavior of the AVN.
</para>
				<informalfigure float="0" id="fig_reaction_diagram">
					<mediaobject>
						<imageobject>
							<objectinfo>
								<title>model diagram</title>
							</objectinfo>
							<imagedata fileref="inada_2009.png"/>
						</imageobject>
					</mediaobject>
					<caption>Schematic diagram of the Inada et al. 2009 cell model.</caption>
				</informalfigure>
				<para>
The original paper reference is cited below:
</para>
				<para>
One-dimensional mathematical model of the atrioventricular node including atrio-nodal, nodal, and nodal-his cells, S. Inada, J.C. Hancox, H. Zhang and M.R. Boyett, 2009, <emphasis>Biophysical Journal</emphasis>, 97, 2117-2127.  <ulink url="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=pubmed&amp;cmd=Retrieve&amp;dopt=AbstractPlus&amp;list_uids=19843444&amp;query_hl=1&amp;itool=pubmed_docsum">PubMed ID: 19843444</ulink>
				</para>
			</sect1>
		</article>
	</documentation>
	<rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#" xmlns:bqs="http://www.cellml.org/bqs/1.0#" xmlns:dc="http://purl.org/dc/elements/1.1/" xmlns:dcterms="http://purl.org/dc/terms/" xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">
		<foo xmlns="http://www.example.org/foo/bar/baz"/>
		<rdf:Description rdf:about="">
			<dc:title>
            One-dimensional mathematical model of the atrioventricular node including atrio-nodal, nodal, and nodal-his cells (Nodal-His Cell)
          </dc:title>
			<dc:creator rdf:parseType="Resource">
				<vCard:N rdf:parseType="Resource">
					<vCard:Family>Noble</vCard:Family>
					<vCard:Given>Penny</vCard:Given>
				</vCard:N>
				<vCard:EMAIL rdf:parseType="Resource">
					<rdf:value>penny.noble@dpag.ox.ac.uk</rdf:value>
					<rdf:type rdf:resource="http://imc.org/vCard/3.0#internet"/>
				</vCard:EMAIL>
				<vCard:ORG rdf:parseType="Resource">
					<vCard:Orgname>The University of Oxford</vCard:Orgname>
				</vCard:ORG>
			</dc:creator>
			<dcterms:created rdf:parseType="Resource">
				<dcterms:W3CDTF>2009-10-26</dcterms:W3CDTF>
			</dcterms:created>
		</rdf:Description>
		<rdf:Description rdf:about="#inada_NH_model_2009">
			<dc:title>
        The Inada et al. 2009 one-dimensional mathematical model of the atrioventricular node including atrio-nodal, nodal, and nodal-his cells
      </dc:title>
			<cmeta:comment xmlns:cmeta="http://www.cellml.org/metadata/1.0#" rdf:parseType="Resource">
				<rdf:value>
          This is the CellML description of Inada et al.'s one-dimensional mathematical model of the atrioventricular node including atrio-nodal, nodal, and nodal-his cells
        </rdf:value>
				<dc:creator rdf:parseType="Resource">
					<vCard:FN>Catherine Lloyd</vCard:FN>
				</dc:creator>
			</cmeta:comment>
			<cmeta:species xmlns:cmeta="http://www.cellml.org/metadata/1.0#">rabbit</cmeta:species>
			<cmeta:bio_entity xmlns:cmeta="http://www.cellml.org/metadata/1.0#">atrioventricular node</cmeta:bio_entity>
			<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>cardiac</rdf:li>
							<rdf:li>atrioventricular node</rdf:li>
							<rdf:li>atrial fibrillation</rdf:li>
						</rdf:Bag>
					</rdf:value>
				</dc:subject>
			</bqs:reference>
			<bqs:reference rdf:parseType="Resource">
				<bqs:Pubmed_id>19843444</bqs:Pubmed_id>
				<bqs:JournalArticle rdf:parseType="Resource">
					<dc:creator>
						<rdf:Seq>
							<rdf:li rdf:parseType="Resource">
								<bqs:Person rdf:parseType="Resource">
									<vCard:N rdf:parseType="Resource">
										<vCard:Family>Inada</vCard:Family>
										<vCard:Given>S</vCard:Given>
									</vCard:N>
								</bqs:Person>
							</rdf:li>
							<rdf:li rdf:parseType="Resource">
								<bqs:Person rdf:parseType="Resource">
									<vCard:N rdf:parseType="Resource">
										<vCard:Family>Hancox</vCard:Family>
										<vCard:Given>J</vCard:Given>
										<vCard:Other>C</vCard:Other>
									</vCard:N>
								</bqs:Person>
							</rdf:li>
							<rdf:li rdf:parseType="Resource">
								<bqs:Person rdf:parseType="Resource">
									<vCard:N rdf:parseType="Resource">
										<vCard:Family>Zhang</vCard:Family>
										<vCard:Given>H</vCard:Given>
									</vCard:N>
								</bqs:Person>
							</rdf:li>
							<rdf:li rdf:parseType="Resource">
								<bqs:Person rdf:parseType="Resource">
									<vCard:N rdf:parseType="Resource">
										<vCard:Family>Boyett</vCard:Family>
										<vCard:Given>M</vCard:Given>
										<vCard:Other>R</vCard:Other>
									</vCard:N>
								</bqs:Person>
							</rdf:li>
						</rdf:Seq>
					</dc:creator>
					<dc:title>
            One-dimensional mathematical model of the atrioventricular node including atrio-nodal, nodal, and nodal-his cells
          </dc:title>
					<dcterms:issued rdf:parseType="Resource">
						<dcterms:W3CDTF>2009</dcterms:W3CDTF>
					</dcterms:issued>
					<bqs:Journal rdf:parseType="Resource">
						<dc:title>Biophysical Journal</dc:title>
					</bqs:Journal>
					<bqs:volume>97</bqs:volume>
					<bqs:first_page>2117</bqs:first_page>
					<bqs:last_page>2127</bqs:last_page>
				</bqs:JournalArticle>
			</bqs:reference>
		</rdf:Description>
	</rdf:RDF>   <units name="per_second">
      <unit units="second" exponent="-1"/>
   </units>
   <units name="millivolt">
      <unit units="volt" prefix="milli"/>
   </units>
   <units name="millivolt2">
      <unit units="millivolt" exponent="2"/>
   </units>
   <units name="per_millivolt">
      <unit units="volt" prefix="milli" exponent="-1"/>
   </units>
   <units name="per_millivolt_second">
      <unit units="millivolt" exponent="-1"/>
      <unit units="second" exponent="-1"/>
   </units>
   <units name="microS">
      <unit units="siemens" prefix="micro"/>
   </units>
   <units name="microF">
      <unit units="farad" prefix="micro"/>
   </units>
   <units name="nanoA">
      <unit units="ampere" prefix="nano"/>
   </units>
   <units name="millimolar">
      <unit units="mole" prefix="milli"/>
      <unit units="litre" exponent="-1"/>
   </units>
   <units name="millimolar_per_second">
      <unit units="millimolar"/>
      <unit units="second" exponent="-1"/>
   </units>
   <units name="per_millimolar_per_second">
      <unit units="millimolar" exponent="-1"/>
      <unit units="second" exponent="-1"/>
   </units>
   <units name="joule_per_kilomole_kelvin">
      <unit units="joule"/>
      <unit units="mole" prefix="kilo" exponent="-1"/>
      <unit units="kelvin" exponent="-1"/>
   </units>
   <units name="coulomb_per_mole">
      <unit units="coulomb"/>
      <unit units="mole" exponent="-1"/>
   </units>
   <units name="micrometre3">
      <unit units="metre" prefix="micro" exponent="3"/>
   </units>
   <units name="microlitre_micrometre3">
      <unit units="litre" prefix="micro"/>
      <unit units="micrometre3" exponent="-1"/>
   </units>
   <units name="microlitre_per_second">
      <unit units="litre" prefix="micro"/>
      <unit units="second" exponent="-1"/>
   </units>
   <component name="environment">
      <variable name="time" units="second" public_interface="out"/>
   </component>
   <component name="membrane">
      <variable name="V" units="millivolt" initial_value="-69.760276376489" public_interface="out"/>
      <variable name="R" units="joule_per_kilomole_kelvin" initial_value="8314.472" public_interface="out"/>
      <variable name="T" units="kelvin" initial_value="310" public_interface="out"/>
      <variable name="F" units="coulomb_per_mole" initial_value="96485.3415" public_interface="out"/>
      <variable name="C" units="microF" initial_value="2.9e-5"/>
      <variable name="RTONF" units="millivolt" public_interface="out"/>
      <variable name="time" units="second" public_interface="in"/>
      <variable name="i_Na" units="nanoA" public_interface="in"/>
      <variable name="i_CaL" units="nanoA" public_interface="in"/>
      <variable name="i_to" units="nanoA" public_interface="in"/>
      <variable name="i_Kr" units="nanoA" public_interface="in"/>
      <variable name="i_f" units="nanoA" public_interface="in"/>
      <variable name="i_st" units="nanoA" public_interface="in"/>
      <variable name="i_K1" units="nanoA" public_interface="in"/>
      <variable name="i_NaCa" units="nanoA" public_interface="in"/>
      <variable name="i_p" units="nanoA" public_interface="in"/>
      <variable name="i_b" units="nanoA" public_interface="in"/>
      <variable name="i_ACh" units="nanoA" public_interface="in"/>
      <variable name="i_Stim" units="nanoA"/>
      <variable name="stim_start" units="second" initial_value="0.1"/>
      <variable name="stim_end" units="second" initial_value="99999"/>
      <variable name="stim_period" units="second" initial_value="1"/>
      <variable name="stim_duration" units="second" initial_value="0.001"/>
      <variable name="stim_amplitude" units="nanoA" initial_value="-2"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_Stim</ci>
            <piecewise>
               <piece>
                  <ci>stim_amplitude</ci>
                  <apply>
                     <and/>
                     <apply>
                        <geq/>
                        <ci>time</ci>
                        <ci>stim_start</ci>
                     </apply>
                     <apply>
                        <leq/>
                        <ci>time</ci>
                        <ci>stim_end</ci>
                     </apply>
                     <apply>
                        <leq/>
                        <apply>
                           <minus/>
                           <apply>
                              <minus/>
                              <ci>time</ci>
                              <ci>stim_start</ci>
                           </apply>
                           <apply>
                              <times/>
                              <apply>
                                 <floor/>
                                 <apply>
                                    <divide/>
                                    <apply>
                                       <minus/>
                                       <ci>time</ci>
                                       <ci>stim_start</ci>
                                    </apply>
                                    <ci>stim_period</ci>
                                 </apply>
                              </apply>
                              <ci>stim_period</ci>
                           </apply>
                        </apply>
                        <ci>stim_duration</ci>
                     </apply>
                  </apply>
               </piece>
               <otherwise>
                  <cn cellml:units="nanoA">0</cn>
               </otherwise>
            </piecewise>
         </apply>
         <apply>
            <eq/>
            <ci>RTONF</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <ci>R</ci>
                  <ci>T</ci>
               </apply>
               <ci>F</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>V</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <apply>
                     <plus/>
                     <ci>i_Na</ci>
                     <ci>i_CaL</ci>
                     <ci>i_to</ci>
                     <ci>i_Kr</ci>
                     <ci>i_f</ci>
                     <ci>i_st</ci>
                     <ci>i_K1</ci>
                     <ci>i_NaCa</ci>
                     <ci>i_p</ci>
                     <ci>i_b</ci>
                     <ci>i_ACh</ci>
                     <ci>i_Stim</ci>
                  </apply>
               </apply>
               <ci>C</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="hyperpolarising_activated_current">
      <variable name="i_f" units="nanoA" public_interface="out"/>
      <variable name="g_f" units="microS" initial_value="0"/>
      <variable name="time" units="second" public_interface="in" private_interface="out"/>
      <variable name="V" units="millivolt" public_interface="in" private_interface="out"/>
      <variable name="ACh" units="millimolar" public_interface="in" private_interface="out"/>
      <variable name="y" units="dimensionless" private_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_f</ci>
            <apply>
               <times/>
               <ci>y</ci>
               <ci>g_f</ci>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <apply>
                     <minus/>
                     <cn cellml:units="millivolt">30</cn>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="hyperpolarising_activated_current_y_gate">
      <variable name="y" units="dimensionless" initial_value="0.19584111039096" public_interface="out"/>
      <variable name="y_inf" units="dimensionless"/>
      <variable name="tau_y" units="second"/>
      <variable name="V" units="millivolt" public_interface="in"/>
      <variable name="time" units="second" public_interface="in"/>
      <variable name="ACh" units="millimolar" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>y_inf</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">83.19</cn>
                           </apply>
                           <apply>
                              <divide/>
                              <apply>
                                 <times/>
                                 <apply>
                                    <minus/>
                                    <cn cellml:units="millivolt">7.2</cn>
                                 </apply>
                                 <apply>
                                    <power/>
                                    <ci>ACh</ci>
                                    <cn cellml:units="dimensionless">0.69</cn>
                                 </apply>
                              </apply>
                              <apply>
                                 <plus/>
                                 <apply>
                                    <power/>
                                    <cn type="e-notation" cellml:units="millimolar">1.26   <sep/>
                                    -5</cn>
                                    <cn cellml:units="dimensionless">0.69</cn>
                                 </apply>
                                 <apply>
                                    <power/>
                                    <ci>ACh</ci>
                                    <cn cellml:units="dimensionless">0.69</cn>
                                 </apply>
                              </apply>
                           </apply>
                        </apply>
                        <cn cellml:units="millivolt">13.56</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>tau_y</ci>
            <apply>
               <plus/>
               <cn cellml:units="second">0.25</cn>
               <apply>
                  <times/>
                  <cn cellml:units="second">2</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <minus/>
                           <apply>
                              <power/>
                              <apply>
                                 <plus/>
                                 <ci>V</ci>
                                 <cn cellml:units="millivolt">70</cn>
                              </apply>
                              <cn cellml:units="dimensionless">2</cn>
                           </apply>
                        </apply>
                        <cn cellml:units="millivolt2">500</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>y</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>y_inf</ci>
                  <ci>y</ci>
               </apply>
               <ci>tau_y</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="rapid_delayed_rectifier_potassium_current">
      <variable name="i_Kr" units="nanoA" public_interface="out"/>
      <variable name="g_Kr" units="microS" initial_value="0.002"/>
      <variable name="time" units="second" public_interface="in" private_interface="out"/>
      <variable name="V" units="millivolt" public_interface="in" private_interface="out"/>
      <variable name="E_K" units="millivolt" public_interface="out"/>
      <variable name="RTONF" units="millivolt" public_interface="in"/>
      <variable name="Ki" units="millimolar" public_interface="in"/>
      <variable name="Kc" units="millimolar" public_interface="in"/>
      <variable name="paf" units="dimensionless" private_interface="in"/>
      <variable name="pas" units="dimensionless" private_interface="in"/>
      <variable name="pik" units="dimensionless" private_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>E_K</ci>
            <apply>
               <times/>
               <ci>RTONF</ci>
               <apply>
                  <ln/>
                  <apply>
                     <divide/>
                     <ci>Kc</ci>
                     <ci>Ki</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>i_Kr</ci>
            <apply>
               <times/>
               <ci>g_Kr</ci>
               <apply>
                  <plus/>
                  <apply>
                     <times/>
                     <cn cellml:units="dimensionless">0.9</cn>
                     <ci>paf</ci>
                  </apply>
                  <apply>
                     <times/>
                     <cn cellml:units="dimensionless">0.1</cn>
                     <ci>pas</ci>
                  </apply>
               </apply>
               <ci>pik</ci>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>E_K</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="rapid_delayed_rectifier_potassium_current_paf_gate">
      <variable name="paf" units="dimensionless" initial_value="0.00141762995766447" public_interface="out"/>
      <variable name="paf_infinity" units="dimensionless"/>
      <variable name="tau_paf" units="second"/>
      <variable name="time" units="second" public_interface="in"/>
      <variable name="V" units="millivolt" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>paf_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">10.22</cn>
                        </apply>
                        <apply>
                           <minus/>
                           <cn cellml:units="millivolt">8.5</cn>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>tau_paf</ci>
            <apply>
               <divide/>
               <cn cellml:units="dimensionless">1</cn>
               <apply>
                  <plus/>
                  <apply>
                     <times/>
                     <cn cellml:units="per_second">17</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <times/>
                           <cn cellml:units="per_millivolt">0.0398</cn>
                           <ci>V</ci>
                        </apply>
                     </apply>
                  </apply>
                  <apply>
                     <times/>
                     <cn cellml:units="per_second">0.211</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <times/>
                           <apply>
                              <minus/>
                              <cn cellml:units="per_millivolt">0.051</cn>
                           </apply>
                           <ci>V</ci>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>paf</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>paf_infinity</ci>
                  <ci>paf</ci>
               </apply>
               <ci>tau_paf</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="rapid_delayed_rectifier_potassium_current_pas_gate">
      <variable name="pas" units="dimensionless" initial_value="0.00539771950846456" public_interface="out"/>
      <variable name="pas_infinity" units="dimensionless"/>
      <variable name="tau_pas" units="second"/>
      <variable name="time" units="second" public_interface="in"/>
      <variable name="V" units="millivolt" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>pas_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">10.22</cn>
                        </apply>
                        <apply>
                           <minus/>
                           <cn cellml:units="millivolt">8.5</cn>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>tau_pas</ci>
            <apply>
               <plus/>
               <cn cellml:units="second">0.33581</cn>
               <apply>
                  <times/>
                  <cn cellml:units="second">0.90673</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <minus/>
                           <apply>
                              <power/>
                              <apply>
                                 <plus/>
                                 <ci>V</ci>
                                 <cn cellml:units="millivolt">10</cn>
                              </apply>
                              <cn cellml:units="dimensionless">2</cn>
                           </apply>
                        </apply>
                        <cn cellml:units="millivolt2">988.05</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>pas</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>pas_infinity</ci>
                  <ci>pas</ci>
               </apply>
               <ci>tau_pas</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="rapid_delayed_rectifier_potassium_current_pik_gate">
      <variable name="pik" units="dimensionless" initial_value="0.98638204514681" public_interface="out"/>
      <variable name="pik_infinity" units="dimensionless"/>
      <variable name="alpha_pik" units="per_second"/>
      <variable name="beta_pik" units="per_second"/>
      <variable name="tau_pik" units="second"/>
      <variable name="time" units="second" public_interface="in"/>
      <variable name="V" units="millivolt" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>pik_infinity</ci>
            <apply>
               <times/>
               <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">4.9</cn>
                           </apply>
                           <cn cellml:units="millivolt">15.14</cn>
                        </apply>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <minus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <times/>
                     <cn cellml:units="dimensionless">0.3</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <minus/>
                              <apply>
                                 <power/>
                                 <ci>V</ci>
                                 <cn cellml:units="dimensionless">2</cn>
                              </apply>
                           </apply>
                           <cn cellml:units="millivolt2">500</cn>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>alpha_pik</ci>
            <apply>
               <times/>
               <cn cellml:units="per_second">92.01</cn>
               <apply>
                  <exp/>
                  <apply>
                     <times/>
                     <apply>
                        <minus/>
                        <cn cellml:units="per_millivolt">0.0183</cn>
                     </apply>
                     <ci>V</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>beta_pik</ci>
            <apply>
               <times/>
               <cn cellml:units="per_second">603.6</cn>
               <apply>
                  <exp/>
                  <apply>
                     <times/>
                     <cn cellml:units="per_millivolt">0.00942</cn>
                     <ci>V</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>tau_pik</ci>
            <apply>
               <divide/>
               <cn cellml:units="dimensionless">1</cn>
               <apply>
                  <plus/>
                  <ci>alpha_pik</ci>
                  <ci>beta_pik</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>pik</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>pik_infinity</ci>
                  <ci>pik</ci>
               </apply>
               <ci>tau_pik</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="time_independent_potassium_current">
      <variable name="i_K1" units="nanoA" public_interface="out"/>
      <variable name="g_K1" units="microS" initial_value="0.015"/>
      <variable name="g_K1_prime" units="microS"/>
      <variable name="V" units="millivolt" public_interface="in"/>
      <variable name="Kc" units="millimolar" public_interface="in"/>
      <variable name="RTONF" units="millivolt" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>g_K1_prime</ci>
            <apply>
               <times/>
               <ci>g_K1</ci>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">0.5</cn>
                  <apply>
                     <divide/>
                     <cn cellml:units="dimensionless">0.5</cn>
                     <apply>
                        <plus/>
                        <cn cellml:units="dimensionless">1</cn>
                        <apply>
                           <exp/>
                           <apply>
                              <divide/>
                              <apply>
                                 <plus/>
                                 <ci>V</ci>
                                 <cn cellml:units="millivolt">30</cn>
                              </apply>
                              <cn cellml:units="millivolt">5</cn>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>i_K1</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <ci>g_K1_prime</ci>
                  <apply>
                     <power/>
                     <apply>
                        <divide/>
                        <ci>Kc</ci>
                        <apply>
                           <plus/>
                           <ci>Kc</ci>
                           <cn cellml:units="millimolar">0.59</cn>
                        </apply>
                     </apply>
                     <cn cellml:units="dimensionless">3</cn>
                  </apply>
                  <apply>
                     <plus/>
                     <ci>V</ci>
                     <cn cellml:units="millivolt">81.9</cn>
                  </apply>
               </apply>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <times/>
                           <cn cellml:units="dimensionless">1.393</cn>
                           <apply>
                              <plus/>
                              <ci>V</ci>
                              <cn cellml:units="millivolt">81.9</cn>
                              <cn cellml:units="millivolt">3.6</cn>
                           </apply>
                        </apply>
                        <ci>RTONF</ci>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="background_current">
      <variable name="i_b" units="nanoA" public_interface="out"/>
      <variable name="g_b" units="microS" initial_value="0.002"/>
      <variable name="V" units="millivolt" public_interface="in"/>
      <variable name="E_b" units="millivolt" initial_value="-40"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_b</ci>
            <apply>
               <times/>
               <ci>g_b</ci>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>E_b</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="sodium_potassium_pump">
      <variable name="i_p" units="nanoA" public_interface="out"/>
      <variable name="I_p" units="nanoA" initial_value="0.1968"/>
      <variable name="V" units="millivolt" public_interface="in"/>
      <variable name="Nai" units="millimolar" public_interface="in"/>
      <variable name="Kc" units="millimolar" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_p</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <ci>I_p</ci>
                  <apply>
                     <power/>
                     <apply>
                        <divide/>
                        <ci>Nai</ci>
                        <apply>
                           <plus/>
                           <cn cellml:units="millimolar">5.64</cn>
                           <ci>Nai</ci>
                        </apply>
                     </apply>
                     <cn cellml:units="dimensionless">3</cn>
                  </apply>
                  <apply>
                     <power/>
                     <apply>
                        <divide/>
                        <ci>Kc</ci>
                        <apply>
                           <plus/>
                           <cn cellml:units="millimolar">0.621</cn>
                           <ci>Kc</ci>
                        </apply>
                     </apply>
                     <cn cellml:units="dimensionless">2</cn>
                  </apply>
                  <cn cellml:units="dimensionless">1.6</cn>
               </apply>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1.5</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <minus/>
                           <apply>
                              <plus/>
                              <ci>V</ci>
                              <cn cellml:units="millivolt">60</cn>
                           </apply>
                        </apply>
                        <cn cellml:units="millivolt">40</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="sodium_calcium_exchange_current">
      <variable name="i_NaCa" units="nanoA" public_interface="out"/>
      <variable name="kNaCa" units="nanoA" initial_value="5.916"/>
      <variable name="x1" units="dimensionless"/>
      <variable name="x2" units="dimensionless"/>
      <variable name="x3" units="dimensionless"/>
      <variable name="x4" units="dimensionless"/>
      <variable name="k41" units="dimensionless"/>
      <variable name="k34" units="dimensionless"/>
      <variable name="k23" units="dimensionless"/>
      <variable name="k21" units="dimensionless"/>
      <variable name="k32" units="dimensionless"/>
      <variable name="k43" units="dimensionless"/>
      <variable name="k12" units="dimensionless"/>
      <variable name="k14" units="dimensionless"/>
      <variable name="Qci" units="dimensionless" initial_value="0.1369"/>
      <variable name="Qn" units="dimensionless" initial_value="0.4315"/>
      <variable name="Qco" units="dimensionless" initial_value="0"/>
      <variable name="Kci" units="millimolar" initial_value="0.0207"/>
      <variable name="K1ni" units="millimolar" initial_value="395.3"/>
      <variable name="K2ni" units="millimolar" initial_value="2.289"/>
      <variable name="K3ni" units="millimolar" initial_value="26.44"/>
      <variable name="Kcni" units="millimolar" initial_value="26.44"/>
      <variable name="K3no" units="millimolar" initial_value="4.663"/>
      <variable name="K1no" units="millimolar" initial_value="1628"/>
      <variable name="K2no" units="millimolar" initial_value="561.4"/>
      <variable name="Kco" units="millimolar" initial_value="3.663"/>
      <variable name="do" units="dimensionless"/>
      <variable name="di" units="dimensionless"/>
      <variable name="Cao" units="millimolar" public_interface="in"/>
      <variable name="Nao" units="millimolar" public_interface="in"/>
      <variable name="Casub" units="millimolar" public_interface="in"/>
      <variable name="Nai" units="millimolar" public_interface="in"/>
      <variable name="RTONF" units="millivolt" public_interface="in"/>
      <variable name="V" units="millivolt" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>do</ci>
            <apply>
               <plus/>
               <cn cellml:units="dimensionless">1</cn>
               <apply>
                  <times/>
                  <apply>
                     <divide/>
                     <ci>Cao</ci>
                     <ci>Kco</ci>
                  </apply>
                  <apply>
                     <plus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <times/>
                              <ci>Qco</ci>
                              <ci>V</ci>
                           </apply>
                           <ci>RTONF</ci>
                        </apply>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <divide/>
                  <ci>Nao</ci>
                  <ci>K1no</ci>
               </apply>
               <apply>
                  <divide/>
                  <apply>
                     <power/>
                     <ci>Nao</ci>
                     <cn cellml:units="dimensionless">2</cn>
                  </apply>
                  <apply>
                     <times/>
                     <ci>K1no</ci>
                     <ci>K2no</ci>
                  </apply>
               </apply>
               <apply>
                  <divide/>
                  <apply>
                     <power/>
                     <ci>Nao</ci>
                     <cn cellml:units="dimensionless">3</cn>
                  </apply>
                  <apply>
                     <times/>
                     <ci>K1no</ci>
                     <ci>K2no</ci>
                     <ci>K3no</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>k32</ci>
            <apply>
               <exp/>
               <apply>
                  <divide/>
                  <apply>
                     <times/>
                     <ci>Qn</ci>
                     <ci>V</ci>
                  </apply>
                  <apply>
                     <times/>
                     <cn cellml:units="dimensionless">2</cn>
                     <ci>RTONF</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>k23</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <apply>
                     <plus/>
                     <apply>
                        <divide/>
                        <apply>
                           <power/>
                           <ci>Nao</ci>
                           <cn cellml:units="dimensionless">2</cn>
                        </apply>
                        <apply>
                           <times/>
                           <ci>K1no</ci>
                           <ci>K2no</ci>
                        </apply>
                     </apply>
                     <apply>
                        <divide/>
                        <apply>
                           <power/>
                           <ci>Nao</ci>
                           <cn cellml:units="dimensionless">3</cn>
                        </apply>
                        <apply>
                           <times/>
                           <ci>K1no</ci>
                           <ci>K2no</ci>
                           <ci>K3no</ci>
                        </apply>
                     </apply>
                  </apply>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <times/>
                           <apply>
                              <minus/>
                              <ci>Qn</ci>
                           </apply>
                           <ci>V</ci>
                        </apply>
                        <apply>
                           <times/>
                           <cn cellml:units="dimensionless">2</cn>
                           <ci>RTONF</ci>
                        </apply>
                     </apply>
                  </apply>
               </apply>
               <ci>do</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>k21</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <apply>
                     <divide/>
                     <ci>Cao</ci>
                     <ci>Kco</ci>
                  </apply>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <times/>
                           <apply>
                              <minus/>
                              <ci>Qco</ci>
                           </apply>
                           <ci>V</ci>
                        </apply>
                        <ci>RTONF</ci>
                     </apply>
                  </apply>
               </apply>
               <ci>do</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>k34</ci>
            <apply>
               <divide/>
               <ci>Nao</ci>
               <apply>
                  <plus/>
                  <ci>K3no</ci>
                  <ci>Nao</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>di</ci>
            <apply>
               <plus/>
               <cn cellml:units="dimensionless">1</cn>
               <apply>
                  <times/>
                  <apply>
                     <divide/>
                     <ci>Casub</ci>
                     <ci>Kci</ci>
                  </apply>
                  <apply>
                     <plus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <times/>
                              <apply>
                                 <minus/>
                                 <ci>Qci</ci>
                              </apply>
                              <ci>V</ci>
                           </apply>
                           <ci>RTONF</ci>
                        </apply>
                     </apply>
                     <apply>
                        <divide/>
                        <ci>Nai</ci>
                        <ci>Kcni</ci>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <divide/>
                  <ci>Nai</ci>
                  <ci>K1ni</ci>
               </apply>
               <apply>
                  <divide/>
                  <apply>
                     <power/>
                     <ci>Nai</ci>
                     <cn cellml:units="dimensionless">2</cn>
                  </apply>
                  <apply>
                     <times/>
                     <ci>K1ni</ci>
                     <ci>K2ni</ci>
                  </apply>
               </apply>
               <apply>
                  <divide/>
                  <apply>
                     <power/>
                     <ci>Nai</ci>
                     <cn cellml:units="dimensionless">3</cn>
                  </apply>
                  <apply>
                     <times/>
                     <ci>K1ni</ci>
                     <ci>K2ni</ci>
                     <ci>K3ni</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>k41</ci>
            <apply>
               <exp/>
               <apply>
                  <divide/>
                  <apply>
                     <times/>
                     <apply>
                        <minus/>
                        <ci>Qn</ci>
                     </apply>
                     <ci>V</ci>
                  </apply>
                  <apply>
                     <times/>
                     <cn cellml:units="dimensionless">2</cn>
                     <ci>RTONF</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>k14</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <apply>
                     <plus/>
                     <apply>
                        <divide/>
                        <apply>
                           <power/>
                           <ci>Nai</ci>
                           <cn cellml:units="dimensionless">2</cn>
                        </apply>
                        <apply>
                           <times/>
                           <ci>K1ni</ci>
                           <ci>K2ni</ci>
                        </apply>
                     </apply>
                     <apply>
                        <divide/>
                        <apply>
                           <power/>
                           <ci>Nai</ci>
                           <cn cellml:units="dimensionless">3</cn>
                        </apply>
                        <apply>
                           <times/>
                           <ci>K1ni</ci>
                           <ci>K2ni</ci>
                           <ci>K3ni</ci>
                        </apply>
                     </apply>
                  </apply>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <times/>
                           <ci>Qn</ci>
                           <ci>V</ci>
                        </apply>
                        <apply>
                           <times/>
                           <cn cellml:units="dimensionless">2</cn>
                           <ci>RTONF</ci>
                        </apply>
                     </apply>
                  </apply>
               </apply>
               <ci>di</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>k12</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <apply>
                     <divide/>
                     <ci>Casub</ci>
                     <ci>Kci</ci>
                  </apply>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <times/>
                           <apply>
                              <minus/>
                              <ci>Qci</ci>
                           </apply>
                           <ci>V</ci>
                        </apply>
                        <ci>RTONF</ci>
                     </apply>
                  </apply>
               </apply>
               <ci>di</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>k43</ci>
            <apply>
               <divide/>
               <ci>Nai</ci>
               <apply>
                  <plus/>
                  <ci>K3ni</ci>
                  <ci>Nai</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>x1</ci>
            <apply>
               <plus/>
               <apply>
                  <times/>
                  <ci>k41</ci>
                  <ci>k34</ci>
                  <apply>
                     <plus/>
                     <ci>k23</ci>
                     <ci>k21</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>k21</ci>
                  <ci>k32</ci>
                  <apply>
                     <plus/>
                     <ci>k43</ci>
                     <ci>k41</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>x2</ci>
            <apply>
               <plus/>
               <apply>
                  <times/>
                  <ci>k32</ci>
                  <ci>k43</ci>
                  <apply>
                     <plus/>
                     <ci>k14</ci>
                     <ci>k12</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>k41</ci>
                  <ci>k12</ci>
                  <apply>
                     <plus/>
                     <ci>k34</ci>
                     <ci>k32</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>x3</ci>
            <apply>
               <plus/>
               <apply>
                  <times/>
                  <ci>k14</ci>
                  <ci>k43</ci>
                  <apply>
                     <plus/>
                     <ci>k23</ci>
                     <ci>k21</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>k12</ci>
                  <ci>k23</ci>
                  <apply>
                     <plus/>
                     <ci>k43</ci>
                     <ci>k41</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>x4</ci>
            <apply>
               <plus/>
               <apply>
                  <times/>
                  <ci>k23</ci>
                  <ci>k34</ci>
                  <apply>
                     <plus/>
                     <ci>k14</ci>
                     <ci>k12</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>k14</ci>
                  <ci>k21</ci>
                  <apply>
                     <plus/>
                     <ci>k34</ci>
                     <ci>k32</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>i_NaCa</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <ci>kNaCa</ci>
                  <apply>
                     <minus/>
                     <apply>
                        <times/>
                        <ci>x2</ci>
                        <ci>k21</ci>
                     </apply>
                     <apply>
                        <times/>
                        <ci>x1</ci>
                        <ci>k12</ci>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <plus/>
                  <ci>x1</ci>
                  <ci>x2</ci>
                  <ci>x3</ci>
                  <ci>x4</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="fast_sodium_current">
      <variable name="i_Na" units="nanoA" public_interface="out"/>
      <variable name="g_Na" units="microlitre_per_second" initial_value="5e-7"/>
      <variable name="E_Na" units="millivolt" public_interface="out"/>
      <variable name="time" units="second" public_interface="in" private_interface="out"/>
      <variable name="V" units="millivolt" public_interface="in" private_interface="out"/>
      <variable name="RTONF" units="millivolt" public_interface="in"/>
      <variable name="F" units="coulomb_per_mole" public_interface="in"/>
      <variable name="Nao" units="millimolar" public_interface="in"/>
      <variable name="Nai" units="millimolar" public_interface="in"/>
      <variable name="Kc" units="millimolar" public_interface="in"/>
      <variable name="Ki" units="millimolar" public_interface="in"/>
      <variable name="m" units="dimensionless" private_interface="in"/>
      <variable name="h1" units="dimensionless" private_interface="in"/>
      <variable name="h2" units="dimensionless" private_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>E_Na</ci>
            <apply>
               <times/>
               <ci>RTONF</ci>
               <apply>
                  <ln/>
                  <apply>
                     <divide/>
                     <ci>Nao</ci>
                     <ci>Nai</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>i_Na</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <ci>g_Na</ci>
                        <apply>
                           <power/>
                           <ci>m</ci>
                           <cn cellml:units="dimensionless">3</cn>
                        </apply>
                        <apply>
                           <plus/>
                           <apply>
                              <times/>
                              <cn cellml:units="dimensionless">0.635</cn>
                              <ci>h1</ci>
                           </apply>
                           <apply>
                              <times/>
                              <cn cellml:units="dimensionless">0.365</cn>
                              <ci>h2</ci>
                           </apply>
                        </apply>
                        <ci>Nao</ci>
                        <ci>V</ci>
                        <ci>F</ci>
                     </apply>
                     <ci>RTONF</ci>
                  </apply>
                  <apply>
                     <minus/>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <minus/>
                              <ci>V</ci>
                              <ci>E_Na</ci>
                           </apply>
                           <ci>RTONF</ci>
                        </apply>
                     </apply>
                     <cn cellml:units="dimensionless">1</cn>
                  </apply>
               </apply>
               <apply>
                  <minus/>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <ci>V</ci>
                        <ci>RTONF</ci>
                     </apply>
                  </apply>
                  <cn cellml:units="dimensionless">1</cn>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="fast_sodium_current_m_gate">
      <variable name="m" units="dimensionless" initial_value="0.0132200747771872" public_interface="out"/>
      <variable name="alpha_m" units="per_second"/>
      <variable name="beta_m" units="per_second"/>
      <variable name="V" units="millivolt" public_interface="in"/>
      <variable name="time" units="second" public_interface="in"/>
      <variable name="delta_m" units="millivolt" initial_value="1e-5"/>
      <variable name="E0_m" units="millivolt"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>E0_m</ci>
            <apply>
               <plus/>
               <ci>V</ci>
               <cn cellml:units="millivolt">44.4</cn>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>alpha_m</ci>
            <piecewise>
               <piece>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <apply>
                           <minus/>
                           <cn cellml:units="per_second">460</cn>
                        </apply>
                        <apply>
                           <minus/>
                           <cn cellml:units="dimensionless">12.673</cn>
                        </apply>
                     </apply>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <ci>E0_m</ci>
                           <apply>
                              <minus/>
                              <cn cellml:units="millivolt">12.673</cn>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
                  <apply>
                     <lt/>
                     <apply>
                        <abs/>
                        <ci>E0_m</ci>
                     </apply>
                     <ci>delta_m</ci>
                  </apply>
               </piece>
               <otherwise>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <apply>
                           <minus/>
                           <cn cellml:units="per_millivolt_second">460</cn>
                        </apply>
                        <ci>E0_m</ci>
                     </apply>
                     <apply>
                        <minus/>
                        <apply>
                           <exp/>
                           <apply>
                              <divide/>
                              <ci>E0_m</ci>
                              <apply>
                                 <minus/>
                                 <cn cellml:units="millivolt">12.673</cn>
                              </apply>
                           </apply>
                        </apply>
                        <cn cellml:units="dimensionless">1</cn>
                     </apply>
                  </apply>
               </otherwise>
            </piecewise>
         </apply>
         <apply>
            <eq/>
            <ci>beta_m</ci>
            <apply>
               <times/>
               <cn cellml:units="per_second">18400</cn>
               <apply>
                  <exp/>
                  <apply>
                     <divide/>
                     <ci>E0_m</ci>
                     <apply>
                        <minus/>
                        <cn cellml:units="millivolt">12.673</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>m</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>alpha_m</ci>
                  <apply>
                     <minus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <ci>m</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>beta_m</ci>
                  <ci>m</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="fast_sodium_current_h1_gate">
      <variable name="h1" units="dimensionless" initial_value="0.706622937059237" public_interface="out"/>
      <variable name="alpha_h1" units="per_second"/>
      <variable name="beta_h1" units="per_second"/>
      <variable name="h1_inf" units="dimensionless"/>
      <variable name="tau_h1" units="second"/>
      <variable name="V" units="millivolt" public_interface="in"/>
      <variable name="time" units="second" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>alpha_h1</ci>
            <apply>
               <times/>
               <cn cellml:units="per_second">44.9</cn>
               <apply>
                  <exp/>
                  <apply>
                     <divide/>
                     <apply>
                        <plus/>
                        <ci>V</ci>
                        <cn cellml:units="millivolt">66.9</cn>
                     </apply>
                     <apply>
                        <minus/>
                        <cn cellml:units="millivolt">5.57</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>beta_h1</ci>
            <apply>
               <divide/>
               <cn cellml:units="per_second">1491</cn>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <times/>
                     <cn cellml:units="dimensionless">323.3</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <plus/>
                              <ci>V</ci>
                              <cn cellml:units="millivolt">94.6</cn>
                           </apply>
                           <apply>
                              <minus/>
                              <cn cellml:units="millivolt">12.9</cn>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>h1_inf</ci>
            <apply>
               <divide/>
               <ci>alpha_h1</ci>
               <apply>
                  <plus/>
                  <ci>alpha_h1</ci>
                  <ci>beta_h1</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>tau_h1</ci>
            <apply>
               <plus/>
               <apply>
                  <divide/>
                  <cn cellml:units="second">0.03</cn>
                  <apply>
                     <plus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <plus/>
                              <ci>V</ci>
                              <cn cellml:units="millivolt">40</cn>
                           </apply>
                           <cn cellml:units="millivolt">6</cn>
                        </apply>
                     </apply>
                  </apply>
               </apply>
               <cn cellml:units="second">0.00035</cn>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>h1</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>h1_inf</ci>
                  <ci>h1</ci>
               </apply>
               <ci>tau_h1</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="fast_sodium_current_h2_gate">
      <variable name="h2" units="dimensionless" initial_value="0.701626826712569" public_interface="out"/>
      <variable name="alpha_h2" units="per_second"/>
      <variable name="beta_h2" units="per_second"/>
      <variable name="h2_inf" units="dimensionless"/>
      <variable name="tau_h2" units="second"/>
      <variable name="V" units="millivolt" public_interface="in"/>
      <variable name="time" units="second" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>alpha_h2</ci>
            <apply>
               <times/>
               <cn cellml:units="per_second">44.9</cn>
               <apply>
                  <exp/>
                  <apply>
                     <divide/>
                     <apply>
                        <plus/>
                        <ci>V</ci>
                        <cn cellml:units="millivolt">66.9</cn>
                     </apply>
                     <apply>
                        <minus/>
                        <cn cellml:units="millivolt">5.57</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>beta_h2</ci>
            <apply>
               <divide/>
               <cn cellml:units="per_second">1491</cn>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <times/>
                     <cn cellml:units="dimensionless">323.3</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <plus/>
                              <ci>V</ci>
                              <cn cellml:units="millivolt">94.6</cn>
                           </apply>
                           <apply>
                              <minus/>
                              <cn cellml:units="millivolt">12.9</cn>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>h2_inf</ci>
            <apply>
               <divide/>
               <ci>alpha_h2</ci>
               <apply>
                  <plus/>
                  <ci>alpha_h2</ci>
                  <ci>beta_h2</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>tau_h2</ci>
            <apply>
               <plus/>
               <apply>
                  <divide/>
                  <cn cellml:units="second">0.12</cn>
                  <apply>
                     <plus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <plus/>
                              <ci>V</ci>
                              <cn cellml:units="millivolt">60</cn>
                           </apply>
                           <cn cellml:units="millivolt">2</cn>
                        </apply>
                     </apply>
                  </apply>
               </apply>
               <cn cellml:units="second">0.00295</cn>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>h2</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>h2_inf</ci>
                  <ci>h2</ci>
               </apply>
               <ci>tau_h2</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="L_type_calcium_current">
      <variable name="i_CaL" units="nanoA" public_interface="out"/>
      <variable name="g_CaL" units="microS" initial_value="0.021"/>
      <variable name="E_CaL" units="millivolt" initial_value="62.1"/>
      <variable name="i_ACh" units="nanoA" public_interface="in"/>
      <variable name="ACh" units="millimolar" public_interface="in"/>
      <variable name="time" units="second" public_interface="in" private_interface="out"/>
      <variable name="V" units="millivolt" public_interface="in" private_interface="out"/>
      <variable name="d" units="dimensionless" private_interface="in"/>
      <variable name="f" units="dimensionless" private_interface="in"/>
      <variable name="f2" units="dimensionless" private_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_CaL</ci>
            <apply>
               <times/>
               <ci>g_CaL</ci>
               <ci>d</ci>
               <apply>
                  <plus/>
                  <apply>
                     <times/>
                     <cn cellml:units="dimensionless">0.675</cn>
                     <ci>f</ci>
                  </apply>
                  <apply>
                     <times/>
                     <cn cellml:units="dimensionless">0.325</cn>
                     <ci>f2</ci>
                  </apply>
               </apply>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>E_CaL</ci>
               </apply>
               <apply>
                  <minus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <divide/>
                     <apply>
                        <divide/>
                        <apply>
                           <times/>
                           <ci>i_ACh</ci>
                           <ci>ACh</ci>
                        </apply>
                        <apply>
                           <plus/>
                           <cn type="e-notation" cellml:units="millimolar">0.9   <sep/>
                           -4</cn>
                           <ci>ACh</ci>
                        </apply>
                     </apply>
                     <cn cellml:units="nanoA">1</cn>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="L_type_calcium_current_d_gate">
      <variable name="d" units="dimensionless" initial_value="4.23500474189711e-5" public_interface="out"/>
      <variable name="alpha_d" units="per_second"/>
      <variable name="beta_d" units="per_second"/>
      <variable name="d_inf" units="dimensionless"/>
      <variable name="tau_d" units="second"/>
      <variable name="V" units="millivolt" public_interface="in"/>
      <variable name="time" units="second" public_interface="in"/>
      <variable name="act_shift" units="millivolt" initial_value="0"/>
      <variable name="slope_factor_act" units="millivolt" initial_value="-6.61"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>alpha_d</ci>
            <apply>
               <plus/>
               <apply>
                  <divide/>
                  <apply>
                     <times/>
                     <apply>
                        <minus/>
                        <cn cellml:units="per_millivolt_second">26.12</cn>
                     </apply>
                     <apply>
                        <plus/>
                        <ci>V</ci>
                        <cn cellml:units="millivolt">35</cn>
                     </apply>
                  </apply>
                  <apply>
                     <minus/>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <plus/>
                              <ci>V</ci>
                              <cn cellml:units="millivolt">35</cn>
                           </apply>
                           <apply>
                              <minus/>
                              <cn cellml:units="millivolt">2.5</cn>
                           </apply>
                        </apply>
                     </apply>
                     <cn cellml:units="dimensionless">1</cn>
                  </apply>
               </apply>
               <apply>
                  <divide/>
                  <apply>
                     <times/>
                     <apply>
                        <minus/>
                        <cn cellml:units="per_millivolt_second">78.11</cn>
                     </apply>
                     <ci>V</ci>
                  </apply>
                  <apply>
                     <minus/>
                     <apply>
                        <exp/>
                        <apply>
                           <times/>
                           <apply>
                              <minus/>
                              <cn cellml:units="per_millivolt">0.208</cn>
                           </apply>
                           <ci>V</ci>
                        </apply>
                     </apply>
                     <cn cellml:units="dimensionless">1</cn>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>beta_d</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <cn cellml:units="per_millivolt_second">10.52</cn>
                  <apply>
                     <minus/>
                     <ci>V</ci>
                     <cn cellml:units="millivolt">5</cn>
                  </apply>
               </apply>
               <apply>
                  <minus/>
                  <apply>
                     <exp/>
                     <apply>
                        <times/>
                        <cn cellml:units="per_millivolt">0.4</cn>
                        <apply>
                           <minus/>
                           <ci>V</ci>
                           <cn cellml:units="millivolt">5</cn>
                        </apply>
                     </apply>
                  </apply>
                  <cn cellml:units="dimensionless">1</cn>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>d_inf</ci>
            <apply>
               <divide/>
               <cn cellml:units="dimensionless">1</cn>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <minus/>
                           <ci>V</ci>
                           <apply>
                              <plus/>
                              <apply>
                                 <minus/>
                                 <cn cellml:units="millivolt">3.2</cn>
                              </apply>
                              <ci>act_shift</ci>
                           </apply>
                        </apply>
                        <ci>slope_factor_act</ci>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>tau_d</ci>
            <apply>
               <divide/>
               <cn cellml:units="dimensionless">1</cn>
               <apply>
                  <plus/>
                  <ci>alpha_d</ci>
                  <ci>beta_d</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>d</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>d_inf</ci>
                  <ci>d</ci>
               </apply>
               <ci>tau_d</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="L_type_calcium_current_f_gate">
      <variable name="f" units="dimensionless" initial_value="0.998435073735753" public_interface="out"/>
      <variable name="f_inf" units="dimensionless"/>
      <variable name="tau_f" units="second"/>
      <variable name="V" units="millivolt" public_interface="in"/>
      <variable name="time" units="second" public_interface="in"/>
      <variable name="inact_shift" units="millivolt" initial_value="-5"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>f_inf</ci>
            <apply>
               <divide/>
               <cn cellml:units="dimensionless">1</cn>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <minus/>
                           <ci>V</ci>
                           <apply>
                              <plus/>
                              <apply>
                                 <minus/>
                                 <cn cellml:units="millivolt">24</cn>
                              </apply>
                              <ci>inact_shift</ci>
                           </apply>
                        </apply>
                        <cn cellml:units="millivolt">6.31</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>tau_f</ci>
            <apply>
               <plus/>
               <cn cellml:units="second">0.01</cn>
               <apply>
                  <times/>
                  <cn cellml:units="second">0.1539</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">185.67</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>f</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>f_inf</ci>
                  <ci>f</ci>
               </apply>
               <ci>tau_f</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="L_type_calcium_current_f2_gate">
      <variable name="f2" units="dimensionless" initial_value="0.998424216938754" public_interface="out"/>
      <variable name="f2_inf" units="dimensionless"/>
      <variable name="tau_f2" units="second"/>
      <variable name="V" units="millivolt" public_interface="in"/>
      <variable name="time" units="second" public_interface="in"/>
      <variable name="inact_shift" units="millivolt" initial_value="-5"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>f2_inf</ci>
            <apply>
               <divide/>
               <cn cellml:units="dimensionless">1</cn>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <minus/>
                           <ci>V</ci>
                           <apply>
                              <plus/>
                              <apply>
                                 <minus/>
                                 <cn cellml:units="millivolt">24</cn>
                              </apply>
                              <ci>inact_shift</ci>
                           </apply>
                        </apply>
                        <cn cellml:units="millivolt">6.31</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>tau_f2</ci>
            <apply>
               <plus/>
               <cn cellml:units="second">0.06</cn>
               <apply>
                  <times/>
                  <cn cellml:units="second">0.48076</cn>
                  <cn cellml:units="dimensionless">2.25</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <minus/>
                           <apply>
                              <power/>
                              <apply>
                                 <minus/>
                                 <ci>V</ci>
                                 <apply>
                                    <minus/>
                                    <cn cellml:units="millivolt">40</cn>
                                 </apply>
                              </apply>
                              <cn cellml:units="dimensionless">2</cn>
                           </apply>
                        </apply>
                        <cn cellml:units="millivolt2">138.04</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>f2</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>f2_inf</ci>
                  <ci>f2</ci>
               </apply>
               <ci>tau_f2</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="transient_outward_potassium_current">
      <variable name="i_to" units="nanoA" public_interface="out"/>
      <variable name="E_K" units="millivolt" public_interface="out"/>
      <variable name="g_to" units="microS" initial_value="0.014"/>
      <variable name="time" units="second" public_interface="in" private_interface="out"/>
      <variable name="V" units="millivolt" public_interface="in" private_interface="out"/>
      <variable name="RTONF" units="millivolt" public_interface="in"/>
      <variable name="Kc" units="millimolar" public_interface="in"/>
      <variable name="Ki" units="millimolar" public_interface="in"/>
      <variable name="r" units="dimensionless" private_interface="in"/>
      <variable name="q_fast" units="dimensionless" private_interface="in"/>
      <variable name="q_slow" units="dimensionless" private_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_to</ci>
            <apply>
               <times/>
               <ci>g_to</ci>
               <ci>r</ci>
               <apply>
                  <plus/>
                  <apply>
                     <times/>
                     <cn cellml:units="dimensionless">0.45</cn>
                     <ci>q_fast</ci>
                  </apply>
                  <apply>
                     <times/>
                     <cn cellml:units="dimensionless">0.55</cn>
                     <ci>q_slow</ci>
                  </apply>
               </apply>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>E_K</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>E_K</ci>
            <apply>
               <times/>
               <ci>RTONF</ci>
               <apply>
                  <ln/>
                  <apply>
                     <divide/>
                     <ci>Kc</ci>
                     <ci>Ki</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="transient_outward_potassium_current_r_gate">
      <variable name="r" units="dimensionless" initial_value="0.00894826428663828" public_interface="out"/>
      <variable name="tau_r" units="second"/>
      <variable name="r_infinity" units="dimensionless"/>
      <variable name="V" units="millivolt" public_interface="in"/>
      <variable name="time" units="second" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>r</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>r_infinity</ci>
                  <ci>r</ci>
               </apply>
               <ci>tau_r</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>r_infinity</ci>
            <apply>
               <divide/>
               <cn cellml:units="dimensionless">1</cn>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <minus/>
                           <ci>V</ci>
                           <cn cellml:units="millivolt">7.44</cn>
                        </apply>
                        <apply>
                           <minus/>
                           <cn cellml:units="millivolt">16.4</cn>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>tau_r</ci>
            <apply>
               <plus/>
               <cn type="e-notation" cellml:units="second">0.596   <sep/>
               -3</cn>
               <apply>
                  <divide/>
                  <cn type="e-notation" cellml:units="second">3.118   <sep/>
                  -3</cn>
                  <apply>
                     <plus/>
                     <apply>
                        <times/>
                        <cn cellml:units="dimensionless">1.037</cn>
                        <apply>
                           <exp/>
                           <apply>
                              <times/>
                              <cn cellml:units="per_millivolt">0.09</cn>
                              <apply>
                                 <plus/>
                                 <ci>V</ci>
                                 <cn cellml:units="millivolt">30.61</cn>
                              </apply>
                           </apply>
                        </apply>
                     </apply>
                     <apply>
                        <times/>
                        <cn cellml:units="dimensionless">0.396</cn>
                        <apply>
                           <exp/>
                           <apply>
                              <times/>
                              <apply>
                                 <minus/>
                                 <cn cellml:units="per_millivolt">0.12</cn>
                              </apply>
                              <apply>
                                 <plus/>
                                 <ci>V</ci>
                                 <cn cellml:units="millivolt">23.84</cn>
                              </apply>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="transient_outward_potassium_current_qfast_gate">
      <variable name="q_fast" units="dimensionless" initial_value="0.994837524153424" public_interface="out"/>
      <variable name="tau_qfast" units="second"/>
      <variable name="qfast_infinity" units="dimensionless"/>
      <variable name="V" units="millivolt" public_interface="in"/>
      <variable name="time" units="second" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>q_fast</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>qfast_infinity</ci>
                  <ci>q_fast</ci>
               </apply>
               <ci>tau_qfast</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>qfast_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">33.8</cn>
                        </apply>
                        <cn cellml:units="millivolt">6.12</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>tau_qfast</ci>
            <apply>
               <plus/>
               <cn type="e-notation" cellml:units="second">12.66   <sep/>
               -3</cn>
               <apply>
                  <divide/>
                  <cn cellml:units="second">4.72716</cn>
                  <apply>
                     <plus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <plus/>
                              <ci>V</ci>
                              <cn cellml:units="millivolt">154.5</cn>
                           </apply>
                           <cn cellml:units="millivolt">23.96</cn>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="transient_outward_potassium_current_qslow_gate">
      <variable name="q_slow" units="dimensionless" initial_value="0.427382372349565" public_interface="out"/>
      <variable name="tau_qslow" units="second"/>
      <variable name="qslow_infinity" units="dimensionless"/>
      <variable name="V" units="millivolt" public_interface="in"/>
      <variable name="time" units="second" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>q_slow</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>qslow_infinity</ci>
                  <ci>q_slow</ci>
               </apply>
               <ci>tau_qslow</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>qslow_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">33.8</cn>
                        </apply>
                        <cn cellml:units="millivolt">6.12</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>tau_qslow</ci>
            <apply>
               <plus/>
               <cn cellml:units="second">0.1</cn>
               <apply>
                  <times/>
                  <cn cellml:units="second">4</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <minus/>
                           <apply>
                              <power/>
                              <apply>
                                 <plus/>
                                 <ci>V</ci>
                                 <cn cellml:units="millivolt">65</cn>
                              </apply>
                              <cn cellml:units="dimensionless">2</cn>
                           </apply>
                        </apply>
                        <cn cellml:units="millivolt2">500</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="sustained_outward_potassium_current">
      <variable name="i_st" units="nanoA" public_interface="out"/>
      <variable name="E_st" units="millivolt" initial_value="-37.4"/>
      <variable name="g_st" units="microS" initial_value="0"/>
      <variable name="time" units="second" public_interface="in" private_interface="out"/>
      <variable name="V" units="millivolt" public_interface="in" private_interface="out"/>
      <variable name="qa" units="dimensionless" private_interface="in"/>
      <variable name="qi" units="dimensionless" private_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_st</ci>
            <apply>
               <times/>
               <ci>g_st</ci>
               <ci>qa</ci>
               <ci>qi</ci>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>E_st</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="sustained_outward_potassium_current_qa_gate">
      <variable name="qa" units="dimensionless" initial_value="0.0910882041816457" public_interface="out"/>
      <variable name="tau_qa" units="second"/>
      <variable name="qa_infinity" units="dimensionless"/>
      <variable name="alpha_qa" units="per_second"/>
      <variable name="beta_qa" units="per_second"/>
      <variable name="V" units="millivolt" public_interface="in"/>
      <variable name="time" units="second" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>qa</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>qa_infinity</ci>
                  <ci>qa</ci>
               </apply>
               <ci>tau_qa</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>qa_infinity</ci>
            <apply>
               <divide/>
               <cn cellml:units="dimensionless">1</cn>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <minus/>
                           <ci>V</ci>
                           <apply>
                              <minus/>
                              <cn cellml:units="millivolt">49.1</cn>
                           </apply>
                        </apply>
                        <apply>
                           <minus/>
                           <cn cellml:units="millivolt">8.98</cn>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>alpha_qa</ci>
            <apply>
               <divide/>
               <cn cellml:units="dimensionless">1</cn>
               <apply>
                  <plus/>
                  <apply>
                     <times/>
                     <cn cellml:units="second">0.15</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <minus/>
                              <ci>V</ci>
                           </apply>
                           <cn cellml:units="millivolt">11</cn>
                        </apply>
                     </apply>
                  </apply>
                  <apply>
                     <times/>
                     <cn cellml:units="second">0.2</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <minus/>
                              <ci>V</ci>
                           </apply>
                           <cn cellml:units="millivolt">700</cn>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>beta_qa</ci>
            <apply>
               <divide/>
               <cn cellml:units="dimensionless">1</cn>
               <apply>
                  <plus/>
                  <apply>
                     <times/>
                     <cn cellml:units="second">16</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <ci>V</ci>
                           <cn cellml:units="millivolt">8</cn>
                        </apply>
                     </apply>
                  </apply>
                  <apply>
                     <times/>
                     <cn cellml:units="second">15</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <ci>V</ci>
                           <cn cellml:units="millivolt">50</cn>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>tau_qa</ci>
            <apply>
               <divide/>
               <cn cellml:units="dimensionless">0.001</cn>
               <apply>
                  <plus/>
                  <ci>alpha_qa</ci>
                  <ci>beta_qa</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="sustained_outward_potassium_current_qi_gate">
      <variable name="qi" units="dimensionless" initial_value="0.890389014175329" public_interface="out"/>
      <variable name="tau_qi" units="second"/>
      <variable name="alpha_qi" units="per_second"/>
      <variable name="beta_qi" units="per_second"/>
      <variable name="qi_infinity" units="dimensionless"/>
      <variable name="V" units="millivolt" public_interface="in"/>
      <variable name="time" units="second" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>qi</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>qi_infinity</ci>
                  <ci>qi</ci>
               </apply>
               <ci>tau_qi</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>alpha_qi</ci>
            <apply>
               <divide/>
               <cn cellml:units="dimensionless">0.1504</cn>
               <apply>
                  <plus/>
                  <apply>
                     <times/>
                     <cn cellml:units="second">3100</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <ci>V</ci>
                           <cn cellml:units="millivolt">13</cn>
                        </apply>
                     </apply>
                  </apply>
                  <apply>
                     <times/>
                     <cn cellml:units="second">700</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <ci>V</ci>
                           <cn cellml:units="millivolt">70</cn>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>beta_qi</ci>
            <apply>
               <plus/>
               <apply>
                  <divide/>
                  <cn cellml:units="dimensionless">0.1504</cn>
                  <apply>
                     <plus/>
                     <apply>
                        <times/>
                        <cn cellml:units="second">95</cn>
                        <apply>
                           <exp/>
                           <apply>
                              <divide/>
                              <apply>
                                 <minus/>
                                 <ci>V</ci>
                              </apply>
                              <cn cellml:units="millivolt">10</cn>
                           </apply>
                        </apply>
                     </apply>
                     <apply>
                        <times/>
                        <cn cellml:units="second">50</cn>
                        <apply>
                           <exp/>
                           <apply>
                              <divide/>
                              <apply>
                                 <minus/>
                                 <ci>V</ci>
                              </apply>
                              <cn cellml:units="millivolt">700</cn>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <divide/>
                  <cn cellml:units="per_second">0.000229</cn>
                  <apply>
                     <plus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <minus/>
                              <ci>V</ci>
                           </apply>
                           <cn cellml:units="millivolt">5</cn>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>qi_infinity</ci>
            <apply>
               <divide/>
               <ci>alpha_qi</ci>
               <apply>
                  <plus/>
                  <ci>alpha_qi</ci>
                  <ci>beta_qi</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>tau_qi</ci>
            <apply>
               <divide/>
               <cn cellml:units="dimensionless">0.001</cn>
               <apply>
                  <plus/>
                  <ci>alpha_qi</ci>
                  <ci>beta_qi</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="acetylcholine_sensitive_current">
      <variable name="i_ACh" units="nanoA" public_interface="out"/>
      <variable name="g_ACh" units="microS"/>
      <variable name="g_ACh_max" units="microS" initial_value="0.0198"/>
      <variable name="time" units="second" public_interface="in" private_interface="out"/>
      <variable name="V" units="millivolt" public_interface="in" private_interface="out"/>
      <variable name="E_K" units="millivolt" public_interface="in"/>
      <variable name="RTONF" units="millivolt" public_interface="in"/>
      <variable name="ACh" units="millimolar" initial_value="0" public_interface="out"/>
      <variable name="K_ACh" units="millimolar" initial_value="0.00035"/>
      <variable name="Kc" units="millimolar" public_interface="in"/>
      <variable name="achf" units="dimensionless" private_interface="in"/>
      <variable name="achs" units="dimensionless" private_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>g_ACh</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <ci>g_ACh_max</ci>
                  <ci>achf</ci>
                  <ci>achs</ci>
                  <apply>
                     <power/>
                     <ci>ACh</ci>
                     <cn cellml:units="dimensionless">1.5</cn>
                  </apply>
               </apply>
               <apply>
                  <plus/>
                  <apply>
                     <power/>
                     <ci>K_ACh</ci>
                     <cn cellml:units="dimensionless">1.5</cn>
                  </apply>
                  <apply>
                     <power/>
                     <ci>ACh</ci>
                     <cn cellml:units="dimensionless">1.5</cn>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>i_ACh</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <ci>g_ACh</ci>
                        <ci>Kc</ci>
                     </apply>
                     <apply>
                        <plus/>
                        <cn cellml:units="millimolar">10</cn>
                        <ci>Kc</ci>
                     </apply>
                  </apply>
                  <apply>
                     <minus/>
                     <ci>V</ci>
                     <ci>E_K</ci>
                  </apply>
               </apply>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <minus/>
                           <apply>
                              <minus/>
                              <ci>V</ci>
                              <ci>E_K</ci>
                           </apply>
                           <cn cellml:units="millivolt">140</cn>
                        </apply>
                        <apply>
                           <times/>
                           <cn cellml:units="dimensionless">2.5</cn>
                           <ci>RTONF</ci>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="acetylcholine_sensitive_current_achf_gate">
      <variable name="achf" units="dimensionless" initial_value="0.74242774587522" public_interface="out"/>
      <variable name="alpha_achf" units="per_second" initial_value="73.1"/>
      <variable name="beta_achf" units="per_second"/>
      <variable name="time" units="second" public_interface="in"/>
      <variable name="V" units="millivolt" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>beta_achf</ci>
            <apply>
               <divide/>
               <cn cellml:units="per_second">120</cn>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <minus/>
                           <apply>
                              <plus/>
                              <ci>V</ci>
                              <cn cellml:units="millivolt">50</cn>
                           </apply>
                        </apply>
                        <cn cellml:units="millivolt">15</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>achf</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>alpha_achf</ci>
                  <apply>
                     <minus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <ci>achf</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>beta_achf</ci>
                  <ci>achf</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="acetylcholine_sensitive_current_achs_gate">
      <variable name="achs" units="dimensionless" initial_value="0.746918323597392" public_interface="out"/>
      <variable name="alpha_achs" units="per_second" initial_value="3.7"/>
      <variable name="beta_achs" units="per_second"/>
      <variable name="time" units="second" public_interface="in"/>
      <variable name="V" units="millivolt" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>beta_achs</ci>
            <apply>
               <divide/>
               <cn cellml:units="per_second">5.82</cn>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <minus/>
                           <apply>
                              <plus/>
                              <ci>V</ci>
                              <cn cellml:units="millivolt">50</cn>
                           </apply>
                        </apply>
                        <cn cellml:units="millivolt">15</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>achs</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>alpha_achs</ci>
                  <apply>
                     <minus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <ci>achs</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>beta_achs</ci>
                  <ci>achs</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="extracellular_sodium_concentration">
      <variable name="Nao" units="millimolar" initial_value="140" public_interface="out"/>
   </component>
   <component name="intracellular_sodium_concentration">
      <variable name="Nai" units="millimolar" initial_value="8" public_interface="out"/>
   </component>
   <component name="extracellular_calcium_concentration">
      <variable name="Cao" units="millimolar" initial_value="2" public_interface="out"/>
   </component>
   <component name="intracellular_calcium_concentration">
      <variable name="Cai" units="millimolar" initial_value="3.73317163732586e-5" public_interface="out"/>
      <variable name="Casub" units="millimolar" initial_value="3.27335718697622e-5" public_interface="out"/>
      <variable name="V_up" units="micrometre3"/>
      <variable name="V_rel" units="micrometre3"/>
      <variable name="V_sub" units="micrometre3"/>
      <variable name="Vi" units="micrometre3"/>
      <variable name="V_cell" units="micrometre3" initial_value="4.39823e-6"/>
      <variable name="i_up" units="millimolar_per_second"/>
      <variable name="i_tr" units="millimolar_per_second"/>
      <variable name="i_rel" units="millimolar_per_second"/>
      <variable name="i_diff" units="millimolar_per_second"/>
      <variable name="Ca_up" units="millimolar" initial_value="0.818671555213184"/>
      <variable name="Ca_rel" units="millimolar" initial_value="0.682159360306652"/>
      <variable name="P_rel" units="per_second" initial_value="1805.6"/>
      <variable name="K_up" units="millimolar" initial_value="0.0006"/>
      <variable name="tau_tr" units="second" initial_value="0.06"/>
      <variable name="time" units="second" public_interface="in"/>
      <variable name="V" units="millivolt" public_interface="in"/>
      <variable name="F" units="coulomb_per_mole" public_interface="in"/>
      <variable name="i_CaL" units="nanoA" public_interface="in"/>
      <variable name="i_NaCa" units="nanoA" public_interface="in"/>
      <variable name="f_TC" units="dimensionless" initial_value="0.00739583869345967"/>
      <variable name="f_TMC" units="dimensionless" initial_value="0.133773505989393"/>
      <variable name="f_TMM" units="dimensionless" initial_value="0.765246381247448"/>
      <variable name="f_CMi" units="dimensionless" initial_value="0.0154714370092264"/>
      <variable name="f_CMs" units="dimensionless" initial_value="0.0135767851016544"/>
      <variable name="f_CQ" units="dimensionless" initial_value="0.449992033196647"/>
      <variable name="f_CSL" units="dimensionless" initial_value="1.21722016147587e-5"/>
      <variable name="diff_f_TC" units="per_second"/>
      <variable name="diff_f_TMC" units="per_second"/>
      <variable name="diff_f_TMM" units="per_second"/>
      <variable name="diff_f_CMi" units="per_second"/>
      <variable name="diff_f_CMs" units="per_second"/>
      <variable name="diff_f_CQ" units="per_second"/>
      <variable name="diff_f_CSL" units="per_second"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_diff</ci>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>Casub</ci>
                  <ci>Cai</ci>
               </apply>
               <cn type="e-notation" cellml:units="second">0.04   <sep/>
               -3</cn>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>i_up</ci>
            <apply>
               <divide/>
               <cn cellml:units="millimolar_per_second">5</cn>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <divide/>
                     <ci>K_up</ci>
                     <ci>Cai</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>i_tr</ci>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>Ca_up</ci>
                  <ci>Ca_rel</ci>
               </apply>
               <ci>tau_tr</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>i_rel</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <ci>P_rel</ci>
                  <apply>
                     <minus/>
                     <ci>Ca_rel</ci>
                     <ci>Casub</ci>
                  </apply>
               </apply>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <power/>
                     <apply>
                        <divide/>
                        <cn cellml:units="millimolar">0.0012</cn>
                        <ci>Casub</ci>
                     </apply>
                     <cn cellml:units="dimensionless">2</cn>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>Ca_up</ci>
            </apply>
            <apply>
               <minus/>
               <ci>i_up</ci>
               <apply>
                  <divide/>
                  <apply>
                     <times/>
                     <ci>i_tr</ci>
                     <ci>V_rel</ci>
                  </apply>
                  <ci>V_up</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>V_up</ci>
            <apply>
               <times/>
               <cn cellml:units="dimensionless">0.0116</cn>
               <ci>V_cell</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>Ca_rel</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <minus/>
                  <ci>i_tr</ci>
                  <ci>i_rel</ci>
               </apply>
               <apply>
                  <times/>
                  <cn cellml:units="millimolar">10</cn>
                  <ci>diff_f_CQ</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>V_rel</ci>
            <apply>
               <times/>
               <cn cellml:units="dimensionless">0.0012</cn>
               <ci>V_cell</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>Cai</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <divide/>
                  <apply>
                     <minus/>
                     <apply>
                        <times/>
                        <ci>i_diff</ci>
                        <ci>V_sub</ci>
                     </apply>
                     <apply>
                        <times/>
                        <ci>i_up</ci>
                        <ci>V_up</ci>
                     </apply>
                  </apply>
                  <ci>Vi</ci>
               </apply>
               <apply>
                  <plus/>
                  <apply>
                     <times/>
                     <cn cellml:units="millimolar">0.045</cn>
                     <ci>diff_f_CMi</ci>
                  </apply>
                  <apply>
                     <times/>
                     <cn cellml:units="millimolar">0.031</cn>
                     <ci>diff_f_TC</ci>
                  </apply>
                  <apply>
                     <times/>
                     <cn cellml:units="millimolar">0.062</cn>
                     <ci>diff_f_TMC</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>Casub</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <minus/>
                  <apply>
                     <minus/>
                     <apply>
                        <divide/>
                        <apply>
                           <plus/>
                           <apply>
                              <divide/>
                              <apply>
                                 <minus/>
                                 <apply>
                                    <minus/>
                                    <ci>i_CaL</ci>
                                    <apply>
                                       <times/>
                                       <cn cellml:units="dimensionless">2</cn>
                                       <ci>i_NaCa</ci>
                                    </apply>
                                 </apply>
                              </apply>
                              <apply>
                                 <times/>
                                 <cn cellml:units="microlitre_micrometre3">2</cn>
                                 <ci>F</ci>
                              </apply>
                           </apply>
                           <apply>
                              <times/>
                              <ci>i_rel</ci>
                              <ci>V_rel</ci>
                           </apply>
                        </apply>
                        <ci>V_sub</ci>
                     </apply>
                     <ci>i_diff</ci>
                  </apply>
                  <apply>
                     <times/>
                     <cn cellml:units="millimolar">0.045</cn>
                     <ci>diff_f_CMs</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <apply>
                     <divide/>
                     <cn cellml:units="millimolar">0.031</cn>
                     <cn cellml:units="dimensionless">1.2</cn>
                  </apply>
                  <ci>diff_f_CSL</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>Vi</ci>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <cn cellml:units="dimensionless">0.46</cn>
                  <ci>V_cell</ci>
               </apply>
               <ci>V_sub</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>V_sub</ci>
            <apply>
               <times/>
               <cn cellml:units="dimensionless">0.01</cn>
               <ci>V_cell</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>diff_f_TC</ci>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <cn type="e-notation" cellml:units="per_millimolar_per_second">88.8   <sep/>
                  3</cn>
                  <ci>Cai</ci>
                  <apply>
                     <minus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <ci>f_TC</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <cn type="e-notation" cellml:units="per_second">0.446   <sep/>
                  3</cn>
                  <ci>f_TC</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>f_TC</ci>
            </apply>
            <ci>diff_f_TC</ci>
         </apply>
         <apply>
            <eq/>
            <ci>diff_f_TMC</ci>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <cn type="e-notation" cellml:units="per_millimolar_per_second">227.7   <sep/>
                  3</cn>
                  <ci>Cai</ci>
                  <apply>
                     <minus/>
                     <apply>
                        <minus/>
                        <cn cellml:units="dimensionless">1</cn>
                        <ci>f_TMC</ci>
                     </apply>
                     <ci>f_TMM</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <cn type="e-notation" cellml:units="per_second">0.00751   <sep/>
                  3</cn>
                  <ci>f_TMC</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>f_TMC</ci>
            </apply>
            <ci>diff_f_TMC</ci>
         </apply>
         <apply>
            <eq/>
            <ci>diff_f_TMM</ci>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <cn type="e-notation" cellml:units="per_millimolar_per_second">2.277   <sep/>
                  3</cn>
                  <cn cellml:units="millimolar">2.5</cn>
                  <apply>
                     <minus/>
                     <apply>
                        <minus/>
                        <cn cellml:units="dimensionless">1</cn>
                        <ci>f_TMC</ci>
                     </apply>
                     <ci>f_TMM</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <cn type="e-notation" cellml:units="per_second">0.751   <sep/>
                  3</cn>
                  <ci>f_TMM</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>f_TMM</ci>
            </apply>
            <ci>diff_f_TMM</ci>
         </apply>
         <apply>
            <eq/>
            <ci>diff_f_CMi</ci>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <cn type="e-notation" cellml:units="per_millimolar_per_second">227.7   <sep/>
                  3</cn>
                  <ci>Cai</ci>
                  <apply>
                     <minus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <ci>f_CMi</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <cn type="e-notation" cellml:units="per_second">0.542   <sep/>
                  3</cn>
                  <ci>f_CMi</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>f_CMi</ci>
            </apply>
            <ci>diff_f_CMi</ci>
         </apply>
         <apply>
            <eq/>
            <ci>diff_f_CMs</ci>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <cn type="e-notation" cellml:units="per_millimolar_per_second">227.7   <sep/>
                  3</cn>
                  <ci>Casub</ci>
                  <apply>
                     <minus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <ci>f_CMs</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <cn type="e-notation" cellml:units="per_second">0.542   <sep/>
                  3</cn>
                  <ci>f_CMs</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>f_CMs</ci>
            </apply>
            <ci>diff_f_CMs</ci>
         </apply>
         <apply>
            <eq/>
            <ci>diff_f_CQ</ci>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <cn type="e-notation" cellml:units="per_millimolar_per_second">0.534   <sep/>
                  3</cn>
                  <ci>Ca_rel</ci>
                  <apply>
                     <minus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <ci>f_CQ</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <cn type="e-notation" cellml:units="per_second">0.445   <sep/>
                  3</cn>
                  <ci>f_CQ</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>f_CQ</ci>
            </apply>
            <ci>diff_f_CQ</ci>
         </apply>
         <apply>
            <eq/>
            <ci>diff_f_CSL</ci>
            <apply>
               <times/>
               <cn cellml:units="dimensionless">0.001</cn>
               <apply>
                  <minus/>
                  <apply>
                     <times/>
                     <cn cellml:units="per_millimolar_per_second">115</cn>
                     <ci>Casub</ci>
                     <apply>
                        <minus/>
                        <cn cellml:units="dimensionless">1</cn>
                        <ci>f_CSL</ci>
                     </apply>
                  </apply>
                  <apply>
                     <times/>
                     <cn cellml:units="per_second">1000</cn>
                     <ci>f_CSL</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>f_CSL</ci>
            </apply>
            <ci>diff_f_CSL</ci>
         </apply>
      </math>
   </component>
   <component name="extracellular_potassium_concentration">
      <variable name="Kc" units="millimolar" initial_value="5.4" public_interface="out"/>
   </component>
   <component name="intracellular_potassium_concentration">
      <variable name="Ki" units="millimolar" initial_value="140" public_interface="out"/>
   </component>
   <group>
      <relationship_ref relationship="encapsulation"/>
      <component_ref component="hyperpolarising_activated_current">
         <component_ref component="hyperpolarising_activated_current_y_gate"/>
      </component_ref>
      <component_ref component="rapid_delayed_rectifier_potassium_current">
         <component_ref component="rapid_delayed_rectifier_potassium_current_paf_gate"/>
         <component_ref component="rapid_delayed_rectifier_potassium_current_pas_gate"/>
         <component_ref component="rapid_delayed_rectifier_potassium_current_pik_gate"/>
      </component_ref>
      <component_ref component="fast_sodium_current">
         <component_ref component="fast_sodium_current_m_gate"/>
         <component_ref component="fast_sodium_current_h1_gate"/>
         <component_ref component="fast_sodium_current_h2_gate"/>
      </component_ref>
      <component_ref component="L_type_calcium_current">
         <component_ref component="L_type_calcium_current_d_gate"/>
         <component_ref component="L_type_calcium_current_f_gate"/>
         <component_ref component="L_type_calcium_current_f2_gate"/>
      </component_ref>
      <component_ref component="transient_outward_potassium_current">
         <component_ref component="transient_outward_potassium_current_r_gate"/>
         <component_ref component="transient_outward_potassium_current_qfast_gate"/>
         <component_ref component="transient_outward_potassium_current_qslow_gate"/>
      </component_ref>
      <component_ref component="sustained_outward_potassium_current">
         <component_ref component="sustained_outward_potassium_current_qa_gate"/>
         <component_ref component="sustained_outward_potassium_current_qi_gate"/>
      </component_ref>
      <component_ref component="acetylcholine_sensitive_current">
         <component_ref component="acetylcholine_sensitive_current_achf_gate"/>
         <component_ref component="acetylcholine_sensitive_current_achs_gate"/>
      </component_ref>
   </group>
   <connection>
      <map_components component_1="membrane" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="hyperpolarising_activated_current" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="rapid_delayed_rectifier_potassium_current" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="fast_sodium_current" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="L_type_calcium_current" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="acetylcholine_sensitive_current" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="transient_outward_potassium_current" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="sustained_outward_potassium_current" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="intracellular_calcium_concentration" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="hyperpolarising_activated_current" component_2="membrane"/>
      <map_variables variable_1="V" variable_2="V"/>
      <map_variables variable_1="i_f" variable_2="i_f"/>
   </connection>
   <connection>
      <map_components component_1="rapid_delayed_rectifier_potassium_current" component_2="membrane"/>
      <map_variables variable_1="V" variable_2="V"/>
      <map_variables variable_1="i_Kr" variable_2="i_Kr"/>
      <map_variables variable_1="RTONF" variable_2="RTONF"/>
   </connection>
   <connection>
      <map_components component_1="time_independent_potassium_current" component_2="membrane"/>
      <map_variables variable_1="V" variable_2="V"/>
      <map_variables variable_1="RTONF" variable_2="RTONF"/>
      <map_variables variable_1="i_K1" variable_2="i_K1"/>
   </connection>
   <connection>
      <map_components component_1="sodium_potassium_pump" component_2="membrane"/>
      <map_variables variable_1="V" variable_2="V"/>
      <map_variables variable_1="i_p" variable_2="i_p"/>
   </connection>
   <connection>
      <map_components component_1="sodium_calcium_exchange_current" component_2="membrane"/>
      <map_variables variable_1="V" variable_2="V"/>
      <map_variables variable_1="i_NaCa" variable_2="i_NaCa"/>
      <map_variables variable_1="RTONF" variable_2="RTONF"/>
   </connection>
   <connection>
      <map_components component_1="fast_sodium_current" component_2="membrane"/>
      <map_variables variable_1="V" variable_2="V"/>
      <map_variables variable_1="i_Na" variable_2="i_Na"/>
      <map_variables variable_1="RTONF" variable_2="RTONF"/>
      <map_variables variable_1="F" variable_2="F"/>
   </connection>
   <connection>
      <map_components component_1="L_type_calcium_current" component_2="membrane"/>
      <map_variables variable_1="V" variable_2="V"/>
      <map_variables variable_1="i_CaL" variable_2="i_CaL"/>
   </connection>
   <connection>
      <map_components component_1="transient_outward_potassium_current" component_2="membrane"/>
      <map_variables variable_1="V" variable_2="V"/>
      <map_variables variable_1="RTONF" variable_2="RTONF"/>
      <map_variables variable_1="i_to" variable_2="i_to"/>
   </connection>
   <connection>
      <map_components component_1="sustained_outward_potassium_current" component_2="membrane"/>
      <map_variables variable_1="V" variable_2="V"/>
      <map_variables variable_1="i_st" variable_2="i_st"/>
   </connection>
   <connection>
      <map_components component_1="background_current" component_2="membrane"/>
      <map_variables variable_1="V" variable_2="V"/>
      <map_variables variable_1="i_b" variable_2="i_b"/>
   </connection>
   <connection>
      <map_components component_1="acetylcholine_sensitive_current" component_2="membrane"/>
      <map_variables variable_1="V" variable_2="V"/>
      <map_variables variable_1="RTONF" variable_2="RTONF"/>
      <map_variables variable_1="i_ACh" variable_2="i_ACh"/>
   </connection>
   <connection>
      <map_components component_1="intracellular_calcium_concentration" component_2="membrane"/>
      <map_variables variable_1="V" variable_2="V"/>
      <map_variables variable_1="F" variable_2="F"/>
   </connection>
   <connection>
      <map_components component_1="rapid_delayed_rectifier_potassium_current" component_2="extracellular_potassium_concentration"/>
      <map_variables variable_1="Kc" variable_2="Kc"/>
   </connection>
   <connection>
      <map_components component_1="rapid_delayed_rectifier_potassium_current" component_2="intracellular_potassium_concentration"/>
      <map_variables variable_1="Ki" variable_2="Ki"/>
   </connection>
   <connection>
      <map_components component_1="transient_outward_potassium_current" component_2="extracellular_potassium_concentration"/>
      <map_variables variable_1="Kc" variable_2="Kc"/>
   </connection>
   <connection>
      <map_components component_1="transient_outward_potassium_current" component_2="intracellular_potassium_concentration"/>
      <map_variables variable_1="Ki" variable_2="Ki"/>
   </connection>
   <connection>
      <map_components component_1="time_independent_potassium_current" component_2="extracellular_potassium_concentration"/>
      <map_variables variable_1="Kc" variable_2="Kc"/>
   </connection>
   <connection>
      <map_components component_1="sodium_potassium_pump" component_2="intracellular_sodium_concentration"/>
      <map_variables variable_1="Nai" variable_2="Nai"/>
   </connection>
   <connection>
      <map_components component_1="sodium_potassium_pump" component_2="extracellular_potassium_concentration"/>
      <map_variables variable_1="Kc" variable_2="Kc"/>
   </connection>
   <connection>
      <map_components component_1="sodium_calcium_exchange_current" component_2="extracellular_sodium_concentration"/>
      <map_variables variable_1="Nao" variable_2="Nao"/>
   </connection>
   <connection>
      <map_components component_1="sodium_calcium_exchange_current" component_2="extracellular_calcium_concentration"/>
      <map_variables variable_1="Cao" variable_2="Cao"/>
   </connection>
   <connection>
      <map_components component_1="sodium_calcium_exchange_current" component_2="intracellular_sodium_concentration"/>
      <map_variables variable_1="Nai" variable_2="Nai"/>
   </connection>
   <connection>
      <map_components component_1="sodium_calcium_exchange_current" component_2="intracellular_calcium_concentration"/>
      <map_variables variable_1="i_NaCa" variable_2="i_NaCa"/>
      <map_variables variable_1="Casub" variable_2="Casub"/>
   </connection>
   <connection>
      <map_components component_1="fast_sodium_current" component_2="intracellular_sodium_concentration"/>
      <map_variables variable_1="Nai" variable_2="Nai"/>
   </connection>
   <connection>
      <map_components component_1="fast_sodium_current" component_2="intracellular_potassium_concentration"/>
      <map_variables variable_1="Ki" variable_2="Ki"/>
   </connection>
   <connection>
      <map_components component_1="fast_sodium_current" component_2="extracellular_potassium_concentration"/>
      <map_variables variable_1="Kc" variable_2="Kc"/>
   </connection>
   <connection>
      <map_components component_1="fast_sodium_current" component_2="extracellular_sodium_concentration"/>
      <map_variables variable_1="Nao" variable_2="Nao"/>
   </connection>
   <connection>
      <map_components component_1="L_type_calcium_current" component_2="intracellular_calcium_concentration"/>
      <map_variables variable_1="i_CaL" variable_2="i_CaL"/>
   </connection>
   <connection>
      <map_components component_1="rapid_delayed_rectifier_potassium_current" component_2="acetylcholine_sensitive_current"/>
      <map_variables variable_1="E_K" variable_2="E_K"/>
   </connection>
   <connection>
      <map_components component_1="acetylcholine_sensitive_current" component_2="extracellular_potassium_concentration"/>
      <map_variables variable_1="Kc" variable_2="Kc"/>
   </connection>
   <connection>
      <map_components component_1="L_type_calcium_current" component_2="acetylcholine_sensitive_current"/>
      <map_variables variable_1="i_ACh" variable_2="i_ACh"/>
      <map_variables variable_1="ACh" variable_2="ACh"/>
   </connection>
   <connection>
      <map_components component_1="hyperpolarising_activated_current" component_2="acetylcholine_sensitive_current"/>
      <map_variables variable_1="ACh" variable_2="ACh"/>
   </connection>
   <connection>
      <map_components component_1="hyperpolarising_activated_current" component_2="hyperpolarising_activated_current_y_gate"/>
      <map_variables variable_1="y" variable_2="y"/>
      <map_variables variable_1="time" variable_2="time"/>
      <map_variables variable_1="V" variable_2="V"/>
      <map_variables variable_1="ACh" variable_2="ACh"/>
   </connection>
   <connection>
      <map_components component_1="rapid_delayed_rectifier_potassium_current" component_2="rapid_delayed_rectifier_potassium_current_paf_gate"/>
      <map_variables variable_1="paf" variable_2="paf"/>
      <map_variables variable_1="time" variable_2="time"/>
      <map_variables variable_1="V" variable_2="V"/>
   </connection>
   <connection>
      <map_components component_1="rapid_delayed_rectifier_potassium_current" component_2="rapid_delayed_rectifier_potassium_current_pas_gate"/>
      <map_variables variable_1="pas" variable_2="pas"/>
      <map_variables variable_1="time" variable_2="time"/>
      <map_variables variable_1="V" variable_2="V"/>
   </connection>
   <connection>
      <map_components component_1="rapid_delayed_rectifier_potassium_current" component_2="rapid_delayed_rectifier_potassium_current_pik_gate"/>
      <map_variables variable_1="pik" variable_2="pik"/>
      <map_variables variable_1="time" variable_2="time"/>
      <map_variables variable_1="V" variable_2="V"/>
   </connection>
   <connection>
      <map_components component_1="fast_sodium_current" component_2="fast_sodium_current_m_gate"/>
      <map_variables variable_1="m" variable_2="m"/>
      <map_variables variable_1="time" variable_2="time"/>
      <map_variables variable_1="V" variable_2="V"/>
   </connection>
   <connection>
      <map_components component_1="fast_sodium_current" component_2="fast_sodium_current_h1_gate"/>
      <map_variables variable_1="h1" variable_2="h1"/>
      <map_variables variable_1="time" variable_2="time"/>
      <map_variables variable_1="V" variable_2="V"/>
   </connection>
   <connection>
      <map_components component_1="fast_sodium_current" component_2="fast_sodium_current_h2_gate"/>
      <map_variables variable_1="h2" variable_2="h2"/>
      <map_variables variable_1="time" variable_2="time"/>
      <map_variables variable_1="V" variable_2="V"/>
   </connection>
   <connection>
      <map_components component_1="L_type_calcium_current" component_2="L_type_calcium_current_d_gate"/>
      <map_variables variable_1="d" variable_2="d"/>
      <map_variables variable_1="time" variable_2="time"/>
      <map_variables variable_1="V" variable_2="V"/>
   </connection>
   <connection>
      <map_components component_1="L_type_calcium_current" component_2="L_type_calcium_current_f_gate"/>
      <map_variables variable_1="f" variable_2="f"/>
      <map_variables variable_1="time" variable_2="time"/>
      <map_variables variable_1="V" variable_2="V"/>
   </connection>
   <connection>
      <map_components component_1="L_type_calcium_current" component_2="L_type_calcium_current_f2_gate"/>
      <map_variables variable_1="f2" variable_2="f2"/>
      <map_variables variable_1="time" variable_2="time"/>
      <map_variables variable_1="V" variable_2="V"/>
   </connection>
   <connection>
      <map_components component_1="transient_outward_potassium_current" component_2="transient_outward_potassium_current_r_gate"/>
      <map_variables variable_1="r" variable_2="r"/>
      <map_variables variable_1="time" variable_2="time"/>
      <map_variables variable_1="V" variable_2="V"/>
   </connection>
   <connection>
      <map_components component_1="transient_outward_potassium_current" component_2="transient_outward_potassium_current_qfast_gate"/>
      <map_variables variable_1="q_fast" variable_2="q_fast"/>
      <map_variables variable_1="time" variable_2="time"/>
      <map_variables variable_1="V" variable_2="V"/>
   </connection>
   <connection>
      <map_components component_1="transient_outward_potassium_current" component_2="transient_outward_potassium_current_qslow_gate"/>
      <map_variables variable_1="q_slow" variable_2="q_slow"/>
      <map_variables variable_1="time" variable_2="time"/>
      <map_variables variable_1="V" variable_2="V"/>
   </connection>
   <connection>
      <map_components component_1="sustained_outward_potassium_current" component_2="sustained_outward_potassium_current_qa_gate"/>
      <map_variables variable_1="qa" variable_2="qa"/>
      <map_variables variable_1="time" variable_2="time"/>
      <map_variables variable_1="V" variable_2="V"/>
   </connection>
   <connection>
      <map_components component_1="sustained_outward_potassium_current" component_2="sustained_outward_potassium_current_qi_gate"/>
      <map_variables variable_1="qi" variable_2="qi"/>
      <map_variables variable_1="time" variable_2="time"/>
      <map_variables variable_1="V" variable_2="V"/>
   </connection>
   <connection>
      <map_components component_1="acetylcholine_sensitive_current" component_2="acetylcholine_sensitive_current_achf_gate"/>
      <map_variables variable_1="achf" variable_2="achf"/>
      <map_variables variable_1="time" variable_2="time"/>
      <map_variables variable_1="V" variable_2="V"/>
   </connection>
   <connection>
      <map_components component_1="acetylcholine_sensitive_current" component_2="acetylcholine_sensitive_current_achs_gate"/>
      <map_variables variable_1="achs" variable_2="achs"/>
      <map_variables variable_1="time" variable_2="time"/>
      <map_variables variable_1="V" variable_2="V"/>
   </connection>
</model>