Rendering of the source text

<?xml version="1.0" encoding="utf-8"?>
<!--
This CellML file was generated on 11/03/2008 at 16:22:14 using:

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

CellML 1.0 was used to generate this cellular model
http://www.CellML.org/
-->
<model xmlns="http://www.cellml.org/cellml/1.0#" xmlns:cmeta="http://www.cellml.org/metadata/1.0#" xmlns:dc="http://purl.org/dc/elements/1.1/" xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#" xmlns:bqs="http://www.cellml.org/bqs/1.0#" xmlns:cellml="http://www.cellml.org/cellml/1.0#" xmlns:dcterms="http://purl.org/dc/terms/" xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" xmlns:ns7="http://www.cellml.org/metadata/simulation/1.0#" cmeta:id="zhang_2000" name="zhang_2000">
   
<documentation xmlns="http://cellml.org/tmp-documentation">
<article>
  <articleinfo>
  <title>Mathematical Model Of Action Potentials In The Periphery Of The Sinoatrial Node, 2000</title>
  <author>
    <firstname>Penny</firstname>
          <surname>Noble</surname>
    <affiliation>
      <shortaffil>Oxford University</shortaffil>
    </affiliation>
  </author>
</articleinfo>
  <section id="sec_status">
    <title>Model Status</title>
    <para>
            The CellML version of this model has been curated and unit checked by Penny Noble and it is known to run in COR to recreate the published results. The Zhang model has since been updated by Alan Garny et al. (2003) and this model has also been described in CellML and can be found in the model repository.
</para>
<para>
ValidateCellML confirms this model as valid CellML with full unit consistency.
          </para>
  </section>
  <sect1 id="sec_structure">
<title>Model Structure</title>

<para>
ABSTRACT:  Mathematical models of the action potential in the periphery and center of the rabbit sinoatrial (SA) node have been developed on the basis of published experimental data. Simulated action potentials are consistent with those recorded experimentally: the model-generated peripheral action potential has a more negative takeoff potential, faster upstroke, more positive peak value, prominent phase 1 repolarization, greater amplitude, shorter duration, and more negative maximum diastolic potential than the model-generated central action potential. In addition, the model peripheral cell shows faster pacemaking. The models behave qualitatively the same as tissue from the periphery and center of the SA node in response to block of tetrodotoxin-sensitive Na(+) current, L- and T-type Ca(2+) currents, 4-aminopyridine-sensitive transient outward current, rapid and slow delayed rectifying K(+) currents, and hyperpolarization-activated current. A one-dimensional model of a string of SA node tissue, incorporating regional heterogeneity, coupled to a string of atrial tissue has been constructed to simulate the behavior of the intact SA node. In the one-dimensional model, the spontaneous action potential initiated in the center propagates to the periphery at approximately 0.06 m/s and then into the atrial muscle at 0.62 m/s.</para>

<para>The original paper is cited below.</para> 

<para>
Mathematical models of action potentials in the periphery and center of the rabbit sinoatrial node, H. Zhang, A.V. Holden, I. Honjo, M. Lei, T. Varghese and M.R. Boyett, 2000, <emphasis>American Journal of Physiology</emphasis>, 279, H397-H421.  <ulink url="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=pubmed&amp;cmd=Retrieve&amp;dopt=AbstractPlus&amp;list_uids=10899081&amp;query_hl=1&amp;itool=pubmed_docsum">PubMed ID: 10899081</ulink>
</para>

<informalfigure float="0" id="fig_cell_diagram">
<mediaobject>
  <imageobject>
    <objectinfo>
      <title>cell diagram of the Zhang et al SAN model showing ionic currents, pumps and exchangers within the sarcolemma and the sarcoplasmic reticulum</title>
    </objectinfo>
    <imagedata fileref="zhang_2000.png"/>
  </imageobject>
</mediaobject>
<caption>A schematic diagram describing the current flows across the cell membrane that are captured in the Zhang <emphasis>et al</emphasis> 2000 model of the action potentials in the peripheral cells of the SA node.</caption>
</informalfigure>

