Location: Noble, Difrancesco, Denyer, 1989 @ f571c7b99aec / noble_difrancesco_denyer_1989.cellml

Author:
Catherine Lloyd <c.lloyd@auckland.ac.nz>
Date:
2009-11-09 15:13:51+13:00
Desc:
Added the curated model and the separate list of corrections to the workspace.
Permanent Source URI:
https://models.cellml.org/workspace/noble_difrancesco_denyer_1989/rawfile/f571c7b99aec652ffc8b748c4e1d3d7cdfaaf0b4/noble_difrancesco_denyer_1989.cellml

<?xml version="1.0" encoding="utf-8"?>
<!--
This CellML file was generated on 19/03/2009 at 12:08:14 using:

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

CellML 1.0 was used to generate this model
http://www.CellML.org/
-->
<model xmlns="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#" name="noble_difrancesco_denyer_1989" cmeta:id="noble_difrancesco_denyer_1989">
   
    <documentation xmlns="http://cellml.org/tmp-documentation">
<article>

  <section id="sec_status">
    <title>Model Status</title>
    <para>
This model has been curated by Ivan Cenci from Bologna University.  It runs in PCEnv and COR to recreate the published results.  The units have been checked and they are consistent.  There is a separate <ulink url="model_curation.doc">curation document</ulink> which lists the changes made to the model. 
          </para>
  </section>
  <sect1 id="sec_structure">
<title>Model Structure</title>

<para>
<informalfigure float="0" id="fig_reaction_diagram">
<mediaobject>
  <imageobject>
    <objectinfo>
      <title>model diagram</title>
    </objectinfo>
    <imagedata fileref="noble_1989.png"/>
  </imageobject>
</mediaobject>
<caption>Schematic diagram of the Noble et al. 1989 model.</caption>
</informalfigure>
</para>

<para>
This model was based on the article: "Ionic Mechanisms in Normal and Abnormal Cardiac Pacemaker Activity" which can be found in the book "Neuronal and Cellular Oscillators", pp. 59-85, edited by Jon W. Jacklet as part of the "Cellular Clocks" series, published by Marcel Dekker Inc.
</para> 

