- Author:
- pmr2.import <nobody@models.cellml.org>
- Date:
- 2009-06-17 15:08:07+12:00
- Desc:
- committing version04 of matsuoka_sarai_kuratomi_ono_noma_2003
- Permanent Source URI:
- https://models.cellml.org/workspace/matsuoka_sarai_kuratomi_ono_noma_2003/rawfile/fc16f608b40f05c69f4fe013aedb6a159ddf6aa7/matsuoka_sarai_kuratomi_ono_noma_2003.cellml
<?xml version='1.0' encoding='utf-8'?>
<!--
This CellML file was generated on 18/04/2008 at 09:15:12 using:
COR (0.9.31.915)
Copyright 2002-2008 Dr Alan Garny
http://COR.physiol.ox.ac.uk/ - COR@physiol.ox.ac.uk
CellML 1.0 was used to generate this cellular model
http://www.CellML.org/
--><model xmlns="http://www.cellml.org/cellml/1.0#" xmlns:cmeta="http://www.cellml.org/metadata/1.0#" xmlns:dc="http://purl.org/dc/elements/1.1/" xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#" xmlns:bqs="http://www.cellml.org/bqs/1.0#" xmlns:cellml="http://www.cellml.org/cellml/1.0#" xmlns:dcterms="http://purl.org/dc/terms/" xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" xmlns:ns7="http://www.cellml.org/metadata/simulation/1.0#" cmeta:id="matsuoka_sarai_kuratomi_ono_noma_2003" name="matsuoka_sarai_kuratomi_ono_noma_2003">
<documentation xmlns="http://cellml.org/tmp-documentation">
<article>
<articleinfo>
<title>Modelling the Roles of Individual Current Systems in Ventricular Cells and in the SA Node</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 model has been curated and unit checked by Alan Garny in collaboration with the model authors and it is known to run in COR to reproduce the published results. It also works in JSim as long as the stimulus protocol is altered.
</para>
</section>
<sect1 id="sec_structure">
<title>Model Structure</title>
<para>
Different types of ion channel, with their distinct gating and conductance properties, play a unique role in generating membrane excitation. Although the behaviours of single ion channels have been extensively studied experimentally, ionic current interactions make it difficult to deduce their individual roles in membrane excitability. In order to better understand the roles of the individual ion currents, Satoshi Matsuoka <emphasis>et al.</emphasis> have developed a mathematical model of membrane excitation for the ventricular cell.
</para>
<para>
The Matsuoka <emphasis>et al.</emphasis> mathematical model is a development of previously published cardiac cellular models. The group's general approach to modelling was based on the methods used in the Difrancesco-Noble, 1985 Model and the Luo-Rudy Ventricular Model II (dynamic), 1994. In addition, the compound <emphasis>Kyoto Model</emphasis> of Matsuoka <emphasis>et al.</emphasis> incorporates the Negroni and Lascano (1996) model of cardiac muscle contraction, allowing simulation of sarcomere shortening, and also, the gating model of Shirokov <emphasis>et al.</emphasis> is used for both the fast Na<superscript>+</superscript> and L-type Ca<superscript>2+</superscript> channels.
</para>
<para>
New experimental data are included, such as new kinetics of the inward rectifier K<superscript>+</superscript> channel, the delayed rectifier K<superscript>+</superscript> channel, and the sustained inward current (for more details on the reaction kinetics, see figure 3 below). Model simulations showed that ion conductance by the L-type Ca<superscript>2+</superscript> channel (I<subscript>CaL</subscript>) dominates the membrane conductance during the action potential. Repolarisation occurs as the following currents are sequentially activated: I<subscript>Ks</subscript>; I<subscript>Kr</subscript>; and I<subscript>K1</subscript>. Agreement between model simulation results and experimental data was taken as supporting evidence for the validity of the model.
</para>
<para>
The complete original paper reference is cited below:
</para>
<para>
<ulink url="http://www.jstage.jst.go.jp/article/jjphysiol/53/2/53_105/_article">Role of Individual Ionic Current Systems in Ventricular Cells Hypothesized by a Model Study</ulink>, Matsuoka S, Sarai N, Kuratomi S, Ono K, and Noma A, 2003, <ulink url="http://wwwsoc.nii.ac.jp/psj/jjp/jjphome.html">
<emphasis>The Japanese Journal of Physiology</emphasis>
</ulink>, 53, 105-123. (A PDF version of the article is available to subscribers on the <emphasis>Japanese Journal of Physiology</emphasis> website.) <ulink url="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?cmd=Retrieve&db=PubMed&list_uids=12877767&dopt=Abstract">PubMed ID: 12877767</ulink>
</para>
<informalfigure float="0" id="fig_ventricular_cell_diagram">
<mediaobject>
<imageobject>
<objectinfo>
<title>ventricular cell diagram</title>
</objectinfo>
<imagedata fileref="matsuoka_2003.png"/>
</imageobject>
</mediaobject>
<caption>A schematic diagram describing the ionic components of the Matsuoka <emphasis>et al.</emphasis> 2003 ventricular cell mathematical model. Ions are exchanged between the intracellular and extracellular environments, and between the cytosol and the sarcoplasmic reticulum (SR) via channels and pumps. Unique to the ventricular cells are the transient outward current, I<subscript>to</subscript>, and the slow component of the delayed rectifier K<superscript>+</superscript> current, I<subscript>Ks</subscript>.</caption>
</informalfigure>
<informalfigure float="0" id="fig_reaction_kinetics">
<mediaobject>
<imageobject>
<objectinfo>
<title>reaction kinetics diagram</title>
</objectinfo>
<imagedata fileref="matsuoka_3_2003.png"/>
</imageobject>
</mediaobject>
<caption>A schematic diagram describing the reaction kinetics of the ion channel gates in the Matsuoka <emphasis>et al.</emphasis> 2003 mathematical models. A) The kinetics of the sodium channel voltage-dependent gate. B) The kinetics of the Na-Ca exchange pump. C) The kinetics of the Na-K pump. D) The kinetics of the inward rectifier K<superscript>+</superscript> current. E) The kinetics of the Ca<superscript>2+</superscript>-dependent gate of the L-type Ca<superscript>2+</superscript> channel. F) The kinetics of the ultra-slow gates of several ion channels. G) The kinetics of the SR Ca<superscript>2+</superscript> pump. H) The kinetics of the hyperpolarisation-activated cation current. I) The kinetics of the RyR channel in the SR.</caption>
</informalfigure>
</sect1>
</article>
</documentation>
<units name="millivolt">
<unit units="volt" prefix="milli"/>
</units>
<units name="per_millivolt">
<unit units="millivolt" exponent="-1"/>
</units>
<units name="millimole">
<unit units="mole" prefix="milli"/>
</units>
<units name="millimolar">
<unit units="mole" prefix="milli"/>
<unit units="litre" exponent="-1"/>
</units>
<units name="millimolar2">
<unit units="millimolar" exponent="2"/>
</units>
<units name="micromolar">
<unit units="mole" prefix="micro"/>
<unit units="litre" exponent="-1"/>
</units>
<units name="millimolar_per_millisecond">
<unit units="millimolar"/>
<unit units="millisecond" exponent="-1"/>
</units>
<units name="millimolar_per_picoA">
<unit units="millimolar"/>
<unit units="picoA" exponent="-1"/>
</units>
<units name="per_millimolar">
<unit units="millimolar" exponent="-1"/>
</units>
<units name="per_micromolar2">
<unit units="micromolar" exponent="-2"/>
</units>
<units name="second_per_micromolar2">
<unit units="second"/>
<unit units="micromolar" exponent="-2"/>
</units>
<units name="millisecond_per_micrometre2">
<unit units="millisecond"/>
<unit units="micrometre" exponent="-2"/>
</units>
<units name="mN_per_mm2">
<unit units="newton" prefix="milli"/>
<unit units="metre" prefix="milli" exponent="-2"/>
</units>
<units name="mN_per_mm2_micrometre">
<unit units="mN_per_mm2"/>
<unit units="micrometre" exponent="-1"/>
</units>
<units name="mN_per_mm2_micrometre5">
<unit units="mN_per_mm2"/>
<unit units="micrometre" exponent="-5"/>
</units>
<units name="mN_per_mm2_micrometre_millimolar">
<unit units="newton" prefix="milli"/>
<unit units="metre" prefix="milli" exponent="-2"/>
<unit units="micrometre" exponent="-1"/>
<unit units="millimolar" exponent="-1"/>
</units>
<units name="per_micrometre2">
<unit units="micrometre" exponent="-2"/>
</units>
<units name="millisecond">
<unit units="second" prefix="milli"/>
</units>
<units name="micromolar_per_second">
<unit units="micromolar"/>
<unit units="second" exponent="-1"/>
</units>
<units name="per_millisecond">
<unit units="millisecond" exponent="-1"/>
</units>
<units name="per_millimolar_millisecond">
<unit units="millimolar" exponent="-1"/>
<unit units="millisecond" exponent="-1"/>
</units>
<units name="picoF">
<unit units="farad" prefix="pico"/>
</units>
<units name="picoF_per_millisecond">
<unit units="picoF"/>
<unit units="millisecond" exponent="-1"/>
</units>
<units name="picoA">
<unit units="ampere" prefix="pico"/>
</units>
<units name="per_picoA_millisecond">
<unit units="picoA" exponent="-1"/>
<unit units="millisecond" exponent="-1"/>
</units>
<units name="micrometre">
<unit units="metre" prefix="micro"/>
</units>
<units name="micrometre_per_millisecond">
<unit units="micrometre"/>
<unit units="millisecond" exponent="-1"/>
</units>
<units name="micrometre3">
<unit units="metre" prefix="micro" exponent="3"/>
</units>
<units name="nanoS">
<unit units="siemens" prefix="nano"/>
</units>
<units name="nanoS_per_picoF">
<unit units="nanoS"/>
<unit units="picoF" exponent="-1"/>
</units>
<units name="picoA_per_picoF">
<unit units="picoA"/>
<unit units="picoF" exponent="-1"/>
</units>
<units name="nanoS_per_millimolar">
<unit units="nanoS"/>
<unit units="millimolar" exponent="-1"/>
</units>
<units name="picoA_per_millimolar">
<unit units="ampere" prefix="pico"/>
<unit units="millimolar" exponent="-1"/>
</units>
<units name="coulomb_millivolt_per_kelvin_millimole">
<unit units="coulomb"/>
<unit units="millivolt"/>
<unit units="kelvin" exponent="-1"/>
<unit units="millimole" exponent="-1"/>
</units>
<units name="coulomb_per_millimole">
<unit units="coulomb"/>
<unit units="mole" prefix="milli" exponent="-1"/>
</units>
<component name="environment">
<variable name="time" units="millisecond" public_interface="out"/>
</component>
<component name="membrane">
<variable name="Vm" units="millivolt" initial_value="-85.95752434460744" public_interface="out"/>
<variable name="R" units="coulomb_millivolt_per_kelvin_millimole" initial_value="8.3143" public_interface="out"/>
<variable name="T" units="kelvin" initial_value="310" public_interface="out"/>
<variable name="F" units="coulomb_per_millimole" initial_value="96.4867" public_interface="out"/>
<variable name="Cm" units="picoF" initial_value="132" public_interface="out"/>
<variable name="i_ext" units="picoA" public_interface="out"/>
<variable name="i_tot" units="picoA"/>
<variable name="i_I" units="picoA"/>
<variable name="time" units="millisecond" public_interface="in"/>
<variable name="i_Na" units="picoA" public_interface="in"/>
<variable name="i_Ca_L" units="picoA" public_interface="in"/>
<variable name="i_Ca_T" units="picoA" public_interface="in"/>
<variable name="i_K1" units="picoA" public_interface="in"/>
<variable name="i_Kr" units="picoA" public_interface="in"/>
<variable name="i_Ks" units="picoA" public_interface="in"/>
<variable name="i_to" units="picoA" public_interface="in"/>
<variable name="i_NaK" units="picoA" public_interface="in"/>
<variable name="i_NaCa" units="picoA" public_interface="in"/>
<variable name="i_bNSC" units="picoA" public_interface="in"/>
<variable name="i_Cab" units="picoA" public_interface="in"/>
<variable name="i_Kpl" units="picoA" public_interface="in"/>
<variable name="i_lCa" units="picoA" public_interface="in"/>
<variable name="i_KATP" units="picoA" public_interface="in"/>
<variable name="stim_start" units="millisecond" initial_value="50"/>
<variable name="stim_end" units="millisecond" initial_value="1000000"/>
<variable name="stim_period" units="millisecond" initial_value="400"/>
<variable name="stim_duration" units="millisecond" initial_value="2"/>
<variable name="stim_amplitude" units="picoA" initial_value="-4000"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>i_ext</ci>
<piecewise>
<piece>
<ci>stim_amplitude</ci>
<apply>
<and/>
<apply>
<geq/>
<ci>time</ci>
<ci>stim_start</ci>
</apply>
<apply>
<leq/>
<ci>time</ci>
<ci>stim_end</ci>
</apply>
<apply>
<leq/>
<apply>
<minus/>
<apply>
<minus/>
<ci>time</ci>
<ci>stim_start</ci>
</apply>
<apply>
<times/>
<apply>
<floor/>
<apply>
<divide/>
<apply>
<minus/>
<ci>time</ci>
<ci>stim_start</ci>
</apply>
<ci>stim_period</ci>
</apply>
</apply>
<ci>stim_period</ci>
</apply>
</apply>
<ci>stim_duration</ci>
</apply>
</apply>
</piece>
<otherwise>
<cn cellml:units="picoA">0</cn>
</otherwise>
</piecewise>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>Vm</ci>
</apply>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<plus/>
<ci>i_tot</ci>
<ci>i_ext</ci>
</apply>
</apply>
<ci>Cm</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>i_tot</ci>
<apply>
<plus/>
<ci>i_Na</ci>
<ci>i_Ca_L</ci>
<ci>i_Ca_T</ci>
<ci>i_K1</ci>
<ci>i_Kr</ci>
<ci>i_Ks</ci>
<ci>i_to</ci>
<ci>i_I</ci>
<ci>i_NaK</ci>
<ci>i_NaCa</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>i_I</ci>
<apply>
<plus/>
<ci>i_bNSC</ci>
<ci>i_Cab</ci>
<ci>i_Kpl</ci>
<ci>i_lCa</ci>
<ci>i_KATP</ci>
</apply>
</apply>
</math>
</component>
<component name="external_ion_concentrations">
<variable name="Nao" units="millimolar" initial_value="140" public_interface="out"/>
<variable name="Cao" units="millimolar" initial_value="1.8" public_interface="out"/>
<variable name="Ko" units="millimolar" initial_value="5.4" public_interface="out"/>
</component>
<component name="internal_ion_concentrations">
<variable name="Nai" units="millimolar" initial_value="4.925761439682025" public_interface="out"/>
<variable name="Cai" units="millimolar" public_interface="out"/>
<variable name="Ki" units="millimolar" initial_value="143.1837333000449" public_interface="out"/>
<variable name="Vi" units="micrometre3" initial_value="8000" public_interface="out"/>
<variable name="i_net_Na" units="picoA"/>
<variable name="i_net_K" units="picoA"/>
<variable name="i_net_Ca" units="picoA"/>
<variable name="time" units="millisecond" public_interface="in"/>
<variable name="F" units="coulomb_per_millimole" public_interface="in"/>
<variable name="i_Na_Na" units="picoA" public_interface="in"/>
<variable name="i_CaL_Na" units="picoA" public_interface="in"/>
<variable name="i_bNSC_Na" units="picoA" public_interface="in"/>
<variable name="i_lCa_Na" units="picoA" public_interface="in"/>
<variable name="i_NaK" units="picoA" public_interface="in"/>
<variable name="i_K1" units="picoA" public_interface="in"/>
<variable name="i_Kr" units="picoA" public_interface="in"/>
<variable name="i_to_K" units="picoA" public_interface="in"/>
<variable name="i_to_Na" units="picoA" public_interface="in"/>
<variable name="i_Ks_K" units="picoA" public_interface="in"/>
<variable name="i_Ks_Na" units="picoA" public_interface="in"/>
<variable name="i_KATP" units="picoA" public_interface="in"/>
<variable name="i_Na_K" units="picoA" public_interface="in"/>
<variable name="i_CaL_K" units="picoA" public_interface="in"/>
<variable name="i_bNSC_K" units="picoA" public_interface="in"/>
<variable name="i_lCa_K" units="picoA" public_interface="in"/>
<variable name="i_Kpl" units="picoA" public_interface="in"/>
<variable name="i_CaL_Ca" units="picoA" public_interface="in"/>
<variable name="i_Ca_T" units="picoA" public_interface="in"/>
<variable name="i_NaCa" units="picoA" public_interface="in"/>
<variable name="i_bNSC" units="picoA" public_interface="in"/>
<variable name="i_Cab" units="picoA" public_interface="in"/>
<variable name="i_ext" units="picoA" public_interface="in"/>
<variable name="i_RyR" units="picoA" public_interface="in"/>
<variable name="i_SR_U" units="picoA" public_interface="in"/>
<variable name="i_SR_L" units="picoA" public_interface="in"/>
<variable name="dCaidt_NL" units="millimolar_per_millisecond" public_interface="in"/>
<variable name="CMDN_max" units="millimolar" initial_value="0.05"/>
<variable name="K_mCMDN" units="millimolar" initial_value="0.00238"/>
<variable name="Ca_Total" units="millimolar" initial_value="4.0180173572968586e-4"/>
<variable name="b1" units="millimolar"/>
<variable name="c1" units="millimolar2"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>b1</ci>
<apply>
<plus/>
<apply>
<minus/>
<ci>CMDN_max</ci>
<ci>Ca_Total</ci>
</apply>
<ci>K_mCMDN</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>c1</ci>
<apply>
<times/>
<ci>K_mCMDN</ci>
<ci>Ca_Total</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>Cai</ci>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<root/>
<apply>
<plus/>
<apply>
<power/>
<ci>b1</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">4</cn>
<ci>c1</ci>
</apply>
</apply>
</apply>
<ci>b1</ci>
</apply>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>Nai</ci>
</apply>
<apply>
<divide/>
<apply>
<minus/>
<ci>i_net_Na</ci>
</apply>
<apply>
<times/>
<ci>F</ci>
<ci>Vi</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>Ki</ci>
</apply>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<plus/>
<ci>i_net_K</ci>
<ci>i_ext</ci>
</apply>
</apply>
<apply>
<times/>
<ci>F</ci>
<ci>Vi</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>Ca_Total</ci>
</apply>
<apply>
<plus/>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<minus/>
<apply>
<minus/>
<apply>
<minus/>
<ci>i_net_Ca</ci>
<ci>i_SR_U</ci>
</apply>
<ci>i_RyR</ci>
</apply>
<ci>i_SR_L</ci>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">2</cn>
<ci>F</ci>
<ci>Vi</ci>
</apply>
</apply>
<ci>dCaidt_NL</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>i_net_Na</ci>
<apply>
<plus/>
<ci>i_Na_Na</ci>
<ci>i_Ks_Na</ci>
<ci>i_to_Na</ci>
<ci>i_CaL_Na</ci>
<ci>i_bNSC_Na</ci>
<ci>i_lCa_Na</ci>
<apply>
<times/>
<cn cellml:units="dimensionless">3</cn>
<ci>i_NaK</ci>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">3</cn>
<ci>i_NaCa</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>i_net_K</ci>
<apply>
<minus/>
<apply>
<plus/>
<ci>i_K1</ci>
<ci>i_Kr</ci>
<ci>i_to_K</ci>
<ci>i_KATP</ci>
<ci>i_Ks_K</ci>
<ci>i_Na_K</ci>
<ci>i_CaL_K</ci>
<ci>i_bNSC_K</ci>
<ci>i_lCa_K</ci>
<ci>i_Kpl</ci>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">2</cn>
<ci>i_NaK</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>i_net_Ca</ci>
<apply>
<minus/>
<apply>
<plus/>
<ci>i_CaL_Ca</ci>
<ci>i_Ca_T</ci>
<ci>i_Cab</ci>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">2</cn>
<ci>i_NaCa</ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="constant_field_equations">
<variable name="CF_Na" units="millimolar" public_interface="out"/>
<variable name="CF_Ca" units="millimolar" public_interface="out"/>
<variable name="CF_K" units="millimolar" public_interface="out"/>
<variable name="Nai" units="millimolar" public_interface="in"/>
<variable name="Nao" units="millimolar" public_interface="in"/>
<variable name="Cai" units="millimolar" public_interface="in"/>
<variable name="Cao" units="millimolar" public_interface="in"/>
<variable name="Ki" units="millimolar" public_interface="in"/>
<variable name="Ko" units="millimolar" public_interface="in"/>
<variable name="R" units="coulomb_millivolt_per_kelvin_millimole" public_interface="in"/>
<variable name="T" units="kelvin" public_interface="in"/>
<variable name="F" units="coulomb_per_millimole" public_interface="in"/>
<variable name="Vm" units="millivolt" public_interface="in"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>CF_Na</ci>
<piecewise>
<piece>
<apply>
<minus/>
<ci>Nao</ci>
</apply>
<apply>
<eq/>
<ci>Vm</ci>
<cn cellml:units="millivolt">0</cn>
</apply>
</piece>
<otherwise>
<apply>
<divide/>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<ci>F</ci>
<ci>Vm</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
<apply>
<minus/>
<ci>Nai</ci>
<apply>
<times/>
<ci>Nao</ci>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<ci>F</ci>
</apply>
<ci>Vm</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<ci>F</ci>
</apply>
<ci>Vm</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</otherwise>
</piecewise>
</apply>
<apply>
<eq/>
<ci>CF_Ca</ci>
<piecewise>
<piece>
<apply>
<minus/>
<ci>Cao</ci>
</apply>
<apply>
<eq/>
<ci>Vm</ci>
<cn cellml:units="millivolt">0</cn>
</apply>
</piece>
<otherwise>
<apply>
<divide/>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<cn cellml:units="dimensionless">2</cn>
<ci>F</ci>
<ci>Vm</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
<apply>
<minus/>
<ci>Cai</ci>
<apply>
<times/>
<ci>Cao</ci>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="dimensionless">2</cn>
</apply>
<ci>F</ci>
<ci>Vm</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="dimensionless">2</cn>
</apply>
<ci>F</ci>
<ci>Vm</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</otherwise>
</piecewise>
</apply>
<apply>
<eq/>
<ci>CF_K</ci>
<piecewise>
<piece>
<ci>Ki</ci>
<apply>
<eq/>
<ci>Vm</ci>
<cn cellml:units="millivolt">0</cn>
</apply>
</piece>
<otherwise>
<apply>
<divide/>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<ci>F</ci>
<ci>Vm</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
<apply>
<minus/>
<ci>Ki</ci>
<apply>
<times/>
<ci>Ko</ci>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<ci>F</ci>
</apply>
<ci>Vm</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<ci>F</ci>
</apply>
<ci>Vm</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</otherwise>
</piecewise>
</apply>
</math>
</component>
<component name="ATP_production">
<variable name="ATPi" units="millimolar" initial_value="4.657102729020499" public_interface="out"/>
<variable name="time" units="millisecond" public_interface="in"/>
<variable name="Vm" units="millivolt" public_interface="in"/>
<variable name="dATPdt_NL" units="millimolar_per_millisecond" public_interface="in"/>
<variable name="F" units="coulomb_per_millimole" public_interface="in"/>
<variable name="Vi" units="micrometre3" public_interface="in"/>
<variable name="i_NaK" units="picoA" public_interface="in"/>
<variable name="i_SR_U" units="picoA" public_interface="in"/>
<variable name="ProducingRate_Max" units="per_millisecond" initial_value="0.003"/>
<variable name="Adenosine_Total" units="millimolar" initial_value="5"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>ATPi</ci>
</apply>
<apply>
<plus/>
<apply>
<minus/>
<apply>
<plus/>
<apply>
<times/>
<ci>ProducingRate_Max</ci>
<apply>
<minus/>
<ci>Adenosine_Total</ci>
<ci>ATPi</ci>
</apply>
</apply>
<ci>dATPdt_NL</ci>
</apply>
<apply>
<divide/>
<ci>i_NaK</ci>
<apply>
<times/>
<ci>F</ci>
<ci>Vi</ci>
</apply>
</apply>
</apply>
<apply>
<divide/>
<ci>i_SR_U</ci>
<apply>
<times/>
<cn cellml:units="dimensionless">4</cn>
<ci>F</ci>
<ci>Vi</ci>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<component name="sodium_current">
<variable name="i_Na" units="picoA" public_interface="out"/>
<variable name="i_Na_Na" units="picoA" public_interface="out"/>
<variable name="i_Na_K" units="picoA" public_interface="out"/>
<variable name="P_Na" units="picoA_per_millimolar" initial_value="2860"/>
<variable name="time" units="millisecond" public_interface="in" private_interface="out"/>
<variable name="Vm" units="millivolt" public_interface="in" private_interface="out"/>
<variable name="CF_Na" units="millimolar" public_interface="in"/>
<variable name="CF_K" units="millimolar" public_interface="in"/>
<variable name="p_AP_Na" units="dimensionless" private_interface="in"/>
<variable name="y" units="dimensionless" private_interface="in"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>i_Na</ci>
<apply>
<plus/>
<ci>i_Na_Na</ci>
<ci>i_Na_K</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>i_Na_Na</ci>
<apply>
<times/>
<ci>P_Na</ci>
<ci>CF_Na</ci>
<ci>p_AP_Na</ci>
<ci>y</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>i_Na_K</ci>
<apply>
<times/>
<cn cellml:units="dimensionless">0.1</cn>
<ci>P_Na</ci>
<ci>CF_K</ci>
<ci>p_AP_Na</ci>
<ci>y</ci>
</apply>
</apply>
</math>
</component>
<component name="sodium_current_voltage_dependent_gate">
<variable name="p_AP_Na" units="dimensionless" initial_value="1.779648367445368e-5" public_interface="out"/>
<variable name="p_RI_Na" units="dimensionless"/>
<variable name="p_RP_Na" units="dimensionless" initial_value="0.3556412697995689"/>
<variable name="p_AI_Na" units="dimensionless" initial_value="0.40285968661346977"/>
<variable name="k_RP_AP" units="per_millisecond"/>
<variable name="k_AP_RP" units="per_millisecond"/>
<variable name="k_RI_AI" units="per_millisecond"/>
<variable name="k_AI_RI" units="per_millisecond"/>
<variable name="k_AP_AI" units="per_millisecond"/>
<variable name="k_AI_AP" units="per_millisecond" initial_value="0.0000875"/>
<variable name="k_RP_RI" units="per_millisecond"/>
<variable name="k_RI_RP" units="per_millisecond"/>
<variable name="time" units="millisecond" public_interface="in"/>
<variable name="Vm" units="millivolt" public_interface="in"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>p_RP_Na</ci>
</apply>
<apply>
<minus/>
<apply>
<plus/>
<apply>
<times/>
<ci>p_AP_Na</ci>
<ci>k_AP_RP</ci>
</apply>
<apply>
<times/>
<ci>p_RI_Na</ci>
<ci>k_RI_RP</ci>
</apply>
</apply>
<apply>
<times/>
<ci>p_RP_Na</ci>
<apply>
<plus/>
<ci>k_RP_RI</ci>
<ci>k_RP_AP</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>p_AP_Na</ci>
</apply>
<apply>
<minus/>
<apply>
<plus/>
<apply>
<times/>
<ci>p_RP_Na</ci>
<ci>k_RP_AP</ci>
</apply>
<apply>
<times/>
<ci>p_AI_Na</ci>
<ci>k_AI_AP</ci>
</apply>
</apply>
<apply>
<times/>
<ci>p_AP_Na</ci>
<apply>
<plus/>
<ci>k_AP_RP</ci>
<ci>k_AP_AI</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>p_AI_Na</ci>
</apply>
<apply>
<minus/>
<apply>
<plus/>
<apply>
<times/>
<ci>p_RI_Na</ci>
<ci>k_RI_AI</ci>
</apply>
<apply>
<times/>
<ci>p_AP_Na</ci>
<ci>k_AP_AI</ci>
</apply>
</apply>
<apply>
<times/>
<ci>p_AI_Na</ci>
<apply>
<plus/>
<ci>k_AI_RI</ci>
<ci>k_AI_AP</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>p_RI_Na</ci>
<apply>
<minus/>
<apply>
<minus/>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<ci>p_RP_Na</ci>
</apply>
<ci>p_AP_Na</ci>
</apply>
<ci>p_AI_Na</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>k_RP_AP</ci>
<apply>
<divide/>
<cn cellml:units="per_millisecond">1</cn>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="dimensionless">0.1027</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<ci>Vm</ci>
</apply>
<cn cellml:units="millivolt">8</cn>
</apply>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">0.25</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<ci>Vm</ci>
</apply>
<cn cellml:units="millivolt">50</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>k_AP_RP</ci>
<apply>
<divide/>
<cn cellml:units="per_millisecond">1</cn>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="dimensionless">26</cn>
<apply>
<exp/>
<apply>
<divide/>
<ci>Vm</ci>
<cn cellml:units="millivolt">17</cn>
</apply>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">0.02</cn>
<apply>
<exp/>
<apply>
<divide/>
<ci>Vm</ci>
<cn cellml:units="millivolt">800</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>k_AP_AI</ci>
<apply>
<divide/>
<cn cellml:units="per_millisecond">1</cn>
<apply>
<times/>
<cn cellml:units="dimensionless">0.8</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<ci>Vm</ci>
</apply>
<cn cellml:units="millivolt">400</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>k_AI_RI</ci>
<apply>
<divide/>
<cn cellml:units="per_millisecond">1</cn>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="dimensionless">1300</cn>
<apply>
<exp/>
<apply>
<divide/>
<ci>Vm</ci>
<cn cellml:units="millivolt">20</cn>
</apply>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">0.04</cn>
<apply>
<exp/>
<apply>
<divide/>
<ci>Vm</ci>
<cn cellml:units="millivolt">800</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>k_RI_AI</ci>
<apply>
<divide/>
<cn cellml:units="per_millisecond">1</cn>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="dimensionless">0.0001027</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<ci>Vm</ci>
</apply>
<cn cellml:units="millivolt">8</cn>
</apply>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">5</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<ci>Vm</ci>
</apply>
<cn cellml:units="millivolt">400</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>k_RP_RI</ci>
<apply>
<divide/>
<cn cellml:units="per_millisecond">0.01</cn>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<divide/>
<apply>
<times/>
<ci>k_AI_AP</ci>
<ci>k_AP_RP</ci>
<ci>k_RI_AI</ci>
</apply>
<apply>
<times/>
<ci>k_AP_AI</ci>
<ci>k_RP_AP</ci>
<ci>k_AI_RI</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>k_RI_RP</ci>
<apply>
<minus/>
<cn cellml:units="per_millisecond">0.01</cn>
<ci>k_RP_RI</ci>
</apply>
</apply>
</math>
</component>
<component name="sodium_current_ultra_slow_gate">
<variable name="y" units="dimensionless" initial_value="0.5861887862983165" public_interface="out"/>
<variable name="alpha_y" units="per_millisecond"/>
<variable name="beta_y" units="per_millisecond"/>
<variable name="Vm" units="millivolt" public_interface="in"/>
<variable name="time" units="millisecond" public_interface="in"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>alpha_y</ci>
<apply>
<divide/>
<cn cellml:units="per_millisecond">1</cn>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="dimensionless">9000000000</cn>
<apply>
<exp/>
<apply>
<divide/>
<ci>Vm</ci>
<cn cellml:units="millivolt">5</cn>
</apply>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">8000</cn>
<apply>
<exp/>
<apply>
<divide/>
<ci>Vm</ci>
<cn cellml:units="millivolt">100</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>beta_y</ci>
<apply>
<divide/>
<cn cellml:units="per_millisecond">1</cn>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="dimensionless">0.014</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<ci>Vm</ci>
</apply>
<cn cellml:units="millivolt">5</cn>
</apply>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">4000</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<ci>Vm</ci>
</apply>
<cn cellml:units="millivolt">100</cn>
</apply>
</apply>
</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="L_type_Ca_channel">
<variable name="i_Ca_L" units="picoA" public_interface="out"/>
<variable name="i_CaL_Na" units="picoA" public_interface="out"/>
<variable name="i_CaL_Ca" units="picoA" public_interface="out"/>
<variable name="i_CaL_K" units="picoA" public_interface="out"/>
<variable name="p_open_CaL" units="dimensionless" public_interface="out" private_interface="out"/>
<variable name="CaDiadic" units="picoA" public_interface="out" private_interface="in"/>
<variable name="P_CaL" units="picoA_per_millimolar" initial_value="8712"/>
<variable name="time" units="millisecond" public_interface="in" private_interface="out"/>
<variable name="Vm" units="millivolt" public_interface="in" private_interface="out"/>
<variable name="CF_Na" units="millimolar" public_interface="in"/>
<variable name="CF_Ca" units="millimolar" public_interface="in" private_interface="out"/>
<variable name="Cai" units="millimolar" public_interface="in" private_interface="out"/>
<variable name="CF_K" units="millimolar" public_interface="in"/>
<variable name="ATPi" units="millimolar" public_interface="in"/>
<variable name="p_AP_CaL" units="dimensionless" private_interface="in"/>
<variable name="p_U" units="dimensionless" private_interface="in"/>
<variable name="p_UCa" units="dimensionless" private_interface="in"/>
<variable name="y" units="dimensionless" private_interface="in"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>i_Ca_L</ci>
<apply>
<plus/>
<ci>i_CaL_Na</ci>
<ci>i_CaL_Ca</ci>
<ci>i_CaL_K</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>i_CaL_Ca</ci>
<apply>
<times/>
<ci>P_CaL</ci>
<ci>CF_Ca</ci>
<ci>p_open_CaL</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>i_CaL_Na</ci>
<apply>
<times/>
<cn cellml:units="dimensionless">0.0000185</cn>
<ci>P_CaL</ci>
<ci>CF_Na</ci>
<ci>p_open_CaL</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>i_CaL_K</ci>
<apply>
<times/>
<cn cellml:units="dimensionless">0.000365</cn>
<ci>P_CaL</ci>
<ci>CF_K</ci>
<ci>p_open_CaL</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>p_open_CaL</ci>
<apply>
<divide/>
<apply>
<times/>
<ci>p_AP_CaL</ci>
<apply>
<plus/>
<ci>p_U</ci>
<ci>p_UCa</ci>
</apply>
<ci>y</ci>
</apply>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<power/>
<apply>
<divide/>
<cn cellml:units="millimolar">1.4</cn>
<ci>ATPi</ci>
</apply>
<cn cellml:units="dimensionless">3</cn>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<component name="L_type_Ca_channel_voltage_dependent_gate">
<variable name="p_AP_CaL" units="dimensionless" initial_value="1.5445004166497696e-6" public_interface="out"/>
<variable name="p_RI_CaL" units="dimensionless"/>
<variable name="p_RP_CaL" units="dimensionless" initial_value="0.9968480629364956"/>
<variable name="p_AI_CaL" units="dimensionless" initial_value="8.77325391245903e-4"/>
<variable name="k_RP_AP" units="per_millisecond"/>
<variable name="k_AP_RP" units="per_millisecond"/>
<variable name="k_RI_AI" units="per_millisecond"/>
<variable name="k_AI_RI" units="per_millisecond"/>
<variable name="k_AP_AI" units="per_millisecond" initial_value="0.004"/>
<variable name="k_AI_AP" units="per_millisecond" initial_value="0.001"/>
<variable name="k_RP_RI" units="per_millisecond"/>
<variable name="k_RI_RP" units="per_millisecond"/>
<variable name="time" units="millisecond" public_interface="in"/>
<variable name="Vm" units="millivolt" public_interface="in"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>p_RP_CaL</ci>
</apply>
<apply>
<minus/>
<apply>
<plus/>
<apply>
<times/>
<ci>p_AP_CaL</ci>
<ci>k_AP_RP</ci>
</apply>
<apply>
<times/>
<ci>p_RI_CaL</ci>
<ci>k_RI_RP</ci>
</apply>
</apply>
<apply>
<times/>
<ci>p_RP_CaL</ci>
<apply>
<plus/>
<ci>k_RP_RI</ci>
<ci>k_RP_AP</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>p_AP_CaL</ci>
</apply>
<apply>
<minus/>
<apply>
<plus/>
<apply>
<times/>
<ci>p_RP_CaL</ci>
<ci>k_RP_AP</ci>
</apply>
<apply>
<times/>
<ci>p_AI_CaL</ci>
<ci>k_AI_AP</ci>
</apply>
</apply>
<apply>
<times/>
<ci>p_AP_CaL</ci>
<apply>
<plus/>
<ci>k_AP_RP</ci>
<ci>k_AP_AI</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>p_AI_CaL</ci>
</apply>
<apply>
<minus/>
<apply>
<plus/>
<apply>
<times/>
<ci>p_RI_CaL</ci>
<ci>k_RI_AI</ci>
</apply>
<apply>
<times/>
<ci>p_AP_CaL</ci>
<ci>k_AP_AI</ci>
</apply>
</apply>
<apply>
<times/>
<ci>p_AI_CaL</ci>
<apply>
<plus/>
<ci>k_AI_RI</ci>
<ci>k_AI_AP</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>p_RI_CaL</ci>
<apply>
<minus/>
<apply>
<minus/>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<ci>p_AP_CaL</ci>
</apply>
<ci>p_RP_CaL</ci>
</apply>
<ci>p_AI_CaL</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>k_RP_AP</ci>
<apply>
<divide/>
<cn cellml:units="per_millisecond">1</cn>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="dimensionless">0.27</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<ci>Vm</ci>
</apply>
<cn cellml:units="millivolt">5.9</cn>
</apply>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">1.5</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<ci>Vm</ci>
</apply>
<cn cellml:units="millivolt">65</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>k_AP_RP</ci>
<apply>
<divide/>
<cn cellml:units="per_millisecond">1</cn>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="dimensionless">480</cn>
<apply>
<exp/>
<apply>
<divide/>
<ci>Vm</ci>
<cn cellml:units="millivolt">7</cn>
</apply>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">2.2</cn>
<apply>
<exp/>
<apply>
<divide/>
<ci>Vm</ci>
<cn cellml:units="millivolt">65</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>k_RI_AI</ci>
<apply>
<divide/>
<cn cellml:units="per_millisecond">1</cn>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="dimensionless">0.0018</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<ci>Vm</ci>
</apply>
<cn cellml:units="millivolt">7.4</cn>
</apply>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">2</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<ci>Vm</ci>
</apply>
<cn cellml:units="millivolt">100</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>k_AI_RI</ci>
<apply>
<divide/>
<cn cellml:units="per_millisecond">1</cn>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="dimensionless">2200000</cn>
<apply>
<exp/>
<apply>
<divide/>
<ci>Vm</ci>
<cn cellml:units="millivolt">7.4</cn>
</apply>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">11</cn>
<apply>
<exp/>
<apply>
<divide/>
<ci>Vm</ci>
<cn cellml:units="millivolt">100</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>k_RP_RI</ci>
<apply>
<divide/>
<cn cellml:units="per_millisecond">0.04</cn>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<divide/>
<apply>
<times/>
<ci>k_AI_AP</ci>
<ci>k_AP_RP</ci>
<ci>k_RI_AI</ci>
</apply>
<apply>
<times/>
<ci>k_AP_AI</ci>
<ci>k_RP_AP</ci>
<ci>k_AI_RI</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>k_RI_RP</ci>
<apply>
<minus/>
<cn cellml:units="per_millisecond">0.04</cn>
<ci>k_RP_RI</ci>
</apply>
</apply>
</math>
</component>
<component name="L_type_Ca_channel_Ca_dependent_gate">
<variable name="p_U" units="dimensionless" initial_value="0.17246483915629204" public_interface="out"/>
<variable name="p_UCa" units="dimensionless" initial_value="6.098246017787626e-5" public_interface="out"/>
<variable name="iCaL" units="picoA"/>
<variable name="Cacm" units="millimolar"/>
<variable name="p_CCa" units="dimensionless"/>
<variable name="p_C" units="dimensionless" initial_value="0.4250747299372254"/>
<variable name="k_CCa_UCa" units="per_millisecond" initial_value="0.0003"/>
<variable name="k_UCa_CCa" units="per_millisecond" initial_value="0.35"/>
<variable name="k_C_U" units="per_millisecond" initial_value="0.143"/>
<variable name="k_U_C" units="per_millisecond" initial_value="0.35"/>
<variable name="k_UCa_U" units="per_millisecond"/>
<variable name="k_U_UCa" units="per_millimolar_millisecond" initial_value="6.954"/>
<variable name="k_CCa_C" units="per_millisecond" initial_value="0.0042"/>
<variable name="k_C_CCa" units="per_millimolar_millisecond" initial_value="6.954"/>
<variable name="time" units="millisecond" public_interface="in"/>
<variable name="Vm" units="millivolt" public_interface="in"/>
<variable name="Cai" units="millimolar" public_interface="in"/>
<variable name="CF_Ca" units="millimolar" public_interface="in"/>
<variable name="p_AP_CaL" units="dimensionless" public_interface="in"/>
<variable name="CaEffC" units="millimolar"/>
<variable name="CaEffU" units="millimolar"/>
<variable name="k_UUCa_Ca" units="per_millisecond"/>
<variable name="k_CCCa_Ca" units="per_millisecond"/>
<variable name="CaDiadic" units="picoA" public_interface="out"/>
<variable name="p_open_CaL" units="dimensionless" public_interface="in"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>iCaL</ci>
<apply>
<times/>
<cn cellml:units="picoA_per_millimolar">0.0676</cn>
<ci>CF_Ca</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>CaDiadic</ci>
<apply>
<times/>
<ci>iCaL</ci>
<ci>p_open_CaL</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>Cacm</ci>
<apply>
<minus/>
<ci>Cai</ci>
<apply>
<times/>
<cn cellml:units="millimolar_per_picoA">0.3</cn>
<ci>iCaL</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>CaEffC</ci>
<apply>
<times/>
<ci>Cacm</ci>
<ci>p_AP_CaL</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>CaEffU</ci>
<apply>
<plus/>
<ci>CaEffC</ci>
<apply>
<times/>
<ci>Cai</ci>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<ci>p_AP_CaL</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>k_UUCa_Ca</ci>
<apply>
<times/>
<ci>k_U_UCa</ci>
<ci>CaEffU</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>k_CCCa_Ca</ci>
<apply>
<times/>
<ci>k_C_CCa</ci>
<ci>CaEffC</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>p_U</ci>
</apply>
<apply>
<minus/>
<apply>
<plus/>
<apply>
<times/>
<ci>p_C</ci>
<ci>k_C_U</ci>
</apply>
<apply>
<times/>
<ci>p_UCa</ci>
<ci>k_UCa_U</ci>
</apply>
</apply>
<apply>
<times/>
<ci>p_U</ci>
<apply>
<plus/>
<ci>k_UUCa_Ca</ci>
<ci>k_U_C</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>p_UCa</ci>
</apply>
<apply>
<minus/>
<apply>
<plus/>
<apply>
<times/>
<ci>p_U</ci>
<ci>k_UUCa_Ca</ci>
</apply>
<apply>
<times/>
<ci>p_CCa</ci>
<ci>k_CCa_UCa</ci>
</apply>
</apply>
<apply>
<times/>
<ci>p_UCa</ci>
<apply>
<plus/>
<ci>k_UCa_CCa</ci>
<ci>k_UCa_U</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>p_C</ci>
</apply>
<apply>
<minus/>
<apply>
<plus/>
<apply>
<times/>
<ci>p_CCa</ci>
<ci>k_CCa_C</ci>
</apply>
<apply>
<times/>
<ci>p_U</ci>
<ci>k_U_C</ci>
</apply>
</apply>
<apply>
<times/>
<ci>p_C</ci>
<apply>
<plus/>
<ci>k_C_U</ci>
<apply>
<times/>
<ci>k_C_CCa</ci>
<ci>Cacm</ci>
<ci>p_AP_CaL</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>p_CCa</ci>
<apply>
<minus/>
<apply>
<minus/>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<ci>p_C</ci>
</apply>
<ci>p_U</ci>
</apply>
<ci>p_UCa</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>k_UCa_U</ci>
<apply>
<divide/>
<apply>
<times/>
<ci>k_CCa_C</ci>
<ci>k_C_U</ci>
<ci>k_U_UCa</ci>
<ci>k_UCa_CCa</ci>
</apply>
<apply>
<times/>
<ci>k_U_C</ci>
<ci>k_C_CCa</ci>
<ci>k_CCa_UCa</ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="L_type_Ca_channel_ultra_slow_gate">
<variable name="y" units="dimensionless" initial_value="0.9985266538252986" public_interface="out"/>
<variable name="alpha_y" units="per_millisecond"/>
<variable name="beta_y" units="per_millisecond"/>
<variable name="Vm" units="millivolt" public_interface="in"/>
<variable name="time" units="millisecond" public_interface="in"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>alpha_y</ci>
<apply>
<divide/>
<cn cellml:units="per_millisecond">1</cn>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="dimensionless">250000</cn>
<apply>
<exp/>
<apply>
<divide/>
<ci>Vm</ci>
<cn cellml:units="millivolt">9</cn>
</apply>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">58</cn>
<apply>
<exp/>
<apply>
<divide/>
<ci>Vm</ci>
<cn cellml:units="millivolt">65</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>beta_y</ci>
<apply>
<divide/>
<cn cellml:units="per_millisecond">1</cn>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="dimensionless">1800</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<ci>Vm</ci>
</apply>
<cn cellml:units="millivolt">14</cn>
</apply>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">66</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<ci>Vm</ci>
</apply>
<cn cellml:units="millivolt">65</cn>
</apply>
</apply>
</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="T_type_Ca_channel">
<variable name="i_Ca_T" units="picoA" public_interface="out"/>
<variable name="P_CaT" units="picoA_per_millimolar" initial_value="612"/>
<variable name="time" units="millisecond" public_interface="in" private_interface="out"/>
<variable name="Vm" units="millivolt" public_interface="in" private_interface="out"/>
<variable name="CF_Ca" units="millimolar" public_interface="in"/>
<variable name="y1" units="dimensionless" private_interface="in"/>
<variable name="y2" units="dimensionless" private_interface="in"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>i_Ca_T</ci>
<apply>
<times/>
<ci>P_CaT</ci>
<ci>CF_Ca</ci>
<ci>y1</ci>
<ci>y2</ci>
</apply>
</apply>
</math>
</component>
<component name="T_type_Ca_channel_y1_gate">
<variable name="y1" units="dimensionless" initial_value="1.6882718240109127e-5" public_interface="out"/>
<variable name="alpha_y1" units="per_millisecond"/>
<variable name="beta_y1" units="per_millisecond"/>
<variable name="Vm" units="millivolt" public_interface="in"/>
<variable name="time" units="millisecond" public_interface="in"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>alpha_y1</ci>
<apply>
<divide/>
<cn cellml:units="per_millisecond">1</cn>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="dimensionless">0.019</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<ci>Vm</ci>
</apply>
<cn cellml:units="millivolt">5.6</cn>
</apply>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">0.82</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<ci>Vm</ci>
</apply>
<cn cellml:units="millivolt">250</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>beta_y1</ci>
<apply>
<divide/>
<cn cellml:units="per_millisecond">1</cn>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="dimensionless">40</cn>
<apply>
<exp/>
<apply>
<divide/>
<ci>Vm</ci>
<cn cellml:units="millivolt">6.3</cn>
</apply>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">1.5</cn>
<apply>
<exp/>
<apply>
<divide/>
<ci>Vm</ci>
<cn cellml:units="millivolt">10000</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>y1</ci>
</apply>
<apply>
<minus/>
<apply>
<times/>
<ci>alpha_y1</ci>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<ci>y1</ci>
</apply>
</apply>
<apply>
<times/>
<ci>beta_y1</ci>
<ci>y1</ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="T_type_Ca_channel_y2_gate">
<variable name="y2" units="dimensionless" initial_value="0.8585352091865849" public_interface="out"/>
<variable name="alpha_y2" units="per_millisecond"/>
<variable name="beta_y2" units="per_millisecond"/>
<variable name="Vm" units="millivolt" public_interface="in"/>
<variable name="time" units="millisecond" public_interface="in"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>alpha_y2</ci>
<apply>
<divide/>
<cn cellml:units="per_millisecond">1</cn>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="dimensionless">62000</cn>
<apply>
<exp/>
<apply>
<divide/>
<ci>Vm</ci>
<cn cellml:units="millivolt">10.1</cn>
</apply>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">30</cn>
<apply>
<exp/>
<apply>
<divide/>
<ci>Vm</ci>
<cn cellml:units="millivolt">3000</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>beta_y2</ci>
<apply>
<divide/>
<cn cellml:units="per_millisecond">1</cn>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="dimensionless">0.0006</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<ci>Vm</ci>
</apply>
<cn cellml:units="millivolt">6.7</cn>
</apply>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">1.2</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<ci>Vm</ci>
</apply>
<cn cellml:units="millivolt">25</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>y2</ci>
</apply>
<apply>
<minus/>
<apply>
<times/>
<ci>alpha_y2</ci>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<ci>y2</ci>
</apply>
</apply>
<apply>
<times/>
<ci>beta_y2</ci>
<ci>y2</ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="time_independent_potassium_current">
<variable name="i_K1" units="picoA" public_interface="out"/>
<variable name="E_K" units="millivolt" public_interface="out" private_interface="out"/>
<variable name="g_K1" units="nanoS"/>
<variable name="P_K1_0" units="nanoS_per_picoF" initial_value="1.146"/>
<variable name="fO" units="dimensionless" private_interface="out"/>
<variable name="fO2" units="dimensionless"/>
<variable name="fO3" units="dimensionless"/>
<variable name="fO4" units="dimensionless"/>
<variable name="fB" units="dimensionless"/>
<variable name="mu" units="per_millisecond"/>
<variable name="lambda" units="per_millisecond"/>
<variable name="time" units="millisecond" public_interface="in" private_interface="out"/>
<variable name="Vm" units="millivolt" public_interface="in" private_interface="out"/>
<variable name="Cm" units="picoF" public_interface="in"/>
<variable name="Ko" units="millimolar" public_interface="in"/>
<variable name="Ki" units="millimolar" public_interface="in"/>
<variable name="R" units="coulomb_millivolt_per_kelvin_millimole" public_interface="in"/>
<variable name="T" units="kelvin" public_interface="in"/>
<variable name="F" units="coulomb_per_millimole" public_interface="in"/>
<variable name="y" units="dimensionless" private_interface="in"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>i_K1</ci>
<apply>
<times/>
<ci>g_K1</ci>
<apply>
<minus/>
<ci>Vm</ci>
<ci>E_K</ci>
</apply>
<apply>
<plus/>
<ci>fO4</ci>
<ci>fO3</ci>
<ci>fO2</ci>
</apply>
<ci>y</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>E_K</ci>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
<ci>F</ci>
</apply>
<apply>
<ln/>
<apply>
<divide/>
<ci>Ko</ci>
<ci>Ki</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>g_K1</ci>
<apply>
<times/>
<ci>P_K1_0</ci>
<ci>Cm</ci>
<apply>
<power/>
<apply>
<divide/>
<ci>Ko</ci>
<cn cellml:units="millimolar">5.4</cn>
</apply>
<cn cellml:units="dimensionless">0.4</cn>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>fB</ci>
<apply>
<divide/>
<ci>mu</ci>
<apply>
<plus/>
<ci>mu</ci>
<ci>lambda</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>fO</ci>
<apply>
<divide/>
<ci>lambda</ci>
<apply>
<plus/>
<ci>mu</ci>
<ci>lambda</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>fO2</ci>
<apply>
<times/>
<cn cellml:units="dimensionless">2</cn>
<apply>
<power/>
<ci>fO</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
<apply>
<power/>
<ci>fB</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>fO3</ci>
<apply>
<times/>
<apply>
<divide/>
<cn cellml:units="dimensionless">8</cn>
<cn cellml:units="dimensionless">3</cn>
</apply>
<apply>
<power/>
<ci>fO</ci>
<cn cellml:units="dimensionless">3</cn>
</apply>
<ci>fB</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>fO4</ci>
<apply>
<power/>
<ci>fO</ci>
<cn cellml:units="dimensionless">4</cn>
</apply>
</apply>
<apply>
<eq/>
<ci>mu</ci>
<apply>
<divide/>
<apply>
<times/>
<cn cellml:units="per_millisecond">0.75</cn>
<apply>
<exp/>
<apply>
<times/>
<cn cellml:units="per_millivolt">0.035</cn>
<apply>
<minus/>
<apply>
<minus/>
<ci>Vm</ci>
<ci>E_K</ci>
</apply>
<cn cellml:units="millivolt">10</cn>
</apply>
</apply>
</apply>
</apply>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<times/>
<cn cellml:units="per_millivolt">0.015</cn>
<apply>
<minus/>
<apply>
<minus/>
<ci>Vm</ci>
<ci>E_K</ci>
</apply>
<cn cellml:units="millivolt">140</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>lambda</ci>
<apply>
<divide/>
<apply>
<times/>
<cn cellml:units="per_millisecond">3</cn>
<apply>
<exp/>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="per_millivolt">0.048</cn>
</apply>
<apply>
<minus/>
<apply>
<minus/>
<ci>Vm</ci>
<ci>E_K</ci>
</apply>
<cn cellml:units="millivolt">10</cn>
</apply>
</apply>
</apply>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<times/>
<cn cellml:units="per_millivolt">0.064</cn>
<apply>
<minus/>
<apply>
<minus/>
<ci>Vm</ci>
<ci>E_K</ci>
</apply>
<cn cellml:units="millivolt">38</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<times/>
<cn cellml:units="per_millivolt">0.03</cn>
<apply>
<minus/>
<apply>
<minus/>
<ci>Vm</ci>
<ci>E_K</ci>
</apply>
<cn cellml:units="millivolt">70</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<component name="time_independent_potassium_current_y_gate">
<variable name="y" units="dimensionless" initial_value="0.6080573900752752" public_interface="out"/>
<variable name="alpha_y" units="per_millisecond"/>
<variable name="beta_y" units="per_millisecond"/>
<variable name="fO" units="dimensionless" public_interface="in"/>
<variable name="Vm" units="millivolt" public_interface="in"/>
<variable name="E_K" units="millivolt" public_interface="in"/>
<variable name="time" units="millisecond" public_interface="in"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>alpha_y</ci>
<apply>
<divide/>
<cn cellml:units="per_millisecond">1</cn>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="dimensionless">8000</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<minus/>
<ci>Vm</ci>
<ci>E_K</ci>
</apply>
<cn cellml:units="millivolt">97</cn>
</apply>
<cn cellml:units="millivolt">8.5</cn>
</apply>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">7</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<minus/>
<ci>Vm</ci>
<ci>E_K</ci>
</apply>
<cn cellml:units="millivolt">97</cn>
</apply>
<cn cellml:units="millivolt">300</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>beta_y</ci>
<apply>
<divide/>
<apply>
<times/>
<apply>
<power/>
<ci>fO</ci>
<cn cellml:units="dimensionless">4</cn>
</apply>
<cn cellml:units="per_millisecond">1</cn>
</apply>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="dimensionless">0.00014</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<minus/>
<apply>
<minus/>
<ci>Vm</ci>
<ci>E_K</ci>
</apply>
<cn cellml:units="millivolt">97</cn>
</apply>
</apply>
<cn cellml:units="millivolt">9.1</cn>
</apply>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">0.2</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<minus/>
<apply>
<minus/>
<ci>Vm</ci>
<ci>E_K</ci>
</apply>
<cn cellml:units="millivolt">97</cn>
</apply>
</apply>
<cn cellml:units="millivolt">500</cn>
</apply>
</apply>
</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="rapid_time_dependent_potassium_current">
<variable name="i_Kr" units="picoA" public_interface="out"/>
<variable name="g_Kr" units="nanoS"/>
<variable name="P_Kr" units="nanoS_per_picoF" initial_value="0.00864"/>
<variable name="time" units="millisecond" public_interface="in" private_interface="out"/>
<variable name="Vm" units="millivolt" public_interface="in" private_interface="out"/>
<variable name="Cm" units="picoF" public_interface="in"/>
<variable name="E_K" units="millivolt" public_interface="in" private_interface="out"/>
<variable name="Ko" units="millimolar" public_interface="in"/>
<variable name="y1" units="dimensionless" private_interface="in"/>
<variable name="y2" units="dimensionless" private_interface="in"/>
<variable name="y3" units="dimensionless" private_interface="in"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>i_Kr</ci>
<apply>
<times/>
<ci>g_Kr</ci>
<apply>
<minus/>
<ci>Vm</ci>
<ci>E_K</ci>
</apply>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="dimensionless">0.6</cn>
<ci>y1</ci>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">0.4</cn>
<ci>y2</ci>
</apply>
</apply>
<ci>y3</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>g_Kr</ci>
<apply>
<times/>
<ci>P_Kr</ci>
<ci>Cm</ci>
<apply>
<power/>
<apply>
<divide/>
<ci>Ko</ci>
<cn cellml:units="millimolar">5.4</cn>
</apply>
<cn cellml:units="dimensionless">0.2</cn>
</apply>
</apply>
</apply>
</math>
</component>
<component name="rapid_time_dependent_potassium_current_y1_gate">
<variable name="y1" units="dimensionless" initial_value="0.0018339931180983765" public_interface="out"/>
<variable name="alpha_y1" units="per_millisecond"/>
<variable name="beta_y1" units="per_millisecond"/>
<variable name="Vm" units="millivolt" public_interface="in"/>
<variable name="E_K" units="millivolt" public_interface="in"/>
<variable name="time" units="millisecond" public_interface="in"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>alpha_y1</ci>
<apply>
<divide/>
<cn cellml:units="per_millisecond">1</cn>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="dimensionless">20</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<ci>Vm</ci>
</apply>
<cn cellml:units="millivolt">11.5</cn>
</apply>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">5</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<ci>Vm</ci>
</apply>
<cn cellml:units="millivolt">300</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>beta_y1</ci>
<apply>
<plus/>
<apply>
<divide/>
<cn cellml:units="per_millisecond">1</cn>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="dimensionless">160</cn>
<apply>
<exp/>
<apply>
<divide/>
<ci>Vm</ci>
<cn cellml:units="millivolt">28</cn>
</apply>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">200</cn>
<apply>
<exp/>
<apply>
<divide/>
<ci>Vm</ci>
<cn cellml:units="millivolt">1000</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<divide/>
<cn cellml:units="per_millisecond">1</cn>
<apply>
<times/>
<cn cellml:units="dimensionless">2500</cn>
<apply>
<exp/>
<apply>
<divide/>
<ci>Vm</ci>
<cn cellml:units="millivolt">20</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>y1</ci>
</apply>
<apply>
<minus/>
<apply>
<times/>
<ci>alpha_y1</ci>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<ci>y1</ci>
</apply>
</apply>
<apply>
<times/>
<ci>beta_y1</ci>
<ci>y1</ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="rapid_time_dependent_potassium_current_y2_gate">
<variable name="y2" units="dimensionless" initial_value="0.20443083454225305" public_interface="out"/>
<variable name="alpha_y2" units="per_millisecond"/>
<variable name="beta_y2" units="per_millisecond"/>
<variable name="Vm" units="millivolt" public_interface="in"/>
<variable name="E_K" units="millivolt" public_interface="in"/>
<variable name="time" units="millisecond" public_interface="in"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>alpha_y2</ci>
<apply>
<divide/>
<cn cellml:units="per_millisecond">1</cn>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="dimensionless">200</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<ci>Vm</ci>
</apply>
<cn cellml:units="millivolt">13</cn>
</apply>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">20</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<ci>Vm</ci>
</apply>
<cn cellml:units="millivolt">300</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>beta_y2</ci>
<apply>
<plus/>
<apply>
<divide/>
<cn cellml:units="per_millisecond">1</cn>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="dimensionless">1600</cn>
<apply>
<exp/>
<apply>
<divide/>
<ci>Vm</ci>
<cn cellml:units="millivolt">28</cn>
</apply>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">2000</cn>
<apply>
<exp/>
<apply>
<divide/>
<ci>Vm</ci>
<cn cellml:units="millivolt">1000</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<divide/>
<cn cellml:units="per_millisecond">1</cn>
<apply>
<times/>
<cn cellml:units="dimensionless">10000</cn>
<apply>
<exp/>
<apply>
<divide/>
<ci>Vm</ci>
<cn cellml:units="millivolt">20</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>y2</ci>
</apply>
<apply>
<minus/>
<apply>
<times/>
<ci>alpha_y2</ci>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<ci>y2</ci>
</apply>
</apply>
<apply>
<times/>
<ci>beta_y2</ci>
<ci>y2</ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="rapid_time_dependent_potassium_current_y3_gate">
<variable name="y3" units="dimensionless" initial_value="0.967887666264921" public_interface="out"/>
<variable name="alpha_y3" units="per_millisecond"/>
<variable name="beta_y3" units="per_millisecond"/>
<variable name="Vm" units="millivolt" public_interface="in"/>
<variable name="E_K" units="millivolt" public_interface="in"/>
<variable name="time" units="millisecond" public_interface="in"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>alpha_y3</ci>
<apply>
<divide/>
<cn cellml:units="per_millisecond">1</cn>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="dimensionless">10</cn>
<apply>
<exp/>
<apply>
<divide/>
<ci>Vm</ci>
<cn cellml:units="millivolt">17</cn>
</apply>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">2.5</cn>
<apply>
<exp/>
<apply>
<divide/>
<ci>Vm</ci>
<cn cellml:units="millivolt">300</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>beta_y3</ci>
<apply>
<divide/>
<cn cellml:units="per_millisecond">1</cn>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="dimensionless">0.35</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<ci>Vm</ci>
</apply>
<cn cellml:units="millivolt">17</cn>
</apply>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">2</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<ci>Vm</ci>
</apply>
<cn cellml:units="millivolt">150</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>y3</ci>
</apply>
<apply>
<minus/>
<apply>
<times/>
<ci>alpha_y3</ci>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<ci>y3</ci>
</apply>
</apply>
<apply>
<times/>
<ci>beta_y3</ci>
<ci>y3</ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="slow_time_dependent_potassium_current">
<variable name="i_Ks" units="picoA" public_interface="out"/>
<variable name="i_Ks_Na" units="picoA" public_interface="out"/>
<variable name="i_Ks_K" units="picoA" public_interface="out"/>
<variable name="time" units="millisecond" public_interface="in" private_interface="out"/>
<variable name="Vm" units="millivolt" public_interface="in" private_interface="out"/>
<variable name="Cai" units="millimolar" public_interface="in" private_interface="out"/>
<variable name="CF_K" units="millimolar" public_interface="in"/>
<variable name="CF_Na" units="millimolar" public_interface="in"/>
<variable name="y1" units="dimensionless" private_interface="in"/>
<variable name="y2" units="dimensionless" private_interface="in"/>
<variable name="P_Ks_K" units="picoA_per_millimolar" initial_value="5.04"/>
<variable name="P_Ks_Na" units="picoA_per_millimolar" initial_value="0.2016"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>i_Ks</ci>
<apply>
<plus/>
<ci>i_Ks_Na</ci>
<ci>i_Ks_K</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>i_Ks_K</ci>
<apply>
<times/>
<ci>P_Ks_K</ci>
<ci>CF_K</ci>
<apply>
<power/>
<ci>y1</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="dimensionless">0.9</cn>
<ci>y2</ci>
</apply>
<cn cellml:units="dimensionless">0.1</cn>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>i_Ks_Na</ci>
<apply>
<times/>
<ci>P_Ks_Na</ci>
<ci>CF_Na</ci>
<apply>
<power/>
<ci>y1</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="dimensionless">0.9</cn>
<ci>y2</ci>
</apply>
<cn cellml:units="dimensionless">0.1</cn>
</apply>
</apply>
</apply>
</math>
</component>
<component name="slow_time_dependent_potassium_current_y1_gate">
<variable name="y1" units="dimensionless" initial_value="0.09738789658609195" public_interface="out"/>
<variable name="alpha_y1" units="per_millisecond"/>
<variable name="beta_y1" units="per_millisecond"/>
<variable name="Vm" units="millivolt" public_interface="in"/>
<variable name="time" units="millisecond" public_interface="in"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>alpha_y1</ci>
<apply>
<divide/>
<cn cellml:units="per_millisecond">1</cn>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="dimensionless">85</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<ci>Vm</ci>
</apply>
<cn cellml:units="millivolt">10.5</cn>
</apply>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">370</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<ci>Vm</ci>
</apply>
<cn cellml:units="millivolt">62</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>beta_y1</ci>
<apply>
<divide/>
<cn cellml:units="per_millisecond">1</cn>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="dimensionless">1450</cn>
<apply>
<exp/>
<apply>
<divide/>
<ci>Vm</ci>
<cn cellml:units="millivolt">20</cn>
</apply>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">260</cn>
<apply>
<exp/>
<apply>
<divide/>
<ci>Vm</ci>
<cn cellml:units="millivolt">100</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>y1</ci>
</apply>
<apply>
<minus/>
<apply>
<times/>
<ci>alpha_y1</ci>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<ci>y1</ci>
</apply>
</apply>
<apply>
<times/>
<ci>beta_y1</ci>
<ci>y1</ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="slow_time_dependent_potassium_current_y2_gate">
<variable name="y2" units="dimensionless" initial_value="0.09745345578743213" public_interface="out"/>
<variable name="alpha_y2" units="per_millisecond"/>
<variable name="beta_y2" units="per_millisecond" initial_value="0.004444"/>
<variable name="Cai" units="millimolar" public_interface="in"/>
<variable name="Vm" units="millivolt" public_interface="in"/>
<variable name="time" units="millisecond" public_interface="in"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>alpha_y2</ci>
<apply>
<times/>
<cn cellml:units="per_millimolar_millisecond">3.7</cn>
<ci>Cai</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>y2</ci>
</apply>
<apply>
<minus/>
<apply>
<times/>
<ci>alpha_y2</ci>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<ci>y2</ci>
</apply>
</apply>
<apply>
<times/>
<ci>beta_y2</ci>
<ci>y2</ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="transient_outward_current">
<variable name="i_to" units="picoA" public_interface="out"/>
<variable name="i_to_Na" units="picoA" public_interface="out"/>
<variable name="i_to_K" units="picoA" public_interface="out"/>
<variable name="time" units="millisecond" public_interface="in" private_interface="out"/>
<variable name="Vm" units="millivolt" public_interface="in" private_interface="out"/>
<variable name="CF_K" units="millimolar" public_interface="in"/>
<variable name="CF_Na" units="millimolar" public_interface="in"/>
<variable name="y1" units="dimensionless" private_interface="in"/>
<variable name="y2" units="dimensionless" private_interface="in"/>
<variable name="P_to_K" units="picoA_per_millimolar" initial_value="0.033"/>
<variable name="P_to_Na" units="picoA_per_millimolar" initial_value="0.00297"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>i_to</ci>
<apply>
<plus/>
<ci>i_to_Na</ci>
<ci>i_to_K</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>i_to_K</ci>
<apply>
<times/>
<ci>P_to_K</ci>
<ci>CF_K</ci>
<apply>
<power/>
<ci>y1</ci>
<cn cellml:units="dimensionless">3</cn>
</apply>
<ci>y2</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>i_to_Na</ci>
<apply>
<times/>
<ci>P_to_Na</ci>
<ci>CF_Na</ci>
<apply>
<power/>
<ci>y1</ci>
<cn cellml:units="dimensionless">3</cn>
</apply>
<ci>y2</ci>
</apply>
</apply>
</math>
</component>
<component name="transient_outward_current_y1_gate">
<variable name="y1" units="dimensionless" initial_value="7.956883250874798e-4" public_interface="out"/>
<variable name="alpha_y1" units="per_millisecond"/>
<variable name="beta_y1" units="per_millisecond"/>
<variable name="Vm" units="millivolt" public_interface="in"/>
<variable name="time" units="millisecond" public_interface="in"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>alpha_y1</ci>
<apply>
<divide/>
<cn cellml:units="per_millisecond">1</cn>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="dimensionless">11</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<ci>Vm</ci>
</apply>
<cn cellml:units="millivolt">28</cn>
</apply>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">0.2</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<ci>Vm</ci>
</apply>
<cn cellml:units="millivolt">400</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>beta_y1</ci>
<apply>
<divide/>
<cn cellml:units="per_millisecond">1</cn>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="dimensionless">4.4</cn>
<apply>
<exp/>
<apply>
<divide/>
<ci>Vm</ci>
<cn cellml:units="millivolt">16</cn>
</apply>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">0.2</cn>
<apply>
<exp/>
<apply>
<divide/>
<ci>Vm</ci>
<cn cellml:units="millivolt">500</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>y1</ci>
</apply>
<apply>
<minus/>
<apply>
<times/>
<ci>alpha_y1</ci>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<ci>y1</ci>
</apply>
</apply>
<apply>
<times/>
<ci>beta_y1</ci>
<ci>y1</ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="transient_outward_current_y2_gate">
<variable name="y2" units="dimensionless" initial_value="0.9999125083105881" public_interface="out"/>
<variable name="alpha_y2" units="per_millisecond"/>
<variable name="beta_y2" units="per_millisecond"/>
<variable name="Vm" units="millivolt" public_interface="in"/>
<variable name="time" units="millisecond" public_interface="in"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>alpha_y2</ci>
<apply>
<divide/>
<apply>
<times/>
<cn cellml:units="per_millisecond">0.0038</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<plus/>
<ci>Vm</ci>
<cn cellml:units="millivolt">13.5</cn>
</apply>
</apply>
<cn cellml:units="millivolt">11.3</cn>
</apply>
</apply>
</apply>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<times/>
<cn cellml:units="dimensionless">0.051335</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<plus/>
<ci>Vm</ci>
<cn cellml:units="millivolt">13.5</cn>
</apply>
</apply>
<cn cellml:units="millivolt">11.3</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>beta_y2</ci>
<apply>
<divide/>
<apply>
<times/>
<cn cellml:units="per_millisecond">0.0038</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<plus/>
<ci>Vm</ci>
<cn cellml:units="millivolt">13.5</cn>
</apply>
<cn cellml:units="millivolt">11.3</cn>
</apply>
</apply>
</apply>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<times/>
<cn cellml:units="dimensionless">0.067083</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<plus/>
<ci>Vm</ci>
<cn cellml:units="millivolt">13.5</cn>
</apply>
<cn cellml:units="millivolt">11.3</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>y2</ci>
</apply>
<apply>
<minus/>
<apply>
<times/>
<ci>alpha_y2</ci>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<ci>y2</ci>
</apply>
</apply>
<apply>
<times/>
<ci>beta_y2</ci>
<ci>y2</ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="background_NSC_current">
<variable name="i_bNSC" units="picoA" public_interface="out"/>
<variable name="i_bNSC_Na" units="picoA" public_interface="out"/>
<variable name="i_bNSC_K" units="picoA" public_interface="out"/>
<variable name="P_bNSC" units="picoA_per_millimolar" initial_value="0.385"/>
<variable name="time" units="millisecond" public_interface="in"/>
<variable name="CF_K" units="millimolar" public_interface="in"/>
<variable name="CF_Na" units="millimolar" public_interface="in"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>i_bNSC</ci>
<apply>
<plus/>
<ci>i_bNSC_K</ci>
<ci>i_bNSC_Na</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>i_bNSC_K</ci>
<apply>
<times/>
<cn cellml:units="dimensionless">0.4</cn>
<ci>P_bNSC</ci>
<ci>CF_K</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>i_bNSC_Na</ci>
<apply>
<times/>
<ci>P_bNSC</ci>
<ci>CF_Na</ci>
</apply>
</apply>
</math>
</component>
<component name="background_Kpl_current">
<variable name="i_Kpl" units="picoA" public_interface="out"/>
<variable name="P_Kpl" units="nanoS_per_millimolar"/>
<variable name="time" units="millisecond" public_interface="in"/>
<variable name="Vm" units="millivolt" public_interface="in"/>
<variable name="CF_K" units="millimolar" public_interface="in"/>
<variable name="Ko" units="millimolar" public_interface="in"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>i_Kpl</ci>
<piecewise>
<piece>
<apply>
<times/>
<ci>P_Kpl</ci>
<ci>CF_K</ci>
<cn cellml:units="millivolt">13.0077</cn>
</apply>
<apply>
<eq/>
<ci>Vm</ci>
<apply>
<minus/>
<cn cellml:units="millivolt">3</cn>
</apply>
</apply>
</piece>
<otherwise>
<apply>
<divide/>
<apply>
<times/>
<ci>P_Kpl</ci>
<ci>CF_K</ci>
<apply>
<plus/>
<ci>Vm</ci>
<cn cellml:units="millivolt">3</cn>
</apply>
</apply>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<plus/>
<ci>Vm</ci>
<cn cellml:units="millivolt">3</cn>
</apply>
</apply>
<cn cellml:units="millivolt">13</cn>
</apply>
</apply>
</apply>
</apply>
</otherwise>
</piecewise>
</apply>
<apply>
<eq/>
<ci>P_Kpl</ci>
<apply>
<times/>
<cn cellml:units="nanoS_per_millimolar">0.00011</cn>
<apply>
<power/>
<apply>
<divide/>
<ci>Ko</ci>
<cn cellml:units="millimolar">5.4</cn>
</apply>
<cn cellml:units="dimensionless">0.16</cn>
</apply>
</apply>
</apply>
</math>
</component>
<component name="background_lCa_current">
<variable name="i_lCa" units="picoA" public_interface="out"/>
<variable name="i_lCa_Na" units="picoA" public_interface="out"/>
<variable name="i_lCa_K" units="picoA" public_interface="out"/>
<variable name="P_lCa" units="picoA_per_millimolar" initial_value="0.11"/>
<variable name="p_open" units="dimensionless"/>
<variable name="time" units="millisecond" public_interface="in"/>
<variable name="CF_K" units="millimolar" public_interface="in"/>
<variable name="CF_Na" units="millimolar" public_interface="in"/>
<variable name="Cai" units="millimolar" public_interface="in"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>i_lCa</ci>
<apply>
<plus/>
<ci>i_lCa_K</ci>
<ci>i_lCa_Na</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>i_lCa_K</ci>
<apply>
<times/>
<ci>P_lCa</ci>
<ci>CF_K</ci>
<ci>p_open</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>i_lCa_Na</ci>
<apply>
<times/>
<ci>P_lCa</ci>
<ci>CF_Na</ci>
<ci>p_open</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>p_open</ci>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<power/>
<apply>
<divide/>
<cn cellml:units="millimolar">0.0012</cn>
<ci>Cai</ci>
</apply>
<cn cellml:units="dimensionless">3</cn>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<component name="background_KATP_current">
<variable name="i_KATP" units="picoA" public_interface="out"/>
<variable name="p_open" units="dimensionless"/>
<variable name="gamma" units="nanoS"/>
<variable name="time" units="millisecond" public_interface="in"/>
<variable name="Vm" units="millivolt" public_interface="in"/>
<variable name="E_K" units="millivolt" public_interface="in"/>
<variable name="Ko" units="millimolar" public_interface="in"/>
<variable name="ATPi" units="millimolar" public_interface="in"/>
<variable name="P_KATP" units="nanoS_per_picoF" initial_value="0.0236"/>
<variable name="N" units="picoF" initial_value="2333"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>i_KATP</ci>
<apply>
<times/>
<ci>gamma</ci>
<apply>
<minus/>
<ci>Vm</ci>
<ci>E_K</ci>
</apply>
<ci>p_open</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>gamma</ci>
<apply>
<times/>
<ci>P_KATP</ci>
<ci>N</ci>
<apply>
<power/>
<apply>
<divide/>
<ci>Ko</ci>
<cn cellml:units="millimolar">1</cn>
</apply>
<cn cellml:units="dimensionless">0.24</cn>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>p_open</ci>
<apply>
<divide/>
<cn cellml:units="dimensionless">0.8</cn>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<power/>
<apply>
<divide/>
<ci>ATPi</ci>
<cn cellml:units="millimolar">0.1</cn>
</apply>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<component name="background_Cab_current">
<variable name="i_Cab" units="picoA" public_interface="out"/>
<variable name="P_Cab" units="picoA_per_millimolar" initial_value="0.04"/>
<variable name="time" units="millisecond" public_interface="in"/>
<variable name="CF_Ca" units="millimolar" public_interface="in"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>i_Cab</ci>
<apply>
<times/>
<ci>P_Cab</ci>
<ci>CF_Ca</ci>
</apply>
</apply>
</math>
</component>
<component name="sodium_calcium_exchanger">
<variable name="i_NaCa" units="picoA" public_interface="out"/>
<variable name="p_E2Na" units="dimensionless" private_interface="out"/>
<variable name="p_E1Na" units="dimensionless" private_interface="out"/>
<variable name="p_E1Ca" units="dimensionless" private_interface="out"/>
<variable name="p_E2Ca" units="dimensionless" private_interface="out"/>
<variable name="k1" units="per_millisecond" private_interface="out"/>
<variable name="k2" units="per_millisecond" private_interface="out"/>
<variable name="k3" units="per_millisecond" initial_value="1" private_interface="out"/>
<variable name="k4" units="per_millisecond" initial_value="1" private_interface="out"/>
<variable name="Km_Nai" units="millimolar" initial_value="8.75"/>
<variable name="Km_Nao" units="millimolar" initial_value="87.5"/>
<variable name="Km_Cai" units="millimolar" initial_value="0.00138"/>
<variable name="Km_Cao" units="millimolar" initial_value="1.38"/>
<variable name="Nai" units="millimolar" public_interface="in"/>
<variable name="Nao" units="millimolar" public_interface="in"/>
<variable name="Cai" units="millimolar" public_interface="in"/>
<variable name="Cao" units="millimolar" public_interface="in"/>
<variable name="R" units="coulomb_millivolt_per_kelvin_millimole" public_interface="in"/>
<variable name="F" units="coulomb_per_millimole" public_interface="in"/>
<variable name="T" units="kelvin" public_interface="in"/>
<variable name="Vm" units="millivolt" public_interface="in"/>
<variable name="Cm" units="picoF" public_interface="in"/>
<variable name="time" units="millisecond" public_interface="in" private_interface="out"/>
<variable name="y" units="dimensionless" private_interface="in"/>
<variable name="P_NaCa" units="picoA_per_picoF" initial_value="6.81"/>
<variable name="Partition" units="dimensionless" initial_value="0.32"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>i_NaCa</ci>
<apply>
<times/>
<ci>P_NaCa</ci>
<ci>Cm</ci>
<cn cellml:units="millisecond">1</cn>
<apply>
<minus/>
<apply>
<times/>
<ci>k1</ci>
<ci>p_E1Na</ci>
<ci>y</ci>
</apply>
<apply>
<times/>
<ci>k2</ci>
<ci>p_E2Na</ci>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<ci>y</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>p_E1Na</ci>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<times/>
<apply>
<power/>
<apply>
<divide/>
<ci>Km_Nai</ci>
<ci>Nai</ci>
</apply>
<cn cellml:units="dimensionless">3</cn>
</apply>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<divide/>
<ci>Cai</ci>
<ci>Km_Cai</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>p_E2Na</ci>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<times/>
<apply>
<power/>
<apply>
<divide/>
<ci>Km_Nao</ci>
<ci>Nao</ci>
</apply>
<cn cellml:units="dimensionless">3</cn>
</apply>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<divide/>
<ci>Cao</ci>
<ci>Km_Cao</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>p_E1Ca</ci>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<times/>
<apply>
<divide/>
<ci>Km_Cai</ci>
<ci>Cai</ci>
</apply>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<power/>
<apply>
<divide/>
<ci>Nai</ci>
<ci>Km_Nai</ci>
</apply>
<cn cellml:units="dimensionless">3</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>p_E2Ca</ci>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<times/>
<apply>
<divide/>
<ci>Km_Cao</ci>
<ci>Cao</ci>
</apply>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<power/>
<apply>
<divide/>
<ci>Nao</ci>
<ci>Km_Nao</ci>
</apply>
<cn cellml:units="dimensionless">3</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>k1</ci>
<apply>
<times/>
<cn cellml:units="per_millisecond">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<ci>Partition</ci>
<ci>F</ci>
<ci>Vm</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>k2</ci>
<apply>
<times/>
<cn cellml:units="per_millisecond">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<ci>Partition</ci>
<cn cellml:units="dimensionless">1</cn>
</apply>
<ci>F</ci>
<ci>Vm</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<component name="sodium_calcium_exchanger_y_gate">
<variable name="y" units="dimensionless" initial_value="0.9891789193465331" public_interface="out"/>
<variable name="alpha_y" units="per_millisecond"/>
<variable name="beta_y" units="per_millisecond"/>
<variable name="time" units="millisecond" public_interface="in"/>
<variable name="p_E1Na" units="dimensionless" public_interface="in"/>
<variable name="p_E2Na" units="dimensionless" public_interface="in"/>
<variable name="p_E1Ca" units="dimensionless" public_interface="in"/>
<variable name="p_E2Ca" units="dimensionless" public_interface="in"/>
<variable name="k1" units="per_millisecond" public_interface="in"/>
<variable name="k2" units="per_millisecond" public_interface="in"/>
<variable name="k3" units="per_millisecond" public_interface="in"/>
<variable name="k4" units="per_millisecond" public_interface="in"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>alpha_y</ci>
<apply>
<plus/>
<apply>
<times/>
<ci>k2</ci>
<ci>p_E2Na</ci>
</apply>
<apply>
<times/>
<ci>k4</ci>
<ci>p_E2Ca</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>beta_y</ci>
<apply>
<plus/>
<apply>
<times/>
<ci>k1</ci>
<ci>p_E1Na</ci>
</apply>
<apply>
<times/>
<ci>k3</ci>
<ci>p_E1Ca</ci>
</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="sodium_potassium_pump">
<variable name="i_NaK" units="picoA" public_interface="out"/>
<variable name="p_E2Na" units="dimensionless" private_interface="out"/>
<variable name="p_E1Na" units="dimensionless" private_interface="out"/>
<variable name="p_E1K" units="dimensionless" private_interface="out"/>
<variable name="p_E2K" units="dimensionless" private_interface="out"/>
<variable name="k1" units="per_millisecond" private_interface="out"/>
<variable name="k2" units="per_millisecond" initial_value="0.04" private_interface="out"/>
<variable name="k3" units="per_millisecond" initial_value="0.01" private_interface="out"/>
<variable name="k4" units="per_millisecond" initial_value="0.165" private_interface="out"/>
<variable name="Km_Nai" units="millimolar" initial_value="4.05"/>
<variable name="Km_Nao" units="millimolar" initial_value="69.8"/>
<variable name="Km_Ki" units="millimolar" initial_value="32.88"/>
<variable name="Km_Ko" units="millimolar" initial_value="0.258"/>
<variable name="Km_ATP" units="millimolar" initial_value="0.094"/>
<variable name="Nao_Eff" units="millimolar"/>
<variable name="ATPi" units="millimolar" public_interface="in"/>
<variable name="Nai" units="millimolar" public_interface="in"/>
<variable name="Nao" units="millimolar" public_interface="in"/>
<variable name="Ki" units="millimolar" public_interface="in"/>
<variable name="Ko" units="millimolar" public_interface="in"/>
<variable name="R" units="coulomb_millivolt_per_kelvin_millimole" public_interface="in"/>
<variable name="F" units="coulomb_per_millimole" public_interface="in"/>
<variable name="T" units="kelvin" public_interface="in"/>
<variable name="Vm" units="millivolt" public_interface="in"/>
<variable name="Cm" units="picoF" public_interface="in"/>
<variable name="time" units="millisecond" public_interface="in" private_interface="out"/>
<variable name="y" units="dimensionless" private_interface="in"/>
<variable name="P_NaK" units="picoA_per_picoF" initial_value="21"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>i_NaK</ci>
<apply>
<times/>
<ci>P_NaK</ci>
<ci>Cm</ci>
<cn cellml:units="millisecond">1</cn>
<apply>
<minus/>
<apply>
<times/>
<ci>k1</ci>
<ci>p_E1Na</ci>
<ci>y</ci>
</apply>
<apply>
<times/>
<ci>k2</ci>
<ci>p_E2Na</ci>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<ci>y</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>p_E1Na</ci>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<times/>
<apply>
<power/>
<apply>
<divide/>
<ci>Km_Nai</ci>
<ci>Nai</ci>
</apply>
<cn cellml:units="dimensionless">1.06</cn>
</apply>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<power/>
<apply>
<divide/>
<ci>Ki</ci>
<ci>Km_Ki</ci>
</apply>
<cn cellml:units="dimensionless">1.12</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>p_E2Na</ci>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<times/>
<apply>
<power/>
<apply>
<divide/>
<ci>Km_Nao</ci>
<ci>Nao_Eff</ci>
</apply>
<cn cellml:units="dimensionless">1.06</cn>
</apply>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<power/>
<apply>
<divide/>
<ci>Ko</ci>
<ci>Km_Ko</ci>
</apply>
<cn cellml:units="dimensionless">1.12</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>p_E1K</ci>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<times/>
<apply>
<power/>
<apply>
<divide/>
<ci>Km_Ki</ci>
<ci>Ki</ci>
</apply>
<cn cellml:units="dimensionless">1.12</cn>
</apply>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<power/>
<apply>
<divide/>
<ci>Nai</ci>
<ci>Km_Nai</ci>
</apply>
<cn cellml:units="dimensionless">1.06</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>p_E2K</ci>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<times/>
<apply>
<power/>
<apply>
<divide/>
<ci>Km_Ko</ci>
<ci>Ko</ci>
</apply>
<cn cellml:units="dimensionless">1.12</cn>
</apply>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<power/>
<apply>
<divide/>
<ci>Nao_Eff</ci>
<ci>Km_Nao</ci>
</apply>
<cn cellml:units="dimensionless">1.06</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>k1</ci>
<apply>
<divide/>
<cn cellml:units="per_millisecond">0.37</cn>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<divide/>
<ci>Km_ATP</ci>
<ci>ATPi</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>Nao_Eff</ci>
<apply>
<times/>
<ci>Nao</ci>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="dimensionless">0.82</cn>
</apply>
<ci>F</ci>
<ci>Vm</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<component name="sodium_potassium_pump_y_gate">
<variable name="y" units="dimensionless" initial_value="0.5910747147428818" public_interface="out"/>
<variable name="alpha_y" units="per_millisecond"/>
<variable name="beta_y" units="per_millisecond"/>
<variable name="time" units="millisecond" public_interface="in"/>
<variable name="p_E1Na" units="dimensionless" public_interface="in"/>
<variable name="p_E2Na" units="dimensionless" public_interface="in"/>
<variable name="p_E1K" units="dimensionless" public_interface="in"/>
<variable name="p_E2K" units="dimensionless" public_interface="in"/>
<variable name="k1" units="per_millisecond" public_interface="in"/>
<variable name="k2" units="per_millisecond" public_interface="in"/>
<variable name="k3" units="per_millisecond" public_interface="in"/>
<variable name="k4" units="per_millisecond" public_interface="in"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>alpha_y</ci>
<apply>
<plus/>
<apply>
<times/>
<ci>k2</ci>
<ci>p_E2Na</ci>
</apply>
<apply>
<times/>
<ci>k4</ci>
<ci>p_E2K</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>beta_y</ci>
<apply>
<plus/>
<apply>
<times/>
<ci>k1</ci>
<ci>p_E1Na</ci>
</apply>
<apply>
<times/>
<ci>k3</ci>
<ci>p_E1K</ci>
</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="SR_calcium_pump">
<variable name="i_SR_U" units="picoA" public_interface="out"/>
<variable name="p_E2Ca" units="dimensionless" private_interface="out"/>
<variable name="p_E1Ca" units="dimensionless" private_interface="out"/>
<variable name="p_E1" units="dimensionless" private_interface="out"/>
<variable name="p_E2" units="dimensionless" private_interface="out"/>
<variable name="k1" units="per_millisecond" initial_value="0.01" private_interface="out"/>
<variable name="k2" units="per_millisecond" private_interface="out"/>
<variable name="k3" units="per_millisecond" initial_value="1" private_interface="out"/>
<variable name="k4" units="per_millisecond" initial_value="0.01" private_interface="out"/>
<variable name="Km_CaSR" units="millimolar" initial_value="0.08"/>
<variable name="Km_CaCyto" units="millimolar" initial_value="0.0008"/>
<variable name="Km_ATP" units="millimolar" initial_value="0.1"/>
<variable name="i_max" units="picoA" initial_value="162500"/>
<variable name="ATPi" units="millimolar" public_interface="in"/>
<variable name="Cai" units="millimolar" public_interface="in"/>
<variable name="Caup" units="millimolar" public_interface="in"/>
<variable name="time" units="millisecond" 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_SR_U</ci>
<apply>
<times/>
<ci>i_max</ci>
<cn cellml:units="millisecond">1</cn>
<apply>
<minus/>
<apply>
<times/>
<ci>k1</ci>
<ci>p_E1Ca</ci>
<ci>y</ci>
</apply>
<apply>
<times/>
<ci>k2</ci>
<ci>p_E2Ca</ci>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<ci>y</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>p_E1Ca</ci>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<divide/>
<ci>Km_CaSR</ci>
<ci>Caup</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>p_E2Ca</ci>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<divide/>
<ci>Km_CaCyto</ci>
<ci>Cai</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>p_E1</ci>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<ci>p_E1Ca</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>p_E2</ci>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<ci>p_E2Ca</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>k2</ci>
<apply>
<divide/>
<cn cellml:units="per_millisecond">1</cn>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<divide/>
<ci>Km_ATP</ci>
<ci>ATPi</ci>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<component name="SR_calcium_pump_y_gate">
<variable name="y" units="dimensionless" initial_value="0.46108441538480216" public_interface="out"/>
<variable name="alpha_y" units="per_millisecond"/>
<variable name="beta_y" units="per_millisecond"/>
<variable name="time" units="millisecond" public_interface="in"/>
<variable name="p_E1Ca" units="dimensionless" public_interface="in"/>
<variable name="p_E2Ca" units="dimensionless" public_interface="in"/>
<variable name="p_E1" units="dimensionless" public_interface="in"/>
<variable name="p_E2" units="dimensionless" public_interface="in"/>
<variable name="k1" units="per_millisecond" public_interface="in"/>
<variable name="k2" units="per_millisecond" public_interface="in"/>
<variable name="k3" units="per_millisecond" public_interface="in"/>
<variable name="k4" units="per_millisecond" public_interface="in"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>alpha_y</ci>
<apply>
<plus/>
<apply>
<times/>
<ci>k2</ci>
<ci>p_E2Ca</ci>
</apply>
<apply>
<times/>
<ci>k4</ci>
<ci>p_E2</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>beta_y</ci>
<apply>
<plus/>
<apply>
<times/>
<ci>k1</ci>
<ci>p_E1Ca</ci>
</apply>
<apply>
<times/>
<ci>k3</ci>
<ci>p_E1</ci>
</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="RyR_channel">
<variable name="i_RyR" units="picoA" public_interface="out"/>
<variable name="P_RyR" units="picoA_per_millimolar" initial_value="62000"/>
<variable name="k1" units="per_millisecond"/>
<variable name="k2" units="per_millisecond"/>
<variable name="k3" units="per_millisecond"/>
<variable name="k4" units="per_millisecond" initial_value="0.000849"/>
<variable name="p_open_RyR" units="dimensionless" initial_value="3.4314360001543243e-4"/>
<variable name="p_close_RyR" units="dimensionless" initial_value="0.19135178123107768"/>
<variable name="p_open_CaL" units="dimensionless" public_interface="in"/>
<variable name="Cai" units="millimolar" public_interface="in"/>
<variable name="Carel" units="millimolar" public_interface="in"/>
<variable name="time" units="millisecond" public_interface="in"/>
<variable name="Diadid_Factor" units="per_picoA_millisecond" initial_value="-150"/>
<variable name="CaDiadic" units="picoA" public_interface="in"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>i_RyR</ci>
<apply>
<times/>
<ci>P_RyR</ci>
<apply>
<minus/>
<ci>Carel</ci>
<ci>Cai</ci>
</apply>
<ci>p_open_RyR</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>p_open_RyR</ci>
</apply>
<apply>
<minus/>
<apply>
<times/>
<ci>p_close_RyR</ci>
<ci>k1</ci>
</apply>
<apply>
<times/>
<ci>p_open_RyR</ci>
<ci>k2</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>p_close_RyR</ci>
</apply>
<apply>
<minus/>
<apply>
<times/>
<ci>k3</ci>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<plus/>
<ci>p_open_RyR</ci>
<ci>p_close_RyR</ci>
</apply>
</apply>
</apply>
<apply>
<times/>
<apply>
<plus/>
<ci>k1</ci>
<ci>k4</ci>
</apply>
<ci>p_close_RyR</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>k1</ci>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="per_millisecond">280000</cn>
<apply>
<power/>
<apply>
<divide/>
<ci>Cai</ci>
<cn cellml:units="millimolar">1</cn>
</apply>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
<apply>
<times/>
<ci>Diadid_Factor</ci>
<ci>CaDiadic</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>k2</ci>
<apply>
<divide/>
<cn cellml:units="per_millisecond">0.08</cn>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<divide/>
<cn cellml:units="millimolar">0.36</cn>
<ci>Carel</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>k3</ci>
<apply>
<times/>
<cn cellml:units="per_millisecond">0.000377</cn>
<apply>
<power/>
<apply>
<divide/>
<ci>Carel</ci>
<cn cellml:units="millimolar">1</cn>
</apply>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
</apply>
</math>
</component>
<component name="SR_T_current">
<variable name="i_SR_T" units="picoA" public_interface="out"/>
<variable name="P_SR_T" units="picoA_per_millimolar" initial_value="386"/>
<variable name="Carel" units="millimolar" public_interface="in"/>
<variable name="Caup" units="millimolar" public_interface="in"/>
<variable name="time" units="millisecond" public_interface="in"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>i_SR_T</ci>
<apply>
<times/>
<ci>P_SR_T</ci>
<apply>
<minus/>
<ci>Caup</ci>
<ci>Carel</ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="SR_L_current">
<variable name="i_SR_L" units="picoA" public_interface="out"/>
<variable name="P_SR_L" units="picoA_per_millimolar" initial_value="459"/>
<variable name="Cai" units="millimolar" public_interface="in"/>
<variable name="Caup" units="millimolar" public_interface="in"/>
<variable name="time" units="millisecond" public_interface="in"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>i_SR_L</ci>
<apply>
<times/>
<ci>P_SR_L</ci>
<apply>
<minus/>
<ci>Caup</ci>
<ci>Cai</ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="Ca_concentrations_in_SR">
<variable name="Ca_Total" units="millimolar" initial_value="9.455741736977666"/>
<variable name="Carel" units="millimolar" public_interface="out"/>
<variable name="Caup" units="millimolar" initial_value="2.611712901567567" public_interface="out"/>
<variable name="V_rel" units="micrometre3" initial_value="160"/>
<variable name="V_up" units="micrometre3" initial_value="400"/>
<variable name="time" units="millisecond" public_interface="in"/>
<variable name="F" units="coulomb_per_millimole" public_interface="in"/>
<variable name="i_SR_L" units="picoA" public_interface="in"/>
<variable name="i_SR_T" units="picoA" public_interface="in"/>
<variable name="i_RyR" units="picoA" public_interface="in"/>
<variable name="i_SR_U" units="picoA" public_interface="in"/>
<variable name="CSQN_max" units="millimolar" initial_value="10"/>
<variable name="K_mCSQN" units="millimolar" initial_value="0.8"/>
<variable name="b1" units="millimolar"/>
<variable name="c1" units="millimolar2"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>b1</ci>
<apply>
<plus/>
<apply>
<minus/>
<ci>CSQN_max</ci>
<ci>Ca_Total</ci>
</apply>
<ci>K_mCSQN</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>c1</ci>
<apply>
<times/>
<ci>K_mCSQN</ci>
<ci>Ca_Total</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>Carel</ci>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<root/>
<apply>
<plus/>
<apply>
<power/>
<ci>b1</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">4</cn>
<ci>c1</ci>
</apply>
</apply>
</apply>
<ci>b1</ci>
</apply>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>Ca_Total</ci>
</apply>
<apply>
<divide/>
<apply>
<minus/>
<ci>i_SR_T</ci>
<ci>i_RyR</ci>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">2</cn>
<ci>F</ci>
<ci>V_rel</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>Caup</ci>
</apply>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<minus/>
<apply>
<minus/>
<ci>i_SR_U</ci>
</apply>
<ci>i_SR_T</ci>
</apply>
<ci>i_SR_L</ci>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">2</cn>
<ci>F</ci>
<ci>V_up</ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="NL_model">
<variable name="dCaidt" units="millimolar_per_millisecond" public_interface="out"/>
<variable name="dATPdt" units="millimolar_per_millisecond" public_interface="out"/>
<variable name="EffFraction" units="dimensionless"/>
<variable name="pCa" units="dimensionless" initial_value="0.02490898775497523"/>
<variable name="pCaCB" units="dimensionless" initial_value="0.001990153835322864"/>
<variable name="pCB" units="dimensionless" initial_value="4.2941813853474524e-4"/>
<variable name="p" units="dimensionless"/>
<variable name="T_t" units="millimolar" initial_value="0.07"/>
<variable name="Q_a" units="per_millisecond"/>
<variable name="Q_b" units="per_millisecond"/>
<variable name="Q_r" units="per_millisecond"/>
<variable name="Q_d" units="per_millisecond"/>
<variable name="Q_d1" units="per_millisecond"/>
<variable name="Q_d2" units="per_millisecond"/>
<variable name="Y_1" units="per_millimolar_millisecond" initial_value="39"/>
<variable name="Y_2" units="per_millisecond" initial_value="0.0039"/>
<variable name="Y_3" units="per_millisecond" initial_value="0.03"/>
<variable name="Y_4" units="per_millisecond" initial_value="0.12"/>
<variable name="Y_d" units="millisecond_per_micrometre2" initial_value="0.027"/>
<variable name="Z_1" units="per_millisecond" initial_value="0.03"/>
<variable name="Z_2" units="per_millisecond" initial_value="0.0039"/>
<variable name="Z_3" units="per_millimolar_millisecond" initial_value="1560"/>
<variable name="h" units="micrometre"/>
<variable name="L_a" units="micrometre" initial_value="1.17"/>
<variable name="L" units="micrometre" initial_value="0.9623799975411884"/>
<variable name="ForceCB" units="mN_per_mm2"/>
<variable name="X" units="micrometre" initial_value="0.9573749975411884"/>
<variable name="time" units="millisecond" public_interface="in"/>
<variable name="Cai" units="millimolar" public_interface="in"/>
<variable name="NewCBF" units="mN_per_mm2_micrometre"/>
<variable name="CBBound" units="millimolar"/>
<variable name="KForceEC" units="mN_per_mm2_micrometre5" initial_value="140000"/>
<variable name="ZeroForceEL" units="micrometre" initial_value="0.97"/>
<variable name="KForceLinearEc" units="mN_per_mm2_micrometre" initial_value="200"/>
<variable name="ForceFactor" units="mN_per_mm2_micrometre_millimolar" initial_value="1800000"/>
<variable name="ForceEcomp" units="mN_per_mm2"/>
<variable name="B" units="per_millisecond" initial_value="1.2"/>
<variable name="h_c" units="micrometre" initial_value="0.005"/>
<variable name="ForceExt" units="mN_per_mm2"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>h</ci>
<apply>
<minus/>
<ci>L</ci>
<ci>X</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>p</ci>
<apply>
<minus/>
<apply>
<minus/>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<ci>pCa</ci>
</apply>
<ci>pCaCB</ci>
</apply>
<ci>pCB</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>Q_b</ci>
<apply>
<minus/>
<apply>
<times/>
<ci>Y_1</ci>
<ci>Cai</ci>
<ci>p</ci>
</apply>
<apply>
<times/>
<ci>Z_1</ci>
<ci>pCa</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>EffFraction</ci>
<apply>
<exp/>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="per_micrometre2">20</cn>
</apply>
<apply>
<power/>
<apply>
<minus/>
<ci>L</ci>
<ci>L_a</ci>
</apply>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>Q_a</ci>
<apply>
<minus/>
<apply>
<times/>
<ci>Y_2</ci>
<ci>pCa</ci>
<ci>EffFraction</ci>
</apply>
<apply>
<times/>
<ci>Z_2</ci>
<ci>pCaCB</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>Q_r</ci>
<apply>
<minus/>
<apply>
<times/>
<ci>Y_3</ci>
<ci>pCaCB</ci>
</apply>
<apply>
<times/>
<ci>Z_3</ci>
<ci>pCB</ci>
<ci>Cai</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>Q_d</ci>
<apply>
<times/>
<ci>Y_4</ci>
<ci>pCB</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>Q_d1</ci>
<apply>
<times/>
<ci>Y_d</ci>
<apply>
<power/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>X</ci>
</apply>
<cn cellml:units="dimensionless">2</cn>
</apply>
<ci>pCB</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>Q_d2</ci>
<apply>
<times/>
<ci>Y_d</ci>
<apply>
<power/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>X</ci>
</apply>
<cn cellml:units="dimensionless">2</cn>
</apply>
<ci>pCaCB</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>pCa</ci>
</apply>
<apply>
<minus/>
<ci>Q_b</ci>
<ci>Q_a</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>pCaCB</ci>
</apply>
<apply>
<minus/>
<apply>
<minus/>
<ci>Q_a</ci>
<ci>Q_r</ci>
</apply>
<ci>Q_d2</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>pCB</ci>
</apply>
<apply>
<minus/>
<apply>
<minus/>
<ci>Q_r</ci>
<ci>Q_d</ci>
</apply>
<ci>Q_d1</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>dCaidt</ci>
<apply>
<times/>
<ci>T_t</ci>
<apply>
<minus/>
<apply>
<plus/>
<ci>Q_d2</ci>
<ci>Q_r</ci>
</apply>
<ci>Q_b</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>dATPdt</ci>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="per_millisecond">0.4</cn>
</apply>
<ci>pCaCB</ci>
<ci>T_t</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>CBBound</ci>
<apply>
<times/>
<ci>T_t</ci>
<apply>
<plus/>
<ci>pCaCB</ci>
<ci>pCB</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>NewCBF</ci>
<apply>
<times/>
<ci>ForceFactor</ci>
<ci>CBBound</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>ForceEcomp</ci>
<apply>
<plus/>
<apply>
<times/>
<ci>KForceEC</ci>
<apply>
<power/>
<apply>
<minus/>
<ci>ZeroForceEL</ci>
<ci>L</ci>
</apply>
<cn cellml:units="dimensionless">5</cn>
</apply>
</apply>
<apply>
<times/>
<ci>KForceLinearEc</ci>
<apply>
<minus/>
<ci>ZeroForceEL</ci>
<ci>L</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>ForceCB</ci>
<apply>
<times/>
<ci>NewCBF</ci>
<ci>h</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>ForceExt</ci>
<apply>
<plus/>
<apply>
<minus/>
<ci>ForceEcomp</ci>
</apply>
<ci>ForceCB</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>X</ci>
</apply>
<apply>
<times/>
<ci>B</ci>
<apply>
<minus/>
<ci>h</ci>
<ci>h_c</ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="internal_concentrations"/>
<component name="intracellular_currents"/>
<component name="transmembrane_currents"/>
<group>
<relationship_ref relationship="containment"/>
<component_ref component="membrane">
<component_ref component="internal_concentrations">
<component_ref component="internal_ion_concentrations"/>
<component_ref component="Ca_concentrations_in_SR"/>
<component_ref component="constant_field_equations"/>
<component_ref component="ATP_production"/>
<component_ref component="NL_model"/>
</component_ref>
<component_ref component="transmembrane_currents">
<component_ref component="sodium_current">
<component_ref component="sodium_current_voltage_dependent_gate"/>
<component_ref component="sodium_current_ultra_slow_gate"/>
</component_ref>
<component_ref component="L_type_Ca_channel">
<component_ref component="L_type_Ca_channel_voltage_dependent_gate"/>
<component_ref component="L_type_Ca_channel_Ca_dependent_gate"/>
<component_ref component="L_type_Ca_channel_ultra_slow_gate"/>
</component_ref>
<component_ref component="T_type_Ca_channel">
<component_ref component="T_type_Ca_channel_y1_gate"/>
<component_ref component="T_type_Ca_channel_y2_gate"/>
</component_ref>
<component_ref component="time_independent_potassium_current">
<component_ref component="time_independent_potassium_current_y_gate"/>
</component_ref>
<component_ref component="rapid_time_dependent_potassium_current">
<component_ref component="rapid_time_dependent_potassium_current_y1_gate"/>
<component_ref component="rapid_time_dependent_potassium_current_y2_gate"/>
<component_ref component="rapid_time_dependent_potassium_current_y3_gate"/>
</component_ref>
<component_ref component="slow_time_dependent_potassium_current">
<component_ref component="slow_time_dependent_potassium_current_y1_gate"/>
<component_ref component="slow_time_dependent_potassium_current_y2_gate"/>
</component_ref>
<component_ref component="transient_outward_current">
<component_ref component="transient_outward_current_y1_gate"/>
<component_ref component="transient_outward_current_y2_gate"/>
</component_ref>
<component_ref component="background_NSC_current"/>
<component_ref component="background_Kpl_current"/>
<component_ref component="background_lCa_current"/>
<component_ref component="background_Cab_current"/>
<component_ref component="background_KATP_current"/>
<component_ref component="sodium_calcium_exchanger">
<component_ref component="sodium_calcium_exchanger_y_gate"/>
</component_ref>
<component_ref component="sodium_potassium_pump">
<component_ref component="sodium_potassium_pump_y_gate"/>
</component_ref>
</component_ref>
<component_ref component="intracellular_currents">
<component_ref component="SR_calcium_pump">
<component_ref component="SR_calcium_pump_y_gate"/>
</component_ref>
<component_ref component="RyR_channel"/>
<component_ref component="SR_T_current"/>
<component_ref component="SR_L_current"/>
</component_ref>
</component_ref>
</group>
<group>
<relationship_ref relationship="encapsulation"/>
<component_ref component="sodium_current">
<component_ref component="sodium_current_voltage_dependent_gate"/>
<component_ref component="sodium_current_ultra_slow_gate"/>
</component_ref>
<component_ref component="L_type_Ca_channel">
<component_ref component="L_type_Ca_channel_voltage_dependent_gate"/>
<component_ref component="L_type_Ca_channel_Ca_dependent_gate"/>
<component_ref component="L_type_Ca_channel_ultra_slow_gate"/>
</component_ref>
<component_ref component="T_type_Ca_channel">
<component_ref component="T_type_Ca_channel_y1_gate"/>
<component_ref component="T_type_Ca_channel_y2_gate"/>
</component_ref>
<component_ref component="time_independent_potassium_current">
<component_ref component="time_independent_potassium_current_y_gate"/>
</component_ref>
<component_ref component="rapid_time_dependent_potassium_current">
<component_ref component="rapid_time_dependent_potassium_current_y1_gate"/>
<component_ref component="rapid_time_dependent_potassium_current_y2_gate"/>
<component_ref component="rapid_time_dependent_potassium_current_y3_gate"/>
</component_ref>
<component_ref component="slow_time_dependent_potassium_current">
<component_ref component="slow_time_dependent_potassium_current_y1_gate"/>
<component_ref component="slow_time_dependent_potassium_current_y2_gate"/>
</component_ref>
<component_ref component="transient_outward_current">
<component_ref component="transient_outward_current_y1_gate"/>
<component_ref component="transient_outward_current_y2_gate"/>
</component_ref>
<component_ref component="sodium_calcium_exchanger">
<component_ref component="sodium_calcium_exchanger_y_gate"/>
</component_ref>
<component_ref component="sodium_potassium_pump">
<component_ref component="sodium_potassium_pump_y_gate"/>
</component_ref>
<component_ref component="SR_calcium_pump">
<component_ref component="SR_calcium_pump_y_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="internal_ion_concentrations" component_2="environment"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection>
<map_components component_1="ATP_production" component_2="environment"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection>
<map_components component_1="sodium_current" component_2="environment"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection>
<map_components component_1="L_type_Ca_channel" component_2="environment"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection>
<map_components component_1="T_type_Ca_channel" component_2="environment"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection>
<map_components component_1="SR_T_current" component_2="environment"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection>
<map_components component_1="SR_L_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="rapid_time_dependent_potassium_current" component_2="environment"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection>
<map_components component_1="slow_time_dependent_potassium_current" component_2="environment"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection>
<map_components component_1="transient_outward_current" component_2="environment"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection>
<map_components component_1="background_NSC_current" component_2="environment"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection>
<map_components component_1="background_Kpl_current" component_2="environment"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection>
<map_components component_1="background_lCa_current" component_2="environment"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection>
<map_components component_1="background_KATP_current" component_2="environment"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection>
<map_components component_1="background_Cab_current" component_2="environment"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection>
<map_components component_1="sodium_calcium_exchanger" 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="SR_calcium_pump" component_2="environment"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection>
<map_components component_1="RyR_channel" component_2="environment"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection>
<map_components component_1="Ca_concentrations_in_SR" component_2="environment"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection>
<map_components component_1="NL_model" component_2="environment"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection>
<map_components component_1="membrane" component_2="internal_ion_concentrations"/>
<map_variables variable_1="F" variable_2="F"/>
<map_variables variable_1="i_ext" variable_2="i_ext"/>
</connection>
<connection>
<map_components component_1="membrane" component_2="constant_field_equations"/>
<map_variables variable_1="Vm" variable_2="Vm"/>
<map_variables variable_1="R" variable_2="R"/>
<map_variables variable_1="T" variable_2="T"/>
<map_variables variable_1="F" variable_2="F"/>
</connection>
<connection>
<map_components component_1="membrane" component_2="ATP_production"/>
<map_variables variable_1="Vm" variable_2="Vm"/>
<map_variables variable_1="F" variable_2="F"/>
</connection>
<connection>
<map_components component_1="membrane" component_2="sodium_current"/>
<map_variables variable_1="Vm" variable_2="Vm"/>
<map_variables variable_1="i_Na" variable_2="i_Na"/>
</connection>
<connection>
<map_components component_1="membrane" component_2="L_type_Ca_channel"/>
<map_variables variable_1="Vm" variable_2="Vm"/>
<map_variables variable_1="i_Ca_L" variable_2="i_Ca_L"/>
</connection>
<connection>
<map_components component_1="membrane" component_2="T_type_Ca_channel"/>
<map_variables variable_1="Vm" variable_2="Vm"/>
<map_variables variable_1="i_Ca_T" variable_2="i_Ca_T"/>
</connection>
<connection>
<map_components component_1="membrane" component_2="time_independent_potassium_current"/>
<map_variables variable_1="Vm" variable_2="Vm"/>
<map_variables variable_1="Cm" variable_2="Cm"/>
<map_variables variable_1="R" variable_2="R"/>
<map_variables variable_1="T" variable_2="T"/>
<map_variables variable_1="F" variable_2="F"/>
<map_variables variable_1="i_K1" variable_2="i_K1"/>
</connection>
<connection>
<map_components component_1="membrane" component_2="rapid_time_dependent_potassium_current"/>
<map_variables variable_1="Vm" variable_2="Vm"/>
<map_variables variable_1="Cm" variable_2="Cm"/>
<map_variables variable_1="i_Kr" variable_2="i_Kr"/>
</connection>
<connection>
<map_components component_1="membrane" component_2="slow_time_dependent_potassium_current"/>
<map_variables variable_1="Vm" variable_2="Vm"/>
<map_variables variable_1="i_Ks" variable_2="i_Ks"/>
</connection>
<connection>
<map_components component_1="membrane" component_2="transient_outward_current"/>
<map_variables variable_1="Vm" variable_2="Vm"/>
<map_variables variable_1="i_to" variable_2="i_to"/>
</connection>
<connection>
<map_components component_1="membrane" component_2="background_NSC_current"/>
<map_variables variable_1="i_bNSC" variable_2="i_bNSC"/>
</connection>
<connection>
<map_components component_1="membrane" component_2="background_Kpl_current"/>
<map_variables variable_1="i_Kpl" variable_2="i_Kpl"/>
<map_variables variable_1="Vm" variable_2="Vm"/>
</connection>
<connection>
<map_components component_1="membrane" component_2="background_lCa_current"/>
<map_variables variable_1="i_lCa" variable_2="i_lCa"/>
</connection>
<connection>
<map_components component_1="membrane" component_2="background_KATP_current"/>
<map_variables variable_1="i_KATP" variable_2="i_KATP"/>
<map_variables variable_1="Vm" variable_2="Vm"/>
</connection>
<connection>
<map_components component_1="membrane" component_2="background_Cab_current"/>
<map_variables variable_1="i_Cab" variable_2="i_Cab"/>
</connection>
<connection>
<map_components component_1="membrane" component_2="sodium_calcium_exchanger"/>
<map_variables variable_1="Vm" variable_2="Vm"/>
<map_variables variable_1="Cm" variable_2="Cm"/>
<map_variables variable_1="R" variable_2="R"/>
<map_variables variable_1="T" variable_2="T"/>
<map_variables variable_1="F" variable_2="F"/>
<map_variables variable_1="i_NaCa" variable_2="i_NaCa"/>
</connection>
<connection>
<map_components component_1="membrane" component_2="sodium_potassium_pump"/>
<map_variables variable_1="Vm" variable_2="Vm"/>
<map_variables variable_1="Cm" variable_2="Cm"/>
<map_variables variable_1="R" variable_2="R"/>
<map_variables variable_1="T" variable_2="T"/>
<map_variables variable_1="F" variable_2="F"/>
<map_variables variable_1="i_NaK" variable_2="i_NaK"/>
</connection>
<connection>
<map_components component_1="membrane" component_2="Ca_concentrations_in_SR"/>
<map_variables variable_1="F" variable_2="F"/>
</connection>
<connection>
<map_components component_1="internal_ion_concentrations" component_2="constant_field_equations"/>
<map_variables variable_1="Nai" variable_2="Nai"/>
<map_variables variable_1="Ki" variable_2="Ki"/>
<map_variables variable_1="Cai" variable_2="Cai"/>
</connection>
<connection>
<map_components component_1="internal_ion_concentrations" component_2="NL_model"/>
<map_variables variable_1="Cai" variable_2="Cai"/>
<map_variables variable_1="dCaidt_NL" variable_2="dCaidt"/>
</connection>
<connection>
<map_components component_1="external_ion_concentrations" component_2="constant_field_equations"/>
<map_variables variable_1="Nao" variable_2="Nao"/>
<map_variables variable_1="Ko" variable_2="Ko"/>
<map_variables variable_1="Cao" variable_2="Cao"/>
</connection>
<connection>
<map_components component_1="external_ion_concentrations" component_2="time_independent_potassium_current"/>
<map_variables variable_1="Ko" variable_2="Ko"/>
</connection>
<connection>
<map_components component_1="external_ion_concentrations" component_2="rapid_time_dependent_potassium_current"/>
<map_variables variable_1="Ko" variable_2="Ko"/>
</connection>
<connection>
<map_components component_1="external_ion_concentrations" component_2="background_Kpl_current"/>
<map_variables variable_1="Ko" variable_2="Ko"/>
</connection>
<connection>
<map_components component_1="external_ion_concentrations" component_2="background_KATP_current"/>
<map_variables variable_1="Ko" variable_2="Ko"/>
</connection>
<connection>
<map_components component_1="internal_ion_concentrations" component_2="background_lCa_current"/>
<map_variables variable_1="Cai" variable_2="Cai"/>
<map_variables variable_1="i_lCa_Na" variable_2="i_lCa_Na"/>
<map_variables variable_1="i_lCa_K" variable_2="i_lCa_K"/>
</connection>
<connection>
<map_components component_1="internal_ion_concentrations" component_2="sodium_calcium_exchanger"/>
<map_variables variable_1="Nai" variable_2="Nai"/>
<map_variables variable_1="Cai" variable_2="Cai"/>
<map_variables variable_1="i_NaCa" variable_2="i_NaCa"/>
</connection>
<connection>
<map_components component_1="external_ion_concentrations" component_2="sodium_calcium_exchanger"/>
<map_variables variable_1="Nao" variable_2="Nao"/>
<map_variables variable_1="Cao" variable_2="Cao"/>
</connection>
<connection>
<map_components component_1="internal_ion_concentrations" component_2="sodium_potassium_pump"/>
<map_variables variable_1="Nai" variable_2="Nai"/>
<map_variables variable_1="Ki" variable_2="Ki"/>
<map_variables variable_1="i_NaK" variable_2="i_NaK"/>
</connection>
<connection>
<map_components component_1="external_ion_concentrations" component_2="sodium_potassium_pump"/>
<map_variables variable_1="Nao" variable_2="Nao"/>
<map_variables variable_1="Ko" variable_2="Ko"/>
</connection>
<connection>
<map_components component_1="internal_ion_concentrations" component_2="SR_calcium_pump"/>
<map_variables variable_1="Cai" variable_2="Cai"/>
<map_variables variable_1="i_SR_U" variable_2="i_SR_U"/>
</connection>
<connection>
<map_components component_1="internal_ion_concentrations" component_2="RyR_channel"/>
<map_variables variable_1="i_RyR" variable_2="i_RyR"/>
<map_variables variable_1="Cai" variable_2="Cai"/>
</connection>
<connection>
<map_components component_1="internal_ion_concentrations" component_2="SR_L_current"/>
<map_variables variable_1="Cai" variable_2="Cai"/>
<map_variables variable_1="i_SR_L" variable_2="i_SR_L"/>
</connection>
<connection>
<map_components component_1="internal_ion_concentrations" component_2="sodium_current"/>
<map_variables variable_1="i_Na_Na" variable_2="i_Na_Na"/>
<map_variables variable_1="i_Na_K" variable_2="i_Na_K"/>
</connection>
<connection>
<map_components component_1="internal_ion_concentrations" component_2="L_type_Ca_channel"/>
<map_variables variable_1="i_CaL_Na" variable_2="i_CaL_Na"/>
<map_variables variable_1="i_CaL_Ca" variable_2="i_CaL_Ca"/>
<map_variables variable_1="i_CaL_K" variable_2="i_CaL_K"/>
<map_variables variable_1="Cai" variable_2="Cai"/>
</connection>
<connection>
<map_components component_1="internal_ion_concentrations" component_2="T_type_Ca_channel"/>
<map_variables variable_1="i_Ca_T" variable_2="i_Ca_T"/>
</connection>
<connection>
<map_components component_1="internal_ion_concentrations" component_2="time_independent_potassium_current"/>
<map_variables variable_1="i_K1" variable_2="i_K1"/>
<map_variables variable_1="Ki" variable_2="Ki"/>
</connection>
<connection>
<map_components component_1="internal_ion_concentrations" component_2="rapid_time_dependent_potassium_current"/>
<map_variables variable_1="i_Kr" variable_2="i_Kr"/>
</connection>
<connection>
<map_components component_1="internal_ion_concentrations" component_2="slow_time_dependent_potassium_current"/>
<map_variables variable_1="i_Ks_Na" variable_2="i_Ks_Na"/>
<map_variables variable_1="i_Ks_K" variable_2="i_Ks_K"/>
<map_variables variable_1="Cai" variable_2="Cai"/>
</connection>
<connection>
<map_components component_1="internal_ion_concentrations" component_2="transient_outward_current"/>
<map_variables variable_1="i_to_Na" variable_2="i_to_Na"/>
<map_variables variable_1="i_to_K" variable_2="i_to_K"/>
</connection>
<connection>
<map_components component_1="internal_ion_concentrations" component_2="background_NSC_current"/>
<map_variables variable_1="i_bNSC" variable_2="i_bNSC"/>
<map_variables variable_1="i_bNSC_Na" variable_2="i_bNSC_Na"/>
<map_variables variable_1="i_bNSC_K" variable_2="i_bNSC_K"/>
</connection>
<connection>
<map_components component_1="internal_ion_concentrations" component_2="background_Kpl_current"/>
<map_variables variable_1="i_Kpl" variable_2="i_Kpl"/>
</connection>
<connection>
<map_components component_1="internal_ion_concentrations" component_2="background_KATP_current"/>
<map_variables variable_1="i_KATP" variable_2="i_KATP"/>
</connection>
<connection>
<map_components component_1="internal_ion_concentrations" component_2="background_Cab_current"/>
<map_variables variable_1="i_Cab" variable_2="i_Cab"/>
</connection>
<connection>
<map_components component_1="constant_field_equations" component_2="sodium_current"/>
<map_variables variable_1="CF_Na" variable_2="CF_Na"/>
<map_variables variable_1="CF_K" variable_2="CF_K"/>
</connection>
<connection>
<map_components component_1="constant_field_equations" component_2="L_type_Ca_channel"/>
<map_variables variable_1="CF_Ca" variable_2="CF_Ca"/>
<map_variables variable_1="CF_Na" variable_2="CF_Na"/>
<map_variables variable_1="CF_K" variable_2="CF_K"/>
</connection>
<connection>
<map_components component_1="constant_field_equations" component_2="T_type_Ca_channel"/>
<map_variables variable_1="CF_Ca" variable_2="CF_Ca"/>
</connection>
<connection>
<map_components component_1="constant_field_equations" component_2="slow_time_dependent_potassium_current"/>
<map_variables variable_1="CF_Na" variable_2="CF_Na"/>
<map_variables variable_1="CF_K" variable_2="CF_K"/>
</connection>
<connection>
<map_components component_1="constant_field_equations" component_2="transient_outward_current"/>
<map_variables variable_1="CF_Na" variable_2="CF_Na"/>
<map_variables variable_1="CF_K" variable_2="CF_K"/>
</connection>
<connection>
<map_components component_1="constant_field_equations" component_2="background_NSC_current"/>
<map_variables variable_1="CF_Na" variable_2="CF_Na"/>
<map_variables variable_1="CF_K" variable_2="CF_K"/>
</connection>
<connection>
<map_components component_1="constant_field_equations" component_2="background_Kpl_current"/>
<map_variables variable_1="CF_K" variable_2="CF_K"/>
</connection>
<connection>
<map_components component_1="constant_field_equations" component_2="background_lCa_current"/>
<map_variables variable_1="CF_Na" variable_2="CF_Na"/>
<map_variables variable_1="CF_K" variable_2="CF_K"/>
</connection>
<connection>
<map_components component_1="constant_field_equations" component_2="background_Cab_current"/>
<map_variables variable_1="CF_Ca" variable_2="CF_Ca"/>
</connection>
<connection>
<map_components component_1="ATP_production" component_2="NL_model"/>
<map_variables variable_1="dATPdt_NL" variable_2="dATPdt"/>
</connection>
<connection>
<map_components component_1="ATP_production" component_2="L_type_Ca_channel"/>
<map_variables variable_1="ATPi" variable_2="ATPi"/>
</connection>
<connection>
<map_components component_1="ATP_production" component_2="background_KATP_current"/>
<map_variables variable_1="ATPi" variable_2="ATPi"/>
</connection>
<connection>
<map_components component_1="ATP_production" component_2="sodium_potassium_pump"/>
<map_variables variable_1="i_NaK" variable_2="i_NaK"/>
<map_variables variable_1="ATPi" variable_2="ATPi"/>
</connection>
<connection>
<map_components component_1="ATP_production" component_2="SR_calcium_pump"/>
<map_variables variable_1="ATPi" variable_2="ATPi"/>
<map_variables variable_1="i_SR_U" variable_2="i_SR_U"/>
</connection>
<connection>
<map_components component_1="ATP_production" component_2="internal_ion_concentrations"/>
<map_variables variable_1="Vi" variable_2="Vi"/>
</connection>
<connection>
<map_components component_1="time_independent_potassium_current" component_2="rapid_time_dependent_potassium_current"/>
<map_variables variable_1="E_K" variable_2="E_K"/>
</connection>
<connection>
<map_components component_1="time_independent_potassium_current" component_2="background_KATP_current"/>
<map_variables variable_1="E_K" variable_2="E_K"/>
</connection>
<connection>
<map_components component_1="Ca_concentrations_in_SR" component_2="SR_calcium_pump"/>
<map_variables variable_1="Caup" variable_2="Caup"/>
<map_variables variable_1="i_SR_U" variable_2="i_SR_U"/>
</connection>
<connection>
<map_components component_1="L_type_Ca_channel" component_2="RyR_channel"/>
<map_variables variable_1="p_open_CaL" variable_2="p_open_CaL"/>
<map_variables variable_1="CaDiadic" variable_2="CaDiadic"/>
</connection>
<connection>
<map_components component_1="Ca_concentrations_in_SR" component_2="RyR_channel"/>
<map_variables variable_1="Carel" variable_2="Carel"/>
<map_variables variable_1="i_RyR" variable_2="i_RyR"/>
</connection>
<connection>
<map_components component_1="Ca_concentrations_in_SR" component_2="SR_T_current"/>
<map_variables variable_1="Caup" variable_2="Caup"/>
<map_variables variable_1="Carel" variable_2="Carel"/>
<map_variables variable_1="i_SR_T" variable_2="i_SR_T"/>
</connection>
<connection>
<map_components component_1="Ca_concentrations_in_SR" component_2="SR_L_current"/>
<map_variables variable_1="Caup" variable_2="Caup"/>
<map_variables variable_1="i_SR_L" variable_2="i_SR_L"/>
</connection>
<connection>
<map_components component_1="sodium_current" component_2="sodium_current_voltage_dependent_gate"/>
<map_variables variable_1="p_AP_Na" variable_2="p_AP_Na"/>
<map_variables variable_1="time" variable_2="time"/>
<map_variables variable_1="Vm" variable_2="Vm"/>
</connection>
<connection>
<map_components component_1="sodium_current" component_2="sodium_current_ultra_slow_gate"/>
<map_variables variable_1="y" variable_2="y"/>
<map_variables variable_1="time" variable_2="time"/>
<map_variables variable_1="Vm" variable_2="Vm"/>
</connection>
<connection>
<map_components component_1="L_type_Ca_channel" component_2="L_type_Ca_channel_voltage_dependent_gate"/>
<map_variables variable_1="p_AP_CaL" variable_2="p_AP_CaL"/>
<map_variables variable_1="time" variable_2="time"/>
<map_variables variable_1="Vm" variable_2="Vm"/>
</connection>
<connection>
<map_components component_1="L_type_Ca_channel_Ca_dependent_gate" component_2="L_type_Ca_channel_voltage_dependent_gate"/>
<map_variables variable_1="p_AP_CaL" variable_2="p_AP_CaL"/>
</connection>
<connection>
<map_components component_1="L_type_Ca_channel" component_2="L_type_Ca_channel_Ca_dependent_gate"/>
<map_variables variable_1="p_U" variable_2="p_U"/>
<map_variables variable_1="p_UCa" variable_2="p_UCa"/>
<map_variables variable_1="Cai" variable_2="Cai"/>
<map_variables variable_1="CF_Ca" variable_2="CF_Ca"/>
<map_variables variable_1="time" variable_2="time"/>
<map_variables variable_1="Vm" variable_2="Vm"/>
<map_variables variable_1="CaDiadic" variable_2="CaDiadic"/>
<map_variables variable_1="p_open_CaL" variable_2="p_open_CaL"/>
</connection>
<connection>
<map_components component_1="L_type_Ca_channel" component_2="L_type_Ca_channel_ultra_slow_gate"/>
<map_variables variable_1="y" variable_2="y"/>
<map_variables variable_1="time" variable_2="time"/>
<map_variables variable_1="Vm" variable_2="Vm"/>
</connection>
<connection>
<map_components component_1="T_type_Ca_channel" component_2="T_type_Ca_channel_y1_gate"/>
<map_variables variable_1="y1" variable_2="y1"/>
<map_variables variable_1="time" variable_2="time"/>
<map_variables variable_1="Vm" variable_2="Vm"/>
</connection>
<connection>
<map_components component_1="T_type_Ca_channel" component_2="T_type_Ca_channel_y2_gate"/>
<map_variables variable_1="y2" variable_2="y2"/>
<map_variables variable_1="time" variable_2="time"/>
<map_variables variable_1="Vm" variable_2="Vm"/>
</connection>
<connection>
<map_components component_1="time_independent_potassium_current" component_2="time_independent_potassium_current_y_gate"/>
<map_variables variable_1="fO" variable_2="fO"/>
<map_variables variable_1="y" variable_2="y"/>
<map_variables variable_1="E_K" variable_2="E_K"/>
<map_variables variable_1="time" variable_2="time"/>
<map_variables variable_1="Vm" variable_2="Vm"/>
</connection>
<connection>
<map_components component_1="rapid_time_dependent_potassium_current" component_2="rapid_time_dependent_potassium_current_y1_gate"/>
<map_variables variable_1="y1" variable_2="y1"/>
<map_variables variable_1="E_K" variable_2="E_K"/>
<map_variables variable_1="time" variable_2="time"/>
<map_variables variable_1="Vm" variable_2="Vm"/>
</connection>
<connection>
<map_components component_1="rapid_time_dependent_potassium_current" component_2="rapid_time_dependent_potassium_current_y2_gate"/>
<map_variables variable_1="y2" variable_2="y2"/>
<map_variables variable_1="E_K" variable_2="E_K"/>
<map_variables variable_1="time" variable_2="time"/>
<map_variables variable_1="Vm" variable_2="Vm"/>
</connection>
<connection>
<map_components component_1="rapid_time_dependent_potassium_current" component_2="rapid_time_dependent_potassium_current_y3_gate"/>
<map_variables variable_1="y3" variable_2="y3"/>
<map_variables variable_1="E_K" variable_2="E_K"/>
<map_variables variable_1="time" variable_2="time"/>
<map_variables variable_1="Vm" variable_2="Vm"/>
</connection>
<connection>
<map_components component_1="slow_time_dependent_potassium_current" component_2="slow_time_dependent_potassium_current_y1_gate"/>
<map_variables variable_1="y1" variable_2="y1"/>
<map_variables variable_1="time" variable_2="time"/>
<map_variables variable_1="Vm" variable_2="Vm"/>
</connection>
<connection>
<map_components component_1="slow_time_dependent_potassium_current" component_2="slow_time_dependent_potassium_current_y2_gate"/>
<map_variables variable_1="y2" variable_2="y2"/>
<map_variables variable_1="Cai" variable_2="Cai"/>
<map_variables variable_1="time" variable_2="time"/>
<map_variables variable_1="Vm" variable_2="Vm"/>
</connection>
<connection>
<map_components component_1="transient_outward_current" component_2="transient_outward_current_y1_gate"/>
<map_variables variable_1="y1" variable_2="y1"/>
<map_variables variable_1="time" variable_2="time"/>
<map_variables variable_1="Vm" variable_2="Vm"/>
</connection>
<connection>
<map_components component_1="transient_outward_current" component_2="transient_outward_current_y2_gate"/>
<map_variables variable_1="y2" variable_2="y2"/>
<map_variables variable_1="time" variable_2="time"/>
<map_variables variable_1="Vm" variable_2="Vm"/>
</connection>
<connection>
<map_components component_1="sodium_calcium_exchanger" component_2="sodium_calcium_exchanger_y_gate"/>
<map_variables variable_1="y" variable_2="y"/>
<map_variables variable_1="p_E1Na" variable_2="p_E1Na"/>
<map_variables variable_1="p_E2Na" variable_2="p_E2Na"/>
<map_variables variable_1="p_E1Ca" variable_2="p_E1Ca"/>
<map_variables variable_1="p_E2Ca" variable_2="p_E2Ca"/>
<map_variables variable_1="time" variable_2="time"/>
<map_variables variable_1="k1" variable_2="k1"/>
<map_variables variable_1="k2" variable_2="k2"/>
<map_variables variable_1="k3" variable_2="k3"/>
<map_variables variable_1="k4" variable_2="k4"/>
</connection>
<connection>
<map_components component_1="sodium_potassium_pump" component_2="sodium_potassium_pump_y_gate"/>
<map_variables variable_1="y" variable_2="y"/>
<map_variables variable_1="p_E1Na" variable_2="p_E1Na"/>
<map_variables variable_1="p_E2Na" variable_2="p_E2Na"/>
<map_variables variable_1="p_E1K" variable_2="p_E1K"/>
<map_variables variable_1="p_E2K" variable_2="p_E2K"/>
<map_variables variable_1="time" variable_2="time"/>
<map_variables variable_1="k1" variable_2="k1"/>
<map_variables variable_1="k2" variable_2="k2"/>
<map_variables variable_1="k3" variable_2="k3"/>
<map_variables variable_1="k4" variable_2="k4"/>
</connection>
<connection>
<map_components component_1="SR_calcium_pump" component_2="SR_calcium_pump_y_gate"/>
<map_variables variable_1="y" variable_2="y"/>
<map_variables variable_1="p_E1" variable_2="p_E1"/>
<map_variables variable_1="p_E2" variable_2="p_E2"/>
<map_variables variable_1="p_E1Ca" variable_2="p_E1Ca"/>
<map_variables variable_1="p_E2Ca" variable_2="p_E2Ca"/>
<map_variables variable_1="time" variable_2="time"/>
<map_variables variable_1="k1" variable_2="k1"/>
<map_variables variable_1="k2" variable_2="k2"/>
<map_variables variable_1="k3" variable_2="k3"/>
<map_variables variable_1="k4" variable_2="k4"/>
</connection>
<rdf:RDF>
<rdf:Seq rdf:about="rdf:#77d509e8-c85e-472b-9697-99db457aae2b">
<rdf:li rdf:resource="rdf:#d9a2749b-91d9-4fed-8c06-a19e41e42763"/>
<rdf:li rdf:resource="rdf:#3557b5fb-08c5-46eb-bbc7-cba9439e0aec"/>
<rdf:li rdf:resource="rdf:#a7ae93dd-06a4-418d-ae64-ac6eb647f6e4"/>
<rdf:li rdf:resource="rdf:#30320d41-add4-431b-9cfe-5c8fd0f2a8dc"/>
<rdf:li rdf:resource="rdf:#e9bfdb1a-e6f6-42f6-ab2b-d55a1fbfabf5"/>
</rdf:Seq>
<rdf:Description rdf:about="">
<dc:publisher>Department of Physiology, Anatomy & Genetics, University of Oxford</dc:publisher>
<cmeta:comment rdf:resource="rdf:#3c8a3f1f-4766-4b01-abff-d627841ee975"/>
<dcterms:created rdf:resource="rdf:#b8a300ce-43dd-46a1-abae-5763f064bfe0"/>
<dc:creator rdf:resource="rdf:#868a1680-dad6-4412-8506-903d23fd0847"/>
<cmeta:modification rdf:resource="rdf:#36721e11-4f7d-4850-a8ad-c86528fc45ae"/>
<cmeta:modification rdf:resource="rdf:#e80f340c-dfe8-4e5b-b3e1-748987784ec8"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#e30a7d74-6c25-4508-bc30-5c071af3b919">
<dc:title>Japanese Journal of Physiology</dc:title>
</rdf:Description>
<rdf:Description rdf:about="rdf:#3c8a3f1f-4766-4b01-abff-d627841ee975">
<dc:creator rdf:resource="rdf:#c7cf0f32-c9ca-4603-a8e9-3ada722d3404"/>
<rdf:value>This model has been curated by Penny Noble and is known to run in PCEnv and COR and reproduce the published results. </rdf:value>
</rdf:Description>
<rdf:Description rdf:about="rdf:#322b339d-300e-4ff7-8aa8-80c3f754d3b8">
<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:#df111ee2-7a8f-41e1-b5c5-0a5e2719726c">
<vCard:Given>Satoshi</vCard:Given>
<vCard:Family>Matsuoka</vCard:Family>
</rdf:Description>
<rdf:Description rdf:about="rdf:#16c204bc-1ce5-49fc-b563-9beab946d2e2">
<dcterms:W3CDTF>2003-01-01 00:00</dcterms:W3CDTF>
</rdf:Description>
<rdf:Description rdf:about="rdf:#e80f340c-dfe8-4e5b-b3e1-748987784ec8">
<dcterms:modified rdf:resource="rdf:#3df98825-d29b-45e7-82e2-49ed3e0b268d"/>
<rdf:value>units checked, curated</rdf:value>
<cmeta:modifier rdf:resource="rdf:#5f0e939e-987a-4931-8aae-8bffe1663a72"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#$AGvhb2">
<ns7:boundIntervals rdf:resource="rdf:#$BGvhb2"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#66fafa6d-a88a-4a0d-bb05-183f0ba7182c">
<dc:creator rdf:resource="rdf:#99d964b1-a0f9-4011-8cc6-086400e0bec8"/>
<rdf:value>This file is a CellML description of Matsuoka et al's 2003 mathematical model of the roles of individual ion currents in ventricular cells.</rdf:value>
</rdf:Description>
<rdf:Description rdf:about="rdf:#5f0e939e-987a-4931-8aae-8bffe1663a72">
<vCard:N rdf:resource="rdf:#fe9549d0-aafd-415b-b4eb-11482022369f"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#5e8bc23d-a6b6-4a3c-bb62-01b081f4303e">
<bqs:JournalArticle rdf:resource="rdf:#5bdef94f-458c-4979-9f79-e2097e71f85f"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#5bdef94f-458c-4979-9f79-e2097e71f85f">
<dc:creator rdf:resource="rdf:#77d509e8-c85e-472b-9697-99db457aae2b"/>
<dc:title>Role of individual ionic current systems in ventricular cells hypothesized by a model study</dc:title>
<bqs:volume>53(2)</bqs:volume>
<bqs:first_page>105</bqs:first_page>
<bqs:Journal rdf:resource="rdf:#e30a7d74-6c25-4508-bc30-5c071af3b919"/>
<dcterms:issued rdf:resource="rdf:#16c204bc-1ce5-49fc-b563-9beab946d2e2"/>
<bqs:last_page>123</bqs:last_page>
</rdf:Description>
<rdf:Description rdf:about="rdf:#850be265-4cb6-43ee-9eea-d24455089f98">
<vCard:Given>Nobuaki</vCard:Given>
<vCard:Family>Sarai</vCard:Family>
</rdf:Description>
<rdf:Description rdf:about="rdf:#fe9549d0-aafd-415b-b4eb-11482022369f">
<vCard:Given>Penny</vCard:Given>
<vCard:Family>Noble</vCard:Family>
<vCard:Other>J</vCard:Other>
</rdf:Description>
<rdf:Description rdf:about="rdf:#$BGvhb2">
<rdf:rest rdf:resource="http://www.w3.org/1999/02/22-rdf-syntax-ns#nil"/>
<rdf:first rdf:resource="rdf:#$CGvhb2"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#7c4e2c38-c0c8-4fc1-9a60-f97f7a87d374">
<vCard:Given>James</vCard:Given>
<vCard:Family>Lawson</vCard:Family>
<vCard:Other>Richard</vCard:Other>
</rdf:Description>
<rdf:Description rdf:about="rdf:#3557b5fb-08c5-46eb-bbc7-cba9439e0aec">
<rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
<vCard:N rdf:resource="rdf:#850be265-4cb6-43ee-9eea-d24455089f98"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#e9bfdb1a-e6f6-42f6-ab2b-d55a1fbfabf5">
<rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
<vCard:N rdf:resource="rdf:#b7200991-dcd5-4df6-94d8-108c225b6732"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#b7200991-dcd5-4df6-94d8-108c225b6732">
<vCard:Given>Akinori</vCard:Given>
<vCard:Family>Noma</vCard:Family>
</rdf:Description>
<rdf:Description rdf:about="rdf:#99d964b1-a0f9-4011-8cc6-086400e0bec8">
<vCard:FN>James Lawson</vCard:FN>
</rdf:Description>
<rdf:Description rdf:about="rdf:#30320d41-add4-431b-9cfe-5c8fd0f2a8dc">
<rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
<vCard:N rdf:resource="rdf:#ff6605e2-beb8-4a61-a7d3-fd4388150776"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#3df98825-d29b-45e7-82e2-49ed3e0b268d">
<dcterms:W3CDTF>2008-02-25T11:29:40+13:00</dcterms:W3CDTF>
</rdf:Description>
<rdf:Description rdf:about="rdf:#c7cf0f32-c9ca-4603-a8e9-3ada722d3404">
<vCard:FN>James Lawson</vCard:FN>
</rdf:Description>
<rdf:Description rdf:about="rdf:#36721e11-4f7d-4850-a8ad-c86528fc45ae">
<dcterms:modified rdf:resource="rdf:#45b24a15-d4e2-4142-9d6b-e2cfb4151eaa"/>
<rdf:value>Removed initial value for variable Vm_old, removing model overconstraint problem in PCEnv.</rdf:value>
<cmeta:modifier rdf:resource="rdf:#9c178f5d-ba61-4baa-b5e7-bd9b29eed314"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#9c178f5d-ba61-4baa-b5e7-bd9b29eed314">
<vCard:N rdf:resource="rdf:#7c4e2c38-c0c8-4fc1-9a60-f97f7a87d374"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#72927c29-71c4-40b9-b313-d8114a40ac20">
<vCard:Given>Penny</vCard:Given>
<vCard:Family>Noble</vCard:Family>
<vCard:Other>J</vCard:Other>
</rdf:Description>
<rdf:Description rdf:about="rdf:#45b24a15-d4e2-4142-9d6b-e2cfb4151eaa">
<dcterms:W3CDTF>2007-10-30T11:57:22+13:00</dcterms:W3CDTF>
</rdf:Description>
<rdf:Description rdf:about="rdf:#4f1da377-87f0-4808-97f0-869e1359ec11">
<vCard:Orgname>Oxford University</vCard:Orgname>
<vCard:Orgunit>Department of Physiology, Anatomy & Genetics</vCard:Orgunit>
</rdf:Description>
<rdf:Description rdf:about="rdf:#52bb414e-a884-49bd-ad99-233709dc9cf3">
<vCard:Given>Shinobu</vCard:Given>
<vCard:Family>Kuratomi</vCard:Family>
</rdf:Description>
<rdf:Description rdf:about="rdf:#b8a300ce-43dd-46a1-abae-5763f064bfe0">
<dcterms:W3CDTF>2007-10-23T00:00:00+00:00</dcterms:W3CDTF>
</rdf:Description>
<rdf:Description rdf:about="rdf:#$CGvhb2">
<ns7:endingValue>2000</ns7:endingValue>
<ns7:nonstandard-pointDensity>10000</ns7:nonstandard-pointDensity>
<ns7:maximumStepSize>0.1</ns7:maximumStepSize>
</rdf:Description>
<rdf:Description rdf:about="#matsuoka_sarai_kuratomi_ono_noma_2003">
<ns7:simulation rdf:resource="rdf:#$AGvhb2"/>
<bqs:reference rdf:resource="rdf:#5e8bc23d-a6b6-4a3c-bb62-01b081f4303e"/>
<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>ventricular myocyte</rdf:li>
<rdf:li>electrophysiology</rdf:li>
</rdf:Bag>
</rdf:value>
</dc:subject>
</bqs:reference>
<cmeta:comment rdf:resource="rdf:#66fafa6d-a88a-4a0d-bb05-183f0ba7182c"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#868a1680-dad6-4412-8506-903d23fd0847">
<vCard:ORG rdf:resource="rdf:#4f1da377-87f0-4808-97f0-869e1359ec11"/>
<vCard:EMAIL rdf:resource="rdf:#322b339d-300e-4ff7-8aa8-80c3f754d3b8"/>
<vCard:N rdf:resource="rdf:#72927c29-71c4-40b9-b313-d8114a40ac20"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#ff6605e2-beb8-4a61-a7d3-fd4388150776">
<vCard:Given>Kyoichi</vCard:Given>
<vCard:Family>Ono</vCard:Family>
</rdf:Description>
<rdf:Description rdf:about="rdf:#a7ae93dd-06a4-418d-ae64-ac6eb647f6e4">
<rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
<vCard:N rdf:resource="rdf:#52bb414e-a884-49bd-ad99-233709dc9cf3"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#d9a2749b-91d9-4fed-8c06-a19e41e42763">
<rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
<vCard:N rdf:resource="rdf:#df111ee2-7a8f-41e1-b5c5-0a5e2719726c"/>
</rdf:Description>
</rdf:RDF>
</model>