</sect1>
</article>
</documentation>
   <units name="per_second">
      <unit units="second" exponent="-1"/>
   </units>
   <units name="millivolt">
      <unit units="volt" prefix="milli"/>
   </units>
   <units name="per_millivolt">
      <unit units="volt" prefix="milli" exponent="-1"/>
   </units>
   <units name="per_millivolt_second">
      <unit units="millivolt" exponent="-1"/>
      <unit units="second" exponent="-1"/>
   </units>
   <units name="microS">
      <unit units="siemens" prefix="micro"/>
   </units>
   <units name="microF">
      <unit units="farad" prefix="micro"/>
   </units>
   <units name="nanoA">
      <unit units="ampere" prefix="nano"/>
   </units>
   <units name="millimolar">
      <unit units="mole" prefix="milli"/>
      <unit units="litre" exponent="-1"/>
   </units>
   <units name="millimolar4">
      <unit units="millimolar" exponent="4"/>
   </units>
   <units name="millijoule_per_mole_kelvin">
      <unit units="joule" prefix="milli"/>
      <unit units="mole" exponent="-1"/>
      <unit units="kelvin" exponent="-1"/>
   </units>
   <units name="coulomb_per_mole">
      <unit units="coulomb"/>
      <unit units="mole" exponent="-1"/>
   </units>
   <units name="microlitre_per_second">
      <unit units="litre" prefix="micro"/>
      <unit units="second" exponent="-1"/>
   </units>
   <component name="environment">
      <variable units="second" public_interface="out" name="time"/>
   </component>
   <component name="membrane">
      <variable units="dimensionless" public_interface="out" name="dCell" initial_value="0"/>
      <variable units="dimensionless" public_interface="out" name="FCell"/>
      <variable units="millivolt" public_interface="out" name="V" initial_value="-39.013558536"/>
      <variable units="millijoule_per_mole_kelvin" public_interface="out" name="R" initial_value="8314"/>
      <variable units="kelvin" public_interface="out" name="T" initial_value="310"/>
      <variable units="coulomb_per_mole" public_interface="out" name="F" initial_value="96845"/>
      <variable units="microF" name="Cm"/>
      <variable units="microF" name="CmCentre" initial_value="2e-5"/>
      <variable units="microF" name="CmPeriphery" initial_value="6.5e-5"/>
      <variable units="second" public_interface="in" name="time"/>
      <variable units="nanoA" public_interface="in" name="i_Na"/>
      <variable units="nanoA" public_interface="in" name="i_Ca_L"/>
      <variable units="nanoA" public_interface="in" name="i_Ca_T"/>
      <variable units="nanoA" public_interface="in" name="i_to"/>
      <variable units="nanoA" public_interface="in" name="i_sus"/>
      <variable units="nanoA" public_interface="in" name="i_K_r"/>
      <variable units="nanoA" public_interface="in" name="i_K_s"/>
      <variable units="nanoA" public_interface="in" name="i_f_Na"/>
      <variable units="nanoA" public_interface="in" name="i_f_K"/>
      <variable units="nanoA" public_interface="in" name="i_b_Na"/>
      <variable units="nanoA" public_interface="in" name="i_b_Ca"/>
      <variable units="nanoA" public_interface="in" name="i_b_K"/>
      <variable units="nanoA" public_interface="in" name="i_NaCa"/>
      <variable units="nanoA" public_interface="in" name="i_p"/>
      <variable units="nanoA" public_interface="in" name="i_Ca_p"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>FCell</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <cn cellml:units="dimensionless">1.07</cn>
                  <apply>
                     <minus/>
                     <apply>
                        <times/>
                        <cn cellml:units="dimensionless">3</cn>
                        <ci>dCell</ci>
                     </apply>
                     <cn cellml:units="dimensionless">0.1</cn>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <cn cellml:units="dimensionless">3</cn>
                  <apply>
                     <plus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <apply>
                        <times/>
                        <cn cellml:units="dimensionless">0.7745</cn>
                        <apply>
                           <exp/>
                           <apply>
                              <divide/>
                              <apply>
                                 <minus/>
                                 <apply>
                                    <minus/>
                                    <apply>
                                       <times/>
                                       <cn cellml:units="dimensionless">3</cn>
                                       <ci>dCell</ci>
                                    </apply>
                                    <cn cellml:units="dimensionless">2.05</cn>
                                 </apply>
                              </apply>
                              <cn cellml:units="dimensionless">0.295</cn>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>Cm</ci>
            <apply>
               <plus/>
               <ci>CmCentre</ci>
               <apply>
                  <times/>
                  <ci>FCell</ci>
                  <apply>
                     <minus/>
                     <ci>CmPeriphery</ci>
                     <ci>CmCentre</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>V</ci>
            </apply>
            <apply>
               <times/>
               <apply>
                  <divide/>
                  <apply>
                     <minus/>
                     <cn cellml:units="dimensionless">1</cn>
                  </apply>
                  <ci>Cm</ci>
               </apply>
               <apply>
                  <plus/>
                  <ci>i_Na</ci>
                  <ci>i_Ca_L</ci>
                  <ci>i_Ca_T</ci>
                  <ci>i_to</ci>
                  <ci>i_sus</ci>
                  <ci>i_K_r</ci>
                  <ci>i_K_s</ci>
                  <ci>i_f_Na</ci>
                  <ci>i_f_K</ci>
                  <ci>i_b_Na</ci>
                  <ci>i_b_Ca</ci>
                  <ci>i_b_K</ci>
                  <ci>i_NaCa</ci>
                  <ci>i_p</ci>
                  <ci>i_Ca_p</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="sodium_current">
      <variable units="dimensionless" public_interface="in" name="FCell"/>
      <variable units="nanoA" public_interface="out" name="i_Na"/>
      <variable units="microlitre_per_second" name="g_Na"/>
      <variable units="microlitre_per_second" name="g_Na_Centre" initial_value="0"/>
      <variable units="microlitre_per_second" name="g_Na_Periphery" initial_value="1.2e-6"/>
      <variable units="millivolt" public_interface="in" name="E_Na"/>
      <variable units="millimolar" public_interface="in" name="Na_o"/>
      <variable units="millijoule_per_mole_kelvin" public_interface="in" name="R"/>
      <variable units="coulomb_per_mole" public_interface="in" name="F"/>
      <variable units="kelvin" public_interface="in" name="T"/>
      <variable units="second" public_interface="in" private_interface="out" name="time"/>
      <variable units="millivolt" public_interface="in" private_interface="out" name="V"/>
      <variable units="dimensionless" private_interface="in" name="m"/>
      <variable units="dimensionless" private_interface="in" name="h"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>g_Na</ci>
            <apply>
               <plus/>
               <ci>g_Na_Centre</ci>
               <apply>
                  <times/>
                  <ci>FCell</ci>
                  <apply>
                     <minus/>
                     <ci>g_Na_Periphery</ci>
                     <ci>g_Na_Centre</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>i_Na</ci>
            <apply>
               <times/>
               <apply>
                  <divide/>
                  <apply>
                     <times/>
                     <apply>
                        <divide/>
                        <apply>
                           <times/>
                           <ci>g_Na</ci>
                           <apply>
                              <power/>
                              <ci>m</ci>
                              <cn cellml:units="dimensionless">3</cn>
                           </apply>
                           <ci>h</ci>
                           <ci>Na_o</ci>
                           <apply>
                              <power/>
                              <ci>F</ci>
                              <cn cellml:units="dimensionless">2</cn>
                           </apply>
                        </apply>
                        <apply>
                           <times/>
                           <ci>R</ci>
                           <ci>T</ci>
                        </apply>
                     </apply>
                     <apply>
                        <minus/>
                        <apply>
                           <exp/>
                           <apply>
                              <divide/>
                              <apply>
                                 <times/>
                                 <apply>
                                    <minus/>
                                    <ci>V</ci>
                                    <ci>E_Na</ci>
                                 </apply>
                                 <ci>F</ci>
                              </apply>
                              <apply>
                                 <times/>
                                 <ci>R</ci>
                                 <ci>T</ci>
                              </apply>
                           </apply>
                        </apply>
                        <cn cellml:units="dimensionless">1</cn>
                     </apply>
                  </apply>
                  <apply>
                     <minus/>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <times/>
                              <ci>V</ci>
                              <ci>F</ci>
                           </apply>
                           <apply>
                              <times/>
                              <ci>R</ci>
                              <ci>T</ci>
                           </apply>
                        </apply>
                     </apply>
                     <cn cellml:units="dimensionless">1</cn>
                  </apply>
               </apply>
               <ci>V</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="sodium_current_m_gate">
      <variable units="dimensionless" public_interface="out" name="m" initial_value="0.092361701692"/>
      <variable units="dimensionless" name="m_infinity"/>
      <variable units="second" name="tau_m"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="second" public_interface="in" name="time"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>m</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>m_infinity</ci>
                  <ci>m</ci>
               </apply>
               <ci>tau_m</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>m_infinity</ci>
            <apply>
               <power/>
               <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>
                           <cn cellml:units="millivolt">5.46</cn>
                        </apply>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <divide/>
                  <cn cellml:units="dimensionless">1</cn>
                  <cn cellml:units="dimensionless">3</cn>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>tau_m</ci>
            <apply>
               <plus/>
               <apply>
                  <divide/>
                  <cn cellml:units="second">0.0006247</cn>
                  <apply>
                     <plus/>
                     <apply>
                        <times/>
                        <cn cellml:units="dimensionless">0.832</cn>
                        <apply>
                           <exp/>
                           <apply>
                              <times/>
                              <apply>
                                 <minus/>
                                 <cn cellml:units="per_millivolt">0.335</cn>
                              </apply>
                              <apply>
                                 <plus/>
                                 <ci>V</ci>
                                 <cn cellml:units="millivolt">56.7</cn>
                              </apply>
                           </apply>
                        </apply>
                     </apply>
                     <apply>
                        <times/>
                        <cn cellml:units="dimensionless">0.627</cn>
                        <apply>
                           <exp/>
                           <apply>
                              <times/>
                              <cn cellml:units="per_millivolt">0.082</cn>
                              <apply>
                                 <plus/>
                                 <ci>V</ci>
                                 <cn cellml:units="millivolt">65.01</cn>
                              </apply>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </apply>
               <cn cellml:units="second" type="e-notation">4<sep/>-5</cn>
            </apply>
         </apply>
      </math>
   </component>
   <component name="sodium_current_h_gate">
      <variable units="dimensionless" public_interface="out" name="h"/>
      <variable units="dimensionless" name="F_Na"/>
      <variable units="dimensionless" name="h1" initial_value="0.015905380261"/>
      <variable units="dimensionless" name="h2" initial_value="0.01445216109"/>
      <variable units="dimensionless" name="h1_infinity"/>
      <variable units="dimensionless" name="h2_infinity"/>
      <variable units="second" name="tau_h1"/>
      <variable units="second" name="tau_h2"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="second" public_interface="in" name="time"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>F_Na</ci>
            <apply>
               <plus/>
               <apply>
                  <divide/>
                  <apply>
                     <times/>
                     <cn cellml:units="dimensionless">0.0952</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <times/>
                           <apply>
                              <minus/>
                              <cn cellml:units="per_millivolt">0.063</cn>
                           </apply>
                           <apply>
                              <plus/>
                              <ci>V</ci>
                              <cn cellml:units="millivolt">34.4</cn>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
                  <apply>
                     <plus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <apply>
                        <times/>
                        <cn cellml:units="dimensionless">1.66</cn>
                        <apply>
                           <exp/>
                           <apply>
                              <times/>
                              <apply>
                                 <minus/>
                                 <cn cellml:units="per_millivolt">0.225</cn>
                              </apply>
                              <apply>
                                 <plus/>
                                 <ci>V</ci>
                                 <cn cellml:units="millivolt">63.7</cn>
                              </apply>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </apply>
               <cn cellml:units="dimensionless">0.0869</cn>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>h</ci>
            <apply>
               <plus/>
               <apply>
                  <times/>
                  <apply>
                     <minus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <ci>F_Na</ci>
                  </apply>
                  <ci>h1</ci>
               </apply>
               <apply>
                  <times/>
                  <ci>F_Na</ci>
                  <ci>h2</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>h1</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>h1_infinity</ci>
                  <ci>h1</ci>
               </apply>
               <ci>tau_h1</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>h2</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>h2_infinity</ci>
                  <ci>h2</ci>
               </apply>
               <ci>tau_h2</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>h1_infinity</ci>
            <apply>
               <divide/>
               <cn cellml:units="dimensionless">1</cn>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <plus/>
                           <ci>V</ci>
                           <cn cellml:units="millivolt">66.1</cn>
                        </apply>
                        <cn cellml:units="millivolt">6.4</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>h2_infinity</ci>
            <ci>h1_infinity</ci>
         </apply>
         <apply>
            <eq/>
            <ci>tau_h1</ci>
            <apply>
               <plus/>
               <apply>
                  <divide/>
                  <apply>
                     <times/>
                     <cn cellml:units="second" type="e-notation">3.717<sep/>-6</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <times/>
                           <apply>
                              <minus/>
                              <cn cellml:units="per_millivolt">0.2815</cn>
                           </apply>
                           <apply>
                              <plus/>
                              <ci>V</ci>
                              <cn cellml:units="millivolt">17.11</cn>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
                  <apply>
                     <plus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <apply>
                        <times/>
                        <cn cellml:units="dimensionless">0.003732</cn>
                        <apply>
                           <exp/>
                           <apply>
                              <times/>
                              <apply>
                                 <minus/>
                                 <cn cellml:units="per_millivolt">0.3426</cn>
                              </apply>
                              <apply>
                                 <plus/>
                                 <ci>V</ci>
                                 <cn cellml:units="millivolt">37.76</cn>
                              </apply>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </apply>
               <cn cellml:units="second">0.0005977</cn>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>tau_h2</ci>
            <apply>
               <plus/>
               <apply>
                  <divide/>
                  <apply>
                     <times/>
                     <cn cellml:units="second" type="e-notation">3.186<sep/>-8</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <times/>
                           <apply>
                              <minus/>
                              <cn cellml:units="per_millivolt">0.6219</cn>
                           </apply>
                           <apply>
                              <plus/>
                              <ci>V</ci>
                              <cn cellml:units="millivolt">18.8</cn>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
                  <apply>
                     <plus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <apply>
                        <times/>
                        <cn cellml:units="dimensionless" type="e-notation">7.189<sep/>-5</cn>
                        <apply>
                           <exp/>
                           <apply>
                              <times/>
                              <apply>
                                 <minus/>
                                 <cn cellml:units="per_millivolt">0.6683</cn>
                              </apply>
                              <apply>
                                 <plus/>
                                 <ci>V</ci>
                                 <cn cellml:units="millivolt">34.07</cn>
                              </apply>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </apply>
               <cn cellml:units="second">0.003556</cn>
            </apply>
         </apply>
      </math>
   </component>
   <component name="L_type_Ca_channel">
      <variable units="dimensionless" public_interface="in" private_interface="out" name="FCell"/>
      <variable units="nanoA" public_interface="out" name="i_Ca_L"/>
      <variable units="microS" name="g_Ca_L_Centre" initial_value="0.0058"/>
      <variable units="microS" name="g_Ca_L_Periphery" initial_value="0.0659"/>
      <variable units="microS" name="g_Ca_L"/>
      <variable units="millivolt" name="E_Ca_L" initial_value="46.4"/>
      <variable units="second" public_interface="in" private_interface="out" name="time"/>
      <variable units="millivolt" public_interface="in" private_interface="out" name="V"/>
      <variable units="dimensionless" private_interface="in" name="d_L"/>
      <variable units="dimensionless" private_interface="in" name="f_L"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>g_Ca_L</ci>
            <apply>
               <plus/>
               <ci>g_Ca_L_Centre</ci>
               <apply>
                  <times/>
                  <ci>FCell</ci>
                  <apply>
                     <minus/>
                     <ci>g_Ca_L_Periphery</ci>
                     <ci>g_Ca_L_Centre</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>i_Ca_L</ci>
            <apply>
               <times/>
               <ci>g_Ca_L</ci>
               <apply>
                  <plus/>
                  <apply>
                     <times/>
                     <ci>f_L</ci>
                     <ci>d_L</ci>
                  </apply>
                  <apply>
                     <divide/>
                     <cn cellml:units="dimensionless">0.006</cn>
                     <apply>
                        <plus/>
                        <cn cellml:units="dimensionless">1</cn>
                        <apply>
                           <exp/>
                           <apply>
                              <divide/>
                              <apply>
                                 <minus/>
                                 <apply>
                                    <plus/>
                                    <ci>V</ci>
                                    <cn cellml:units="millivolt">14.1</cn>
                                 </apply>
                              </apply>
                              <cn cellml:units="millivolt">6</cn>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>E_Ca_L</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="L_type_Ca_channel_d_gate">
      <variable units="dimensionless" public_interface="in" name="FCell"/>
      <variable units="dimensionless" public_interface="out" name="d_L" initial_value="0.04804900895"/>
      <variable units="per_second" name="alpha_d_L"/>
      <variable units="per_second" name="beta_d_L"/>
      <variable units="dimensionless" name="d_L_infinity"/>
      <variable units="second" name="tau_d_L"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="second" public_interface="in" name="time"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>d_L</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>d_L_infinity</ci>
                  <ci>d_L</ci>
               </apply>
               <ci>tau_d_L</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>alpha_d_L</ci>
            <apply>
               <minus/>
               <apply>
                  <divide/>
                  <apply>
                     <times/>
                     <apply>
                        <minus/>
                        <cn cellml:units="per_millivolt_second">28.38</cn>
                     </apply>
                     <apply>
                        <plus/>
                        <ci>V</ci>
                        <cn cellml:units="millivolt">35</cn>
                     </apply>
                  </apply>
                  <apply>
                     <minus/>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <minus/>
                              <apply>
                                 <plus/>
                                 <ci>V</ci>
                                 <cn cellml:units="millivolt">35</cn>
                              </apply>
                           </apply>
                           <cn cellml:units="millivolt">2.5</cn>
                        </apply>
                     </apply>
                     <cn cellml:units="dimensionless">1</cn>
                  </apply>
               </apply>
               <apply>
                  <divide/>
                  <apply>
                     <times/>
                     <cn cellml:units="per_millivolt_second">84.9</cn>
                     <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_L</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <cn cellml:units="per_millivolt_second">11.42</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>tau_d_L</ci>
            <apply>
               <divide/>
               <cn cellml:units="dimensionless">2</cn>
               <apply>
                  <plus/>
                  <ci>alpha_d_L</ci>
                  <ci>beta_d_L</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>d_L_infinity</ci>
            <apply>
               <divide/>
               <cn cellml:units="dimensionless">1</cn>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <minus/>
                           <apply>
                              <plus/>
                              <ci>V</ci>
                              <cn cellml:units="millivolt">23.1</cn>
                           </apply>
                        </apply>
                        <cn cellml:units="millivolt">6</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="L_type_Ca_channel_f_gate">
      <variable units="dimensionless" public_interface="in" name="FCell"/>
      <variable units="dimensionless" public_interface="out" name="f_L" initial_value="0.48779845203"/>
      <variable units="per_second" name="alpha_f_L"/>
      <variable units="per_second" name="beta_f_L"/>
      <variable units="dimensionless" name="f_L_infinity"/>
      <variable units="second" name="tau_f_L"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="second" public_interface="in" name="time"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>f_L</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>f_L_infinity</ci>
                  <ci>f_L</ci>
               </apply>
               <ci>tau_f_L</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>alpha_f_L</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <cn cellml:units="per_millivolt_second">3.12</cn>
                  <apply>
                     <plus/>
                     <ci>V</ci>
                     <cn cellml:units="millivolt">28</cn>
                  </apply>
               </apply>
               <apply>
                  <minus/>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <plus/>
                           <ci>V</ci>
                           <cn cellml:units="millivolt">28</cn>
                        </apply>
                        <cn cellml:units="millivolt">4</cn>
                     </apply>
                  </apply>
                  <cn cellml:units="dimensionless">1</cn>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>beta_f_L</ci>
            <apply>
               <divide/>
               <cn cellml:units="per_second">25</cn>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <minus/>
                           <apply>
                              <plus/>
                              <ci>V</ci>
                              <cn cellml:units="millivolt">28</cn>
                           </apply>
                        </apply>
                        <cn cellml:units="millivolt">4</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>tau_f_L</ci>
            <apply>
               <divide/>
               <cn cellml:units="dimensionless">1</cn>
               <apply>
                  <plus/>
                  <ci>alpha_f_L</ci>
                  <ci>beta_f_L</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>f_L_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">45</cn>
                        </apply>
                        <cn cellml:units="millivolt">5</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="T_type_Ca_channel">
      <variable units="dimensionless" public_interface="in" private_interface="out" name="FCell"/>
      <variable units="nanoA" public_interface="out" name="i_Ca_T"/>
      <variable units="microS" name="g_Ca_T_Centre" initial_value="0.0043"/>
      <variable units="microS" name="g_Ca_T_Periphery" initial_value="0.0139"/>
      <variable units="microS" name="g_Ca_T"/>
      <variable units="millivolt" name="E_Ca_T" initial_value="45"/>
      <variable units="second" public_interface="in" private_interface="out" name="time"/>
      <variable units="millivolt" public_interface="in" private_interface="out" name="V"/>
      <variable units="dimensionless" private_interface="in" name="d_T"/>
      <variable units="dimensionless" private_interface="in" name="f_T"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>g_Ca_T</ci>
            <apply>
               <plus/>
               <ci>g_Ca_T_Centre</ci>
               <apply>
                  <times/>
                  <ci>FCell</ci>
                  <apply>
                     <minus/>
                     <ci>g_Ca_T_Periphery</ci>
                     <ci>g_Ca_T_Centre</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>i_Ca_T</ci>
            <apply>
               <times/>
               <ci>g_Ca_T</ci>
               <ci>d_T</ci>
               <ci>f_T</ci>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>E_Ca_T</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="T_type_Ca_channel_d_gate">
      <variable units="dimensionless" public_interface="out" name="d_T" initial_value="0.42074047435"/>
      <variable units="per_second" name="alpha_d_T"/>
      <variable units="per_second" name="beta_d_T"/>
      <variable units="dimensionless" name="d_T_infinity"/>
      <variable units="second" name="tau_d_T"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="second" public_interface="in" name="time"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>d_T</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>d_T_infinity</ci>
                  <ci>d_T</ci>
               </apply>
               <ci>tau_d_T</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>alpha_d_T</ci>
            <apply>
               <times/>
               <cn cellml:units="per_second">1068</cn>
               <apply>
                  <exp/>
                  <apply>
                     <divide/>
                     <apply>
                        <plus/>
                        <ci>V</ci>
                        <cn cellml:units="millivolt">26.3</cn>
                     </apply>
                     <cn cellml:units="millivolt">30</cn>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>beta_d_T</ci>
            <apply>
               <times/>
               <cn cellml:units="per_second">1068</cn>
               <apply>
                  <exp/>
                  <apply>
                     <divide/>
                     <apply>
                        <minus/>
                        <apply>
                           <plus/>
                           <ci>V</ci>
                           <cn cellml:units="millivolt">26.3</cn>
                        </apply>
                     </apply>
                     <cn cellml:units="millivolt">30</cn>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>tau_d_T</ci>
            <apply>
               <divide/>
               <cn cellml:units="dimensionless">1</cn>
               <apply>
                  <plus/>
                  <ci>alpha_d_T</ci>
                  <ci>beta_d_T</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>d_T_infinity</ci>
            <apply>
               <divide/>
               <cn cellml:units="dimensionless">1</cn>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <minus/>
                           <apply>
                              <plus/>
                              <ci>V</ci>
                              <cn cellml:units="millivolt">37</cn>
                           </apply>
                        </apply>
                        <cn cellml:units="millivolt">6.8</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="T_type_Ca_channel_f_gate">
      <variable units="dimensionless" public_interface="in" name="FCell"/>
      <variable units="dimensionless" public_interface="out" name="f_T" initial_value="0.038968420558"/>
      <variable units="per_second" name="alpha_f_T"/>
      <variable units="per_second" name="beta_f_T"/>
      <variable units="dimensionless" name="f_T_infinity"/>
      <variable units="second" name="tau_f_T"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="second" public_interface="in" name="time"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>f_T</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>f_T_infinity</ci>
                  <ci>f_T</ci>
               </apply>
               <ci>tau_f_T</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>alpha_f_T</ci>
            <apply>
               <times/>
               <cn cellml:units="per_second">15.3</cn>
               <apply>
                  <exp/>
                  <apply>
                     <divide/>
                     <apply>
                        <minus/>
                        <apply>
                           <plus/>
                           <ci>V</ci>
                           <cn cellml:units="millivolt">71.7</cn>
                        </apply>
                     </apply>
                     <cn cellml:units="millivolt">83.3</cn>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>beta_f_T</ci>
            <apply>
               <times/>
               <cn cellml:units="per_second">15</cn>
               <apply>
                  <exp/>
                  <apply>
                     <divide/>
                     <apply>
                        <plus/>
                        <ci>V</ci>
                        <cn cellml:units="millivolt">71.7</cn>
                     </apply>
                     <cn cellml:units="millivolt">15.38</cn>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>tau_f_T</ci>
            <apply>
               <divide/>
               <cn cellml:units="dimensionless">1</cn>
               <apply>
                  <plus/>
                  <ci>alpha_f_T</ci>
                  <ci>beta_f_T</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>f_T_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">71</cn>
                        </apply>
                        <cn cellml:units="millivolt">9</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="four_AP_sensitive_currents">
      <variable units="dimensionless" public_interface="in" private_interface="out" name="FCell"/>
      <variable units="nanoA" public_interface="out" name="i_to"/>
      <variable units="nanoA" public_interface="out" name="i_sus"/>
      <variable units="microS" name="g_to_Centre" initial_value="0.00491"/>
      <variable units="microS" name="g_to_Periphery" initial_value="0.03649"/>
      <variable units="microS" name="g_to"/>
      <variable units="microS" name="g_sus_Centre" initial_value="6.65e-5"/>
      <variable units="microS" name="g_sus_Periphery" initial_value="0.0114"/>
      <variable units="microS" name="g_sus"/>
      <variable units="millivolt" public_interface="in" name="E_K"/>
      <variable units="millivolt" public_interface="in" private_interface="out" name="V"/>
      <variable units="second" public_interface="in" private_interface="out" name="time"/>
      <variable units="dimensionless" private_interface="in" name="q"/>
      <variable units="dimensionless" private_interface="in" name="r"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>g_to</ci>
            <apply>
               <plus/>
               <ci>g_to_Centre</ci>
               <apply>
                  <times/>
                  <ci>FCell</ci>
                  <apply>
                     <minus/>
                     <ci>g_to_Periphery</ci>
                     <ci>g_to_Centre</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>g_sus</ci>
            <apply>
               <plus/>
               <ci>g_sus_Centre</ci>
               <apply>
                  <times/>
                  <ci>FCell</ci>
                  <apply>
                     <minus/>
                     <ci>g_sus_Periphery</ci>
                     <ci>g_sus_Centre</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>i_to</ci>
            <apply>
               <times/>
               <ci>g_to</ci>
               <ci>q</ci>
               <ci>r</ci>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>E_K</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>i_sus</ci>
            <apply>
               <times/>
               <ci>g_sus</ci>
               <ci>r</ci>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>E_K</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="four_AP_sensitive_currents_q_gate">
      <variable units="dimensionless" public_interface="in" name="FCell"/>
      <variable units="dimensionless" public_interface="out" name="q" initial_value="0.29760539675"/>
      <variable units="dimensionless" name="q_infinity"/>
      <variable units="second" name="tau_q"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="second" public_interface="in" name="time"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>q</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>q_infinity</ci>
                  <ci>q</ci>
               </apply>
               <ci>tau_q</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>q_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">59.37</cn>
                        </apply>
                        <cn cellml:units="millivolt">13.1</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>tau_q</ci>
            <apply>
               <plus/>
               <cn cellml:units="second">0.0101</cn>
               <apply>
                  <divide/>
                  <cn cellml:units="second">0.06517</cn>
                  <apply>
                     <times/>
                     <cn cellml:units="dimensionless">0.57</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <times/>
                           <apply>
                              <minus/>
                              <cn cellml:units="per_millivolt">0.08</cn>
                           </apply>
                           <apply>
                              <plus/>
                              <ci>V</ci>
                              <cn cellml:units="millivolt">49</cn>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <cn cellml:units="second" type="e-notation">2.4<sep/>-5</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <times/>
                        <cn cellml:units="per_millivolt">0.1</cn>
                        <apply>
                           <plus/>
                           <ci>V</ci>
                           <cn cellml:units="millivolt">50.93</cn>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="four_AP_sensitive_currents_r_gate">
      <variable units="dimensionless" public_interface="out" name="r" initial_value="0.064402950262"/>
      <variable units="dimensionless" name="r_infinity"/>
      <variable units="second" name="tau_r"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="second" public_interface="in" name="time"/>
      <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/>
                           <apply>
                              <minus/>
                              <ci>V</ci>
                              <cn cellml:units="millivolt">10.93</cn>
                           </apply>
                        </apply>
                        <cn cellml:units="millivolt">19.7</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>tau_r</ci>
            <apply>
               <times/>
               <cn cellml:units="dimensionless">0.001</cn>
               <apply>
                  <plus/>
                  <cn cellml:units="second">2.98</cn>
                  <apply>
                     <divide/>
                     <cn cellml:units="second">15.59</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.369</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>
         </apply>
      </math>
   </component>
   <component name="rapid_delayed_rectifying_potassium_current">
      <variable units="dimensionless" public_interface="in" name="FCell"/>
      <variable units="nanoA" public_interface="out" name="i_K_r"/>
      <variable units="microS" name="g_K_r_Centre" initial_value="0.000797"/>
      <variable units="microS" name="g_K_r_Periphery" initial_value="0.016"/>
      <variable units="microS" name="g_K_r"/>
      <variable units="dimensionless" name="P_a"/>
      <variable units="millivolt" public_interface="in" name="E_K"/>
      <variable units="second" public_interface="in" private_interface="out" name="time"/>
      <variable units="millivolt" public_interface="in" private_interface="out" name="V"/>
      <variable units="dimensionless" private_interface="in" name="P_af"/>
      <variable units="dimensionless" private_interface="in" name="P_as"/>
      <variable units="dimensionless" private_interface="in" name="P_i"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>g_K_r</ci>
            <apply>
               <plus/>
               <ci>g_K_r_Centre</ci>
               <apply>
                  <times/>
                  <ci>FCell</ci>
                  <apply>
                     <minus/>
                     <ci>g_K_r_Periphery</ci>
                     <ci>g_K_r_Centre</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>i_K_r</ci>
            <apply>
               <times/>
               <ci>g_K_r</ci>
               <ci>P_a</ci>
               <ci>P_i</ci>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>E_K</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>P_a</ci>
            <apply>
               <plus/>
               <apply>
                  <times/>
                  <cn cellml:units="dimensionless">0.6</cn>
                  <ci>P_af</ci>
               </apply>
               <apply>
                  <times/>
                  <cn cellml:units="dimensionless">0.4</cn>
                  <ci>P_as</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="rapid_delayed_rectifying_potassium_current_P_af_gate">
      <variable units="dimensionless" public_interface="out" name="P_af" initial_value="0.13034201158"/>
      <variable units="dimensionless" public_interface="out" name="P_af_infinity"/>
      <variable units="second" name="tau_P_af"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="second" public_interface="in" name="time"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>P_af</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>P_af_infinity</ci>
                  <ci>P_af</ci>
               </apply>
               <ci>tau_P_af</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>P_af_infinity</ci>
            <apply>
               <divide/>
               <cn cellml:units="dimensionless">1</cn>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <minus/>
                           <apply>
                              <plus/>
                              <ci>V</ci>
                              <cn cellml:units="millivolt">14.2</cn>
                           </apply>
                        </apply>
                        <cn cellml:units="millivolt">10.6</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>tau_P_af</ci>
            <apply>
               <divide/>
               <cn cellml:units="second">1</cn>
               <apply>
                  <plus/>
                  <apply>
                     <times/>
                     <cn cellml:units="dimensionless">37.2</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <minus/>
                              <ci>V</ci>
                              <cn cellml:units="millivolt">9</cn>
                           </apply>
                           <cn cellml:units="millivolt">15.9</cn>
                        </apply>
                     </apply>
                  </apply>
                  <apply>
                     <times/>
                     <cn cellml:units="dimensionless">0.96</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <minus/>
                              <apply>
                                 <minus/>
                                 <ci>V</ci>
                                 <cn cellml:units="millivolt">9</cn>
                              </apply>
                           </apply>
                           <cn cellml:units="millivolt">22.5</cn>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="rapid_delayed_rectifying_potassium_current_P_as_gate">
      <variable units="dimensionless" public_interface="out" name="P_as" initial_value="0.46960956028"/>
      <variable units="dimensionless" name="P_as_infinity"/>
      <variable units="second" name="tau_P_as"/>
      <variable units="dimensionless" public_interface="in" name="P_af_infinity"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="second" public_interface="in" name="time"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>P_as</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>P_as_infinity</ci>
                  <ci>P_as</ci>
               </apply>
               <ci>tau_P_as</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>P_as_infinity</ci>
            <ci>P_af_infinity</ci>
         </apply>
         <apply>
            <eq/>
            <ci>tau_P_as</ci>
            <apply>
               <divide/>
               <cn cellml:units="second">1</cn>
               <apply>
                  <plus/>
                  <apply>
                     <times/>
                     <cn cellml:units="dimensionless">4.2</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <minus/>
                              <ci>V</ci>
                              <cn cellml:units="millivolt">9</cn>
                           </apply>
                           <cn cellml:units="millivolt">17</cn>
                        </apply>
                     </apply>
                  </apply>
                  <apply>
                     <times/>
                     <cn cellml:units="dimensionless">0.15</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <minus/>
                              <apply>
                                 <minus/>
                                 <ci>V</ci>
                                 <cn cellml:units="millivolt">9</cn>
                              </apply>
                           </apply>
                           <cn cellml:units="millivolt">21.6</cn>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="rapid_delayed_rectifying_potassium_current_P_i_gate">
      <variable units="dimensionless" public_interface="out" name="P_i" initial_value="0.87993375273"/>
      <variable units="dimensionless" name="P_i_infinity"/>
      <variable units="second" name="tau_P_i" initial_value="0.002"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="second" public_interface="in" name="time"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>P_i</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>P_i_infinity</ci>
                  <ci>P_i</ci>
               </apply>
               <ci>tau_P_i</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>P_i_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">18.6</cn>
                        </apply>
                        <cn cellml:units="millivolt">10.1</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="slow_delayed_rectifying_potassium_current">
      <variable units="dimensionless" public_interface="in" name="FCell"/>
      <variable units="nanoA" public_interface="out" name="i_K_s"/>
      <variable units="microS" name="g_K_s_Centre" initial_value="0.000518"/>
      <variable units="microS" name="g_K_s_Periphery" initial_value="0.0104"/>
      <variable units="microS" name="g_K_s"/>
      <variable units="millivolt" public_interface="in" name="E_K_s"/>
      <variable units="second" public_interface="in" private_interface="out" name="time"/>
      <variable units="millivolt" public_interface="in" private_interface="out" name="V"/>
      <variable units="dimensionless" private_interface="in" name="xs"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>g_K_s</ci>
            <apply>
               <plus/>
               <ci>g_K_s_Centre</ci>
               <apply>
                  <times/>
                  <ci>FCell</ci>
                  <apply>
                     <minus/>
                     <ci>g_K_s_Periphery</ci>
                     <ci>g_K_s_Centre</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>i_K_s</ci>
            <apply>
               <times/>
               <ci>g_K_s</ci>
               <apply>
                  <power/>
                  <ci>xs</ci>
                  <cn cellml:units="dimensionless">2</cn>
               </apply>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>E_K_s</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="slow_delayed_rectifying_potassium_current_xs_gate">
      <variable units="dimensionless" public_interface="out" name="xs" initial_value="0.082293827208"/>
      <variable units="per_second" name="alpha_xs"/>
      <variable units="per_second" name="beta_xs"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="second" public_interface="in" name="time"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>xs</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>alpha_xs</ci>
                  <apply>
                     <minus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <ci>xs</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>beta_xs</ci>
                  <ci>xs</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>alpha_xs</ci>
            <apply>
               <divide/>
               <cn cellml:units="per_second">14</cn>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <minus/>
                           <apply>
                              <minus/>
                              <ci>V</ci>
                              <cn cellml:units="millivolt">40</cn>
                           </apply>
                        </apply>
                        <cn cellml:units="millivolt">9</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>beta_xs</ci>
            <apply>
               <times/>
               <cn cellml:units="per_second">1</cn>
               <apply>
                  <exp/>
                  <apply>
                     <divide/>
                     <apply>
                        <minus/>
                        <ci>V</ci>
                     </apply>
                     <cn cellml:units="millivolt">45</cn>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="hyperpolarisation_activated_current">
      <variable units="dimensionless" public_interface="in" name="FCell"/>
      <variable units="nanoA" public_interface="out" name="i_f_Na"/>
      <variable units="nanoA" public_interface="out" name="i_f_K"/>
      <variable units="microS" name="g_f_Na_Centre" initial_value="0.000548"/>
      <variable units="microS" name="g_f_Na_Periphery" initial_value="0.0069"/>
      <variable units="microS" name="g_f_Na"/>
      <variable units="microS" name="g_f_K_Centre" initial_value="0.000548"/>
      <variable units="microS" name="g_f_K_Periphery" initial_value="0.0069"/>
      <variable units="microS" name="g_f_K"/>
      <variable units="millivolt" public_interface="in" name="E_K"/>
      <variable units="millivolt" public_interface="in" name="E_Na"/>
      <variable units="second" public_interface="in" private_interface="out" name="time"/>
      <variable units="millivolt" public_interface="in" private_interface="out" name="V"/>
      <variable units="dimensionless" private_interface="in" name="y"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>g_f_Na</ci>
            <apply>
               <plus/>
               <ci>g_f_Na_Centre</ci>
               <apply>
                  <times/>
                  <ci>FCell</ci>
                  <apply>
                     <minus/>
                     <ci>g_f_Na_Periphery</ci>
                     <ci>g_f_Na_Centre</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>i_f_Na</ci>
            <apply>
               <times/>
               <ci>g_f_Na</ci>
               <ci>y</ci>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>E_Na</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>g_f_K</ci>
            <apply>
               <plus/>
               <ci>g_f_K_Centre</ci>
               <apply>
                  <times/>
                  <ci>FCell</ci>
                  <apply>
                     <minus/>
                     <ci>g_f_K_Periphery</ci>
                     <ci>g_f_K_Centre</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>i_f_K</ci>
            <apply>
               <times/>
               <ci>g_f_K</ci>
               <ci>y</ci>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>E_K</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="hyperpolarisation_activated_current_y_gate">
      <variable units="dimensionless" public_interface="out" name="y" initial_value="0.03889291759"/>
      <variable units="per_second" name="alpha_y"/>
      <variable units="per_second" name="beta_y"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="second" public_interface="in" name="time"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>y</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>alpha_y</ci>
                  <apply>
                     <minus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <ci>y</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>beta_y</ci>
                  <ci>y</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>alpha_y</ci>
            <apply>
               <times/>
               <cn cellml:units="per_second">1</cn>
               <apply>
                  <exp/>
                  <apply>
                     <divide/>
                     <apply>
                        <minus/>
                        <apply>
                           <plus/>
                           <ci>V</ci>
                           <cn cellml:units="millivolt">78.91</cn>
                        </apply>
                     </apply>
                     <cn cellml:units="millivolt">26.62</cn>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>beta_y</ci>
            <apply>
               <times/>
               <cn cellml:units="per_second">1</cn>
               <apply>
                  <exp/>
                  <apply>
                     <divide/>
                     <apply>
                        <plus/>
                        <ci>V</ci>
                        <cn cellml:units="millivolt">75.13</cn>
                     </apply>
                     <cn cellml:units="millivolt">21.25</cn>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="sodium_background_current">
      <variable units="dimensionless" public_interface="in" name="FCell"/>
      <variable units="nanoA" public_interface="out" name="i_b_Na"/>
      <variable units="microS" name="g_b_Na_Centre" initial_value="5.8e-5"/>
      <variable units="microS" name="g_b_Na_Periphery" initial_value="0.000189"/>
      <variable units="microS" name="g_b_Na"/>
      <variable units="millivolt" public_interface="in" name="E_Na"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>g_b_Na</ci>
            <apply>
               <plus/>
               <ci>g_b_Na_Centre</ci>
               <apply>
                  <times/>
                  <ci>FCell</ci>
                  <apply>
                     <minus/>
                     <ci>g_b_Na_Periphery</ci>
                     <ci>g_b_Na_Centre</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>i_b_Na</ci>
            <apply>
               <times/>
               <ci>g_b_Na</ci>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>E_Na</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="potassium_background_current">
      <variable units="dimensionless" public_interface="in" name="FCell"/>
      <variable units="nanoA" public_interface="out" name="i_b_K"/>
      <variable units="microS" name="g_b_K_Centre" initial_value="2.52e-5"/>
      <variable units="microS" name="g_b_K_Periphery" initial_value="8.19e-5"/>
      <variable units="microS" name="g_b_K"/>
      <variable units="millivolt" public_interface="in" name="E_K"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>g_b_K</ci>
            <apply>
               <plus/>
               <ci>g_b_K_Centre</ci>
               <apply>
                  <times/>
                  <ci>FCell</ci>
                  <apply>
                     <minus/>
                     <ci>g_b_K_Periphery</ci>
                     <ci>g_b_K_Centre</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>i_b_K</ci>
            <apply>
               <times/>
               <ci>g_b_K</ci>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>E_K</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="calcium_background_current">
      <variable units="dimensionless" public_interface="in" name="FCell"/>
      <variable units="nanoA" public_interface="out" name="i_b_Ca"/>
      <variable units="microS" name="g_b_Ca_Centre" initial_value="1.32e-5"/>
      <variable units="microS" name="g_b_Ca_Periphery" initial_value="4.3e-5"/>
      <variable units="microS" name="g_b_Ca"/>
      <variable units="millivolt" public_interface="in" name="E_Ca"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>g_b_Ca</ci>
            <apply>
               <plus/>
               <ci>g_b_Ca_Centre</ci>
               <apply>
                  <times/>
                  <ci>FCell</ci>
                  <apply>
                     <minus/>
                     <ci>g_b_Ca_Periphery</ci>
                     <ci>g_b_Ca_Centre</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>i_b_Ca</ci>
            <apply>
               <times/>
               <ci>g_b_Ca</ci>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>E_Ca</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="sodium_calcium_exchanger">
      <variable units="dimensionless" public_interface="in" name="FCell"/>
      <variable units="nanoA" public_interface="out" name="i_NaCa"/>
      <variable units="nanoA" name="k_NaCa_Centre" initial_value="2.7e-6"/>
      <variable units="nanoA" name="k_NaCa_Periphery" initial_value="8.8e-6"/>
      <variable units="nanoA" name="k_NaCa"/>
      <variable units="dimensionless" name="d_NaCa" initial_value="0.0001"/>
      <variable units="dimensionless" name="gamma_NaCa" initial_value="0.5"/>
      <variable units="millimolar" public_interface="in" name="Na_i"/>
      <variable units="millimolar" public_interface="in" name="Na_o"/>
      <variable units="millimolar" public_interface="in" name="Ca_i"/>
      <variable units="millimolar" public_interface="in" name="Ca_o"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>k_NaCa</ci>
            <apply>
               <plus/>
               <ci>k_NaCa_Centre</ci>
               <apply>
                  <times/>
                  <ci>FCell</ci>
                  <apply>
                     <minus/>
                     <ci>k_NaCa_Periphery</ci>
                     <ci>k_NaCa_Centre</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>i_NaCa</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <ci>k_NaCa</ci>
                  <apply>
                     <minus/>
                     <apply>
                        <times/>
                        <apply>
                           <power/>
                           <ci>Na_i</ci>
                           <cn cellml:units="dimensionless">3</cn>
                        </apply>
                        <ci>Ca_o</ci>
                        <apply>
                           <exp/>
                           <apply>
                              <times/>
                              <cn cellml:units="per_millivolt">0.03743</cn>
                              <ci>V</ci>
                              <ci>gamma_NaCa</ci>
                           </apply>
                        </apply>
                     </apply>
                     <apply>
                        <times/>
                        <apply>
                           <power/>
                           <ci>Na_o</ci>
                           <cn cellml:units="dimensionless">3</cn>
                        </apply>
                        <ci>Ca_i</ci>
                        <apply>
                           <exp/>
                           <apply>
                              <times/>
                              <cn cellml:units="per_millivolt">0.0374</cn>
                              <ci>V</ci>
                              <apply>
                                 <minus/>
                                 <ci>gamma_NaCa</ci>
                                 <cn cellml:units="dimensionless">1</cn>
                              </apply>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <plus/>
                  <cn cellml:units="millimolar4">1</cn>
                  <apply>
                     <times/>
                     <ci>d_NaCa</ci>
                     <apply>
                        <plus/>
                        <apply>
                           <times/>
                           <ci>Ca_i</ci>
                           <apply>
                              <power/>
                              <ci>Na_o</ci>
                              <cn cellml:units="dimensionless">3</cn>
                           </apply>
                        </apply>
                        <apply>
                           <times/>
                           <ci>Ca_o</ci>
                           <apply>
                              <power/>
                              <ci>Na_i</ci>
                              <cn cellml:units="dimensionless">3</cn>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="sodium_potassium_pump">
      <variable units="dimensionless" public_interface="in" name="FCell"/>
      <variable units="nanoA" public_interface="out" name="i_p"/>
      <variable units="millimolar" name="K_m_Na" initial_value="5.64"/>
      <variable units="millimolar" name="K_m_K" initial_value="0.621"/>
      <variable units="nanoA" name="i_p_max_Centre" initial_value="0.0478"/>
      <variable units="nanoA" name="i_p_max_Periphery" initial_value="0.16"/>
      <variable units="nanoA" name="i_p_max"/>
      <variable units="millimolar" public_interface="in" name="Na_i"/>
      <variable units="millimolar" public_interface="in" name="K_o"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_p_max</ci>
            <apply>
               <plus/>
               <ci>i_p_max_Centre</ci>
               <apply>
                  <times/>
                  <ci>FCell</ci>
                  <apply>
                     <minus/>
                     <ci>i_p_max_Periphery</ci>
                     <ci>i_p_max_Centre</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>i_p</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <ci>i_p_max</ci>
                  <apply>
                     <power/>
                     <apply>
                        <divide/>
                        <ci>Na_i</ci>
                        <apply>
                           <plus/>
                           <ci>K_m_Na</ci>
                           <ci>Na_i</ci>
                        </apply>
                     </apply>
                     <cn cellml:units="dimensionless">3</cn>
                  </apply>
                  <apply>
                     <power/>
                     <apply>
                        <divide/>
                        <ci>K_o</ci>
                        <apply>
                           <plus/>
                           <ci>K_m_K</ci>
                           <ci>K_o</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="persistent_calcium_current">
      <variable units="dimensionless" public_interface="in" name="FCell"/>
      <variable units="nanoA" public_interface="out" name="i_Ca_p"/>
      <variable units="nanoA" name="i_Ca_p_max_Centre" initial_value="0"/>
      <variable units="nanoA" name="i_Ca_p_max_Periphery" initial_value="0"/>
      <variable units="nanoA" name="i_Ca_p_max"/>
      <variable units="millimolar" public_interface="in" name="Ca_i"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_Ca_p_max</ci>
            <apply>
               <plus/>
               <ci>i_Ca_p_max_Centre</ci>
               <apply>
                  <times/>
                  <ci>FCell</ci>
                  <apply>
                     <minus/>
                     <ci>i_Ca_p_max_Periphery</ci>
                     <ci>i_Ca_p_max_Centre</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>i_Ca_p</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <ci>i_Ca_p_max</ci>
                  <ci>Ca_i</ci>
               </apply>
               <apply>
                  <plus/>
                  <ci>Ca_i</ci>
                  <cn cellml:units="millimolar">0.0004</cn>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="ionic_concentrations">
      <variable units="millimolar" public_interface="out" name="Na_i" initial_value="8"/>
      <variable units="millimolar" public_interface="out" name="Na_o" initial_value="140"/>
      <variable units="millimolar" public_interface="out" name="Ca_i" initial_value="0.0001"/>
      <variable units="millimolar" public_interface="out" name="Ca_o" initial_value="2"/>
      <variable units="millimolar" public_interface="out" name="K_i" initial_value="140"/>
      <variable units="millimolar" public_interface="out" name="K_o" initial_value="5.4"/>
   </component>
   <component name="reversal_and_equilibrium_potentials">
      <variable units="millivolt" public_interface="out" name="E_Na"/>
      <variable units="millivolt" public_interface="out" name="E_K"/>
      <variable units="millivolt" public_interface="out" name="E_K_s"/>
      <variable units="millivolt" public_interface="out" name="E_Ca"/>
      <variable units="millimolar" public_interface="in" name="K_o"/>
      <variable units="millimolar" public_interface="in" name="Na_o"/>
      <variable units="millimolar" public_interface="in" name="K_i"/>
      <variable units="millimolar" public_interface="in" name="Na_i"/>
      <variable units="millimolar" public_interface="in" name="Ca_o"/>
      <variable units="millimolar" public_interface="in" name="Ca_i"/>
      <variable units="millijoule_per_mole_kelvin" public_interface="in" name="R"/>
      <variable units="coulomb_per_mole" public_interface="in" name="F"/>
      <variable units="kelvin" public_interface="in" name="T"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>E_Na</ci>
            <apply>
               <times/>
               <apply>
                  <divide/>
                  <apply>
                     <times/>
                     <ci>R</ci>
                     <ci>T</ci>
                  </apply>
                  <ci>F</ci>
               </apply>
               <apply>
                  <ln/>
                  <apply>
                     <divide/>
                     <ci>Na_o</ci>
                     <ci>Na_i</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>E_K</ci>
            <apply>
               <times/>
               <apply>
                  <divide/>
                  <apply>
                     <times/>
                     <ci>R</ci>
                     <ci>T</ci>
                  </apply>
                  <ci>F</ci>
               </apply>
               <apply>
                  <ln/>
                  <apply>
                     <divide/>
                     <ci>K_o</ci>
                     <ci>K_i</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>E_Ca</ci>
            <apply>
               <times/>
               <apply>
                  <divide/>
                  <apply>
                     <times/>
                     <ci>R</ci>
                     <ci>T</ci>
                  </apply>
                  <apply>
                     <times/>
                     <cn cellml:units="dimensionless">2</cn>
                     <ci>F</ci>
                  </apply>
               </apply>
               <apply>
                  <ln/>
                  <apply>
                     <divide/>
                     <ci>Ca_o</ci>
                     <ci>Ca_i</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>E_K_s</ci>
            <apply>
               <times/>
               <apply>
                  <divide/>
                  <apply>
                     <times/>
                     <ci>R</ci>
                     <ci>T</ci>
                  </apply>
                  <ci>F</ci>
               </apply>
               <apply>
                  <ln/>
                  <apply>
                     <divide/>
                     <apply>
                        <plus/>
                        <ci>K_o</ci>
                        <apply>
                           <times/>
                           <cn cellml:units="dimensionless">0.12</cn>
                           <ci>Na_o</ci>
                        </apply>
                     </apply>
                     <apply>
                        <plus/>
                        <ci>K_i</ci>
                        <apply>
                           <times/>
                           <cn cellml:units="dimensionless">0.12</cn>
                           <ci>Na_i</ci>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <group>
      <relationship_ref relationship="containment"/>
      <component_ref component="membrane">
         <component_ref component="sodium_current">
            <component_ref component="sodium_current_m_gate"/>
            <component_ref component="sodium_current_h_gate"/>
         </component_ref>
         <component_ref component="L_type_Ca_channel">
            <component_ref component="L_type_Ca_channel_d_gate"/>
            <component_ref component="L_type_Ca_channel_f_gate"/>
         </component_ref>
         <component_ref component="T_type_Ca_channel">
            <component_ref component="T_type_Ca_channel_d_gate"/>
            <component_ref component="T_type_Ca_channel_f_gate"/>
         </component_ref>
         <component_ref component="four_AP_sensitive_currents">
            <component_ref component="four_AP_sensitive_currents_q_gate"/>
            <component_ref component="four_AP_sensitive_currents_r_gate"/>
         </component_ref>
         <component_ref component="rapid_delayed_rectifying_potassium_current">
            <component_ref component="rapid_delayed_rectifying_potassium_current_P_af_gate"/>
            <component_ref component="rapid_delayed_rectifying_potassium_current_P_as_gate"/>
            <component_ref component="rapid_delayed_rectifying_potassium_current_P_i_gate"/>
         </component_ref>
         <component_ref component="slow_delayed_rectifying_potassium_current">
            <component_ref component="slow_delayed_rectifying_potassium_current_xs_gate"/>
         </component_ref>
         <component_ref component="hyperpolarisation_activated_current">
            <component_ref component="hyperpolarisation_activated_current_y_gate"/>
         </component_ref>
         <component_ref component="sodium_background_current"/>
         <component_ref component="calcium_background_current"/>
         <component_ref component="potassium_background_current"/>
         <component_ref component="sodium_calcium_exchanger"/>
         <component_ref component="sodium_potassium_pump"/>
         <component_ref component="persistent_calcium_current"/>
         <component_ref component="ionic_concentrations"/>
         <component_ref component="reversal_and_equilibrium_potentials"/>
      </component_ref>
   </group>
   <group>
      <relationship_ref relationship="encapsulation"/>
      <component_ref component="sodium_current">
         <component_ref component="sodium_current_m_gate"/>
         <component_ref component="sodium_current_h_gate"/>
      </component_ref>
      <component_ref component="L_type_Ca_channel">
         <component_ref component="L_type_Ca_channel_d_gate"/>
         <component_ref component="L_type_Ca_channel_f_gate"/>
      </component_ref>
      <component_ref component="T_type_Ca_channel">
         <component_ref component="T_type_Ca_channel_d_gate"/>
         <component_ref component="T_type_Ca_channel_f_gate"/>
      </component_ref>
      <component_ref component="four_AP_sensitive_currents">
         <component_ref component="four_AP_sensitive_currents_q_gate"/>
         <component_ref component="four_AP_sensitive_currents_r_gate"/>
      </component_ref>
      <component_ref component="rapid_delayed_rectifying_potassium_current">
         <component_ref component="rapid_delayed_rectifying_potassium_current_P_af_gate"/>
         <component_ref component="rapid_delayed_rectifying_potassium_current_P_as_gate"/>
         <component_ref component="rapid_delayed_rectifying_potassium_current_P_i_gate"/>
      </component_ref>
      <component_ref component="slow_delayed_rectifying_potassium_current">
         <component_ref component="slow_delayed_rectifying_potassium_current_xs_gate"/>
      </component_ref>
      <component_ref component="hyperpolarisation_activated_current">
         <component_ref component="hyperpolarisation_activated_current_y_gate"/>
      </component_ref>
   </group>
   <connection>
      <map_components component_2="environment" component_1="membrane"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="environment" component_1="sodium_current"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="environment" component_1="L_type_Ca_channel"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="environment" component_1="T_type_Ca_channel"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="environment" component_1="four_AP_sensitive_currents"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="environment" component_1="rapid_delayed_rectifying_potassium_current"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="environment" component_1="slow_delayed_rectifying_potassium_current"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="environment" component_1="hyperpolarisation_activated_current"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="sodium_current" component_1="membrane"/>
      <map_variables variable_2="FCell" variable_1="FCell"/>
      <map_variables variable_2="i_Na" variable_1="i_Na"/>
      <map_variables variable_2="V" variable_1="V"/>
      <map_variables variable_2="R" variable_1="R"/>
      <map_variables variable_2="T" variable_1="T"/>
      <map_variables variable_2="F" variable_1="F"/>
   </connection>
   <connection>
      <map_components component_2="L_type_Ca_channel" component_1="membrane"/>
      <map_variables variable_2="FCell" variable_1="FCell"/>
      <map_variables variable_2="i_Ca_L" variable_1="i_Ca_L"/>
      <map_variables variable_2="V" variable_1="V"/>
   </connection>
   <connection>
      <map_components component_2="T_type_Ca_channel" component_1="membrane"/>
      <map_variables variable_2="FCell" variable_1="FCell"/>
      <map_variables variable_2="i_Ca_T" variable_1="i_Ca_T"/>
      <map_variables variable_2="V" variable_1="V"/>
   </connection>
   <connection>
      <map_components component_2="four_AP_sensitive_currents" component_1="membrane"/>
      <map_variables variable_2="FCell" variable_1="FCell"/>
      <map_variables variable_2="i_to" variable_1="i_to"/>
      <map_variables variable_2="i_sus" variable_1="i_sus"/>
      <map_variables variable_2="V" variable_1="V"/>
   </connection>
   <connection>
      <map_components component_2="rapid_delayed_rectifying_potassium_current" component_1="membrane"/>
      <map_variables variable_2="FCell" variable_1="FCell"/>
      <map_variables variable_2="i_K_r" variable_1="i_K_r"/>
      <map_variables variable_2="V" variable_1="V"/>
   </connection>
   <connection>
      <map_components component_2="slow_delayed_rectifying_potassium_current" component_1="membrane"/>
      <map_variables variable_2="FCell" variable_1="FCell"/>
      <map_variables variable_2="i_K_s" variable_1="i_K_s"/>
      <map_variables variable_2="V" variable_1="V"/>
   </connection>
   <connection>
      <map_components component_2="hyperpolarisation_activated_current" component_1="membrane"/>
      <map_variables variable_2="FCell" variable_1="FCell"/>
      <map_variables variable_2="i_f_Na" variable_1="i_f_Na"/>
      <map_variables variable_2="i_f_K" variable_1="i_f_K"/>
      <map_variables variable_2="V" variable_1="V"/>
   </connection>
   <connection>
      <map_components component_2="sodium_background_current" component_1="membrane"/>
      <map_variables variable_2="FCell" variable_1="FCell"/>
      <map_variables variable_2="V" variable_1="V"/>
      <map_variables variable_2="i_b_Na" variable_1="i_b_Na"/>
   </connection>
   <connection>
      <map_components component_2="potassium_background_current" component_1="membrane"/>
      <map_variables variable_2="FCell" variable_1="FCell"/>
      <map_variables variable_2="V" variable_1="V"/>
      <map_variables variable_2="i_b_K" variable_1="i_b_K"/>
   </connection>
   <connection>
      <map_components component_2="calcium_background_current" component_1="membrane"/>
      <map_variables variable_2="FCell" variable_1="FCell"/>
      <map_variables variable_2="V" variable_1="V"/>
      <map_variables variable_2="i_b_Ca" variable_1="i_b_Ca"/>
   </connection>
   <connection>
      <map_components component_2="sodium_calcium_exchanger" component_1="membrane"/>
      <map_variables variable_2="FCell" variable_1="FCell"/>
      <map_variables variable_2="i_NaCa" variable_1="i_NaCa"/>
      <map_variables variable_2="V" variable_1="V"/>
   </connection>
   <connection>
      <map_components component_2="sodium_potassium_pump" component_1="membrane"/>
      <map_variables variable_2="FCell" variable_1="FCell"/>
      <map_variables variable_2="i_p" variable_1="i_p"/>
      <map_variables variable_2="V" variable_1="V"/>
   </connection>
   <connection>
      <map_components component_2="persistent_calcium_current" component_1="membrane"/>
      <map_variables variable_2="FCell" variable_1="FCell"/>
      <map_variables variable_2="i_Ca_p" variable_1="i_Ca_p"/>
      <map_variables variable_2="V" variable_1="V"/>
   </connection>
   <connection>
      <map_components component_2="reversal_and_equilibrium_potentials" component_1="membrane"/>
      <map_variables variable_2="R" variable_1="R"/>
      <map_variables variable_2="T" variable_1="T"/>
      <map_variables variable_2="F" variable_1="F"/>
   </connection>
   <connection>
      <map_components component_2="sodium_current" component_1="reversal_and_equilibrium_potentials"/>
      <map_variables variable_2="E_Na" variable_1="E_Na"/>
   </connection>
   <connection>
      <map_components component_2="sodium_current" component_1="ionic_concentrations"/>
      <map_variables variable_2="Na_o" variable_1="Na_o"/>
   </connection>
   <connection>
      <map_components component_2="four_AP_sensitive_currents" component_1="reversal_and_equilibrium_potentials"/>
      <map_variables variable_2="E_K" variable_1="E_K"/>
   </connection>
   <connection>
      <map_components component_2="rapid_delayed_rectifying_potassium_current" component_1="reversal_and_equilibrium_potentials"/>
      <map_variables variable_2="E_K" variable_1="E_K"/>
   </connection>
   <connection>
      <map_components component_2="slow_delayed_rectifying_potassium_current" component_1="reversal_and_equilibrium_potentials"/>
      <map_variables variable_2="E_K_s" variable_1="E_K_s"/>
   </connection>
   <connection>
      <map_components component_2="hyperpolarisation_activated_current" component_1="reversal_and_equilibrium_potentials"/>
      <map_variables variable_2="E_K" variable_1="E_K"/>
      <map_variables variable_2="E_Na" variable_1="E_Na"/>
   </connection>
   <connection>
      <map_components component_2="sodium_background_current" component_1="reversal_and_equilibrium_potentials"/>
      <map_variables variable_2="E_Na" variable_1="E_Na"/>
   </connection>
   <connection>
      <map_components component_2="potassium_background_current" component_1="reversal_and_equilibrium_potentials"/>
      <map_variables variable_2="E_K" variable_1="E_K"/>
   </connection>
   <connection>
      <map_components component_2="calcium_background_current" component_1="reversal_and_equilibrium_potentials"/>
      <map_variables variable_2="E_Ca" variable_1="E_Ca"/>
   </connection>
   <connection>
      <map_components component_2="ionic_concentrations" component_1="sodium_calcium_exchanger"/>
      <map_variables variable_2="Ca_i" variable_1="Ca_i"/>
      <map_variables variable_2="Na_i" variable_1="Na_i"/>
      <map_variables variable_2="Ca_o" variable_1="Ca_o"/>
      <map_variables variable_2="Na_o" variable_1="Na_o"/>
   </connection>
   <connection>
      <map_components component_2="ionic_concentrations" component_1="sodium_potassium_pump"/>
      <map_variables variable_2="Na_i" variable_1="Na_i"/>
      <map_variables variable_2="K_o" variable_1="K_o"/>
   </connection>
   <connection>
      <map_components component_2="ionic_concentrations" component_1="reversal_and_equilibrium_potentials"/>
      <map_variables variable_2="K_o" variable_1="K_o"/>
      <map_variables variable_2="Na_o" variable_1="Na_o"/>
      <map_variables variable_2="K_i" variable_1="K_i"/>
      <map_variables variable_2="Na_i" variable_1="Na_i"/>
      <map_variables variable_2="Ca_o" variable_1="Ca_o"/>
      <map_variables variable_2="Ca_i" variable_1="Ca_i"/>
   </connection>
   <connection>
      <map_components component_2="ionic_concentrations" component_1="persistent_calcium_current"/>
      <map_variables variable_2="Ca_i" variable_1="Ca_i"/>
   </connection>
   <connection>
      <map_components component_2="sodium_current_m_gate" component_1="sodium_current"/>
      <map_variables variable_2="m" variable_1="m"/>
      <map_variables variable_2="time" variable_1="time"/>
      <map_variables variable_2="V" variable_1="V"/>
   </connection>
   <connection>
      <map_components component_2="sodium_current_h_gate" component_1="sodium_current"/>
      <map_variables variable_2="h" variable_1="h"/>
      <map_variables variable_2="time" variable_1="time"/>
      <map_variables variable_2="V" variable_1="V"/>
   </connection>
   <connection>
      <map_components component_2="L_type_Ca_channel_d_gate" component_1="L_type_Ca_channel"/>
      <map_variables variable_2="FCell" variable_1="FCell"/>
      <map_variables variable_2="d_L" variable_1="d_L"/>
      <map_variables variable_2="time" variable_1="time"/>
      <map_variables variable_2="V" variable_1="V"/>
   </connection>
   <connection>
      <map_components component_2="L_type_Ca_channel_f_gate" component_1="L_type_Ca_channel"/>
      <map_variables variable_2="FCell" variable_1="FCell"/>
      <map_variables variable_2="f_L" variable_1="f_L"/>
      <map_variables variable_2="time" variable_1="time"/>
      <map_variables variable_2="V" variable_1="V"/>
   </connection>
   <connection>
      <map_components component_2="T_type_Ca_channel_d_gate" component_1="T_type_Ca_channel"/>
      <map_variables variable_2="d_T" variable_1="d_T"/>
      <map_variables variable_2="time" variable_1="time"/>
      <map_variables variable_2="V" variable_1="V"/>
   </connection>
   <connection>
      <map_components component_2="T_type_Ca_channel_f_gate" component_1="T_type_Ca_channel"/>
      <map_variables variable_2="FCell" variable_1="FCell"/>
      <map_variables variable_2="f_T" variable_1="f_T"/>
      <map_variables variable_2="time" variable_1="time"/>
      <map_variables variable_2="V" variable_1="V"/>
   </connection>
   <connection>
      <map_components component_2="four_AP_sensitive_currents_q_gate" component_1="four_AP_sensitive_currents"/>
      <map_variables variable_2="FCell" variable_1="FCell"/>
      <map_variables variable_2="q" variable_1="q"/>
      <map_variables variable_2="time" variable_1="time"/>
      <map_variables variable_2="V" variable_1="V"/>
   </connection>
   <connection>
      <map_components component_2="four_AP_sensitive_currents_r_gate" component_1="four_AP_sensitive_currents"/>
      <map_variables variable_2="r" variable_1="r"/>
      <map_variables variable_2="time" variable_1="time"/>
      <map_variables variable_2="V" variable_1="V"/>
   </connection>
   <connection>
      <map_components component_2="rapid_delayed_rectifying_potassium_current_P_af_gate" component_1="rapid_delayed_rectifying_potassium_current"/>
      <map_variables variable_2="P_af" variable_1="P_af"/>
      <map_variables variable_2="time" variable_1="time"/>
      <map_variables variable_2="V" variable_1="V"/>
   </connection>
   <connection>
      <map_components component_2="rapid_delayed_rectifying_potassium_current_P_as_gate" component_1="rapid_delayed_rectifying_potassium_current"/>
      <map_variables variable_2="P_as" variable_1="P_as"/>
      <map_variables variable_2="time" variable_1="time"/>
      <map_variables variable_2="V" variable_1="V"/>
   </connection>
   <connection>
      <map_components component_2="rapid_delayed_rectifying_potassium_current_P_as_gate" component_1="rapid_delayed_rectifying_potassium_current_P_af_gate"/>
      <map_variables variable_2="P_af_infinity" variable_1="P_af_infinity"/>
   </connection>
   <connection>
      <map_components component_2="rapid_delayed_rectifying_potassium_current_P_i_gate" component_1="rapid_delayed_rectifying_potassium_current"/>
      <map_variables variable_2="P_i" variable_1="P_i"/>
      <map_variables variable_2="time" variable_1="time"/>
      <map_variables variable_2="V" variable_1="V"/>
   </connection>
   <connection>
      <map_components component_2="slow_delayed_rectifying_potassium_current_xs_gate" component_1="slow_delayed_rectifying_potassium_current"/>
      <map_variables variable_2="xs" variable_1="xs"/>
      <map_variables variable_2="time" variable_1="time"/>
      <map_variables variable_2="V" variable_1="V"/>
   </connection>
   <connection>
      <map_components component_2="hyperpolarisation_activated_current_y_gate" component_1="hyperpolarisation_activated_current"/>
      <map_variables variable_2="y" variable_1="y"/>
      <map_variables variable_2="time" variable_1="time"/>
      <map_variables variable_2="V" variable_1="V"/>
   </connection>