</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="millivolt_second">
      <unit units="millivolt"/>
      <unit units="second"/>
   </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="mA_nA">
      <unit units="ampere" prefix="milli"/>
      <unit units="nanoA" exponent="-1"/>
   </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="nanoA_per_millimolar">
      <unit units="nanoA"/>
      <unit units="millimolar" exponent="-1"/>
   </units>
   <units name="joule_per_kilomole_kelvin">
      <unit units="joule"/>
      <unit units="mole" prefix="kilo" exponent="-1"/>
      <unit units="kelvin" exponent="-1"/>
   </units>
   <units name="coulomb_per_mole">
      <unit units="coulomb"/>
      <unit units="mole" exponent="-1"/>
   </units>
   <units name="millimetre">
      <unit units="metre" prefix="milli"/>
   </units>
   <units name="millimetre3">
      <unit units="metre" prefix="milli" exponent="3"/>
   </units>
   <units name="litre_millimetre3">
      <unit units="litre"/>
      <unit units="millimetre3" exponent="-1"/>
   </units>
   <component name="environment">
      <variable name="time" units="second" public_interface="out"/>
   </component>
   <component name="membrane">
      <variable name="V" units="millivolt" initial_value="-50" public_interface="out"/>
      <variable name="R" units="joule_per_kilomole_kelvin" initial_value="8314.472" public_interface="out"/>
      <variable name="T" units="kelvin" initial_value="310" public_interface="out"/>
      <variable name="F" units="coulomb_per_mole" initial_value="96485.3415" public_interface="out"/>
      <variable name="C" units="microF" initial_value="0.000027"/>
      <variable name="RTONF" units="millivolt" public_interface="out"/>
      <variable name="time" units="second" public_interface="in"/>
      <variable name="i_f" units="nanoA" public_interface="in"/>
      <variable name="i_K" units="nanoA" public_interface="in"/>
      <variable name="i_K1" units="nanoA" public_interface="in"/>
      <variable name="i_Na_b" units="nanoA" public_interface="in"/>
      <variable name="i_Ca_b" units="nanoA" public_interface="in"/>
      <variable name="i_p" units="nanoA" public_interface="in"/>
      <variable name="i_NaCa" units="nanoA" public_interface="in"/>
      <variable name="i_Na" units="nanoA" public_interface="in"/>
      <variable name="i_Caf" units="nanoA" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>RTONF</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <ci>R</ci>
                  <ci>T</ci>
               </apply>
               <ci>F</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>V</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <apply>
                     <plus/>
                     <ci>i_f</ci>
                     <ci>i_K</ci>
                     <ci>i_K1</ci>
                     <ci>i_Na_b</ci>
                     <ci>i_Ca_b</ci>
                     <ci>i_p</ci>
                     <ci>i_NaCa</ci>
                     <ci>i_Na</ci>
                     <ci>i_Caf</ci>
                  </apply>
               </apply>
               <ci>C</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="hyperpolarising_activated_current">
      <variable name="i_f" units="nanoA" public_interface="out"/>
      <variable name="i_fNa" units="nanoA" public_interface="out"/>
      <variable name="E_Na" units="millivolt" public_interface="out"/>
      <variable name="E_K" units="millivolt" public_interface="out"/>
      <variable name="i_fK" units="nanoA" public_interface="out"/>
      <variable name="g_f_Na" units="microS" initial_value="0.06"/>
      <variable name="g_f_K" units="microS" initial_value="0.06"/>
      <variable name="Km_f" units="millimolar" initial_value="45"/>
      <variable name="time" units="second" public_interface="in" private_interface="out"/>
      <variable name="V" units="millivolt" public_interface="in" private_interface="out"/>
      <variable name="RTONF" units="millivolt" public_interface="in"/>
      <variable name="Kc" units="millimolar" public_interface="in"/>
      <variable name="Ki" units="millimolar" public_interface="in"/>
      <variable name="Nai" units="millimolar" public_interface="in"/>
      <variable name="Nao" units="millimolar" public_interface="in"/>
      <variable name="y" units="dimensionless" private_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>E_Na</ci>
            <apply>
               <times/>
               <ci>RTONF</ci>
               <apply>
                  <ln/>
                  <apply>
                     <divide/>
                     <ci>Nao</ci>
                     <ci>Nai</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>E_K</ci>
            <apply>
               <times/>
               <ci>RTONF</ci>
               <apply>
                  <ln/>
                  <apply>
                     <divide/>
                     <ci>Kc</ci>
                     <ci>Ki</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>i_fNa</ci>
            <apply>
               <times/>
               <apply>
                  <divide/>
                  <apply>
                     <times/>
                     <apply>
                        <power/>
                        <ci>y</ci>
                        <cn cellml:units="dimensionless">2</cn>
                     </apply>
                     <ci>Kc</ci>
                  </apply>
                  <apply>
                     <plus/>
                     <ci>Kc</ci>
                     <ci>Km_f</ci>
                  </apply>
               </apply>
               <ci>g_f_Na</ci>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>E_Na</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>i_fK</ci>
            <apply>
               <times/>
               <apply>
                  <divide/>
                  <apply>
                     <times/>
                     <apply>
                        <power/>
                        <ci>y</ci>
                        <cn cellml:units="dimensionless">2</cn>
                     </apply>
                     <ci>Kc</ci>
                  </apply>
                  <apply>
                     <plus/>
                     <ci>Kc</ci>
                     <ci>Km_f</ci>
                  </apply>
               </apply>
               <ci>g_f_K</ci>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>E_K</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>i_f</ci>
            <apply>
               <plus/>
               <ci>i_fNa</ci>
               <ci>i_fK</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="hyperpolarising_activated_current_y_gate">
      <variable name="y" units="dimensionless" initial_value="0.0822" public_interface="out"/>
      <variable name="alpha_y" units="per_second"/>
      <variable name="beta_y" units="per_second"/>
      <variable name="V" units="millivolt" public_interface="in"/>
      <variable name="time" units="second" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>alpha_y</ci>
            <apply>
               <times/>
               <cn cellml:units="per_second">0.02</cn>
               <apply>
                  <exp/>
                  <apply>
                     <divide/>
                     <apply>
                        <minus/>
                        <ci>V</ci>
                     </apply>
                     <cn cellml:units="millivolt">18</cn>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>beta_y</ci>
            <apply>
               <times/>
               <cn cellml:units="per_second">19.5</cn>
               <apply>
                  <exp/>
                  <apply>
                     <divide/>
                     <ci>V</ci>
                     <cn cellml:units="millivolt">19</cn>
                  </apply>
               </apply>
            </apply>
         </apply>
         <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>
      </math>
   </component>
   <component name="time_dependent_potassium_current">
      <variable name="i_K" units="nanoA" public_interface="out"/>
      <variable name="I_K" units="nanoA"/>
      <variable name="i_K_max" units="nanoA" initial_value="0.8"/>
      <variable name="time" units="second" public_interface="in" private_interface="out"/>
      <variable name="V" units="millivolt" public_interface="in" private_interface="out"/>
      <variable name="RTONF" units="millivolt" public_interface="in"/>
      <variable name="Ki" units="millimolar" public_interface="in"/>
      <variable name="Kc" units="millimolar" public_interface="in"/>
      <variable name="x" units="dimensionless" private_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>I_K</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <ci>i_K_max</ci>
                  <apply>
                     <minus/>
                     <ci>Ki</ci>
                     <apply>
                        <times/>
                        <ci>Kc</ci>
                        <apply>
                           <exp/>
                           <apply>
                              <divide/>
                              <apply>
                                 <minus/>
                                 <ci>V</ci>
                              </apply>
                              <cn cellml:units="millivolt">25</cn>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </apply>
               <cn cellml:units="millimolar">140</cn>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>i_K</ci>
            <apply>
               <times/>
               <ci>x</ci>
               <ci>I_K</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="time_dependent_potassium_current_x_gate">
      <variable name="x" units="dimensionless" initial_value="0.225" public_interface="out"/>
      <variable name="alpha_x" units="per_second"/>
      <variable name="beta_x" units="per_second"/>
      <variable name="time" units="second" public_interface="in"/>
      <variable name="V" units="millivolt" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>alpha_x</ci>
            <apply>
               <times/>
               <cn cellml:units="per_second">2.1</cn>
               <apply>
                  <exp/>
                  <apply>
                     <divide/>
                     <ci>V</ci>
                     <cn cellml:units="millivolt">28</cn>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>beta_x</ci>
            <apply>
               <times/>
               <cn cellml:units="per_second">0.96</cn>
               <apply>
                  <exp/>
                  <apply>
                     <divide/>
                     <apply>
                        <minus/>
                        <ci>V</ci>
                     </apply>
                     <cn cellml:units="millivolt">24</cn>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>x</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>alpha_x</ci>
                  <apply>
                     <minus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <ci>x</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>beta_x</ci>
                  <ci>x</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="time_independent_potassium_current">
      <variable name="i_K1" units="nanoA" public_interface="out"/>
      <variable name="g_K1" units="microS" initial_value="0.0075"/>
      <variable name="Km_K1" units="millimolar" initial_value="10"/>
      <variable name="time" units="second" public_interface="in"/>
      <variable name="V" units="millivolt" public_interface="in"/>
      <variable name="E_K" units="millivolt" public_interface="in"/>
      <variable name="Kc" units="millimolar" public_interface="in"/>
      <variable name="RTONF" units="millivolt" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_K1</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <ci>g_K1</ci>
                        <ci>Kc</ci>
                     </apply>
                     <apply>
                        <plus/>
                        <ci>Kc</ci>
                        <ci>Km_K1</ci>
                     </apply>
                  </apply>
                  <apply>
                     <minus/>
                     <ci>V</ci>
                     <ci>E_K</ci>
                  </apply>
               </apply>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <times/>
                           <apply>
                              <minus/>
                              <apply>
                                 <plus/>
                                 <ci>V</ci>
                                 <cn cellml:units="millivolt">10</cn>
                              </apply>
                              <ci>E_K</ci>
                           </apply>
                           <cn cellml:units="dimensionless">2</cn>
                        </apply>
                        <ci>RTONF</ci>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="sodium_background_current">
      <variable name="i_Na_b" units="nanoA" public_interface="out"/>
      <variable name="g_Nab" units="microS" initial_value="0.0007"/>
      <variable name="time" units="second" public_interface="in"/>
      <variable name="V" units="millivolt" public_interface="in"/>
      <variable name="E_Na" units="millivolt" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_Na_b</ci>
            <apply>
               <times/>
               <ci>g_Nab</ci>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>E_Na</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="calcium_background_current">
      <variable name="i_Ca_b" units="nanoA" public_interface="out"/>
      <variable name="E_Ca" units="millivolt" public_interface="out"/>
      <variable name="g_Cab" units="microS" initial_value="0.0001"/>
      <variable name="time" units="second" public_interface="in"/>
      <variable name="V" units="millivolt" public_interface="in"/>
      <variable name="RTONF" units="millivolt" public_interface="in"/>
      <variable name="Cai" units="millimolar" public_interface="in"/>
      <variable name="Cao" units="millimolar" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>E_Ca</ci>
            <apply>
               <times/>
               <cn cellml:units="dimensionless">0.5</cn>
               <ci>RTONF</ci>
               <apply>
                  <ln/>
                  <apply>
                     <divide/>
                     <ci>Cao</ci>
                     <ci>Cai</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>i_Ca_b</ci>
            <apply>
               <times/>
               <ci>g_Cab</ci>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>E_Ca</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="sodium_potassium_pump">
      <variable name="i_p" units="nanoA" public_interface="out"/>
      <variable name="I_p" units="nanoA" initial_value="0.5"/>
      <variable name="K_mK" units="millimolar" initial_value="1"/>
      <variable name="K_mNa" units="millimolar" initial_value="40"/>
      <variable name="time" units="second" public_interface="in"/>
      <variable name="V" units="millivolt" public_interface="in"/>
      <variable name="Nai" units="millimolar" public_interface="in"/>
      <variable name="Kc" units="millimolar" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_p</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <ci>I_p</ci>
                        <ci>Kc</ci>
                     </apply>
                     <apply>
                        <plus/>
                        <ci>K_mK</ci>
                        <ci>Kc</ci>
                     </apply>
                  </apply>
                  <ci>Nai</ci>
               </apply>
               <apply>
                  <plus/>
                  <ci>K_mNa</ci>
                  <ci>Nai</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="Na_Ca_exchanger">
      <variable name="i_NaCa" units="nanoA" public_interface="out"/>
      <variable name="n_NaCa" units="dimensionless" initial_value="3" public_interface="out"/>
      <variable name="K_NaCa" units="nanoA" initial_value="0.00002"/>
      <variable name="d_NaCa" units="dimensionless" initial_value="0.0001"/>
      <variable name="gamma" units="dimensionless" initial_value="0.5"/>
      <variable name="time" units="second" public_interface="in"/>
      <variable name="V" units="millivolt" public_interface="in"/>
      <variable name="RTONF" units="millivolt" public_interface="in"/>
      <variable name="Cao" units="millimolar" public_interface="in"/>
      <variable name="Cai" units="millimolar" public_interface="in"/>
      <variable name="Nai" units="millimolar" public_interface="in"/>
      <variable name="Nao" units="millimolar" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_NaCa</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <ci>K_NaCa</ci>
                  <apply>
                     <minus/>
                     <apply>
                        <times/>
                        <apply>
                           <exp/>
                           <apply>
                              <divide/>
                              <apply>
                                 <times/>
                                 <ci>gamma</ci>
                                 <apply>
                                    <minus/>
                                    <ci>n_NaCa</ci>
                                    <cn cellml:units="dimensionless">2</cn>
                                 </apply>
                                 <ci>V</ci>
                              </apply>
                              <ci>RTONF</ci>
                           </apply>
                        </apply>
                        <apply>
                           <power/>
                           <ci>Nai</ci>
                           <ci>n_NaCa</ci>
                        </apply>
                        <ci>Cao</ci>
                     </apply>
                     <apply>
                        <times/>
                        <apply>
                           <exp/>
                           <apply>
                              <divide/>
                              <apply>
                                 <times/>
                                 <apply>
                                    <minus/>
                                    <ci>gamma</ci>
                                    <cn cellml:units="dimensionless">1</cn>
                                 </apply>
                                 <apply>
                                    <minus/>
                                    <ci>n_NaCa</ci>
                                    <cn cellml:units="dimensionless">2</cn>
                                 </apply>
                                 <ci>V</ci>
                              </apply>
                              <ci>RTONF</ci>
                           </apply>
                        </apply>
                        <apply>
                           <power/>
                           <ci>Nao</ci>
                           <ci>n_NaCa</ci>
                        </apply>
                        <ci>Cai</ci>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <apply>
                     <plus/>
                     <cn cellml:units="millimolar4">1</cn>
                     <apply>
                        <times/>
                        <ci>d_NaCa</ci>
                        <apply>
                           <plus/>
                           <apply>
                              <times/>
                              <ci>Cai</ci>
                              <apply>
                                 <power/>
                                 <ci>Nao</ci>
                                 <ci>n_NaCa</ci>
                              </apply>
                           </apply>
                           <apply>
                              <times/>
                              <ci>Cao</ci>
                              <apply>
                                 <power/>
                                 <ci>Nai</ci>
                                 <ci>n_NaCa</ci>
                              </apply>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
                  <apply>
                     <plus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <apply>
                        <divide/>
                        <ci>Cai</ci>
                        <cn cellml:units="millimolar">0.0069</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="fast_sodium_current">
      <variable name="i_Na" units="nanoA" public_interface="out"/>
      <variable name="g_Na" units="microS" initial_value="0.0125"/>
      <variable name="E_mh" units="millivolt"/>
      <variable name="time" units="second" public_interface="in" private_interface="out"/>
      <variable name="V" units="millivolt" public_interface="in" private_interface="out"/>
      <variable name="RTONF" units="millivolt" public_interface="in"/>
      <variable name="Nao" units="millimolar" public_interface="in"/>
      <variable name="Nai" units="millimolar" public_interface="in"/>
      <variable name="Kc" units="millimolar" public_interface="in"/>
      <variable name="Ki" units="millimolar" public_interface="in"/>
      <variable name="m" units="dimensionless" private_interface="in"/>
      <variable name="h" units="dimensionless" private_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>E_mh</ci>
            <apply>
               <times/>
               <ci>RTONF</ci>
               <apply>
                  <ln/>
                  <apply>
                     <divide/>
                     <apply>
                        <plus/>
                        <ci>Nao</ci>
                        <apply>
                           <times/>
                           <cn cellml:units="dimensionless">0.12</cn>
                           <ci>Kc</ci>
                        </apply>
                     </apply>
                     <apply>
                        <plus/>
                        <ci>Nai</ci>
                        <apply>
                           <times/>
                           <cn cellml:units="dimensionless">0.12</cn>
                           <ci>Ki</ci>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>i_Na</ci>
            <apply>
               <times/>
               <ci>g_Na</ci>
               <apply>
                  <power/>
                  <ci>m</ci>
                  <cn cellml:units="dimensionless">3</cn>
               </apply>
               <ci>h</ci>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>E_mh</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="fast_sodium_current_m_gate">
      <variable name="m" units="dimensionless" initial_value="0.0365" public_interface="out"/>
      <variable name="alpha_m" units="per_second"/>
      <variable name="beta_m" units="per_second"/>
      <variable name="V" units="millivolt" public_interface="in"/>
      <variable name="time" units="second" public_interface="in"/>
      <variable name="delta_m" units="millivolt" initial_value="1e-5"/>
      <variable name="E0_m" units="millivolt"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>E0_m</ci>
            <apply>
               <plus/>
               <ci>V</ci>
               <cn cellml:units="millivolt">41</cn>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>alpha_m</ci>
            <piecewise>
               <piece>
                  <cn cellml:units="per_second">2000</cn>
                  <apply>
                     <lt/>
                     <apply>
                        <abs/>
                        <ci>E0_m</ci>
                     </apply>
                     <ci>delta_m</ci>
                  </apply>
               </piece>
               <otherwise>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <cn cellml:units="per_millivolt_second">200</cn>
                        <ci>E0_m</ci>
                     </apply>
                     <apply>
                        <minus/>
                        <cn cellml:units="dimensionless">1</cn>
                        <apply>
                           <exp/>
                           <apply>
                              <times/>
                              <apply>
                                 <minus/>
                                 <cn cellml:units="per_millivolt">0.1</cn>
                              </apply>
                              <ci>E0_m</ci>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </otherwise>
            </piecewise>
         </apply>
         <apply>
            <eq/>
            <ci>beta_m</ci>
            <apply>
               <times/>
               <cn cellml:units="per_second">8000</cn>
               <apply>
                  <exp/>
                  <apply>
                     <times/>
                     <apply>
                        <minus/>
                        <cn cellml:units="per_millivolt">0.056</cn>
                     </apply>
                     <apply>
                        <plus/>
                        <ci>V</ci>
                        <cn cellml:units="millivolt">66</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>m</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>alpha_m</ci>
                  <apply>
                     <minus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <ci>m</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>beta_m</ci>
                  <ci>m</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="fast_sodium_current_h_gate">
      <variable name="h" units="dimensionless" initial_value="0.1969" public_interface="out"/>
      <variable name="alpha_h" units="per_second"/>
      <variable name="beta_h" units="per_second"/>
      <variable name="V" units="millivolt" public_interface="in"/>
      <variable name="time" units="second" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>alpha_h</ci>
            <apply>
               <times/>
               <cn cellml:units="per_second">20</cn>
               <apply>
                  <exp/>
                  <apply>
                     <times/>
                     <apply>
                        <minus/>
                        <cn cellml:units="per_millivolt">0.125</cn>
                     </apply>
                     <apply>
                        <plus/>
                        <ci>V</ci>
                        <cn cellml:units="millivolt">75</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>beta_h</ci>
            <apply>
               <divide/>
               <cn cellml:units="per_second">2000</cn>
               <apply>
                  <plus/>
                  <apply>
                     <times/>
                     <cn cellml:units="dimensionless">320</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <times/>
                           <apply>
                              <minus/>
                              <cn cellml:units="per_millivolt">0.1</cn>
                           </apply>
                           <apply>
                              <plus/>
                              <ci>V</ci>
                              <cn cellml:units="millivolt">75</cn>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
                  <cn cellml:units="dimensionless">1</cn>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>h</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>alpha_h</ci>
                  <apply>
                     <minus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <ci>h</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>beta_h</ci>
                  <ci>h</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="second_inward_current">
      <variable name="i_Caf" units="nanoA" public_interface="out"/>
      <variable name="i_siCa" units="nanoA" public_interface="out"/>
      <variable name="i_siK" units="nanoA" public_interface="out"/>
      <variable name="i_siNa" units="nanoA" public_interface="out"/>
      <variable name="P_si" units="nanoA_per_millimolar" initial_value="0.12"/>
      <variable name="time" units="second" public_interface="in" private_interface="out"/>
      <variable name="V" units="millivolt" public_interface="in" private_interface="out"/>
      <variable name="RTONF" units="millivolt" public_interface="in"/>
      <variable name="Nao" units="millimolar" public_interface="in"/>
      <variable name="Nai" units="millimolar" public_interface="in"/>
      <variable name="Kc" units="millimolar" public_interface="in"/>
      <variable name="Ki" units="millimolar" public_interface="in"/>
      <variable name="Cao" units="millimolar" public_interface="in"/>
      <variable name="Cai" units="millimolar" public_interface="in" private_interface="out"/>
      <variable name="d" units="dimensionless" private_interface="in"/>
      <variable name="f" units="dimensionless" private_interface="in"/>
      <variable name="f2" units="dimensionless" private_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_siCa</ci>
            <apply>
               <times/>
               <apply>
                  <divide/>
                  <apply>
                     <times/>
                     <cn cellml:units="dimensionless">4</cn>
                     <ci>P_si</ci>
                     <apply>
                        <minus/>
                        <ci>V</ci>
                        <cn cellml:units="millivolt">50</cn>
                     </apply>
                  </apply>
                  <apply>
                     <times/>
                     <ci>RTONF</ci>
                     <apply>
                        <minus/>
                        <cn cellml:units="dimensionless">1</cn>
                        <apply>
                           <exp/>
                           <apply>
                              <divide/>
                              <apply>
                                 <times/>
                                 <apply>
                                    <minus/>
                                    <cn cellml:units="dimensionless">1</cn>
                                 </apply>
                                 <apply>
                                    <minus/>
                                    <ci>V</ci>
                                    <cn cellml:units="millivolt">50</cn>
                                 </apply>
                                 <cn cellml:units="dimensionless">2</cn>
                              </apply>
                              <ci>RTONF</ci>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <minus/>
                  <apply>
                     <times/>
                     <ci>Cai</ci>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <cn cellml:units="millivolt">100</cn>
                           <ci>RTONF</ci>
                        </apply>
                     </apply>
                  </apply>
                  <apply>
                     <times/>
                     <ci>Cao</ci>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <times/>
                              <apply>
                                 <minus/>
                                 <cn cellml:units="dimensionless">2</cn>
                              </apply>
                              <apply>
                                 <minus/>
                                 <ci>V</ci>
                                 <cn cellml:units="millivolt">50</cn>
                              </apply>
                           </apply>
                           <ci>RTONF</ci>
                        </apply>
                     </apply>
                  </apply>
               </apply>
               <ci>d</ci>
               <ci>f</ci>
               <ci>f2</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>i_siK</ci>
            <apply>
               <times/>
               <apply>
                  <divide/>
                  <apply>
                     <times/>
                     <cn cellml:units="dimensionless">0.01</cn>
                     <ci>P_si</ci>
                     <apply>
                        <minus/>
                        <ci>V</ci>
                        <cn cellml:units="millivolt">50</cn>
                     </apply>
                  </apply>
                  <apply>
                     <times/>
                     <ci>RTONF</ci>
                     <apply>
                        <minus/>
                        <cn cellml:units="dimensionless">1</cn>
                        <apply>
                           <exp/>
                           <apply>
                              <divide/>
                              <apply>
                                 <times/>
                                 <apply>
                                    <minus/>
                                    <cn cellml:units="dimensionless">1</cn>
                                 </apply>
                                 <apply>
                                    <minus/>
                                    <ci>V</ci>
                                    <cn cellml:units="millivolt">50</cn>
                                 </apply>
                              </apply>
                              <ci>RTONF</ci>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <minus/>
                  <apply>
                     <times/>
                     <ci>Ki</ci>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <cn cellml:units="millivolt">50</cn>
                           <ci>RTONF</ci>
                        </apply>
                     </apply>
                  </apply>
                  <apply>
                     <times/>
                     <ci>Kc</ci>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <times/>
                              <apply>
                                 <minus/>
                                 <cn cellml:units="dimensionless">1</cn>
                              </apply>
                              <apply>
                                 <minus/>
                                 <ci>V</ci>
                                 <cn cellml:units="millivolt">50</cn>
                              </apply>
                           </apply>
                           <ci>RTONF</ci>
                        </apply>
                     </apply>
                  </apply>
               </apply>
               <ci>d</ci>
               <ci>f</ci>
               <ci>f2</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>i_siNa</ci>
            <apply>
               <times/>
               <apply>
                  <divide/>
                  <apply>
                     <times/>
                     <cn cellml:units="dimensionless">0.01</cn>
                     <ci>P_si</ci>
                     <apply>
                        <minus/>
                        <ci>V</ci>
                        <cn cellml:units="millivolt">50</cn>
                     </apply>
                  </apply>
                  <apply>
                     <times/>
                     <ci>RTONF</ci>
                     <apply>
                        <minus/>
                        <cn cellml:units="dimensionless">1</cn>
                        <apply>
                           <exp/>
                           <apply>
                              <divide/>
                              <apply>
                                 <times/>
                                 <apply>
                                    <minus/>
                                    <cn cellml:units="dimensionless">1</cn>
                                 </apply>
                                 <apply>
                                    <minus/>
                                    <ci>V</ci>
                                    <cn cellml:units="millivolt">50</cn>
                                 </apply>
                              </apply>
                              <ci>RTONF</ci>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <minus/>
                  <apply>
                     <times/>
                     <ci>Nai</ci>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <cn cellml:units="millivolt">50</cn>
                           <ci>RTONF</ci>
                        </apply>
                     </apply>
                  </apply>
                  <apply>
                     <times/>
                     <ci>Nao</ci>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <times/>
                              <apply>
                                 <minus/>
                                 <cn cellml:units="dimensionless">1</cn>
                              </apply>
                              <apply>
                                 <minus/>
                                 <ci>V</ci>
                                 <cn cellml:units="millivolt">50</cn>
                              </apply>
                           </apply>
                           <ci>RTONF</ci>
                        </apply>
                     </apply>
                  </apply>
               </apply>
               <ci>d</ci>
               <ci>f</ci>
               <ci>f2</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>i_Caf</ci>
            <apply>
               <plus/>
               <ci>i_siCa</ci>
               <ci>i_siK</ci>
               <ci>i_siNa</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="second_inward_current_d_gate">
      <variable name="d" units="dimensionless" initial_value="0" public_interface="out"/>
      <variable name="alpha_d" units="per_second"/>
      <variable name="beta_d" units="per_second"/>
      <variable name="V" units="millivolt" public_interface="in"/>
      <variable name="time" units="second" public_interface="in"/>
      <variable name="delta_d" units="millivolt" initial_value="0.0001"/>
      <variable name="E0_d" units="millivolt"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>E0_d</ci>
            <apply>
               <minus/>
               <apply>
                  <plus/>
                  <ci>V</ci>
                  <cn cellml:units="millivolt">24</cn>
               </apply>
               <cn cellml:units="millivolt">5</cn>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>alpha_d</ci>
            <piecewise>
               <piece>
                  <cn cellml:units="per_second">120</cn>
                  <apply>
                     <lt/>
                     <apply>
                        <abs/>
                        <ci>E0_d</ci>
                     </apply>
                     <ci>delta_d</ci>
                  </apply>
               </piece>
               <otherwise>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <cn cellml:units="per_millivolt_second">30</cn>
                        <ci>E0_d</ci>
                     </apply>
                     <apply>
                        <minus/>
                        <cn cellml:units="dimensionless">1</cn>
                        <apply>
                           <exp/>
                           <apply>
                              <divide/>
                              <apply>
                                 <times/>
                                 <apply>
                                    <minus/>
                                    <cn cellml:units="dimensionless">1</cn>
                                 </apply>
                                 <ci>E0_d</ci>
                              </apply>
                              <cn cellml:units="millivolt">4</cn>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </otherwise>
            </piecewise>
         </apply>
         <apply>
            <eq/>
            <ci>beta_d</ci>
            <piecewise>
               <piece>
                  <cn cellml:units="per_second">120</cn>
                  <apply>
                     <lt/>
                     <apply>
                        <abs/>
                        <ci>E0_d</ci>
                     </apply>
                     <ci>delta_d</ci>
                  </apply>
               </piece>
               <otherwise>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <cn cellml:units="per_millivolt_second">12</cn>
                        <ci>E0_d</ci>
                     </apply>
                     <apply>
                        <minus/>
                        <apply>
                           <exp/>
                           <apply>
                              <divide/>
                              <ci>E0_d</ci>
                              <cn cellml:units="millivolt">10</cn>
                           </apply>
                        </apply>
                        <cn cellml:units="dimensionless">1</cn>
                     </apply>
                  </apply>
               </otherwise>
            </piecewise>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>d</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>alpha_d</ci>
                  <apply>
                     <minus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <ci>d</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>beta_d</ci>
                  <ci>d</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="second_inward_current_f_gate">
      <variable name="f" units="dimensionless" initial_value="0.9997" public_interface="out"/>
      <variable name="alpha_f" units="per_second"/>
      <variable name="beta_f" units="per_second"/>
      <variable name="V" units="millivolt" public_interface="in"/>
      <variable name="time" units="second" public_interface="in"/>
      <variable name="delta_f" units="millivolt" initial_value="0.0001"/>
      <variable name="E0_f" units="millivolt"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>E0_f</ci>
            <apply>
               <plus/>
               <ci>V</ci>
               <cn cellml:units="millivolt">34</cn>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>alpha_f</ci>
            <piecewise>
               <piece>
                  <cn cellml:units="per_second">25</cn>
                  <apply>
                     <lt/>
                     <apply>
                        <abs/>
                        <ci>E0_f</ci>
                     </apply>
                     <ci>delta_f</ci>
                  </apply>
               </piece>
               <otherwise>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <cn cellml:units="per_millivolt_second">6.25</cn>
                        <ci>E0_f</ci>
                     </apply>
                     <apply>
                        <minus/>
                        <apply>
                           <exp/>
                           <apply>
                              <divide/>
                              <ci>E0_f</ci>
                              <cn cellml:units="millivolt">4</cn>
                           </apply>
                        </apply>
                        <cn cellml:units="dimensionless">1</cn>
                     </apply>
                  </apply>
               </otherwise>
            </piecewise>
         </apply>
         <apply>
            <eq/>
            <ci>beta_f</ci>
            <apply>
               <divide/>
               <cn cellml:units="per_second">50</cn>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <times/>
                           <apply>
                              <minus/>
                              <cn cellml:units="dimensionless">1</cn>
                           </apply>
                           <apply>
                              <plus/>
                              <ci>V</ci>
                              <cn cellml:units="millivolt">34</cn>
                           </apply>
                        </apply>
                        <cn cellml:units="millivolt">4</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>f</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>alpha_f</ci>
                  <apply>
                     <minus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <ci>f</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>beta_f</ci>
                  <ci>f</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="second_inward_current_f2_gate">
      <variable name="f2" units="dimensionless" initial_value="0.5765" public_interface="out"/>
      <variable name="alpha_f2" units="per_second" initial_value="10"/>
      <variable name="beta_f2" units="per_second"/>
      <variable name="K_mf2" units="millimolar" initial_value="0.0005"/>
      <variable name="V" units="millivolt" public_interface="in"/>
      <variable name="time" units="second" public_interface="in"/>
      <variable name="Cai" units="millimolar" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>beta_f2</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <ci>Cai</ci>
                  <ci>alpha_f2</ci>
               </apply>
               <ci>K_mf2</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>f2</ci>
            </apply>
            <apply>
               <minus/>
               <ci>alpha_f2</ci>
               <apply>
                  <times/>
                  <ci>f2</ci>
                  <apply>
                     <plus/>
                     <ci>alpha_f2</ci>
                     <ci>beta_f2</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="extracellular_sodium_concentration">
      <variable name="Nao" units="millimolar" initial_value="140" public_interface="out"/>
   </component>
   <component name="intracellular_sodium_concentration">
      <variable name="Nai" units="millimolar" initial_value="7.5" public_interface="out"/>
      <variable name="radius" units="millimetre" initial_value="0.008"/>
      <variable name="length" units="millimetre" initial_value="0.11"/>
      <variable name="V_e_ratio" units="dimensionless" initial_value="0.1" public_interface="out"/>
      <variable name="V_Cell" units="millimetre3" public_interface="out"/>
      <variable name="Vi" units="millimetre3" public_interface="out"/>
      <variable name="time" units="second" public_interface="in"/>
      <variable name="F" units="coulomb_per_mole" public_interface="in"/>
      <variable name="i_Na" units="nanoA" public_interface="in"/>
      <variable name="i_Na_b" units="nanoA" public_interface="in"/>
      <variable name="i_siNa" units="nanoA" public_interface="in"/>
      <variable name="i_p" units="nanoA" public_interface="in"/>
      <variable name="i_fNa" units="nanoA" public_interface="in"/>
      <variable name="i_NaCa" units="nanoA" public_interface="in"/>
      <variable name="n_NaCa" units="dimensionless" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>V_Cell</ci>
            <apply>
               <times/>
               <cn cellml:units="dimensionless">3.141592654</cn>
               <apply>
                  <power/>
                  <ci>radius</ci>
                  <cn cellml:units="dimensionless">2</cn>
               </apply>
               <ci>length</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>Vi</ci>
            <apply>
               <times/>
               <ci>V_Cell</ci>
               <apply>
                  <minus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <ci>V_e_ratio</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>Nai</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <apply>
                     <minus/>
                     <cn cellml:units="mA_nA">1</cn>
                  </apply>
                  <apply>
                     <plus/>
                     <ci>i_Na</ci>
                     <ci>i_Na_b</ci>
                     <ci>i_fNa</ci>
                     <ci>i_siNa</ci>
                     <apply>
                        <times/>
                        <ci>i_p</ci>
                        <cn cellml:units="dimensionless">3</cn>
                     </apply>
                     <apply>
                        <divide/>
                        <apply>
                           <times/>
                           <ci>i_NaCa</ci>
                           <ci>n_NaCa</ci>
                        </apply>
                        <apply>
                           <minus/>
                           <ci>n_NaCa</ci>
                           <cn cellml:units="dimensionless">2</cn>
                        </apply>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <cn cellml:units="litre_millimetre3">1</cn>
                  <ci>Vi</ci>
                  <ci>F</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="extracellular_calcium_concentration">
      <variable name="Cao" units="millimolar" initial_value="2" public_interface="out"/>
   </component>
   <component name="intracellular_calcium_concentration">
      <variable name="Cai" units="millimolar" initial_value="5.6e-5" public_interface="out"/>
      <variable name="V_up" units="millimetre3"/>
      <variable name="V_rel" units="millimetre3"/>
      <variable name="Vi" units="millimetre3" public_interface="in"/>
      <variable name="i_up" units="nanoA"/>
      <variable name="i_tr" units="nanoA"/>
      <variable name="i_rel" units="nanoA"/>
      <variable name="Ca_up" units="millimolar" initial_value="1.98"/>
      <variable name="Ca_rel" units="millimolar" initial_value="0.55"/>
      <variable name="Ca_up_max" units="millimolar" initial_value="5"/>
      <variable name="K_mCa" units="millimolar" initial_value="0.002"/>
      <variable name="p" units="dimensionless" initial_value="0.237"/>
      <variable name="alpha_p" units="per_second"/>
      <variable name="beta_p" units="per_second"/>
      <variable name="E0_p" units="millivolt"/>
      <variable name="tau_up" units="second" initial_value="0.005"/>
      <variable name="tau_rep" units="second" initial_value="0.2"/>
      <variable name="tau_rel" units="second" initial_value="0.01"/>
      <variable name="rCa" units="dimensionless" initial_value="2"/>
      <variable name="time" units="second" public_interface="in"/>
      <variable name="V" units="millivolt" public_interface="in"/>
      <variable name="F" units="coulomb_per_mole" public_interface="in"/>
      <variable name="i_Ca_b" units="nanoA" public_interface="in"/>
      <variable name="i_siCa" units="nanoA" public_interface="in"/>
      <variable name="i_NaCa" units="nanoA" public_interface="in"/>
      <variable name="n_NaCa" units="dimensionless" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_up</ci>
            <apply>
               <times/>
               <apply>
                  <divide/>
                  <apply>
                     <times/>
                     <cn cellml:units="dimensionless">2</cn>
                     <cn cellml:units="litre_millimetre3">1</cn>
                     <ci>Vi</ci>
                     <ci>F</ci>
                  </apply>
                  <apply>
                     <times/>
                     <cn cellml:units="mA_nA">1</cn>
                     <ci>tau_up</ci>
                     <ci>Ca_up_max</ci>
                  </apply>
               </apply>
               <ci>Cai</ci>
               <apply>
                  <minus/>
                  <ci>Ca_up_max</ci>
                  <ci>Ca_up</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>i_tr</ci>
            <apply>
               <times/>
               <apply>
                  <divide/>
                  <apply>
                     <times/>
                     <cn cellml:units="dimensionless">2</cn>
                     <cn cellml:units="litre_millimetre3">1</cn>
                     <ci>V_rel</ci>
                     <ci>F</ci>
                  </apply>
                  <apply>
                     <times/>
                     <cn cellml:units="mA_nA">1</cn>
                     <ci>tau_rep</ci>
                  </apply>
               </apply>
               <ci>p</ci>
               <apply>
                  <minus/>
                  <ci>Ca_up</ci>
                  <ci>Ca_rel</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>E0_p</ci>
            <apply>
               <minus/>
               <apply>
                  <plus/>
                  <ci>V</ci>
                  <cn cellml:units="millivolt">34</cn>
               </apply>
               <apply>
                  <minus/>
                  <cn cellml:units="millivolt">30</cn>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>alpha_p</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <cn cellml:units="per_millivolt_second">0.625</cn>
                  <ci>E0_p</ci>
               </apply>
               <apply>
                  <minus/>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <ci>E0_p</ci>
                        <cn cellml:units="millivolt">4</cn>
                     </apply>
                  </apply>
                  <cn cellml:units="dimensionless">1</cn>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>beta_p</ci>
            <apply>
               <divide/>
               <cn cellml:units="per_second">5</cn>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <times/>
                           <apply>
                              <minus/>
                              <cn cellml:units="dimensionless">1</cn>
                           </apply>
                           <ci>E0_p</ci>
                        </apply>
                        <cn cellml:units="millivolt">4</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>i_rel</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <cn cellml:units="dimensionless">2</cn>
                        <cn cellml:units="litre_millimetre3">1</cn>
                        <ci>V_rel</ci>
                        <ci>F</ci>
                     </apply>
                     <apply>
                        <times/>
                        <cn cellml:units="mA_nA">1</cn>
                        <ci>tau_rel</ci>
                     </apply>
                  </apply>
                  <ci>Ca_rel</ci>
                  <apply>
                     <power/>
                     <ci>Cai</ci>
                     <ci>rCa</ci>
                  </apply>
               </apply>
               <apply>
                  <plus/>
                  <apply>
                     <power/>
                     <ci>Cai</ci>
                     <ci>rCa</ci>
                  </apply>
                  <apply>
                     <power/>
                     <ci>K_mCa</ci>
                     <ci>rCa</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>p</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>alpha_p</ci>
                  <apply>
                     <minus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <ci>p</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>beta_p</ci>
                  <ci>p</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>Ca_up</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <cn cellml:units="mA_nA">1</cn>
                  <apply>
                     <minus/>
                     <ci>i_up</ci>
                     <ci>i_tr</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <cn cellml:units="dimensionless">2</cn>
                  <cn cellml:units="litre_millimetre3">1</cn>
                  <ci>V_up</ci>
                  <ci>F</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>V_up</ci>
            <apply>
               <times/>
               <ci>Vi</ci>
               <cn cellml:units="dimensionless">0.05</cn>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>Ca_rel</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <cn cellml:units="mA_nA">1</cn>
                  <apply>
                     <minus/>
                     <ci>i_tr</ci>
                     <ci>i_rel</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <cn cellml:units="dimensionless">2</cn>
                  <cn cellml:units="litre_millimetre3">1</cn>
                  <ci>V_rel</ci>
                  <ci>F</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>V_rel</ci>
            <apply>
               <times/>
               <ci>Vi</ci>
               <cn cellml:units="dimensionless">0.02</cn>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>Cai</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <apply>
                     <minus/>
                     <cn cellml:units="mA_nA">1</cn>
                  </apply>
                  <apply>
                     <plus/>
                     <apply>
                        <minus/>
                        <apply>
                           <minus/>
                           <apply>
                              <plus/>
                              <ci>i_siCa</ci>
                              <ci>i_Ca_b</ci>
                           </apply>
                           <apply>
                              <divide/>
                              <apply>
                                 <times/>
                                 <cn cellml:units="dimensionless">2</cn>
                                 <ci>i_NaCa</ci>
                              </apply>
                              <apply>
                                 <minus/>
                                 <ci>n_NaCa</ci>
                                 <cn cellml:units="dimensionless">2</cn>
                              </apply>
                           </apply>
                        </apply>
                        <ci>i_rel</ci>
                     </apply>
                     <ci>i_up</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <cn cellml:units="dimensionless">2</cn>
                  <cn cellml:units="litre_millimetre3">1</cn>
                  <ci>Vi</ci>
                  <ci>F</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="extracellular_potassium_concentration">
      <variable name="Kc" units="millimolar" initial_value="3" public_interface="out"/>
      <variable name="V_e" units="millimetre3"/>
      <variable name="V_e_ratio" units="dimensionless" public_interface="in"/>
      <variable name="V_Cell" units="millimetre3" public_interface="in"/>
      <variable name="Kb" units="millimolar" initial_value="3"/>
      <variable name="i_mK" units="nanoA" public_interface="out"/>
      <variable name="pf" units="per_second" initial_value="1"/>
      <variable name="time" units="second" public_interface="in"/>
      <variable name="F" units="coulomb_per_mole" public_interface="in"/>
      <variable name="i_K1" units="nanoA" public_interface="in"/>
      <variable name="i_K" units="nanoA" public_interface="in"/>
      <variable name="i_fK" units="nanoA" public_interface="in"/>
      <variable name="i_siK" units="nanoA" public_interface="in"/>
      <variable name="i_p" units="nanoA" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>V_e</ci>
            <apply>
               <times/>
               <ci>V_Cell</ci>
               <ci>V_e_ratio</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>i_mK</ci>
            <apply>
               <minus/>
               <apply>
                  <plus/>
                  <ci>i_K1</ci>
                  <ci>i_K</ci>
                  <ci>i_fK</ci>
                  <ci>i_siK</ci>
               </apply>
               <apply>
                  <times/>
                  <cn cellml:units="dimensionless">2</cn>
                  <ci>i_p</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>Kc</ci>
            </apply>
            <apply>
               <plus/>
               <apply>
                  <times/>
                  <apply>
                     <minus/>
                     <ci>pf</ci>
                  </apply>
                  <apply>
                     <minus/>
                     <ci>Kc</ci>
                     <ci>Kb</ci>
                  </apply>
               </apply>
               <apply>
                  <divide/>
                  <apply>
                     <times/>
                     <cn cellml:units="mA_nA">1</cn>
                     <ci>i_mK</ci>
                  </apply>
                  <apply>
                     <times/>
                     <cn cellml:units="litre_millimetre3">1</cn>
                     <ci>V_e</ci>
                     <ci>F</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="intracellular_potassium_concentration">
      <variable name="Ki" units="millimolar" initial_value="140" public_interface="out"/>
      <variable name="Vi" units="millimetre3" public_interface="in"/>
      <variable name="i_mK" units="nanoA" public_interface="in"/>
      <variable name="time" units="second" public_interface="in"/>
      <variable name="F" units="coulomb_per_mole" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>Ki</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <apply>
                     <minus/>
                     <cn cellml:units="mA_nA">1</cn>
                  </apply>
                  <ci>i_mK</ci>
               </apply>
               <apply>
                  <times/>
                  <cn cellml:units="litre_millimetre3">1</cn>
                  <ci>Vi</ci>
                  <ci>F</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <group>
      <relationship_ref relationship="containment"/>
      <component_ref component="membrane">
         <component_ref component="hyperpolarising_activated_current">
            <component_ref component="hyperpolarising_activated_current_y_gate"/>
         </component_ref>
         <component_ref component="time_dependent_potassium_current">
            <component_ref component="time_dependent_potassium_current_x_gate"/>
         </component_ref>
         <component_ref component="time_independent_potassium_current"/>
         <component_ref component="sodium_background_current"/>
         <component_ref component="calcium_background_current"/>
         <component_ref component="sodium_potassium_pump"/>
         <component_ref component="Na_Ca_exchanger"/>
         <component_ref component="fast_sodium_current">
            <component_ref component="fast_sodium_current_m_gate"/>
            <component_ref component="fast_sodium_current_h_gate"/>
         </component_ref>
         <component_ref component="second_inward_current">
            <component_ref component="second_inward_current_d_gate"/>
            <component_ref component="second_inward_current_f_gate"/>
            <component_ref component="second_inward_current_f2_gate"/>
         </component_ref>
         <component_ref component="intracellular_sodium_concentration"/>
         <component_ref component="extracellular_sodium_concentration"/>
         <component_ref component="intracellular_calcium_concentration"/>
         <component_ref component="extracellular_calcium_concentration"/>
         <component_ref component="extracellular_potassium_concentration"/>
         <component_ref component="intracellular_potassium_concentration"/>
      </component_ref>
   </group>
   <group>
      <relationship_ref relationship="encapsulation"/>
      <component_ref component="hyperpolarising_activated_current">
         <component_ref component="hyperpolarising_activated_current_y_gate"/>
      </component_ref>
      <component_ref component="time_dependent_potassium_current">
         <component_ref component="time_dependent_potassium_current_x_gate"/>
      </component_ref>
      <component_ref component="fast_sodium_current">
         <component_ref component="fast_sodium_current_m_gate"/>
         <component_ref component="fast_sodium_current_h_gate"/>
      </component_ref>
      <component_ref component="second_inward_current">
         <component_ref component="second_inward_current_d_gate"/>
         <component_ref component="second_inward_current_f_gate"/>
         <component_ref component="second_inward_current_f2_gate"/>
      </component_ref>
   </group>
   <connection>
      <map_components component_1="intracellular_sodium_concentration" component_2="extracellular_potassium_concentration"/>
      <map_variables variable_1="V_e_ratio" variable_2="V_e_ratio"/>
      <map_variables variable_1="V_Cell" variable_2="V_Cell"/>
   </connection>
   <connection>
      <map_components component_1="membrane" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="hyperpolarising_activated_current" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="time_dependent_potassium_current" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="time_independent_potassium_current" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="sodium_background_current" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="calcium_background_current" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="sodium_potassium_pump" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="Na_Ca_exchanger" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="fast_sodium_current" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="second_inward_current" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="intracellular_sodium_concentration" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="intracellular_calcium_concentration" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="extracellular_potassium_concentration" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="intracellular_potassium_concentration" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="hyperpolarising_activated_current" component_2="membrane"/>
      <map_variables variable_1="V" variable_2="V"/>
      <map_variables variable_1="RTONF" variable_2="RTONF"/>
      <map_variables variable_1="i_f" variable_2="i_f"/>
   </connection>
   <connection>
      <map_components component_1="time_dependent_potassium_current" component_2="membrane"/>
      <map_variables variable_1="V" variable_2="V"/>
      <map_variables variable_1="i_K" variable_2="i_K"/>
      <map_variables variable_1="RTONF" variable_2="RTONF"/>
   </connection>
   <connection>
      <map_components component_1="time_independent_potassium_current" component_2="membrane"/>
      <map_variables variable_1="V" variable_2="V"/>
      <map_variables variable_1="RTONF" variable_2="RTONF"/>
      <map_variables variable_1="i_K1" variable_2="i_K1"/>
   </connection>
   <connection>
      <map_components component_1="sodium_background_current" component_2="membrane"/>
      <map_variables variable_1="V" variable_2="V"/>
      <map_variables variable_1="i_Na_b" variable_2="i_Na_b"/>
   </connection>
   <connection>
      <map_components component_1="calcium_background_current" component_2="membrane"/>
      <map_variables variable_1="V" variable_2="V"/>
      <map_variables variable_1="RTONF" variable_2="RTONF"/>
      <map_variables variable_1="i_Ca_b" variable_2="i_Ca_b"/>
   </connection>
   <connection>
      <map_components component_1="sodium_potassium_pump" component_2="membrane"/>
      <map_variables variable_1="V" variable_2="V"/>
      <map_variables variable_1="i_p" variable_2="i_p"/>
   </connection>
   <connection>
      <map_components component_1="Na_Ca_exchanger" component_2="membrane"/>
      <map_variables variable_1="V" variable_2="V"/>
      <map_variables variable_1="i_NaCa" variable_2="i_NaCa"/>
      <map_variables variable_1="RTONF" variable_2="RTONF"/>
   </connection>
   <connection>
      <map_components component_1="fast_sodium_current" component_2="membrane"/>
      <map_variables variable_1="V" variable_2="V"/>
      <map_variables variable_1="i_Na" variable_2="i_Na"/>
      <map_variables variable_1="RTONF" variable_2="RTONF"/>
   </connection>
   <connection>
      <map_components component_1="second_inward_current" component_2="membrane"/>
      <map_variables variable_1="V" variable_2="V"/>
      <map_variables variable_1="i_Caf" variable_2="i_Caf"/>
      <map_variables variable_1="RTONF" variable_2="RTONF"/>
   </connection>
   <connection>
      <map_components component_1="intracellular_sodium_concentration" component_2="membrane"/>
      <map_variables variable_1="F" variable_2="F"/>
   </connection>
   <connection>
      <map_components component_1="intracellular_calcium_concentration" component_2="membrane"/>
      <map_variables variable_1="V" variable_2="V"/>
      <map_variables variable_1="F" variable_2="F"/>
   </connection>
   <connection>
      <map_components component_1="extracellular_potassium_concentration" component_2="membrane"/>
      <map_variables variable_1="F" variable_2="F"/>
   </connection>
   <connection>
      <map_components component_1="intracellular_potassium_concentration" component_2="membrane"/>
      <map_variables variable_1="F" variable_2="F"/>
   </connection>
   <connection>
      <map_components component_1="hyperpolarising_activated_current" component_2="extracellular_potassium_concentration"/>
      <map_variables variable_1="Kc" variable_2="Kc"/>
      <map_variables variable_1="i_fK" variable_2="i_fK"/>
   </connection>
   <connection>
      <map_components component_1="hyperpolarising_activated_current" component_2="intracellular_potassium_concentration"/>
      <map_variables variable_1="Ki" variable_2="Ki"/>
   </connection>
   <connection>
      <map_components component_1="hyperpolarising_activated_current" component_2="time_independent_potassium_current"/>
      <map_variables variable_1="E_K" variable_2="E_K"/>
   </connection>
   <connection>
      <map_components component_1="hyperpolarising_activated_current" component_2="sodium_background_current"/>
      <map_variables variable_1="E_Na" variable_2="E_Na"/>
   </connection>
   <connection>
      <map_components component_1="hyperpolarising_activated_current" component_2="intracellular_sodium_concentration"/>
      <map_variables variable_1="i_fNa" variable_2="i_fNa"/>
      <map_variables variable_1="Nai" variable_2="Nai"/>
   </connection>
   <connection>
      <map_components component_1="hyperpolarising_activated_current" component_2="extracellular_sodium_concentration"/>
      <map_variables variable_1="Nao" variable_2="Nao"/>
   </connection>
   <connection>
      <map_components component_1="calcium_background_current" component_2="extracellular_calcium_concentration"/>
      <map_variables variable_1="Cao" variable_2="Cao"/>
   </connection>
   <connection>
      <map_components component_1="time_dependent_potassium_current" component_2="extracellular_potassium_concentration"/>
      <map_variables variable_1="Kc" variable_2="Kc"/>
      <map_variables variable_1="i_K" variable_2="i_K"/>
   </connection>
   <connection>
      <map_components component_1="time_dependent_potassium_current" component_2="intracellular_potassium_concentration"/>
      <map_variables variable_1="Ki" variable_2="Ki"/>
   </connection>
   <connection>
      <map_components component_1="time_independent_potassium_current" component_2="extracellular_potassium_concentration"/>
      <map_variables variable_1="Kc" variable_2="Kc"/>
      <map_variables variable_1="i_K1" variable_2="i_K1"/>
   </connection>
   <connection>
      <map_components component_1="sodium_background_current" component_2="intracellular_sodium_concentration"/>
      <map_variables variable_1="i_Na_b" variable_2="i_Na_b"/>
   </connection>
   <connection>
      <map_components component_1="calcium_background_current" component_2="intracellular_calcium_concentration"/>
      <map_variables variable_1="i_Ca_b" variable_2="i_Ca_b"/>
      <map_variables variable_1="Cai" variable_2="Cai"/>
   </connection>
   <connection>
      <map_components component_1="sodium_potassium_pump" component_2="extracellular_potassium_concentration"/>
      <map_variables variable_1="i_p" variable_2="i_p"/>
      <map_variables variable_1="Kc" variable_2="Kc"/>
   </connection>
   <connection>
      <map_components component_1="sodium_potassium_pump" component_2="intracellular_sodium_concentration"/>
      <map_variables variable_1="i_p" variable_2="i_p"/>
      <map_variables variable_1="Nai" variable_2="Nai"/>
   </connection>
   <connection>
      <map_components component_1="Na_Ca_exchanger" component_2="extracellular_sodium_concentration"/>
      <map_variables variable_1="Nao" variable_2="Nao"/>
   </connection>
   <connection>
      <map_components component_1="Na_Ca_exchanger" component_2="extracellular_calcium_concentration"/>
      <map_variables variable_1="Cao" variable_2="Cao"/>
   </connection>
   <connection>
      <map_components component_1="Na_Ca_exchanger" component_2="intracellular_sodium_concentration"/>
      <map_variables variable_1="n_NaCa" variable_2="n_NaCa"/>
      <map_variables variable_1="i_NaCa" variable_2="i_NaCa"/>
      <map_variables variable_1="Nai" variable_2="Nai"/>
   </connection>
   <connection>
      <map_components component_1="Na_Ca_exchanger" component_2="intracellular_calcium_concentration"/>
      <map_variables variable_1="n_NaCa" variable_2="n_NaCa"/>
      <map_variables variable_1="i_NaCa" variable_2="i_NaCa"/>
      <map_variables variable_1="Cai" variable_2="Cai"/>
   </connection>
   <connection>
      <map_components component_1="fast_sodium_current" component_2="intracellular_sodium_concentration"/>
      <map_variables variable_1="i_Na" variable_2="i_Na"/>
      <map_variables variable_1="Nai" variable_2="Nai"/>
   </connection>
   <connection>
      <map_components component_1="fast_sodium_current" component_2="intracellular_potassium_concentration"/>
      <map_variables variable_1="Ki" variable_2="Ki"/>
   </connection>
   <connection>
      <map_components component_1="fast_sodium_current" component_2="extracellular_potassium_concentration"/>
      <map_variables variable_1="Kc" variable_2="Kc"/>
   </connection>
   <connection>
      <map_components component_1="fast_sodium_current" component_2="extracellular_sodium_concentration"/>
      <map_variables variable_1="Nao" variable_2="Nao"/>
   </connection>
   <connection>
      <map_components component_1="second_inward_current" component_2="intracellular_calcium_concentration"/>
      <map_variables variable_1="i_siCa" variable_2="i_siCa"/>
      <map_variables variable_1="Cai" variable_2="Cai"/>
   </connection>
   <connection>
      <map_components component_1="second_inward_current" component_2="extracellular_calcium_concentration"/>
      <map_variables variable_1="Cao" variable_2="Cao"/>
   </connection>
   <connection>
      <map_components component_1="second_inward_current" component_2="extracellular_sodium_concentration"/>
      <map_variables variable_1="Nao" variable_2="Nao"/>
   </connection>
   <connection>
      <map_components component_1="second_inward_current" component_2="extracellular_potassium_concentration"/>
      <map_variables variable_1="i_siK" variable_2="i_siK"/>
      <map_variables variable_1="Kc" variable_2="Kc"/>
   </connection>
   <connection>
      <map_components component_1="second_inward_current" component_2="intracellular_sodium_concentration"/>
      <map_variables variable_1="i_siNa" variable_2="i_siNa"/>
      <map_variables variable_1="Nai" variable_2="Nai"/>
   </connection>
   <connection>
      <map_components component_1="second_inward_current" component_2="intracellular_potassium_concentration"/>
      <map_variables variable_1="Ki" variable_2="Ki"/>
   </connection>
   <connection>
      <map_components component_1="hyperpolarising_activated_current" component_2="hyperpolarising_activated_current_y_gate"/>
      <map_variables variable_1="y" variable_2="y"/>
      <map_variables variable_1="time" variable_2="time"/>
      <map_variables variable_1="V" variable_2="V"/>
   </connection>
   <connection>
      <map_components component_1="time_dependent_potassium_current" component_2="time_dependent_potassium_current_x_gate"/>
      <map_variables variable_1="x" variable_2="x"/>
      <map_variables variable_1="time" variable_2="time"/>
      <map_variables variable_1="V" variable_2="V"/>
   </connection>
   <connection>
      <map_components component_1="fast_sodium_current" component_2="fast_sodium_current_m_gate"/>
      <map_variables variable_1="m" variable_2="m"/>
      <map_variables variable_1="time" variable_2="time"/>
      <map_variables variable_1="V" variable_2="V"/>
   </connection>
   <connection>
      <map_components component_1="fast_sodium_current" component_2="fast_sodium_current_h_gate"/>
      <map_variables variable_1="h" variable_2="h"/>
      <map_variables variable_1="time" variable_2="time"/>
      <map_variables variable_1="V" variable_2="V"/>
   </connection>
   <connection>
      <map_components component_1="second_inward_current" component_2="second_inward_current_d_gate"/>
      <map_variables variable_1="d" variable_2="d"/>
      <map_variables variable_1="time" variable_2="time"/>
      <map_variables variable_1="V" variable_2="V"/>
   </connection>
   <connection>
      <map_components component_1="second_inward_current" component_2="second_inward_current_f_gate"/>
      <map_variables variable_1="f" variable_2="f"/>
      <map_variables variable_1="time" variable_2="time"/>
      <map_variables variable_1="V" variable_2="V"/>
   </connection>
   <connection>
      <map_components component_1="second_inward_current" component_2="second_inward_current_f2_gate"/>
      <map_variables variable_1="f2" variable_2="f2"/>
      <map_variables variable_1="time" variable_2="time"/>
      <map_variables variable_1="V" variable_2="V"/>
      <map_variables variable_1="Cai" variable_2="Cai"/>
   </connection>
   <connection>
      <map_components component_1="intracellular_sodium_concentration" component_2="intracellular_calcium_concentration"/>
      <map_variables variable_1="Vi" variable_2="Vi"/>
   </connection>
   <connection>
      <map_components component_1="intracellular_sodium_concentration" component_2="intracellular_potassium_concentration"/>
      <map_variables variable_1="Vi" variable_2="Vi"/>
   </connection>
   <connection>
      <map_components component_1="extracellular_potassium_concentration" component_2="intracellular_potassium_concentration"/>
      <map_variables variable_1="i_mK" variable_2="i_mK"/>
   </connection>
   
   <rdf:RDF>
  <rdf:Seq rdf:about="rdf:#5ce26374-249d-49c3-9af9-7e2627bee5e2">
    <rdf:li rdf:resource="rdf:#bae873c2-9b4e-42e5-9d7b-1b07ef8fbefb"/>
    <rdf:li rdf:resource="rdf:#8aa35288-c71a-45e8-8d33-6c6287ffa6fe"/>
    <rdf:li rdf:resource="rdf:#4f4a857a-8aac-4c3d-b1e6-7004aaad5001"/>
  </rdf:Seq>
  <rdf:Description rdf:about="rdf:#ba0baeba-09d2-412b-8530-4640da1b3226">
    <dc:title>Cellular and Neuronal Oscillators</dc:title>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#335cfb99-a4d2-4d6e-ab7e-4e07850481f0">
    <dcterms:W3CDTF>1989-01-01</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#f52c5ecd-ce40-4663-9397-208a13d92a79">
    <vCard:Orgname>Oxford University</vCard:Orgname>
    <vCard:Orgunit>unknown</vCard:Orgunit>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#68364949-45e3-4d35-a168-e2c469b7454b">
    <vCard:Given>Penny</vCard:Given>
    <vCard:Family>Noble</vCard:Family>
    <vCard:Other>J</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#109a6917-0710-4a06-b023-90c54eebef1d">
    <vCard:Given>Penny</vCard:Given>
    <vCard:Family>Noble</vCard:Family>
    <vCard:Other>unknown</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#a4ecf7ea-da89-4e45-b4e5-100e53ad98d5">
    <vCard:Given>D</vCard:Given>
    <vCard:Family>DiFrancesco</vCard:Family>
  </rdf:Description>
  <rdf:Description rdf:about="#noble_difrancesco_denyer_1989">
    <ns7:simulation rdf:resource="rdf:#$7WTPD2"/>
    <bqs:reference rdf:resource="rdf:#9ed26340-1df5-4a2e-8523-c1eeb8854c55"/>
