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