<rdf:RDF>
  <rdf:Seq rdf:about="rdf:#6da5cf72-8f95-4078-a061-48a14818c86a">
    <rdf:li rdf:resource="rdf:#654cc1fc-ea78-4f12-95d4-abff62ebf81d"/>
    <rdf:li rdf:resource="rdf:#c06cea6d-5ab0-4bd9-960f-2f901e37eee8"/>
    <rdf:li rdf:resource="rdf:#74266e92-0c6d-4f45-bf1d-6cbce647c78b"/>
    <rdf:li rdf:resource="rdf:#718dcd20-947d-480e-8652-55aab920ce89"/>
    <rdf:li rdf:resource="rdf:#4ca79487-b0d3-4583-bfdc-cba930a8dbf6"/>
    <rdf:li rdf:resource="rdf:#a01af73c-1cf8-43ad-9f7b-c47509eecb9c"/>
    <rdf:li rdf:resource="rdf:#b1243c8e-4b18-4a29-81ed-dd8a329f5157"/>
  </rdf:Seq>
  <rdf:Description rdf:about="rdf:#232702d6-f5b1-41be-b744-7fa585f2c94f">
    <vCard:Given>Mark</vCard:Given>
    <vCard:Family>Boyett</vCard:Family>
  </rdf:Description>
  <rdf:Description rdf:about="">
    <dc:publisher/>
    <cmeta:comment rdf:resource="rdf:#f15c907b-e3bb-4504-bf36-4e03787a744f"/>
    <dcterms:created rdf:resource="rdf:#df68909f-e0bc-4906-999f-92067e950e53"/>
    <dc:creator rdf:resource="rdf:#76fdeb6b-b9c9-46a4-b4de-59ce53347ce2"/>
    <cmeta:modification rdf:resource="rdf:#75ff48d2-afc3-4258-8302-4a66dd4280ea"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#f15c907b-e3bb-4504-bf36-4e03787a744f">
    <dc:creator rdf:resource="rdf:#81d88dd4-0667-4eb9-bedb-6d373bf4c706"/>
    <rdf:value>The CellML version of this model has been curated and unit checked by Penny Noble and it is known to run in COR to recreate the published results. The Zhang model has since been updated by Alan Garny et al. (2003) and this model has also been described in CellML and can be found in the model repository.</rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#3e398f48-8f32-41e2-bf98-716a41352180">
    <dc:title>Am J Physiol Heart Circ Physiol. </dc:title>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#ef8aecf0-0565-4668-8ede-6020f2d8fb35">
    <vCard:Given>Alan</vCard:Given>
    <vCard:Family>Garny</vCard:Family>
    <vCard:Other/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#c06cea6d-5ab0-4bd9-960f-2f901e37eee8">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#b622959d-663f-432b-ad2a-cb65fe367b8b"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#6a4a28f6-b611-49e5-8c18-ea59837f0695">
    <vCard:Given>Tony</vCard:Given>
    <vCard:Family>Varghese</vCard:Family>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#b622959d-663f-432b-ad2a-cb65fe367b8b">
    <vCard:Given>Arun</vCard:Given>
    <vCard:Family>Holden</vCard:Family>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#55c4986e-440f-4ab3-9bc9-4960ff2d5122">
    <vCard:Orgname>University of Oxford</vCard:Orgname>
    <vCard:Orgunit/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#ec90b28c-2d26-4cce-a02d-08147e835902">
    <dcterms:W3CDTF>2008-04-29T23:06:39+12:00</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#76fdeb6b-b9c9-46a4-b4de-59ce53347ce2">
    <vCard:ORG rdf:resource="rdf:#55c4986e-440f-4ab3-9bc9-4960ff2d5122"/>
    <vCard:EMAIL rdf:resource="rdf:#1a1f23a0-c270-44c1-a268-28443201b042"/>
    <vCard:N rdf:resource="rdf:#ef8aecf0-0565-4668-8ede-6020f2d8fb35"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#81d88dd4-0667-4eb9-bedb-6d373bf4c706">
    <vCard:FN>Penny Noble</vCard:FN>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#73a734bc-f981-451f-aad0-efe51dd6727c">
    <vCard:Given>H</vCard:Given>
    <vCard:Family>Honjo</vCard:Family>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#4ca79487-b0d3-4583-bfdc-cba930a8dbf6">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#e9892937-e44b-4d4c-826a-688ae1646cbd"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#1a1f23a0-c270-44c1-a268-28443201b042">
    <rdf:type rdf:resource="http://imc.org/vCard/3.0#internet"/>
    <rdf:value>alan.garny@dpag.ox.ac.uk</rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#b1243c8e-4b18-4a29-81ed-dd8a329f5157">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#232702d6-f5b1-41be-b744-7fa585f2c94f"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#df68909f-e0bc-4906-999f-92067e950e53">
    <dcterms:W3CDTF>2008-04-29T00:00:00+00:00</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#0d4c8e2f-7ecb-457b-a679-cf707b1a3492">
    <dc:creator rdf:resource="rdf:#6849be5c-4bd3-40c1-a16a-24d069ddb66e"/>
    <rdf:value/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#15502543-9b81-41c6-a3b7-a471e07ca46c">
    <vCard:N rdf:resource="rdf:#41910c2d-acd4-4d7d-9886-cdb62e91ef02"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#11d4915c-47b8-4605-9684-577d06aed09a">
    <bqs:Pubmed_id>10899081</bqs:Pubmed_id>
    <bqs:JournalArticle rdf:resource="rdf:#a3c37f74-bdb0-4505-9f3e-e9603ff98970"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#718dcd20-947d-480e-8652-55aab920ce89">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#73a734bc-f981-451f-aad0-efe51dd6727c"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#74266e92-0c6d-4f45-bf1d-6cbce647c78b">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#72e89503-8a74-4cc9-87b2-851dbb8dc585"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#654cc1fc-ea78-4f12-95d4-abff62ebf81d">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#7500ad3c-5d14-4d62-bd07-05ac80fc36b2"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#6849be5c-4bd3-40c1-a16a-24d069ddb66e">
    <vCard:FN/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#65c27604-2f1c-4066-9000-14520d67b4fb">
    <dcterms:W3CDTF>2000-07-00 00:00</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#41910c2d-acd4-4d7d-9886-cdb62e91ef02">
    <vCard:Given>Alan</vCard:Given>
    <vCard:Family>Garny</vCard:Family>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#7500ad3c-5d14-4d62-bd07-05ac80fc36b2">
    <vCard:Given>Henngui</vCard:Given>
    <vCard:Family>Zhang</vCard:Family>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#e9892937-e44b-4d4c-826a-688ae1646cbd">
    <vCard:Given>M</vCard:Given>
    <vCard:Family>Lei</vCard:Family>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#a3c37f74-bdb0-4505-9f3e-e9603ff98970">
    <dc:creator rdf:resource="rdf:#6da5cf72-8f95-4078-a061-48a14818c86a"/>
    <dc:title>Mathematical models of action potentials in the periphery and center of the rabbit sinoatrial node.</dc:title>
    <bqs:volume>279</bqs:volume>
    <bqs:first_page>397</bqs:first_page>
    <bqs:Journal rdf:resource="rdf:#3e398f48-8f32-41e2-bf98-716a41352180"/>
    <dcterms:issued rdf:resource="rdf:#65c27604-2f1c-4066-9000-14520d67b4fb"/>
    <bqs:last_page>421</bqs:last_page>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#75ff48d2-afc3-4258-8302-4a66dd4280ea">
    <dcterms:modified rdf:resource="rdf:#ec90b28c-2d26-4cce-a02d-08147e835902"/>
    <rdf:value>Worked with the original model authors to curate the CellML code and produce the published results.</rdf:value>
    <cmeta:modifier rdf:resource="rdf:#15502543-9b81-41c6-a3b7-a471e07ca46c"/>
  </rdf:Description>
  <rdf:Description rdf:about="#zhang_2000">
    <bqs:reference rdf:resource="rdf:#11d4915c-47b8-4605-9684-577d06aed09a"/>
<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>Sinoatrial Node</rdf:li>
      </rdf:Bag>
    </rdf:value>
  </dc:subject>
</bqs:reference>
    <cmeta:comment rdf:resource="rdf:#0d4c8e2f-7ecb-457b-a679-cf707b1a3492"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#72e89503-8a74-4cc9-87b2-851dbb8dc585">
    <vCard:Given>I</vCard:Given>
    <vCard:Family>Kodama</vCard:Family>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#a01af73c-1cf8-43ad-9f7b-c47509eecb9c">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#6a4a28f6-b611-49e5-8c18-ea59837f0695"/>
  </rdf:Description>
</rdf:RDF>
</model>