<bqs:reference rdf:parseType="Resource">
  <dc:subject rdf:parseType="Resource">
    <bqs:subject_type>keyword</bqs:subject_type>
    <rdf:value>
      <rdf:Bag>
        <rdf:li>cardiac</rdf:li>
        <rdf:li>electrophysiology</rdf:li>
        <rdf:li>pacemaker</rdf:li>
      </rdf:Bag>
    </rdf:value>
  </dc:subject>
</bqs:reference>
    <cmeta:comment rdf:resource="rdf:#76c0570b-573f-4c68-88c9-ae1435d6270e"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#76c0570b-573f-4c68-88c9-ae1435d6270e">
    <dc:creator rdf:resource="rdf:#3815eb1e-986f-435b-a46a-d9995b775ee9"/>
    <rdf:value>This file is a CellML description of the model described in the article "Ionic Mechanisms in Normal and Abnormal Cardiac Pacemaker Activity" in the book "Neuronal and Cellular Oscillators," pages 59-85</rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#42aa4574-d98b-4504-8dd3-995a97663916">
    <dcterms:modified rdf:resource="rdf:#24c6b2dd-c914-4b07-b405-b38c3ed537ec"/>
    <rdf:value>Units checked, curated.</rdf:value>
    <cmeta:modifier rdf:resource="rdf:#bcf78a16-f89e-45fc-b8dc-eb3e3f462adc"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#0b666db7-deb8-4b53-b3ab-8fa76d96c7ab">
    <dcterms:W3CDTF>2006-03-31</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#f58ef4d2-71a1-4dde-8112-6f215d814600">
    <vCard:Given>J</vCard:Given>
    <vCard:Family>Denyer</vCard:Family>
    <vCard:Other>C</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#eed511e8-073f-491f-b75b-7731ce47e223">
    <dcterms:W3CDTF>2005-07-05</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#$9WTPD2">
    <ns7:nonstandard-pointDensity>100000</ns7:nonstandard-pointDensity>
    <ns7:maximumStepSize>0.001</ns7:maximumStepSize>
  </rdf:Description>
  <rdf:Description rdf:about="">
    <cmeta:modification rdf:resource="rdf:#17b42d92-6343-4695-8355-b54fbffaf7e6"/>
    <dcterms:created rdf:resource="rdf:#eed511e8-073f-491f-b75b-7731ce47e223"/>
    <dc:creator rdf:resource="rdf:#88f15bab-7991-4d4b-84b1-f20b57db7361"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#8aa35288-c71a-45e8-8d33-6c6287ffa6fe">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#a4ecf7ea-da89-4e45-b4e5-100e53ad98d5"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#16a5cce5-27a8-4c82-aa55-3afea2b9fb71">
    <dc:creator rdf:resource="rdf:#631a6439-7e9e-4cbd-a151-8a00065825ce"/>
    <rdf:value>This model has been curated by Penny Noble of Oxford University and is known to run in PCEnv and COR and reproduce the results published in the paper it is based on.</rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#496f0fe2-8556-46ee-9378-7479511fcdcd">
    <vCard:Orgname>Oxford University</vCard:Orgname>
    <vCard:Orgunit>Department of Physiology, Anatomy &amp; Genetics</vCard:Orgunit>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#$8WTPD2">
    <rdf:rest rdf:resource="http://www.w3.org/1999/02/22-rdf-syntax-ns#nil"/>
    <rdf:first rdf:resource="rdf:#$9WTPD2"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#61b1fc83-d5f0-45fa-b281-4fbeec9cdfbb">
    <vCard:N rdf:resource="rdf:#60f7b64d-5cc1-4b16-8ad3-46d289d2e902"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#d3ddf607-2cd3-440f-acb5-aec927fc62c7">
    <rdf:type rdf:resource="http://imc.org/vCard/3.0#internet"/>
    <rdf:value>penny.noble@physiol.ox.ac.uk</rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#4f4a857a-8aac-4c3d-b1e6-7004aaad5001">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#f58ef4d2-71a1-4dde-8112-6f215d814600"/>
  </rdf:Description>
  <rdf:Description rdf:about="">
    <dc:publisher>Department of Physiology, Anatomy &amp; Genetics, University of Oxford</dc:publisher>
    <cmeta:comment rdf:resource="rdf:#16a5cce5-27a8-4c82-aa55-3afea2b9fb71"/>
    <dcterms:created rdf:resource="rdf:#8e8365a5-6c4f-45b8-b69e-49e319d4e2d1"/>
    <dc:creator rdf:resource="rdf:#59f41a5f-2c1d-4cd3-929c-e393213e5708"/>
    <cmeta:modification rdf:resource="rdf:#42aa4574-d98b-4504-8dd3-995a97663916"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#ca282cb6-45c8-4136-a669-5be3b8dd1dc3">
    <dc:creator rdf:resource="rdf:#5ce26374-249d-49c3-9af9-7e2627bee5e2"/>
    <dc:title>Ionic Mechanisms in Normal and Abnormal Cardiac Pacemaker Activity</dc:title>
    <bqs:volume/>
    <bqs:first_page>59</bqs:first_page>
    <bqs:Journal rdf:resource="rdf:#ba0baeba-09d2-412b-8530-4640da1b3226"/>
    <dcterms:issued rdf:resource="rdf:#335cfb99-a4d2-4d6e-ab7e-4e07850481f0"/>
    <bqs:last_page>85</bqs:last_page>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#9ed26340-1df5-4a2e-8523-c1eeb8854c55">
    <bqs:Pubmed_id/>
    <bqs:JournalArticle rdf:resource="rdf:#ca282cb6-45c8-4136-a669-5be3b8dd1dc3"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#c6e85afb-789f-4ac5-a9ca-ffaa678467f9">
    <rdf:type rdf:resource="http://imc.org/vCard/3.0#internet"/>
    <rdf:value>penny.noble@dpag.ox.ac.uk</rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#254dc983-3007-447a-8c03-aba525176b6f">
    <vCard:Given>Penny</vCard:Given>
    <vCard:Family>Noble</vCard:Family>
    <vCard:Other>J</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#$7WTPD2">
    <ns7:boundIntervals rdf:resource="rdf:#$8WTPD2"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#24c6b2dd-c914-4b07-b405-b38c3ed537ec">
    <dcterms:W3CDTF>2007-12-04T15:37:59+13:00</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#bae873c2-9b4e-42e5-9d7b-1b07ef8fbefb">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#37071d3d-cb37-41db-8edf-53e6606de12b"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#631a6439-7e9e-4cbd-a151-8a00065825ce">
    <vCard:FN>James Lawson</vCard:FN>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#37071d3d-cb37-41db-8edf-53e6606de12b">
    <vCard:Given>D</vCard:Given>
    <vCard:Family>Noble</vCard:Family>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#17b42d92-6343-4695-8355-b54fbffaf7e6">
    <dcterms:modified rdf:resource="rdf:#0b666db7-deb8-4b53-b3ab-8fa76d96c7ab"/>
    <rdf:value>added metadata</rdf:value>
    <cmeta:modifier rdf:resource="rdf:#61b1fc83-d5f0-45fa-b281-4fbeec9cdfbb"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#8e8365a5-6c4f-45b8-b69e-49e319d4e2d1">
    <dcterms:W3CDTF>2007-09-20T00:00:00+00:00</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#88f15bab-7991-4d4b-84b1-f20b57db7361">
    <vCard:ORG rdf:resource="rdf:#f52c5ecd-ce40-4663-9397-208a13d92a79"/>
    <vCard:EMAIL rdf:resource="rdf:#d3ddf607-2cd3-440f-acb5-aec927fc62c7"/>
    <vCard:N rdf:resource="rdf:#109a6917-0710-4a06-b023-90c54eebef1d"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#60f7b64d-5cc1-4b16-8ad3-46d289d2e902">
    <vCard:Given>Peter</vCard:Given>
    <vCard:Family>Villiger</vCard:Family>
    <vCard:Other>J</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#59f41a5f-2c1d-4cd3-929c-e393213e5708">
    <vCard:ORG rdf:resource="rdf:#496f0fe2-8556-46ee-9378-7479511fcdcd"/>
    <vCard:EMAIL rdf:resource="rdf:#c6e85afb-789f-4ac5-a9ca-ffaa678467f9"/>
    <vCard:N rdf:resource="rdf:#68364949-45e3-4d35-a168-e2c469b7454b"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#3815eb1e-986f-435b-a46a-d9995b775ee9">
    <vCard:FN>James Lawson</vCard:FN>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#bcf78a16-f89e-45fc-b8dc-eb3e3f462adc">
    <vCard:N rdf:resource="rdf:#254dc983-3007-447a-8c03-aba525176b6f"/>
  </rdf:Description>
</rdf:RDF>

</model>