- Author:
- Alan Garny <agarny@hellix.com>
- Date:
- 2024-07-03 23:04:47+12:00
- Desc:
- Set the correct unit to a couple of variables.
- Permanent Source URI:
- https://models.cellml.org/workspace/698/rawfile/6ae199148333679e941e0896dc08f96f313f8143/fabbri_et_al_based_composite_SAN_model.cellml
<?xml version='1.0' encoding='UTF-8'?>
<model cmeta:id="fabbri_et_al_based_composite_SAN_model" name="fabbri_et_al_based_composite_SAN_model" xmlns="http://www.cellml.org/cellml/1.1#" xmlns:cellml="http://www.cellml.org/cellml/1.1#" xmlns:cmeta="http://www.cellml.org/metadata/1.0#">
<units name="per_millivolt">
<unit exponent="-1.0" prefix="milli" units="volt"/>
</units>
<units name="coulomb_per_mole">
<unit exponent="-1.0" units="mole"/>
<unit units="coulomb"/>
</units>
<units name="microS">
<unit prefix="micro" units="siemens"/>
</units>
<units name="microF">
<unit prefix="micro" units="farad"/>
</units>
<units name="per_millisecond">
<unit exponent="-1.0" units="millisecond"/>
</units>
<units name="joule_per_kilomole_kelvin">
<unit units="joule"/>
<unit exponent="-1.0" units="kelvin"/>
<unit exponent="-1.0" prefix="kilo" units="mole"/>
</units>
<units name="per_second">
<unit exponent="-1.0" units="second"/>
</units>
<units name="millimetre">
<unit prefix="milli" units="metre"/>
</units>
<units name="millimolar">
<unit prefix="milli" units="mole"/>
<unit exponent="-1.0" units="litre"/>
</units>
<units name="millivolt_second">
<unit units="second"/>
<unit units="millivolt"/>
</units>
<units name="nanoA_per_millimolar">
<unit exponent="-1.0" units="millimolar"/>
<unit units="nanoA"/>
</units>
<units name="mA_nA">
<unit prefix="milli" units="ampere"/>
<unit exponent="-1.0" units="nanoA"/>
</units>
<units name="per_millivolt_second">
<unit exponent="-1.0" units="millivolt"/>
<unit exponent="-1.0" units="second"/>
</units>
<units name="millivolt">
<unit prefix="milli" units="volt"/>
</units>
<units name="per_millimolar_second">
<unit exponent="-1.0" units="millimolar"/>
<unit exponent="-1.0" units="second"/>
</units>
<units name="millisecond">
<unit prefix="milli" units="second"/>
</units>
<units name="litre_millimetre3">
<unit exponent="-1.0" units="millimetre3"/>
<unit units="litre"/>
</units>
<units name="millimetre3">
<unit exponent="3.0" prefix="milli" units="metre"/>
</units>
<units name="nanoA">
<unit prefix="nano" units="ampere"/>
</units>
<units name="micrometre">
<unit prefix="micro" units="metre"/>
</units>
<units name="nanoS">
<unit prefix="nano" units="siemens"/>
</units>
<units name="millimolar_per_second">
<unit exponent="-1.0" units="second"/>
<unit units="millimolar"/>
</units>
<units name="per_millimolar">
<unit exponent="-1.0" units="millimolar"/>
</units>
<units name="millimolar4">
<unit exponent="4.0" units="millimolar"/>
</units>
<units name="per_millimolar2_second">
<unit exponent="-2.0" units="millimolar"/>
<unit exponent="-1.0" units="second"/>
</units>
<units name="micrometre_per_femtometre">
<unit exponent="-1.0" prefix="femto" units="metre"/>
<unit units="micrometre"/>
</units>
<component name="environment">
<variable cmeta:id="environment.time" name="time" public_interface="out" units="second"/>
</component>
<component name="Brain_stem">
<!-- See https://brian2.readthedocs.io/en/stable/user/converting_from_integrated_form.html.-->
<variable name="t" public_interface="in" units="millisecond"/>
<variable initial_value="300" name="t_period" units="millisecond"/>
<variable name="t_relative" units="millisecond"/>
<variable name="t_w_relative" units="millisecond"/>
<variable initial_value="0" name="g" public_interface="out" units="nanoA"/>
<variable initial_value="10" name="tau" units="millisecond"/>
<variable name="w" units="nanoA"/>
<variable initial_value="13" name="w_start" units="millisecond"/>
<variable initial_value="1" name="w_duration" units="millisecond"/>
<variable initial_value="3" name="w_period" units="millisecond"/>
<variable initial_value="10" name="w_n" units="dimensionless"/>
<variable initial_value="0.1" name="w_value" units="nanoA"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>t_relative</ci>
<apply>
<minus/>
<apply>
<rem/>
<ci>t</ci>
<ci>t_period</ci>
</apply>
<ci>w_start</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>t_w_relative</ci>
<apply>
<rem/>
<ci>t_relative</ci>
<ci>w_period</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>w</ci>
<piecewise>
<piece>
<ci>w_value</ci>
<apply>
<and/>
<apply>
<geq/>
<ci>t_relative</ci>
<cn cellml:units="millisecond">0</cn>
</apply>
<apply>
<lt/>
<ci>t_relative</ci>
<apply>
<times/>
<ci>w_n</ci>
<ci>w_period</ci>
</apply>
</apply>
<apply>
<lt/>
<ci>t_w_relative</ci>
<ci>w_duration</ci>
</apply>
</apply>
</piece>
<otherwise>
<cn cellml:units="nanoA">0</cn>
</otherwise>
</piecewise>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>g</ci>
</apply>
<apply>
<divide/>
<apply>
<plus/>
<apply>
<minus/>
<ci>g</ci>
</apply>
<ci>w</ci>
</apply>
<ci>tau</ci>
</apply>
</apply>
</math>
</component>
<component name="Membrane">
<variable cmeta:id="Membrane.R" initial_value="8314.472" name="R" public_interface="out" units="joule_per_kilomole_kelvin"/>
<variable cmeta:id="Membrane.T" initial_value="310" name="T" public_interface="out" units="kelvin"/>
<variable cmeta:id="Membrane.F" initial_value="96485.3415" name="F" public_interface="out" units="coulomb_per_mole"/>
<variable cmeta:id="Membrane.C" initial_value="5.7e-5" name="C" units="microF"/>
<variable cmeta:id="Membrane.RTONF" name="RTONF" public_interface="out" units="millivolt"/>
<variable cmeta:id="Membrane.time" name="time" private_interface="out" public_interface="in" units="second"/>
<variable cmeta:id="Membrane.i_f" name="i_f" public_interface="in" units="nanoA"/>
<variable cmeta:id="Membrane.i_NaK" name="i_NaK" public_interface="in" units="nanoA"/>
<variable cmeta:id="Membrane.i_NaCa" name="i_NaCa" public_interface="in" units="nanoA"/>
<variable cmeta:id="Membrane.i_Na" name="i_Na" public_interface="in" units="nanoA"/>
<variable cmeta:id="Membrane.i_Kr" name="i_Kr" public_interface="in" units="nanoA"/>
<variable cmeta:id="Membrane.i_Ks" name="i_Ks" public_interface="in" units="nanoA"/>
<variable cmeta:id="Membrane.i_to" name="i_to" public_interface="in" units="nanoA"/>
<variable cmeta:id="Membrane.i_CaL" name="i_CaL" public_interface="in" units="nanoA"/>
<variable cmeta:id="Membrane.i_CaT" name="i_CaT" public_interface="in" units="nanoA"/>
<variable cmeta:id="Membrane.i_KACh" name="i_KACh" public_interface="in" units="nanoA"/>
<variable cmeta:id="Membrane.i_Kur" name="i_Kur" public_interface="in" units="nanoA"/>
<variable cmeta:id="Membrane.V" name="V" public_interface="out" units="millivolt"/>
<variable cmeta:id="Membrane.clamp_mode" initial_value="0" name="clamp_mode" units="dimensionless"/>
<variable cmeta:id="Membrane.V_clamp" name="V_clamp" private_interface="in" units="millivolt"/>
<variable cmeta:id="Membrane.V_ode" initial_value="-47.787168" name="V_ode" units="millivolt"/>
<variable cmeta:id="Membrane.i_tot" name="i_tot" units="nanoA"/>
<variable name="g" public_interface="in" units="nanoA"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>i_tot</ci>
<apply>
<minus/>
<apply>
<plus/>
<ci>i_f</ci>
<ci>i_Kr</ci>
<ci>i_Ks</ci>
<ci>i_to</ci>
<ci>i_NaK</ci>
<ci>i_NaCa</ci>
<ci>i_Na</ci>
<ci>i_CaL</ci>
<ci>i_CaT</ci>
<ci>i_KACh</ci>
<ci>i_Kur</ci>
</apply>
<ci>g</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>V_ode</ci>
</apply>
<apply>
<divide/>
<apply>
<minus/>
<ci>i_tot</ci>
</apply>
<ci>C</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>RTONF</ci>
<apply>
<divide/>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
<ci>F</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>V</ci>
<piecewise>
<piece>
<ci>V_clamp</ci>
<apply>
<geq/>
<ci>clamp_mode</ci>
<cn cellml:units="dimensionless">1</cn>
</apply>
</piece>
<otherwise>
<ci>V_ode</ci>
</otherwise>
</piecewise>
</apply>
</math>
</component>
<component name="Voltage_clamp">
<variable cmeta:id="Voltage_clamp.t_holding" initial_value="0.5" name="t_holding" units="second"/>
<variable cmeta:id="Voltage_clamp.t_test" initial_value="0.5" name="t_test" units="second"/>
<variable cmeta:id="Voltage_clamp.V_clamp" name="V_clamp" public_interface="out" units="millivolt"/>
<variable cmeta:id="Voltage_clamp.V_test" initial_value="-35" name="V_test" units="millivolt"/>
<variable cmeta:id="Voltage_clamp.V_holding" initial_value="-45" name="V_holding" units="millivolt"/>
<variable cmeta:id="Voltage_clamp.time" name="time" public_interface="in" units="second"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>V_clamp</ci>
<piecewise>
<piece>
<ci>V_test</ci>
<apply>
<and/>
<apply>
<gt/>
<ci>time</ci>
<ci>t_holding</ci>
</apply>
<apply>
<lt/>
<ci>time</ci>
<apply>
<plus/>
<ci>t_holding</ci>
<ci>t_test</ci>
</apply>
</apply>
</apply>
</piece>
<otherwise>
<ci>V_holding</ci>
</otherwise>
</piecewise>
</apply>
</math>
</component>
<component name="Rate_modulation_experiments">
<variable cmeta:id="Rate_modulation_experiments.ACh" initial_value="0" name="ACh" public_interface="out" units="millimolar"/>
<variable cmeta:id="Rate_modulation_experiments.Iso_1_uM" initial_value="0" name="Iso_1_uM" public_interface="out" units="dimensionless"/>
</component>
<component name="Ionic_values">
<variable cmeta:id="Ionic_values.Nai" name="Nai" public_interface="in" units="millimolar"/>
<variable cmeta:id="Ionic_values.Nao" initial_value="140" name="Nao" public_interface="out" units="millimolar"/>
<variable cmeta:id="Ionic_values.Ki" initial_value="140" name="Ki" public_interface="out" units="millimolar"/>
<variable cmeta:id="Ionic_values.Ko" initial_value="5.4" name="Ko" public_interface="out" units="millimolar"/>
<variable cmeta:id="Ionic_values.Ca_sub" name="Ca_sub" public_interface="in" units="millimolar"/>
<variable cmeta:id="Ionic_values.Cao" initial_value="1.8" name="Cao" public_interface="out" units="millimolar"/>
<variable cmeta:id="Ionic_values.RTONF" name="RTONF" public_interface="in" units="millivolt"/>
<variable cmeta:id="Ionic_values.E_Na" name="E_Na" public_interface="out" units="millivolt"/>
<variable cmeta:id="Ionic_values.E_K" name="E_K" public_interface="out" units="millivolt"/>
<variable cmeta:id="Ionic_values.E_Ca" name="E_Ca" public_interface="out" units="millivolt"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>E_Na</ci>
<apply>
<times/>
<ci>RTONF</ci>
<apply>
<ln/>
<apply>
<divide/>
<ci>Nao</ci>
<ci>Nai</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>E_K</ci>
<apply>
<times/>
<ci>RTONF</ci>
<apply>
<ln/>
<apply>
<divide/>
<ci>Ko</ci>
<ci>Ki</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>E_Ca</ci>
<apply>
<times/>
<cn cellml:units="dimensionless">0.5</cn>
<ci>RTONF</ci>
<apply>
<ln/>
<apply>
<divide/>
<ci>Cao</ci>
<ci>Ca_sub</ci>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<component name="Nai_concentration">
<variable cmeta:id="Nai_concentration.Nai" name="Nai" public_interface="out" units="millimolar"/>
<variable cmeta:id="Nai_concentration.time" name="time" public_interface="in" units="second"/>
<variable cmeta:id="Nai_concentration.F" name="F" public_interface="in" units="coulomb_per_mole"/>
<variable cmeta:id="Nai_concentration.V_sub" name="V_sub" public_interface="in" units="millimetre3"/>
<variable cmeta:id="Nai_concentration.V_i" name="V_i" public_interface="in" units="millimetre3"/>
<variable cmeta:id="Nai_concentration.i_Na" name="i_Na" public_interface="in" units="nanoA"/>
<variable cmeta:id="Nai_concentration.i_NaCa" name="i_NaCa" public_interface="in" units="nanoA"/>
<variable cmeta:id="Nai_concentration.i_NaK" name="i_NaK" public_interface="in" units="nanoA"/>
<variable cmeta:id="Nai_concentration.i_fNa" name="i_fNa" public_interface="in" units="nanoA"/>
<variable cmeta:id="Nai_concentration.i_siNa" name="i_siNa" public_interface="in" units="nanoA"/>
<variable cmeta:id="Nai_concentration.Nai_" initial_value="5" name="Nai_" units="millimolar"/>
<variable cmeta:id="Nai_concentration.Nai_clamp" initial_value="1" name="Nai_clamp" units="dimensionless"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>Nai_</ci>
</apply>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<ci>Nai_clamp</ci>
</apply>
<apply>
<minus/>
<cn cellml:units="mA_nA">1</cn>
</apply>
<apply>
<plus/>
<ci>i_Na</ci>
<ci>i_fNa</ci>
<ci>i_siNa</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>
<times/>
<cn cellml:units="litre_millimetre3">1</cn>
<apply>
<plus/>
<ci>V_i</ci>
<ci>V_sub</ci>
</apply>
<ci>F</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>Nai</ci>
<ci>Nai_</ci>
</apply>
</math>
</component>
<component name="i_f">
<variable cmeta:id="i_f.ACh" name="ACh" private_interface="out" public_interface="in" units="millimolar"/>
<variable cmeta:id="i_f.Iso_1_uM" name="Iso_1_uM" private_interface="out" public_interface="in" units="dimensionless"/>
<variable cmeta:id="i_f.i_f" name="i_f" public_interface="out" units="nanoA"/>
<variable cmeta:id="i_f.i_fNa" name="i_fNa" public_interface="out" units="nanoA"/>
<variable cmeta:id="i_f.E_Na" name="E_Na" public_interface="in" units="millivolt"/>
<variable cmeta:id="i_f.E_K" name="E_K" public_interface="in" units="millivolt"/>
<variable cmeta:id="i_f.i_fK" name="i_fK" public_interface="out" units="nanoA"/>
<variable cmeta:id="i_f.g_f" initial_value="0.00427" name="g_f" units="microS"/>
<variable cmeta:id="i_f.G_f" name="G_f" units="microS"/>
<variable cmeta:id="i_f.g_f_Na" name="g_f_Na" public_interface="out" units="microS"/>
<variable cmeta:id="i_f.G_f_Na" name="G_f_Na" units="microS"/>
<variable cmeta:id="i_f.g_f_K" name="g_f_K" public_interface="out" units="microS"/>
<variable cmeta:id="i_f.G_f_K" name="G_f_K" units="microS"/>
<variable cmeta:id="i_f.Km_f" initial_value="45" name="Km_f" units="millimolar"/>
<variable cmeta:id="i_f.alpha" initial_value="0.5927" name="alpha" units="dimensionless"/>
<variable cmeta:id="i_f.time" name="time" private_interface="out" public_interface="in" units="second"/>
<variable cmeta:id="i_f.V" name="V" private_interface="out" public_interface="in" units="millivolt"/>
<variable cmeta:id="i_f.RTONF" name="RTONF" public_interface="in" units="millivolt"/>
<variable cmeta:id="i_f.Ko" name="Ko" public_interface="in" units="millimolar"/>
<variable cmeta:id="i_f.y" name="y" private_interface="in" units="dimensionless"/>
<variable cmeta:id="i_f.blockade" initial_value="0" name="blockade" units="dimensionless"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>G_f</ci>
<apply>
<divide/>
<ci>g_f</ci>
<apply>
<divide/>
<ci>Ko</ci>
<apply>
<plus/>
<ci>Ko</ci>
<ci>Km_f</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>G_f_K</ci>
<apply>
<divide/>
<ci>G_f</ci>
<apply>
<plus/>
<ci>alpha</ci>
<cn cellml:units="dimensionless">1</cn>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>G_f_Na</ci>
<apply>
<times/>
<ci>alpha</ci>
<ci>G_f_K</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>g_f_Na</ci>
<apply>
<divide/>
<apply>
<times/>
<ci>G_f_Na</ci>
<ci>Ko</ci>
</apply>
<apply>
<plus/>
<ci>Ko</ci>
<ci>Km_f</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>g_f_K</ci>
<apply>
<divide/>
<apply>
<times/>
<ci>G_f_K</ci>
<ci>Ko</ci>
</apply>
<apply>
<plus/>
<ci>Ko</ci>
<ci>Km_f</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>i_fNa</ci>
<apply>
<times/>
<ci>y</ci>
<ci>g_f_Na</ci>
<apply>
<minus/>
<ci>V</ci>
<ci>E_Na</ci>
</apply>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<ci>blockade</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>i_fK</ci>
<apply>
<times/>
<ci>y</ci>
<ci>g_f_K</ci>
<apply>
<minus/>
<ci>V</ci>
<ci>E_K</ci>
</apply>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<ci>blockade</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>i_f</ci>
<apply>
<plus/>
<ci>i_fNa</ci>
<ci>i_fK</ci>
</apply>
</apply>
</math>
</component>
<component name="i_f_y_gate">
<variable cmeta:id="i_f_y_gate.y" initial_value="0.009508" name="y" public_interface="out" units="dimensionless"/>
<variable cmeta:id="i_f_y_gate.tau_y" name="tau_y" units="second"/>
<variable cmeta:id="i_f_y_gate.y_infinity" name="y_infinity" units="dimensionless"/>
<variable cmeta:id="i_f_y_gate.V" name="V" public_interface="in" units="millivolt"/>
<variable cmeta:id="i_f_y_gate.time" name="time" public_interface="in" units="second"/>
<variable cmeta:id="i_f_y_gate.ACh_shift" name="ACh_shift" units="millivolt"/>
<variable cmeta:id="i_f_y_gate.ACh" name="ACh" public_interface="in" units="millimolar"/>
<variable cmeta:id="i_f_y_gate.Iso_1_uM" name="Iso_1_uM" public_interface="in" units="dimensionless"/>
<variable cmeta:id="i_f_y_gate.Iso_shift" name="Iso_shift" units="millivolt"/>
<variable cmeta:id="i_f_y_gate.y_shift" initial_value="0" name="y_shift" units="millivolt"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>ACh_shift</ci>
<piecewise>
<piece>
<apply>
<minus/>
<apply>
<minus/>
<cn cellml:units="millivolt">1</cn>
</apply>
<apply>
<divide/>
<apply>
<times/>
<cn cellml:units="millivolt">9.898</cn>
<apply>
<power/>
<apply>
<times/>
<cn cellml:units="per_millimolar">1</cn>
<ci>ACh</ci>
</apply>
<cn cellml:units="dimensionless">0.618</cn>
</apply>
</apply>
<apply>
<plus/>
<apply>
<power/>
<apply>
<times/>
<cn cellml:units="per_millimolar">1</cn>
<ci>ACh</ci>
</apply>
<cn cellml:units="dimensionless">0.618</cn>
</apply>
<cn cellml:units="dimensionless">0.00122423</cn>
</apply>
</apply>
</apply>
<apply>
<gt/>
<ci>ACh</ci>
<cn cellml:units="millimolar">0</cn>
</apply>
</piece>
<otherwise>
<cn cellml:units="millivolt">0</cn>
</otherwise>
</piecewise>
</apply>
<apply>
<eq/>
<ci>Iso_shift</ci>
<piecewise>
<piece>
<cn cellml:units="millivolt">7.5</cn>
<apply>
<gt/>
<ci>Iso_1_uM</ci>
<cn cellml:units="dimensionless">0</cn>
</apply>
</piece>
<otherwise>
<cn cellml:units="millivolt">0</cn>
</otherwise>
</piecewise>
</apply>
<apply>
<eq/>
<ci>tau_y</ci>
<apply>
<minus/>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<plus/>
<apply>
<divide/>
<apply>
<times/>
<cn cellml:units="per_millivolt_second">0.36</cn>
<apply>
<minus/>
<apply>
<minus/>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="millivolt">148.8</cn>
</apply>
<ci>ACh_shift</ci>
</apply>
<ci>Iso_shift</ci>
</apply>
</apply>
<apply>
<minus/>
<apply>
<exp/>
<apply>
<times/>
<cn cellml:units="per_millivolt">0.066</cn>
<apply>
<minus/>
<apply>
<minus/>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="millivolt">148.8</cn>
</apply>
<ci>ACh_shift</ci>
</apply>
<ci>Iso_shift</ci>
</apply>
</apply>
</apply>
<cn cellml:units="dimensionless">1</cn>
</apply>
</apply>
<apply>
<divide/>
<apply>
<times/>
<cn cellml:units="per_millivolt_second">0.1</cn>
<apply>
<minus/>
<apply>
<minus/>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="millivolt">87.3</cn>
</apply>
<ci>ACh_shift</ci>
</apply>
<ci>Iso_shift</ci>
</apply>
</apply>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="per_millivolt">0.2</cn>
</apply>
<apply>
<minus/>
<apply>
<minus/>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="millivolt">87.3</cn>
</apply>
<ci>ACh_shift</ci>
</apply>
<ci>Iso_shift</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<cn cellml:units="second">0.054</cn>
</apply>
</apply>
<apply>
<eq/>
<ci>y_infinity</ci>
<piecewise>
<piece>
<apply>
<plus/>
<cn cellml:units="dimensionless">0.01329</cn>
<apply>
<divide/>
<cn cellml:units="dimensionless">0.99921</cn>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<minus/>
<apply>
<minus/>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="millivolt">97.134</cn>
</apply>
<ci>ACh_shift</ci>
</apply>
<ci>Iso_shift</ci>
</apply>
<ci>y_shift</ci>
</apply>
<cn cellml:units="millivolt">8.1752</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<lt/>
<ci>V</ci>
<apply>
<minus/>
<apply>
<minus/>
<apply>
<minus/>
<apply>
<minus/>
<cn cellml:units="millivolt">80</cn>
<ci>ACh_shift</ci>
</apply>
<ci>Iso_shift</ci>
</apply>
<ci>y_shift</ci>
</apply>
</apply>
</apply>
</piece>
<otherwise>
<apply>
<times/>
<cn cellml:units="dimensionless">0.0002501</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<minus/>
<apply>
<minus/>
<apply>
<minus/>
<ci>V</ci>
<ci>ACh_shift</ci>
</apply>
<ci>Iso_shift</ci>
</apply>
<ci>y_shift</ci>
</apply>
</apply>
<cn cellml:units="millivolt">12.861</cn>
</apply>
</apply>
</apply>
</otherwise>
</piecewise>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>y</ci>
</apply>
<apply>
<divide/>
<apply>
<minus/>
<ci>y_infinity</ci>
<ci>y</ci>
</apply>
<ci>tau_y</ci>
</apply>
</apply>
</math>
</component>
<component name="i_NaK">
<variable cmeta:id="i_NaK.i_NaK" name="i_NaK" public_interface="out" units="nanoA"/>
<variable cmeta:id="i_NaK.time" name="time" public_interface="in" units="second"/>
<variable cmeta:id="i_NaK.V" name="V" public_interface="in" units="millivolt"/>
<variable cmeta:id="i_NaK.Km_Kp" initial_value="1.4" name="Km_Kp" units="millimolar"/>
<variable cmeta:id="i_NaK.Km_Nap" initial_value="14" name="Km_Nap" units="millimolar"/>
<variable cmeta:id="i_NaK.i_NaK_max" initial_value="0.08105" name="i_NaK_max" units="nanoA"/>
<variable cmeta:id="i_NaK.Nai" name="Nai" public_interface="in" units="millimolar"/>
<variable cmeta:id="i_NaK.Ko" name="Ko" public_interface="in" units="millimolar"/>
<variable cmeta:id="i_NaK.E_Na" name="E_Na" public_interface="in" units="millivolt"/>
<variable cmeta:id="i_NaK.Iso_increase" name="Iso_increase" units="dimensionless"/>
<variable cmeta:id="i_NaK.Iso_1_uM" name="Iso_1_uM" private_interface="out" public_interface="in" units="dimensionless"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>Iso_increase</ci>
<piecewise>
<piece>
<cn cellml:units="dimensionless">1.2</cn>
<apply>
<gt/>
<ci>Iso_1_uM</ci>
<cn cellml:units="dimensionless">0</cn>
</apply>
</piece>
<otherwise>
<cn cellml:units="dimensionless">1</cn>
</otherwise>
</piecewise>
</apply>
<apply>
<eq/>
<ci>i_NaK</ci>
<apply>
<times/>
<ci>Iso_increase</ci>
<ci>i_NaK_max</ci>
<apply>
<power/>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<power/>
<apply>
<divide/>
<ci>Km_Kp</ci>
<ci>Ko</ci>
</apply>
<cn cellml:units="dimensionless">1.2</cn>
</apply>
</apply>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
</apply>
</apply>
<apply>
<power/>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<power/>
<apply>
<divide/>
<ci>Km_Nap</ci>
<ci>Nai</ci>
</apply>
<cn cellml:units="dimensionless">1.3</cn>
</apply>
</apply>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
</apply>
</apply>
<apply>
<power/>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<plus/>
<apply>
<minus/>
<ci>V</ci>
<ci>E_Na</ci>
</apply>
<cn cellml:units="millivolt">110</cn>
</apply>
</apply>
<cn cellml:units="millivolt">20</cn>
</apply>
</apply>
</apply>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<component name="i_NaCa">
<variable cmeta:id="i_NaCa.i_NaCa" name="i_NaCa" public_interface="out" units="nanoA"/>
<variable cmeta:id="i_NaCa.K_NaCa" initial_value="3.343" name="K_NaCa" units="nanoA"/>
<variable cmeta:id="i_NaCa.x1" name="x1" units="dimensionless"/>
<variable cmeta:id="i_NaCa.x2" name="x2" units="dimensionless"/>
<variable cmeta:id="i_NaCa.x3" name="x3" units="dimensionless"/>
<variable cmeta:id="i_NaCa.x4" name="x4" units="dimensionless"/>
<variable cmeta:id="i_NaCa.k41" name="k41" units="dimensionless"/>
<variable cmeta:id="i_NaCa.k34" name="k34" units="dimensionless"/>
<variable cmeta:id="i_NaCa.k23" name="k23" units="dimensionless"/>
<variable cmeta:id="i_NaCa.k21" name="k21" units="dimensionless"/>
<variable cmeta:id="i_NaCa.k32" name="k32" units="dimensionless"/>
<variable cmeta:id="i_NaCa.k43" name="k43" units="dimensionless"/>
<variable cmeta:id="i_NaCa.k12" name="k12" units="dimensionless"/>
<variable cmeta:id="i_NaCa.k14" name="k14" units="dimensionless"/>
<variable cmeta:id="i_NaCa.Qci" initial_value="0.1369" name="Qci" units="dimensionless"/>
<variable cmeta:id="i_NaCa.Qn" initial_value="0.4315" name="Qn" units="dimensionless"/>
<variable cmeta:id="i_NaCa.Qco" initial_value="0" name="Qco" units="dimensionless"/>
<variable cmeta:id="i_NaCa.K3ni" initial_value="26.44" name="K3ni" units="millimolar"/>
<variable cmeta:id="i_NaCa.Kci" initial_value="0.0207" name="Kci" units="millimolar"/>
<variable cmeta:id="i_NaCa.K1ni" initial_value="395.3" name="K1ni" units="millimolar"/>
<variable cmeta:id="i_NaCa.K2ni" initial_value="2.289" name="K2ni" units="millimolar"/>
<variable cmeta:id="i_NaCa.Kcni" initial_value="26.44" name="Kcni" units="millimolar"/>
<variable cmeta:id="i_NaCa.K3no" initial_value="4.663" name="K3no" units="millimolar"/>
<variable cmeta:id="i_NaCa.K1no" initial_value="1628" name="K1no" units="millimolar"/>
<variable cmeta:id="i_NaCa.K2no" initial_value="561.4" name="K2no" units="millimolar"/>
<variable cmeta:id="i_NaCa.Kco" initial_value="3.663" name="Kco" units="millimolar"/>
<variable cmeta:id="i_NaCa.do" name="do" units="dimensionless"/>
<variable cmeta:id="i_NaCa.di" name="di" units="dimensionless"/>
<variable cmeta:id="i_NaCa.Cao" name="Cao" public_interface="in" units="millimolar"/>
<variable cmeta:id="i_NaCa.Nao" name="Nao" public_interface="in" units="millimolar"/>
<variable cmeta:id="i_NaCa.Ca_sub" name="Ca_sub" public_interface="in" units="millimolar"/>
<variable cmeta:id="i_NaCa.Nai" name="Nai" public_interface="in" units="millimolar"/>
<variable cmeta:id="i_NaCa.time" name="time" public_interface="in" units="second"/>
<variable cmeta:id="i_NaCa.V" name="V" public_interface="in" units="millivolt"/>
<variable cmeta:id="i_NaCa.RTONF" name="RTONF" public_interface="in" units="millivolt"/>
<variable cmeta:id="i_NaCa.blockade_NaCa" initial_value="0" name="blockade_NaCa" units="dimensionless"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>i_NaCa</ci>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<ci>blockade_NaCa</ci>
</apply>
<ci>K_NaCa</ci>
<apply>
<minus/>
<apply>
<times/>
<ci>x2</ci>
<ci>k21</ci>
</apply>
<apply>
<times/>
<ci>x1</ci>
<ci>k12</ci>
</apply>
</apply>
</apply>
<apply>
<plus/>
<ci>x1</ci>
<ci>x2</ci>
<ci>x3</ci>
<ci>x4</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>x1</ci>
<apply>
<plus/>
<apply>
<times/>
<ci>k41</ci>
<ci>k34</ci>
<apply>
<plus/>
<ci>k23</ci>
<ci>k21</ci>
</apply>
</apply>
<apply>
<times/>
<ci>k21</ci>
<ci>k32</ci>
<apply>
<plus/>
<ci>k43</ci>
<ci>k41</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>x2</ci>
<apply>
<plus/>
<apply>
<times/>
<ci>k32</ci>
<ci>k43</ci>
<apply>
<plus/>
<ci>k14</ci>
<ci>k12</ci>
</apply>
</apply>
<apply>
<times/>
<ci>k41</ci>
<ci>k12</ci>
<apply>
<plus/>
<ci>k34</ci>
<ci>k32</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>x3</ci>
<apply>
<plus/>
<apply>
<times/>
<ci>k14</ci>
<ci>k43</ci>
<apply>
<plus/>
<ci>k23</ci>
<ci>k21</ci>
</apply>
</apply>
<apply>
<times/>
<ci>k12</ci>
<ci>k23</ci>
<apply>
<plus/>
<ci>k43</ci>
<ci>k41</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>x4</ci>
<apply>
<plus/>
<apply>
<times/>
<ci>k23</ci>
<ci>k34</ci>
<apply>
<plus/>
<ci>k14</ci>
<ci>k12</ci>
</apply>
</apply>
<apply>
<times/>
<ci>k14</ci>
<ci>k21</ci>
<apply>
<plus/>
<ci>k34</ci>
<ci>k32</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>k43</ci>
<apply>
<divide/>
<ci>Nai</ci>
<apply>
<plus/>
<ci>K3ni</ci>
<ci>Nai</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>k12</ci>
<apply>
<divide/>
<apply>
<times/>
<apply>
<divide/>
<ci>Ca_sub</ci>
<ci>Kci</ci>
</apply>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<ci>Qci</ci>
</apply>
<ci>V</ci>
</apply>
<ci>RTONF</ci>
</apply>
</apply>
</apply>
<ci>di</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>k14</ci>
<apply>
<divide/>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<divide/>
<ci>Nai</ci>
<ci>K1ni</ci>
</apply>
<ci>Nai</ci>
</apply>
<ci>K2ni</ci>
</apply>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<divide/>
<ci>Nai</ci>
<ci>K3ni</ci>
</apply>
</apply>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<ci>Qn</ci>
<ci>V</ci>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">2</cn>
<ci>RTONF</ci>
</apply>
</apply>
</apply>
</apply>
<ci>di</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>k41</ci>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<ci>Qn</ci>
</apply>
<ci>V</ci>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">2</cn>
<ci>RTONF</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>di</ci>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<times/>
<apply>
<divide/>
<ci>Ca_sub</ci>
<ci>Kci</ci>
</apply>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<ci>Qci</ci>
</apply>
<ci>V</ci>
</apply>
<ci>RTONF</ci>
</apply>
</apply>
<apply>
<divide/>
<ci>Nai</ci>
<ci>Kcni</ci>
</apply>
</apply>
</apply>
<apply>
<times/>
<apply>
<divide/>
<ci>Nai</ci>
<ci>K1ni</ci>
</apply>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<times/>
<apply>
<divide/>
<ci>Nai</ci>
<ci>K2ni</ci>
</apply>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<divide/>
<ci>Nai</ci>
<ci>K3ni</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>k34</ci>
<apply>
<divide/>
<ci>Nao</ci>
<apply>
<plus/>
<ci>K3no</ci>
<ci>Nao</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>k21</ci>
<apply>
<divide/>
<apply>
<times/>
<apply>
<divide/>
<ci>Cao</ci>
<ci>Kco</ci>
</apply>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<ci>Qco</ci>
<ci>V</ci>
</apply>
<ci>RTONF</ci>
</apply>
</apply>
</apply>
<ci>do</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>k23</ci>
<apply>
<divide/>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<divide/>
<ci>Nao</ci>
<ci>K1no</ci>
</apply>
<ci>Nao</ci>
</apply>
<ci>K2no</ci>
</apply>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<divide/>
<ci>Nao</ci>
<ci>K3no</ci>
</apply>
</apply>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<ci>Qn</ci>
</apply>
<ci>V</ci>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">2</cn>
<ci>RTONF</ci>
</apply>
</apply>
</apply>
</apply>
<ci>do</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>k32</ci>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<ci>Qn</ci>
<ci>V</ci>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">2</cn>
<ci>RTONF</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>do</ci>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<times/>
<apply>
<divide/>
<ci>Cao</ci>
<ci>Kco</ci>
</apply>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<ci>Qco</ci>
<ci>V</ci>
</apply>
<ci>RTONF</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<times/>
<apply>
<divide/>
<ci>Nao</ci>
<ci>K1no</ci>
</apply>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<times/>
<apply>
<divide/>
<ci>Nao</ci>
<ci>K2no</ci>
</apply>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<divide/>
<ci>Nao</ci>
<ci>K3no</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<component name="i_Na">
<variable cmeta:id="i_Na.i_Na_" name="i_Na_" public_interface="out" units="nanoA"/>
<variable cmeta:id="i_Na.i_Na_L" name="i_Na_L" public_interface="out" units="nanoA"/>
<variable cmeta:id="i_Na.i_Na" name="i_Na" public_interface="out" units="nanoA"/>
<variable cmeta:id="i_Na.g_Na" initial_value="0.0223" name="g_Na" units="microS"/>
<variable cmeta:id="i_Na.g_Na_L" initial_value="0" name="g_Na_L" units="microS"/>
<variable cmeta:id="i_Na.E_mh" name="E_mh" units="millivolt"/>
<variable cmeta:id="i_Na.time" name="time" private_interface="out" public_interface="in" units="second"/>
<variable cmeta:id="i_Na.V" name="V" private_interface="out" public_interface="in" units="millivolt"/>
<variable cmeta:id="i_Na.RTONF" name="RTONF" public_interface="in" units="millivolt"/>
<variable cmeta:id="i_Na.Nao" name="Nao" public_interface="in" units="millimolar"/>
<variable cmeta:id="i_Na.Nai" name="Nai" public_interface="in" units="millimolar"/>
<variable cmeta:id="i_Na.Ko" name="Ko" public_interface="in" units="millimolar"/>
<variable cmeta:id="i_Na.Ki" name="Ki" public_interface="in" units="millimolar"/>
<variable cmeta:id="i_Na.m" name="m" private_interface="in" units="dimensionless"/>
<variable cmeta:id="i_Na.h" name="h" private_interface="in" units="dimensionless"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>E_mh</ci>
<apply>
<times/>
<ci>RTONF</ci>
<apply>
<ln/>
<apply>
<divide/>
<apply>
<plus/>
<ci>Nao</ci>
<apply>
<times/>
<cn cellml:units="dimensionless">0.12</cn>
<ci>Ko</ci>
</apply>
</apply>
<apply>
<plus/>
<ci>Nai</ci>
<apply>
<times/>
<cn cellml:units="dimensionless">0.12</cn>
<ci>Ki</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>i_Na_</ci>
<apply>
<times/>
<ci>g_Na</ci>
<apply>
<power/>
<ci>m</ci>
<cn cellml:units="dimensionless">3</cn>
</apply>
<ci>h</ci>
<apply>
<minus/>
<ci>V</ci>
<ci>E_mh</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>i_Na_L</ci>
<apply>
<times/>
<ci>g_Na_L</ci>
<apply>
<power/>
<ci>m</ci>
<cn cellml:units="dimensionless">3</cn>
</apply>
<apply>
<minus/>
<ci>V</ci>
<ci>E_mh</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>i_Na</ci>
<apply>
<plus/>
<ci>i_Na_</ci>
<ci>i_Na_L</ci>
</apply>
</apply>
</math>
</component>
<component name="i_Na_m_gate">
<variable cmeta:id="i_Na_m_gate.m" initial_value="0.447724" name="m" public_interface="out" units="dimensionless"/>
<variable cmeta:id="i_Na_m_gate.alpha_m" name="alpha_m" units="per_second"/>
<variable cmeta:id="i_Na_m_gate.beta_m" name="beta_m" units="per_second"/>
<variable cmeta:id="i_Na_m_gate.m_infinity" name="m_infinity" units="dimensionless"/>
<variable cmeta:id="i_Na_m_gate.tau_m" name="tau_m" units="second"/>
<variable cmeta:id="i_Na_m_gate.V" name="V" public_interface="in" units="millivolt"/>
<variable cmeta:id="i_Na_m_gate.time" name="time" public_interface="in" units="second"/>
<variable cmeta:id="i_Na_m_gate.delta_m" initial_value="1e-5" name="delta_m" units="millivolt"/>
<variable cmeta:id="i_Na_m_gate.E0_m" name="E0_m" units="millivolt"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>m_infinity</ci>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="millivolt">42.0504</cn>
</apply>
</apply>
<cn cellml:units="millivolt">8.3106</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>E0_m</ci>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="millivolt">41</cn>
</apply>
</apply>
<apply>
<eq/>
<ci>alpha_m</ci>
<piecewise>
<piece>
<cn cellml:units="per_second">2000</cn>
<apply>
<lt/>
<apply>
<abs/>
<ci>E0_m</ci>
</apply>
<ci>delta_m</ci>
</apply>
</piece>
<otherwise>
<apply>
<divide/>
<apply>
<times/>
<cn cellml:units="per_millivolt_second">200</cn>
<ci>E0_m</ci>
</apply>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="per_millivolt">0.1</cn>
</apply>
<ci>E0_m</ci>
</apply>
</apply>
</apply>
</apply>
</otherwise>
</piecewise>
</apply>
<apply>
<eq/>
<ci>beta_m</ci>
<apply>
<times/>
<cn cellml:units="per_second">8000</cn>
<apply>
<exp/>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="per_millivolt">0.056</cn>
</apply>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="millivolt">66</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>tau_m</ci>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<plus/>
<ci>alpha_m</ci>
<ci>beta_m</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>m</ci>
</apply>
<apply>
<divide/>
<apply>
<minus/>
<ci>m_infinity</ci>
<ci>m</ci>
</apply>
<ci>tau_m</ci>
</apply>
</apply>
</math>
</component>
<component name="i_Na_h_gate">
<variable cmeta:id="i_Na_h_gate.h" initial_value="0.003058" name="h" public_interface="out" units="dimensionless"/>
<variable cmeta:id="i_Na_h_gate.alpha_h" name="alpha_h" units="per_second"/>
<variable cmeta:id="i_Na_h_gate.beta_h" name="beta_h" units="per_second"/>
<variable cmeta:id="i_Na_h_gate.h_infinity" name="h_infinity" units="dimensionless"/>
<variable cmeta:id="i_Na_h_gate.tau_h" name="tau_h" units="second"/>
<variable cmeta:id="i_Na_h_gate.V" name="V" public_interface="in" units="millivolt"/>
<variable cmeta:id="i_Na_h_gate.time" name="time" public_interface="in" units="second"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>h_infinity</ci>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="millivolt">69.804</cn>
</apply>
<cn cellml:units="millivolt">4.4565</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>alpha_h</ci>
<apply>
<times/>
<cn cellml:units="per_second">20</cn>
<apply>
<exp/>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="per_millivolt">0.125</cn>
</apply>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="millivolt">75</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>beta_h</ci>
<apply>
<divide/>
<cn cellml:units="per_second">2000</cn>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="dimensionless">320</cn>
<apply>
<exp/>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="per_millivolt">0.1</cn>
</apply>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="millivolt">75</cn>
</apply>
</apply>
</apply>
</apply>
<cn cellml:units="dimensionless">1</cn>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>tau_h</ci>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<plus/>
<ci>alpha_h</ci>
<ci>beta_h</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>h</ci>
</apply>
<apply>
<divide/>
<apply>
<minus/>
<ci>h_infinity</ci>
<ci>h</ci>
</apply>
<ci>tau_h</ci>
</apply>
</apply>
</math>
</component>
<component name="i_CaL">
<variable cmeta:id="i_CaL.i_CaL" name="i_CaL" public_interface="out" units="nanoA"/>
<variable cmeta:id="i_CaL.i_siCa" name="i_siCa" public_interface="out" units="nanoA"/>
<variable cmeta:id="i_CaL.i_siK" name="i_siK" public_interface="out" units="nanoA"/>
<variable cmeta:id="i_CaL.i_siNa" name="i_siNa" public_interface="out" units="nanoA"/>
<variable cmeta:id="i_CaL.ACh_block" name="ACh_block" units="dimensionless"/>
<variable cmeta:id="i_CaL.ACh" name="ACh" public_interface="in" units="millimolar"/>
<variable cmeta:id="i_CaL.time" name="time" private_interface="out" public_interface="in" units="second"/>
<variable cmeta:id="i_CaL.V" name="V" private_interface="out" public_interface="in" units="millivolt"/>
<variable cmeta:id="i_CaL.P_CaL" initial_value="0.4578" name="P_CaL" units="nanoA_per_millimolar"/>
<variable cmeta:id="i_CaL.RTONF" name="RTONF" public_interface="in" units="millivolt"/>
<variable cmeta:id="i_CaL.Nao" name="Nao" public_interface="in" units="millimolar"/>
<variable cmeta:id="i_CaL.Nai" name="Nai" public_interface="in" units="millimolar"/>
<variable cmeta:id="i_CaL.Ko" name="Ko" public_interface="in" units="millimolar"/>
<variable cmeta:id="i_CaL.Ki" name="Ki" public_interface="in" units="millimolar"/>
<variable cmeta:id="i_CaL.Cao" name="Cao" public_interface="in" units="millimolar"/>
<variable cmeta:id="i_CaL.Ca_sub" name="Ca_sub" private_interface="out" public_interface="in" units="millimolar"/>
<variable cmeta:id="i_CaL.dL" name="dL" private_interface="in" units="dimensionless"/>
<variable cmeta:id="i_CaL.fL" name="fL" private_interface="in" units="dimensionless"/>
<variable cmeta:id="i_CaL.fCa" name="fCa" private_interface="in" units="dimensionless"/>
<variable cmeta:id="i_CaL.Iso_increase" name="Iso_increase" units="dimensionless"/>
<variable cmeta:id="i_CaL.Iso_1_uM" name="Iso_1_uM" private_interface="out" public_interface="in" units="dimensionless"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>Iso_increase</ci>
<piecewise>
<piece>
<cn cellml:units="dimensionless">1.23</cn>
<apply>
<gt/>
<ci>Iso_1_uM</ci>
<cn cellml:units="dimensionless">0</cn>
</apply>
</piece>
<otherwise>
<cn cellml:units="dimensionless">1</cn>
</otherwise>
</piecewise>
</apply>
<apply>
<eq/>
<ci>i_siCa</ci>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<cn cellml:units="dimensionless">2</cn>
<ci>P_CaL</ci>
<apply>
<minus/>
<ci>V</ci>
<cn cellml:units="millivolt">0</cn>
</apply>
</apply>
<apply>
<times/>
<ci>RTONF</ci>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
</apply>
<apply>
<minus/>
<ci>V</ci>
<cn cellml:units="millivolt">0</cn>
</apply>
<cn cellml:units="dimensionless">2</cn>
</apply>
<ci>RTONF</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<minus/>
<ci>Ca_sub</ci>
<apply>
<times/>
<ci>Cao</ci>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="dimensionless">2</cn>
</apply>
<apply>
<minus/>
<ci>V</ci>
<cn cellml:units="millivolt">0</cn>
</apply>
</apply>
<ci>RTONF</ci>
</apply>
</apply>
</apply>
</apply>
<ci>dL</ci>
<ci>fL</ci>
<ci>fCa</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>i_siK</ci>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<cn cellml:units="dimensionless">0.000365</cn>
<ci>P_CaL</ci>
<apply>
<minus/>
<ci>V</ci>
<cn cellml:units="millivolt">0</cn>
</apply>
</apply>
<apply>
<times/>
<ci>RTONF</ci>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
</apply>
<apply>
<minus/>
<ci>V</ci>
<cn cellml:units="millivolt">0</cn>
</apply>
</apply>
<ci>RTONF</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<minus/>
<ci>Ki</ci>
<apply>
<times/>
<ci>Ko</ci>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
</apply>
<apply>
<minus/>
<ci>V</ci>
<cn cellml:units="millivolt">0</cn>
</apply>
</apply>
<ci>RTONF</ci>
</apply>
</apply>
</apply>
</apply>
<ci>dL</ci>
<ci>fL</ci>
<ci>fCa</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>i_siNa</ci>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<cn cellml:units="dimensionless">0.0000185</cn>
<ci>P_CaL</ci>
<apply>
<minus/>
<ci>V</ci>
<cn cellml:units="millivolt">0</cn>
</apply>
</apply>
<apply>
<times/>
<ci>RTONF</ci>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
</apply>
<apply>
<minus/>
<ci>V</ci>
<cn cellml:units="millivolt">0</cn>
</apply>
</apply>
<ci>RTONF</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<minus/>
<ci>Nai</ci>
<apply>
<times/>
<ci>Nao</ci>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
</apply>
<apply>
<minus/>
<ci>V</ci>
<cn cellml:units="millivolt">0</cn>
</apply>
</apply>
<ci>RTONF</ci>
</apply>
</apply>
</apply>
</apply>
<ci>dL</ci>
<ci>fL</ci>
<ci>fCa</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>ACh_block</ci>
<apply>
<divide/>
<apply>
<times/>
<cn cellml:units="dimensionless">0.31</cn>
<ci>ACh</ci>
</apply>
<apply>
<plus/>
<ci>ACh</ci>
<cn cellml:units="millimolar">0.00009</cn>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>i_CaL</ci>
<apply>
<times/>
<apply>
<plus/>
<ci>i_siCa</ci>
<ci>i_siK</ci>
<ci>i_siNa</ci>
</apply>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<ci>ACh_block</ci>
</apply>
<cn cellml:units="dimensionless">1</cn>
<ci>Iso_increase</ci>
</apply>
</apply>
</math>
</component>
<component name="i_CaL_dL_gate">
<variable cmeta:id="i_CaL_dL_gate.dL" initial_value="0.001921" name="dL" public_interface="out" units="dimensionless"/>
<variable cmeta:id="i_CaL_dL_gate.dL_infinity" name="dL_infinity" units="dimensionless"/>
<variable cmeta:id="i_CaL_dL_gate.tau_dL" name="tau_dL" units="second"/>
<variable cmeta:id="i_CaL_dL_gate.alpha_dL" name="alpha_dL" units="per_second"/>
<variable cmeta:id="i_CaL_dL_gate.beta_dL" name="beta_dL" units="per_second"/>
<variable cmeta:id="i_CaL_dL_gate.adVm" name="adVm" units="millivolt"/>
<variable cmeta:id="i_CaL_dL_gate.bdVm" name="bdVm" units="millivolt"/>
<variable cmeta:id="i_CaL_dL_gate.V" name="V" public_interface="in" units="millivolt"/>
<variable cmeta:id="i_CaL_dL_gate.time" name="time" public_interface="in" units="second"/>
<variable cmeta:id="i_CaL_dL_gate.Iso_1_uM" name="Iso_1_uM" public_interface="in" units="dimensionless"/>
<variable cmeta:id="i_CaL_dL_gate.k_dL" initial_value="4.3371" name="k_dL" units="millivolt"/>
<variable cmeta:id="i_CaL_dL_gate.V_dL" initial_value="-16.4508" name="V_dL" units="millivolt"/>
<variable cmeta:id="i_CaL_dL_gate.Iso_shift_dL" name="Iso_shift_dL" units="millivolt"/>
<variable cmeta:id="i_CaL_dL_gate.Iso_slope_dL" name="Iso_slope_dL" units="dimensionless"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>Iso_shift_dL</ci>
<piecewise>
<piece>
<apply>
<minus/>
<cn cellml:units="millivolt">8</cn>
</apply>
<apply>
<gt/>
<ci>Iso_1_uM</ci>
<cn cellml:units="dimensionless">0</cn>
</apply>
</piece>
<otherwise>
<cn cellml:units="millivolt">0</cn>
</otherwise>
</piecewise>
</apply>
<apply>
<eq/>
<ci>Iso_slope_dL</ci>
<piecewise>
<piece>
<apply>
<minus/>
<cn cellml:units="dimensionless">27</cn>
</apply>
<apply>
<gt/>
<ci>Iso_1_uM</ci>
<cn cellml:units="dimensionless">0</cn>
</apply>
</piece>
<otherwise>
<cn cellml:units="dimensionless">0</cn>
</otherwise>
</piecewise>
</apply>
<apply>
<eq/>
<ci>dL_infinity</ci>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<minus/>
<apply>
<minus/>
<ci>V</ci>
<ci>V_dL</ci>
</apply>
<ci>Iso_shift_dL</ci>
</apply>
</apply>
<apply>
<times/>
<ci>k_dL</ci>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<divide/>
<ci>Iso_slope_dL</ci>
<cn cellml:units="dimensionless">100</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>tau_dL</ci>
<apply>
<divide/>
<cn cellml:units="dimensionless">0.001</cn>
<apply>
<plus/>
<ci>alpha_dL</ci>
<ci>beta_dL</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>alpha_dL</ci>
<apply>
<minus/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="per_millivolt_second">0.02839</cn>
</apply>
<apply>
<plus/>
<ci>adVm</ci>
<cn cellml:units="millivolt">41.8</cn>
</apply>
</apply>
<apply>
<minus/>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<plus/>
<ci>adVm</ci>
<cn cellml:units="millivolt">41.8</cn>
</apply>
</apply>
<cn cellml:units="millivolt">2.5</cn>
</apply>
</apply>
<cn cellml:units="dimensionless">1</cn>
</apply>
</apply>
<apply>
<divide/>
<apply>
<times/>
<cn cellml:units="per_millivolt_second">0.0849</cn>
<apply>
<plus/>
<ci>adVm</ci>
<cn cellml:units="millivolt">6.8</cn>
</apply>
</apply>
<apply>
<minus/>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<plus/>
<ci>adVm</ci>
<cn cellml:units="millivolt">6.8</cn>
</apply>
</apply>
<cn cellml:units="millivolt">4.8</cn>
</apply>
</apply>
<cn cellml:units="dimensionless">1</cn>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>adVm</ci>
<piecewise>
<piece>
<apply>
<minus/>
<cn cellml:units="millivolt">41.80001</cn>
</apply>
<apply>
<eq/>
<ci>V</ci>
<apply>
<minus/>
<cn cellml:units="millivolt">41.8</cn>
</apply>
</apply>
</piece>
<piece>
<cn cellml:units="millivolt">0</cn>
<apply>
<eq/>
<ci>V</ci>
<cn cellml:units="millivolt">0</cn>
</apply>
</piece>
<piece>
<apply>
<minus/>
<cn cellml:units="millivolt">6.80001</cn>
</apply>
<apply>
<eq/>
<ci>V</ci>
<apply>
<minus/>
<cn cellml:units="millivolt">6.8</cn>
</apply>
</apply>
</piece>
<otherwise>
<ci>V</ci>
</otherwise>
</piecewise>
</apply>
<apply>
<eq/>
<ci>beta_dL</ci>
<apply>
<divide/>
<apply>
<times/>
<cn cellml:units="per_millivolt_second">0.01143</cn>
<apply>
<plus/>
<ci>bdVm</ci>
<cn cellml:units="millivolt">1.8</cn>
</apply>
</apply>
<apply>
<minus/>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<plus/>
<ci>bdVm</ci>
<cn cellml:units="millivolt">1.8</cn>
</apply>
<cn cellml:units="millivolt">2.5</cn>
</apply>
</apply>
<cn cellml:units="dimensionless">1</cn>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>bdVm</ci>
<piecewise>
<piece>
<apply>
<minus/>
<cn cellml:units="millivolt">1.80001</cn>
</apply>
<apply>
<eq/>
<ci>V</ci>
<apply>
<minus/>
<cn cellml:units="millivolt">1.8</cn>
</apply>
</apply>
</piece>
<otherwise>
<ci>V</ci>
</otherwise>
</piecewise>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>dL</ci>
</apply>
<apply>
<divide/>
<apply>
<minus/>
<ci>dL_infinity</ci>
<ci>dL</ci>
</apply>
<ci>tau_dL</ci>
</apply>
</apply>
</math>
</component>
<component name="i_CaL_fL_gate">
<variable cmeta:id="i_CaL_fL_gate.fL" initial_value="0.846702" name="fL" public_interface="out" units="dimensionless"/>
<variable cmeta:id="i_CaL_fL_gate.fL_infinity" name="fL_infinity" units="dimensionless"/>
<variable cmeta:id="i_CaL_fL_gate.tau_fL" name="tau_fL" units="second"/>
<variable cmeta:id="i_CaL_fL_gate.V" name="V" public_interface="in" units="millivolt"/>
<variable cmeta:id="i_CaL_fL_gate.time" name="time" public_interface="in" units="second"/>
<variable cmeta:id="i_CaL_fL_gate.shift_fL" initial_value="0" name="shift_fL" units="millivolt"/>
<variable cmeta:id="i_CaL_fL_gate.k_fL" initial_value="0" name="k_fL" units="millivolt"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>fL_infinity</ci>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="millivolt">37.4</cn>
<ci>shift_fL</ci>
</apply>
<apply>
<plus/>
<cn cellml:units="millivolt">5.3</cn>
<ci>k_fL</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>tau_fL</ci>
<apply>
<times/>
<cn cellml:units="dimensionless">0.001</cn>
<apply>
<plus/>
<cn cellml:units="second">44.3</cn>
<apply>
<times/>
<cn cellml:units="second">230</cn>
<apply>
<exp/>
<apply>
<minus/>
<apply>
<power/>
<apply>
<divide/>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="millivolt">36</cn>
</apply>
<cn cellml:units="millivolt">10</cn>
</apply>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>fL</ci>
</apply>
<apply>
<divide/>
<apply>
<minus/>
<ci>fL_infinity</ci>
<ci>fL</ci>
</apply>
<ci>tau_fL</ci>
</apply>
</apply>
</math>
</component>
<component name="i_CaL_fCa_gate">
<variable cmeta:id="i_CaL_fCa_gate.fCa" initial_value="0.844449" name="fCa" public_interface="out" units="dimensionless"/>
<variable cmeta:id="i_CaL_fCa_gate.alpha_fCa" initial_value="0.0075" name="alpha_fCa" units="per_second"/>
<variable cmeta:id="i_CaL_fCa_gate.fCa_infinity" name="fCa_infinity" units="dimensionless"/>
<variable cmeta:id="i_CaL_fCa_gate.tau_fCa" name="tau_fCa" units="second"/>
<variable cmeta:id="i_CaL_fCa_gate.Km_fCa" initial_value="0.000338" name="Km_fCa" units="millimolar"/>
<variable cmeta:id="i_CaL_fCa_gate.V" name="V" public_interface="in" units="millivolt"/>
<variable cmeta:id="i_CaL_fCa_gate.Ca_sub" name="Ca_sub" public_interface="in" units="millimolar"/>
<variable cmeta:id="i_CaL_fCa_gate.time" name="time" public_interface="in" units="second"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>fCa_infinity</ci>
<apply>
<divide/>
<ci>Km_fCa</ci>
<apply>
<plus/>
<ci>Km_fCa</ci>
<ci>Ca_sub</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>tau_fCa</ci>
<apply>
<divide/>
<apply>
<times/>
<cn cellml:units="dimensionless">0.001</cn>
<ci>fCa_infinity</ci>
</apply>
<ci>alpha_fCa</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>fCa</ci>
</apply>
<apply>
<divide/>
<apply>
<minus/>
<ci>fCa_infinity</ci>
<ci>fCa</ci>
</apply>
<ci>tau_fCa</ci>
</apply>
</apply>
</math>
</component>
<component name="i_CaT">
<variable cmeta:id="i_CaT.i_CaT" name="i_CaT" public_interface="out" units="nanoA"/>
<variable cmeta:id="i_CaT.time" name="time" private_interface="out" public_interface="in" units="second"/>
<variable cmeta:id="i_CaT.V" name="V" private_interface="out" public_interface="in" units="millivolt"/>
<variable cmeta:id="i_CaT.P_CaT" initial_value="0.04132" name="P_CaT" units="nanoA_per_millimolar"/>
<variable cmeta:id="i_CaT.RTONF" name="RTONF" public_interface="in" units="millivolt"/>
<variable cmeta:id="i_CaT.Cao" name="Cao" public_interface="in" units="millimolar"/>
<variable cmeta:id="i_CaT.Ca_sub" name="Ca_sub" private_interface="out" public_interface="in" units="millimolar"/>
<variable cmeta:id="i_CaT.dT" name="dT" private_interface="in" units="dimensionless"/>
<variable cmeta:id="i_CaT.fT" name="fT" private_interface="in" units="dimensionless"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>i_CaT</ci>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<cn cellml:units="dimensionless">2</cn>
<ci>P_CaT</ci>
<ci>V</ci>
</apply>
<apply>
<times/>
<ci>RTONF</ci>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
</apply>
<ci>V</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
<ci>RTONF</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<minus/>
<ci>Ca_sub</ci>
<apply>
<times/>
<ci>Cao</ci>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="dimensionless">2</cn>
</apply>
<ci>V</ci>
</apply>
<ci>RTONF</ci>
</apply>
</apply>
</apply>
</apply>
<ci>dT</ci>
<ci>fT</ci>
</apply>
</apply>
</math>
</component>
<component name="i_CaT_dT_gate">
<variable cmeta:id="i_CaT_dT_gate.dT" initial_value="0.268909" name="dT" public_interface="out" units="dimensionless"/>
<variable cmeta:id="i_CaT_dT_gate.dT_infinity" name="dT_infinity" units="dimensionless"/>
<variable cmeta:id="i_CaT_dT_gate.tau_dT" name="tau_dT" units="second"/>
<variable cmeta:id="i_CaT_dT_gate.V" name="V" public_interface="in" units="millivolt"/>
<variable cmeta:id="i_CaT_dT_gate.time" name="time" public_interface="in" units="second"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>dT_infinity</ci>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="millivolt">38.3</cn>
</apply>
</apply>
<cn cellml:units="millivolt">5.5</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>tau_dT</ci>
<apply>
<divide/>
<cn cellml:units="second">0.001</cn>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="dimensionless">1.068</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="millivolt">38.3</cn>
</apply>
<cn cellml:units="millivolt">30</cn>
</apply>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">1.068</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="millivolt">38.3</cn>
</apply>
</apply>
<cn cellml:units="millivolt">30</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>dT</ci>
</apply>
<apply>
<divide/>
<apply>
<minus/>
<ci>dT_infinity</ci>
<ci>dT</ci>
</apply>
<ci>tau_dT</ci>
</apply>
</apply>
</math>
</component>
<component name="i_CaT_fT_gate">
<variable cmeta:id="i_CaT_fT_gate.fT" initial_value="0.020484" name="fT" public_interface="out" units="dimensionless"/>
<variable cmeta:id="i_CaT_fT_gate.fT_infinity" name="fT_infinity" units="dimensionless"/>
<variable cmeta:id="i_CaT_fT_gate.tau_fT" name="tau_fT" units="second"/>
<variable cmeta:id="i_CaT_fT_gate.V" name="V" public_interface="in" units="millivolt"/>
<variable cmeta:id="i_CaT_fT_gate.time" name="time" public_interface="in" units="second"/>
<variable cmeta:id="i_CaT_fT_gate.offset_fT" initial_value="0" name="offset_fT" units="second"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>fT_infinity</ci>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="millivolt">58.7</cn>
</apply>
<cn cellml:units="millivolt">3.8</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>tau_fT</ci>
<apply>
<plus/>
<apply>
<divide/>
<cn cellml:units="second">1</cn>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="dimensionless">16.67</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="millivolt">75</cn>
</apply>
</apply>
<cn cellml:units="millivolt">83.3</cn>
</apply>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">16.67</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="millivolt">75</cn>
</apply>
<cn cellml:units="millivolt">15.38</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
<ci>offset_fT</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>fT</ci>
</apply>
<apply>
<divide/>
<apply>
<minus/>
<ci>fT_infinity</ci>
<ci>fT</ci>
</apply>
<ci>tau_fT</ci>
</apply>
</apply>
</math>
</component>
<component name="Ca_SR_release">
<variable cmeta:id="Ca_SR_release.j_SRCarel" name="j_SRCarel" public_interface="out" units="millimolar_per_second"/>
<variable cmeta:id="Ca_SR_release.R" initial_value="0.9308" name="R" units="dimensionless"/>
<variable cmeta:id="Ca_SR_release.O" initial_value="6.181512e-9" name="O" units="dimensionless"/>
<variable cmeta:id="Ca_SR_release.I" initial_value="4.595622e-10" name="I" units="dimensionless"/>
<variable cmeta:id="Ca_SR_release.RI" initial_value="0.069199" name="RI" units="dimensionless"/>
<variable cmeta:id="Ca_SR_release.ks" initial_value="148041085.1" name="ks" units="per_second"/>
<variable cmeta:id="Ca_SR_release.MaxSR" initial_value="15" name="MaxSR" units="dimensionless"/>
<variable cmeta:id="Ca_SR_release.MinSR" initial_value="1" name="MinSR" units="dimensionless"/>
<variable cmeta:id="Ca_SR_release.EC50_SR" initial_value="0.45" name="EC50_SR" units="millimolar"/>
<variable cmeta:id="Ca_SR_release.HSR" initial_value="2.5" name="HSR" units="dimensionless"/>
<variable cmeta:id="Ca_SR_release.koSRCa" name="koSRCa" units="per_millimolar2_second"/>
<variable cmeta:id="Ca_SR_release.kiSRCa" name="kiSRCa" units="per_millimolar_second"/>
<variable cmeta:id="Ca_SR_release.koCa" initial_value="10000" name="koCa" units="per_millimolar2_second"/>
<variable cmeta:id="Ca_SR_release.kiCa" initial_value="500" name="kiCa" units="per_millimolar_second"/>
<variable cmeta:id="Ca_SR_release.kCaSR" name="kCaSR" units="dimensionless"/>
<variable cmeta:id="Ca_SR_release.kim" initial_value="5" name="kim" units="per_second"/>
<variable cmeta:id="Ca_SR_release.kom" initial_value="660" name="kom" units="per_second"/>
<variable cmeta:id="Ca_SR_release.time" name="time" public_interface="in" units="second"/>
<variable cmeta:id="Ca_SR_release.Ca_sub" name="Ca_sub" public_interface="in" units="millimolar"/>
<variable cmeta:id="Ca_SR_release.Ca_jsr" name="Ca_jsr" public_interface="in" units="millimolar"/>
<variable cmeta:id="Ca_SR_release.diff" name="diff" units="millimolar"/>
<variable cmeta:id="Ca_SR_release.P_tot" name="P_tot" units="dimensionless"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>j_SRCarel</ci>
<apply>
<times/>
<ci>ks</ci>
<ci>O</ci>
<apply>
<minus/>
<ci>Ca_jsr</ci>
<ci>Ca_sub</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>diff</ci>
<apply>
<minus/>
<ci>Ca_jsr</ci>
<ci>Ca_sub</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>kCaSR</ci>
<apply>
<minus/>
<ci>MaxSR</ci>
<apply>
<divide/>
<apply>
<minus/>
<ci>MaxSR</ci>
<ci>MinSR</ci>
</apply>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<power/>
<apply>
<divide/>
<ci>EC50_SR</ci>
<ci>Ca_jsr</ci>
</apply>
<ci>HSR</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>koSRCa</ci>
<apply>
<divide/>
<ci>koCa</ci>
<ci>kCaSR</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>kiSRCa</ci>
<apply>
<times/>
<ci>kiCa</ci>
<ci>kCaSR</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>R</ci>
</apply>
<apply>
<minus/>
<apply>
<minus/>
<apply>
<times/>
<ci>kim</ci>
<ci>RI</ci>
</apply>
<apply>
<times/>
<ci>kiSRCa</ci>
<ci>Ca_sub</ci>
<ci>R</ci>
</apply>
</apply>
<apply>
<minus/>
<apply>
<times/>
<ci>koSRCa</ci>
<apply>
<power/>
<ci>Ca_sub</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
<ci>R</ci>
</apply>
<apply>
<times/>
<ci>kom</ci>
<ci>O</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>O</ci>
</apply>
<apply>
<minus/>
<apply>
<minus/>
<apply>
<times/>
<ci>koSRCa</ci>
<apply>
<power/>
<ci>Ca_sub</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
<ci>R</ci>
</apply>
<apply>
<times/>
<ci>kom</ci>
<ci>O</ci>
</apply>
</apply>
<apply>
<minus/>
<apply>
<times/>
<ci>kiSRCa</ci>
<ci>Ca_sub</ci>
<ci>O</ci>
</apply>
<apply>
<times/>
<ci>kim</ci>
<ci>I</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>I</ci>
</apply>
<apply>
<minus/>
<apply>
<minus/>
<apply>
<times/>
<ci>kiSRCa</ci>
<ci>Ca_sub</ci>
<ci>O</ci>
</apply>
<apply>
<times/>
<ci>kim</ci>
<ci>I</ci>
</apply>
</apply>
<apply>
<minus/>
<apply>
<times/>
<ci>kom</ci>
<ci>I</ci>
</apply>
<apply>
<times/>
<ci>koSRCa</ci>
<apply>
<power/>
<ci>Ca_sub</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
<ci>RI</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>RI</ci>
</apply>
<apply>
<minus/>
<apply>
<minus/>
<apply>
<times/>
<ci>kom</ci>
<ci>I</ci>
</apply>
<apply>
<times/>
<ci>koSRCa</ci>
<apply>
<power/>
<ci>Ca_sub</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
<ci>RI</ci>
</apply>
</apply>
<apply>
<minus/>
<apply>
<times/>
<ci>kim</ci>
<ci>RI</ci>
</apply>
<apply>
<times/>
<ci>kiSRCa</ci>
<ci>Ca_sub</ci>
<ci>R</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>P_tot</ci>
<apply>
<plus/>
<ci>R</ci>
<ci>O</ci>
<ci>I</ci>
<ci>RI</ci>
</apply>
</apply>
</math>
</component>
<component name="Ca_intracellular_fluxes">
<variable cmeta:id="Ca_intracellular_fluxes.j_Ca_dif" name="j_Ca_dif" public_interface="out" units="millimolar_per_second"/>
<variable cmeta:id="Ca_intracellular_fluxes.j_up" name="j_up" public_interface="out" units="millimolar_per_second"/>
<variable cmeta:id="Ca_intracellular_fluxes.j_tr" name="j_tr" public_interface="out" units="millimolar_per_second"/>
<variable cmeta:id="Ca_intracellular_fluxes.tau_dif_Ca" initial_value="5.469e-5" name="tau_dif_Ca" units="second"/>
<variable cmeta:id="Ca_intracellular_fluxes.tau_tr" initial_value="0.04" name="tau_tr" units="second"/>
<variable cmeta:id="Ca_intracellular_fluxes.P_up" name="P_up" units="millimolar_per_second"/>
<variable cmeta:id="Ca_intracellular_fluxes.P_up_basal" initial_value="5" name="P_up_basal" units="millimolar_per_second"/>
<variable cmeta:id="Ca_intracellular_fluxes.b_up" name="b_up" units="dimensionless"/>
<variable cmeta:id="Ca_intracellular_fluxes.ACh" name="ACh" public_interface="in" units="millimolar"/>
<variable cmeta:id="Ca_intracellular_fluxes.Iso_1_uM" name="Iso_1_uM" public_interface="in" units="dimensionless"/>
<variable cmeta:id="Ca_intracellular_fluxes.K_up" initial_value="0.000286113" name="K_up" units="millimolar"/>
<variable cmeta:id="Ca_intracellular_fluxes.Ca_nsr" name="Ca_nsr" public_interface="in" units="millimolar"/>
<variable cmeta:id="Ca_intracellular_fluxes.Cai" name="Cai" public_interface="in" units="millimolar"/>
<variable cmeta:id="Ca_intracellular_fluxes.Ca_sub" name="Ca_sub" public_interface="in" units="millimolar"/>
<variable cmeta:id="Ca_intracellular_fluxes.Ca_jsr" name="Ca_jsr" public_interface="in" units="millimolar"/>
<variable cmeta:id="Ca_intracellular_fluxes.slope_up" initial_value="5e-5" name="slope_up" units="millimolar"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>b_up</ci>
<piecewise>
<piece>
<apply>
<minus/>
<cn cellml:units="dimensionless">0.25</cn>
</apply>
<apply>
<gt/>
<ci>Iso_1_uM</ci>
<cn cellml:units="dimensionless">0</cn>
</apply>
</piece>
<piece>
<apply>
<divide/>
<apply>
<times/>
<cn cellml:units="dimensionless">0.7</cn>
<ci>ACh</ci>
</apply>
<apply>
<plus/>
<cn cellml:units="millimolar">0.00009</cn>
<ci>ACh</ci>
</apply>
</apply>
<apply>
<gt/>
<ci>ACh</ci>
<cn cellml:units="millimolar">0</cn>
</apply>
</piece>
<otherwise>
<cn cellml:units="dimensionless">0</cn>
</otherwise>
</piecewise>
</apply>
<apply>
<eq/>
<ci>P_up</ci>
<apply>
<times/>
<ci>P_up_basal</ci>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<ci>b_up</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>j_Ca_dif</ci>
<apply>
<divide/>
<apply>
<minus/>
<ci>Ca_sub</ci>
<ci>Cai</ci>
</apply>
<ci>tau_dif_Ca</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>j_up</ci>
<apply>
<divide/>
<ci>P_up</ci>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<plus/>
<apply>
<minus/>
<ci>Cai</ci>
</apply>
<ci>K_up</ci>
</apply>
<ci>slope_up</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>j_tr</ci>
<apply>
<divide/>
<apply>
<minus/>
<ci>Ca_nsr</ci>
<ci>Ca_jsr</ci>
</apply>
<ci>tau_tr</ci>
</apply>
</apply>
</math>
</component>
<component name="Ca_buffering">
<variable cmeta:id="Ca_buffering.TC_tot" initial_value="0.031" name="TC_tot" public_interface="out" units="millimolar"/>
<variable cmeta:id="Ca_buffering.TMC_tot" initial_value="0.062" name="TMC_tot" public_interface="out" units="millimolar"/>
<variable cmeta:id="Ca_buffering.CM_tot" initial_value="0.045" name="CM_tot" public_interface="out" units="millimolar"/>
<variable cmeta:id="Ca_buffering.CQ_tot" initial_value="10" name="CQ_tot" public_interface="out" units="millimolar"/>
<variable cmeta:id="Ca_buffering.delta_fTC" name="delta_fTC" public_interface="out" units="per_second"/>
<variable cmeta:id="Ca_buffering.delta_fTMC" name="delta_fTMC" public_interface="out" units="per_second"/>
<variable cmeta:id="Ca_buffering.delta_fCMs" name="delta_fCMs" public_interface="out" units="per_second"/>
<variable cmeta:id="Ca_buffering.delta_fCMi" name="delta_fCMi" public_interface="out" units="per_second"/>
<variable cmeta:id="Ca_buffering.delta_fCQ" name="delta_fCQ" public_interface="out" units="per_second"/>
<variable cmeta:id="Ca_buffering.delta_fTMM" name="delta_fTMM" units="per_second"/>
<variable cmeta:id="Ca_buffering.fTMM" initial_value="0.653777" name="fTMM" units="dimensionless"/>
<variable cmeta:id="Ca_buffering.fCMi" initial_value="0.217311" name="fCMi" units="dimensionless"/>
<variable cmeta:id="Ca_buffering.fCMs" initial_value="0.158521" name="fCMs" units="dimensionless"/>
<variable cmeta:id="Ca_buffering.fTC" initial_value="0.017929" name="fTC" units="dimensionless"/>
<variable cmeta:id="Ca_buffering.fTMC" initial_value="0.259947" name="fTMC" units="dimensionless"/>
<variable cmeta:id="Ca_buffering.fCQ" initial_value="0.138975" name="fCQ" units="dimensionless"/>
<variable cmeta:id="Ca_buffering.kf_TC" initial_value="88800" name="kf_TC" units="per_millimolar_second"/>
<variable cmeta:id="Ca_buffering.kf_TMM" initial_value="2277" name="kf_TMM" units="per_millimolar_second"/>
<variable cmeta:id="Ca_buffering.kf_TMC" initial_value="227700" name="kf_TMC" units="per_millimolar_second"/>
<variable cmeta:id="Ca_buffering.kf_CM" initial_value="1.642e6" name="kf_CM" units="per_millimolar_second"/>
<variable cmeta:id="Ca_buffering.kf_CQ" initial_value="175.4" name="kf_CQ" units="per_millimolar_second"/>
<variable cmeta:id="Ca_buffering.kb_TC" initial_value="446" name="kb_TC" units="per_second"/>
<variable cmeta:id="Ca_buffering.kb_TMC" initial_value="7.51" name="kb_TMC" units="per_second"/>
<variable cmeta:id="Ca_buffering.kb_TMM" initial_value="751" name="kb_TMM" units="per_second"/>
<variable cmeta:id="Ca_buffering.kb_CM" initial_value="542" name="kb_CM" units="per_second"/>
<variable cmeta:id="Ca_buffering.kb_CQ" initial_value="445" name="kb_CQ" units="per_second"/>
<variable cmeta:id="Ca_buffering.time" name="time" public_interface="in" units="second"/>
<variable cmeta:id="Ca_buffering.Mgi" initial_value="2.5" name="Mgi" units="millimolar"/>
<variable cmeta:id="Ca_buffering.Cai" name="Cai" public_interface="in" units="millimolar"/>
<variable cmeta:id="Ca_buffering.Ca_sub" name="Ca_sub" public_interface="in" units="millimolar"/>
<variable cmeta:id="Ca_buffering.Ca_jsr" name="Ca_jsr" public_interface="in" units="millimolar"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>fTC</ci>
</apply>
<ci>delta_fTC</ci>
</apply>
<apply>
<eq/>
<ci>delta_fTC</ci>
<apply>
<minus/>
<apply>
<times/>
<ci>kf_TC</ci>
<ci>Cai</ci>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<ci>fTC</ci>
</apply>
</apply>
<apply>
<times/>
<ci>kb_TC</ci>
<ci>fTC</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>fTMC</ci>
</apply>
<ci>delta_fTMC</ci>
</apply>
<apply>
<eq/>
<ci>delta_fTMC</ci>
<apply>
<minus/>
<apply>
<times/>
<ci>kf_TMC</ci>
<ci>Cai</ci>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<plus/>
<ci>fTMC</ci>
<ci>fTMM</ci>
</apply>
</apply>
</apply>
<apply>
<times/>
<ci>kb_TMC</ci>
<ci>fTMC</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>fTMM</ci>
</apply>
<ci>delta_fTMM</ci>
</apply>
<apply>
<eq/>
<ci>delta_fTMM</ci>
<apply>
<minus/>
<apply>
<times/>
<ci>kf_TMM</ci>
<ci>Mgi</ci>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<plus/>
<ci>fTMC</ci>
<ci>fTMM</ci>
</apply>
</apply>
</apply>
<apply>
<times/>
<ci>kb_TMM</ci>
<ci>fTMM</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>fCMi</ci>
</apply>
<ci>delta_fCMi</ci>
</apply>
<apply>
<eq/>
<ci>delta_fCMi</ci>
<apply>
<minus/>
<apply>
<times/>
<ci>kf_CM</ci>
<ci>Cai</ci>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<ci>fCMi</ci>
</apply>
</apply>
<apply>
<times/>
<ci>kb_CM</ci>
<ci>fCMi</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>fCMs</ci>
</apply>
<ci>delta_fCMs</ci>
</apply>
<apply>
<eq/>
<ci>delta_fCMs</ci>
<apply>
<minus/>
<apply>
<times/>
<ci>kf_CM</ci>
<ci>Ca_sub</ci>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<ci>fCMs</ci>
</apply>
</apply>
<apply>
<times/>
<ci>kb_CM</ci>
<ci>fCMs</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>fCQ</ci>
</apply>
<ci>delta_fCQ</ci>
</apply>
<apply>
<eq/>
<ci>delta_fCQ</ci>
<apply>
<minus/>
<apply>
<times/>
<ci>kf_CQ</ci>
<ci>Ca_jsr</ci>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<ci>fCQ</ci>
</apply>
</apply>
<apply>
<times/>
<ci>kb_CQ</ci>
<ci>fCQ</ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="Ca_dynamics">
<variable cmeta:id="Ca_dynamics.Cai" initial_value="9.15641e-6" name="Cai" public_interface="out" units="millimolar"/>
<variable cmeta:id="Ca_dynamics.Ca_nsr" initial_value="0.435148" name="Ca_nsr" public_interface="out" units="millimolar"/>
<variable cmeta:id="Ca_dynamics.Ca_jsr" initial_value="0.409551" name="Ca_jsr" public_interface="out" units="millimolar"/>
<variable cmeta:id="Ca_dynamics.Ca_sub" initial_value="6.226104e-5" name="Ca_sub" public_interface="out" units="millimolar"/>
<variable cmeta:id="Ca_dynamics.V_i" name="V_i" public_interface="in" units="millimetre3"/>
<variable cmeta:id="Ca_dynamics.V_jsr" name="V_jsr" public_interface="in" units="millimetre3"/>
<variable cmeta:id="Ca_dynamics.V_nsr" name="V_nsr" public_interface="in" units="millimetre3"/>
<variable cmeta:id="Ca_dynamics.V_sub" name="V_sub" public_interface="in" units="millimetre3"/>
<variable cmeta:id="Ca_dynamics.time" name="time" public_interface="in" units="second"/>
<variable cmeta:id="Ca_dynamics.F" name="F" public_interface="in" units="coulomb_per_mole"/>
<variable cmeta:id="Ca_dynamics.i_siCa" name="i_siCa" public_interface="in" units="nanoA"/>
<variable cmeta:id="Ca_dynamics.i_CaT" name="i_CaT" public_interface="in" units="nanoA"/>
<variable cmeta:id="Ca_dynamics.i_NaCa" name="i_NaCa" public_interface="in" units="nanoA"/>
<variable cmeta:id="Ca_dynamics.j_SRCarel" name="j_SRCarel" public_interface="in" units="millimolar_per_second"/>
<variable cmeta:id="Ca_dynamics.j_Ca_dif" name="j_Ca_dif" public_interface="in" units="millimolar_per_second"/>
<variable cmeta:id="Ca_dynamics.j_tr" name="j_tr" public_interface="in" units="millimolar_per_second"/>
<variable cmeta:id="Ca_dynamics.j_up" name="j_up" public_interface="in" units="millimolar_per_second"/>
<variable cmeta:id="Ca_dynamics.TMC_tot" name="TMC_tot" public_interface="in" units="millimolar"/>
<variable cmeta:id="Ca_dynamics.CM_tot" name="CM_tot" public_interface="in" units="millimolar"/>
<variable cmeta:id="Ca_dynamics.TC_tot" name="TC_tot" public_interface="in" units="millimolar"/>
<variable cmeta:id="Ca_dynamics.CQ_tot" name="CQ_tot" public_interface="in" units="millimolar"/>
<variable cmeta:id="Ca_dynamics.delta_fTMC" name="delta_fTMC" public_interface="in" units="per_second"/>
<variable cmeta:id="Ca_dynamics.delta_fCMi" name="delta_fCMi" public_interface="in" units="per_second"/>
<variable cmeta:id="Ca_dynamics.delta_fCMs" name="delta_fCMs" public_interface="in" units="per_second"/>
<variable cmeta:id="Ca_dynamics.delta_fTC" name="delta_fTC" public_interface="in" units="per_second"/>
<variable cmeta:id="Ca_dynamics.delta_fCQ" name="delta_fCQ" public_interface="in" units="per_second"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>Cai</ci>
</apply>
<apply>
<minus/>
<apply>
<divide/>
<apply>
<times/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<minus/>
<apply>
<times/>
<ci>j_Ca_dif</ci>
<ci>V_sub</ci>
</apply>
<apply>
<times/>
<ci>j_up</ci>
<ci>V_nsr</ci>
</apply>
</apply>
</apply>
<ci>V_i</ci>
</apply>
<apply>
<plus/>
<apply>
<times/>
<ci>CM_tot</ci>
<ci>delta_fCMi</ci>
</apply>
<apply>
<times/>
<ci>TC_tot</ci>
<ci>delta_fTC</ci>
</apply>
<apply>
<times/>
<ci>TMC_tot</ci>
<ci>delta_fTMC</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>Ca_sub</ci>
</apply>
<apply>
<minus/>
<apply>
<divide/>
<apply>
<times/>
<ci>j_SRCarel</ci>
<ci>V_jsr</ci>
</apply>
<ci>V_sub</ci>
</apply>
<apply>
<plus/>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<plus/>
<ci>i_siCa</ci>
<ci>i_CaT</ci>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">2</cn>
<ci>i_NaCa</ci>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">2</cn>
<ci>F</ci>
<ci>V_sub</ci>
</apply>
</apply>
<ci>j_Ca_dif</ci>
<apply>
<times/>
<ci>CM_tot</ci>
<ci>delta_fCMs</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>Ca_nsr</ci>
</apply>
<apply>
<minus/>
<ci>j_up</ci>
<apply>
<divide/>
<apply>
<times/>
<ci>j_tr</ci>
<ci>V_jsr</ci>
</apply>
<ci>V_nsr</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>Ca_jsr</ci>
</apply>
<apply>
<minus/>
<ci>j_tr</ci>
<apply>
<plus/>
<ci>j_SRCarel</ci>
<apply>
<times/>
<ci>CQ_tot</ci>
<ci>delta_fCQ</ci>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<component name="Cell_parameters">
<variable cmeta:id="Cell_parameters.V_sub" name="V_sub" public_interface="out" units="millimetre3"/>
<variable cmeta:id="Cell_parameters.V_jsr" name="V_jsr" public_interface="out" units="millimetre3"/>
<variable cmeta:id="Cell_parameters.V_nsr" name="V_nsr" public_interface="out" units="millimetre3"/>
<variable cmeta:id="Cell_parameters.V_i" name="V_i" public_interface="out" units="millimetre3"/>
<variable cmeta:id="Cell_parameters.V_cell" name="V_cell" units="millimetre3"/>
<variable cmeta:id="Cell_parameters.V_jsr_part" initial_value="0.0012" name="V_jsr_part" units="dimensionless"/>
<variable cmeta:id="Cell_parameters.V_i_part" initial_value="0.46" name="V_i_part" units="dimensionless"/>
<variable cmeta:id="Cell_parameters.V_nsr_part" initial_value="0.0116" name="V_nsr_part" units="dimensionless"/>
<variable cmeta:id="Cell_parameters.R_cell" initial_value="3.9" name="R_cell" units="micrometre"/>
<variable cmeta:id="Cell_parameters.L_cell" initial_value="67" name="L_cell" units="micrometre"/>
<variable cmeta:id="Cell_parameters.L_sub" initial_value="0.02" name="L_sub" units="micrometre"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>V_cell</ci>
<apply>
<times/>
<cn cellml:units="micrometre_per_femtometre">0.000000001</cn>
<pi/>
<apply>
<power/>
<ci>R_cell</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
<ci>L_cell</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>V_sub</ci>
<apply>
<times/>
<cn cellml:units="micrometre_per_femtometre">0.000000001</cn>
<cn cellml:units="dimensionless">2</cn>
<pi/>
<ci>L_sub</ci>
<apply>
<minus/>
<ci>R_cell</ci>
<apply>
<divide/>
<ci>L_sub</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
<ci>L_cell</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>V_jsr</ci>
<apply>
<times/>
<ci>V_jsr_part</ci>
<ci>V_cell</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>V_i</ci>
<apply>
<minus/>
<apply>
<times/>
<ci>V_i_part</ci>
<ci>V_cell</ci>
</apply>
<ci>V_sub</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>V_nsr</ci>
<apply>
<times/>
<ci>V_nsr_part</ci>
<ci>V_cell</ci>
</apply>
</apply>
</math>
</component>
<component name="i_Kur">
<variable cmeta:id="i_Kur.i_Kur" name="i_Kur" public_interface="out" units="nanoA"/>
<variable cmeta:id="i_Kur.g_Kur" initial_value="0.1539e-3" name="g_Kur" units="microS"/>
<variable cmeta:id="i_Kur.time" name="time" private_interface="out" public_interface="in" units="second"/>
<variable cmeta:id="i_Kur.V" name="V" private_interface="out" public_interface="in" units="millivolt"/>
<variable cmeta:id="i_Kur.E_K" name="E_K" public_interface="in" units="millivolt"/>
<variable cmeta:id="i_Kur.r_Kur" name="r_Kur" private_interface="in" units="dimensionless"/>
<variable cmeta:id="i_Kur.s_Kur" name="s_Kur" private_interface="in" units="dimensionless"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>i_Kur</ci>
<apply>
<times/>
<ci>g_Kur</ci>
<ci>r_Kur</ci>
<ci>s_Kur</ci>
<apply>
<minus/>
<ci>V</ci>
<ci>E_K</ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="i_Kur_rKur_gate">
<variable cmeta:id="i_Kur_rKur_gate.r_Kur" initial_value="0.011845" name="r_Kur" public_interface="out" units="dimensionless"/>
<variable cmeta:id="i_Kur_rKur_gate.tau_r_Kur" name="tau_r_Kur" units="second"/>
<variable cmeta:id="i_Kur_rKur_gate.r_Kur_infinity" name="r_Kur_infinity" units="dimensionless"/>
<variable cmeta:id="i_Kur_rKur_gate.V" name="V" public_interface="in" units="millivolt"/>
<variable cmeta:id="i_Kur_rKur_gate.time" name="time" public_interface="in" units="second"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>r_Kur</ci>
</apply>
<apply>
<divide/>
<apply>
<minus/>
<ci>r_Kur_infinity</ci>
<ci>r_Kur</ci>
</apply>
<ci>tau_r_Kur</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>r_Kur_infinity</ci>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="millivolt">6</cn>
</apply>
<apply>
<minus/>
<cn cellml:units="millivolt">8.6</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>tau_r_Kur</ci>
<apply>
<plus/>
<apply>
<divide/>
<cn cellml:units="second">0.009</cn>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="millivolt">5</cn>
</apply>
<cn cellml:units="millivolt">12</cn>
</apply>
</apply>
</apply>
</apply>
<cn cellml:units="second">0.0005</cn>
</apply>
</apply>
</math>
</component>
<component name="i_Kur_sKur_gate">
<variable cmeta:id="i_Kur_sKur_gate.s_Kur" initial_value="0.845304" name="s_Kur" public_interface="out" units="dimensionless"/>
<variable cmeta:id="i_Kur_sKur_gate.tau_s_Kur" name="tau_s_Kur" units="second"/>
<variable cmeta:id="i_Kur_sKur_gate.s_Kur_infinity" name="s_Kur_infinity" units="dimensionless"/>
<variable cmeta:id="i_Kur_sKur_gate.V" name="V" public_interface="in" units="millivolt"/>
<variable cmeta:id="i_Kur_sKur_gate.time" name="time" public_interface="in" units="second"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>s_Kur</ci>
</apply>
<apply>
<divide/>
<apply>
<minus/>
<ci>s_Kur_infinity</ci>
<ci>s_Kur</ci>
</apply>
<ci>tau_s_Kur</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>s_Kur_infinity</ci>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="millivolt">7.5</cn>
</apply>
<cn cellml:units="millivolt">10</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>tau_s_Kur</ci>
<apply>
<plus/>
<apply>
<divide/>
<cn cellml:units="second">0.59</cn>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="millivolt">60</cn>
</apply>
<cn cellml:units="millivolt">10</cn>
</apply>
</apply>
</apply>
</apply>
<cn cellml:units="second">3.05</cn>
</apply>
</apply>
</math>
</component>
<component name="i_to">
<variable cmeta:id="i_to.i_to" name="i_to" public_interface="out" units="nanoA"/>
<variable cmeta:id="i_to.g_to" initial_value="3.5e-3" name="g_to" units="microS"/>
<variable cmeta:id="i_to.E_K" name="E_K" public_interface="in" units="millivolt"/>
<variable cmeta:id="i_to.time" name="time" private_interface="out" public_interface="in" units="second"/>
<variable cmeta:id="i_to.V" name="V" private_interface="out" public_interface="in" units="millivolt"/>
<variable cmeta:id="i_to.q" name="q" private_interface="in" units="dimensionless"/>
<variable cmeta:id="i_to.r" name="r" private_interface="in" units="dimensionless"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>i_to</ci>
<apply>
<times/>
<ci>g_to</ci>
<apply>
<minus/>
<ci>V</ci>
<ci>E_K</ci>
</apply>
<ci>q</ci>
<ci>r</ci>
</apply>
</apply>
</math>
</component>
<component name="i_to_q_gate">
<variable cmeta:id="i_to_q_gate.q" initial_value="0.430836" name="q" public_interface="out" units="dimensionless"/>
<variable cmeta:id="i_to_q_gate.q_infinity" name="q_infinity" units="dimensionless"/>
<variable cmeta:id="i_to_q_gate.tau_q" name="tau_q" units="second"/>
<variable cmeta:id="i_to_q_gate.V" name="V" public_interface="in" units="millivolt"/>
<variable cmeta:id="i_to_q_gate.time" name="time" public_interface="in" units="second"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>q_infinity</ci>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="millivolt">49</cn>
</apply>
<cn cellml:units="millivolt">13</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>tau_q</ci>
<apply>
<times/>
<cn cellml:units="dimensionless">0.001</cn>
<cn cellml:units="dimensionless">0.6</cn>
<apply>
<plus/>
<apply>
<divide/>
<cn cellml:units="dimensionless">65.17</cn>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="per_second">0.57</cn>
<apply>
<exp/>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="per_millivolt">0.08</cn>
</apply>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="millivolt">44</cn>
</apply>
</apply>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="per_second">0.065</cn>
<apply>
<exp/>
<apply>
<times/>
<cn cellml:units="per_millivolt">0.1</cn>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="millivolt">45.93</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<cn cellml:units="second">10.1</cn>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>q</ci>
</apply>
<apply>
<divide/>
<apply>
<minus/>
<ci>q_infinity</ci>
<ci>q</ci>
</apply>
<ci>tau_q</ci>
</apply>
</apply>
</math>
</component>
<component name="i_to_r_gate">
<variable cmeta:id="i_to_r_gate.r" initial_value="0.014523" name="r" public_interface="out" units="dimensionless"/>
<variable cmeta:id="i_to_r_gate.r_infinity" name="r_infinity" units="dimensionless"/>
<variable cmeta:id="i_to_r_gate.tau_r" name="tau_r" units="second"/>
<variable cmeta:id="i_to_r_gate.V" name="V" public_interface="in" units="millivolt"/>
<variable cmeta:id="i_to_r_gate.time" name="time" public_interface="in" units="second"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>r_infinity</ci>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<minus/>
<ci>V</ci>
<cn cellml:units="millivolt">19.3</cn>
</apply>
</apply>
<cn cellml:units="millivolt">15</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>tau_r</ci>
<apply>
<times/>
<cn cellml:units="dimensionless">0.001</cn>
<cn cellml:units="dimensionless">0.66</cn>
<cn cellml:units="dimensionless">1.4</cn>
<apply>
<plus/>
<apply>
<divide/>
<cn cellml:units="second">15.59</cn>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="dimensionless">1.037</cn>
<apply>
<exp/>
<apply>
<times/>
<cn cellml:units="per_millivolt">0.09</cn>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="millivolt">30.61</cn>
</apply>
</apply>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">0.369</cn>
<apply>
<exp/>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="per_millivolt">0.12</cn>
</apply>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="millivolt">23.84</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<cn cellml:units="second">2.98</cn>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>r</ci>
</apply>
<apply>
<divide/>
<apply>
<minus/>
<ci>r_infinity</ci>
<ci>r</ci>
</apply>
<ci>tau_r</ci>
</apply>
</apply>
</math>
</component>
<component name="i_Kr">
<variable cmeta:id="i_Kr.i_Kr" name="i_Kr" public_interface="out" units="nanoA"/>
<variable cmeta:id="i_Kr.g_Kr" initial_value="0.00424" name="g_Kr" units="microS"/>
<variable cmeta:id="i_Kr.time" name="time" private_interface="out" public_interface="in" units="second"/>
<variable cmeta:id="i_Kr.V" name="V" private_interface="out" public_interface="in" units="millivolt"/>
<variable cmeta:id="i_Kr.E_K" name="E_K" public_interface="in" units="millivolt"/>
<variable cmeta:id="i_Kr.Ko" name="Ko" public_interface="in" units="millimolar"/>
<variable cmeta:id="i_Kr.paS" name="paS" private_interface="in" units="dimensionless"/>
<variable cmeta:id="i_Kr.paF" name="paF" private_interface="in" units="dimensionless"/>
<variable cmeta:id="i_Kr.piy" name="piy" private_interface="in" units="dimensionless"/>
<variable cmeta:id="i_Kr.RTONF" name="RTONF" public_interface="in" units="millivolt"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>i_Kr</ci>
<apply>
<times/>
<ci>g_Kr</ci>
<apply>
<minus/>
<ci>V</ci>
<ci>E_K</ci>
</apply>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="dimensionless">0.9</cn>
<ci>paF</ci>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">0.1</cn>
<ci>paS</ci>
</apply>
</apply>
<ci>piy</ci>
</apply>
</apply>
</math>
</component>
<component name="i_Kr_pa_gate">
<variable cmeta:id="i_Kr_pa_gate.paS" initial_value="0.283185" name="paS" public_interface="out" units="dimensionless"/>
<variable cmeta:id="i_Kr_pa_gate.paF" initial_value="0.011068" name="paF" public_interface="out" units="dimensionless"/>
<variable cmeta:id="i_Kr_pa_gate.pa_infinity" name="pa_infinity" units="dimensionless"/>
<variable cmeta:id="i_Kr_pa_gate.alfapaF" name="alfapaF" units="per_second"/>
<variable cmeta:id="i_Kr_pa_gate.betapaF" name="betapaF" units="per_second"/>
<variable cmeta:id="i_Kr_pa_gate.tau_paS" name="tau_paS" units="second"/>
<variable cmeta:id="i_Kr_pa_gate.tau_paF" name="tau_paF" units="second"/>
<variable cmeta:id="i_Kr_pa_gate.V" name="V" public_interface="in" units="millivolt"/>
<variable cmeta:id="i_Kr_pa_gate.time" name="time" public_interface="in" units="second"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>alfapaF</ci>
<apply>
<divide/>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="millivolt">23.2</cn>
</apply>
</apply>
<cn cellml:units="millivolt">6.6</cn>
</apply>
</apply>
</apply>
</apply>
<apply>
<divide/>
<cn cellml:units="second">0.84655354</cn>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="dimensionless">37.2</cn>
<apply>
<exp/>
<apply>
<divide/>
<ci>V</ci>
<cn cellml:units="millivolt">11.9</cn>
</apply>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">0.96</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<ci>V</ci>
</apply>
<cn cellml:units="millivolt">18.5</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>betapaF</ci>
<apply>
<times/>
<cn cellml:units="dimensionless">4</cn>
<apply>
<minus/>
<apply>
<divide/>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="dimensionless">37.2</cn>
<apply>
<exp/>
<apply>
<divide/>
<ci>V</ci>
<cn cellml:units="millivolt">15.9</cn>
</apply>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">0.96</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<ci>V</ci>
</apply>
<cn cellml:units="millivolt">22.5</cn>
</apply>
</apply>
</apply>
</apply>
<cn cellml:units="second">0.84655354</cn>
</apply>
<apply>
<divide/>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="millivolt">23.2</cn>
</apply>
</apply>
<cn cellml:units="millivolt">10.6</cn>
</apply>
</apply>
</apply>
</apply>
<apply>
<divide/>
<cn cellml:units="second">0.84655354</cn>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="dimensionless">37.2</cn>
<apply>
<exp/>
<apply>
<divide/>
<ci>V</ci>
<cn cellml:units="millivolt">15.9</cn>
</apply>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">0.96</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<ci>V</ci>
</apply>
<cn cellml:units="millivolt">22.5</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>pa_infinity</ci>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="millivolt">10.0144</cn>
</apply>
</apply>
<cn cellml:units="millivolt">7.6607</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>tau_paS</ci>
<apply>
<divide/>
<cn cellml:units="second">0.84655354</cn>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="dimensionless">4.2</cn>
<apply>
<exp/>
<apply>
<divide/>
<ci>V</ci>
<cn cellml:units="millivolt">17</cn>
</apply>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">0.15</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<ci>V</ci>
</apply>
<cn cellml:units="millivolt">21.6</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>tau_paF</ci>
<apply>
<divide/>
<cn cellml:units="second">1</cn>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="dimensionless">30</cn>
<apply>
<exp/>
<apply>
<divide/>
<ci>V</ci>
<cn cellml:units="millivolt">10</cn>
</apply>
</apply>
</apply>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<ci>V</ci>
</apply>
<cn cellml:units="millivolt">12</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>paS</ci>
</apply>
<apply>
<divide/>
<apply>
<minus/>
<ci>pa_infinity</ci>
<ci>paS</ci>
</apply>
<ci>tau_paS</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>paF</ci>
</apply>
<apply>
<divide/>
<apply>
<minus/>
<ci>pa_infinity</ci>
<ci>paF</ci>
</apply>
<ci>tau_paF</ci>
</apply>
</apply>
</math>
</component>
<component name="i_Kr_pi_gate">
<variable cmeta:id="i_Kr_pi_gate.piy" initial_value="0.709051" name="piy" public_interface="out" units="dimensionless"/>
<variable cmeta:id="i_Kr_pi_gate.pi_infinity" name="pi_infinity" units="dimensionless"/>
<variable cmeta:id="i_Kr_pi_gate.tau_pi" name="tau_pi" units="second"/>
<variable cmeta:id="i_Kr_pi_gate.V" name="V" public_interface="in" units="millivolt"/>
<variable cmeta:id="i_Kr_pi_gate.time" name="time" public_interface="in" units="second"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>tau_pi</ci>
<apply>
<divide/>
<cn cellml:units="second">1</cn>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="dimensionless">100</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<ci>V</ci>
</apply>
<cn cellml:units="millivolt">54.645</cn>
</apply>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">656</cn>
<apply>
<exp/>
<apply>
<divide/>
<ci>V</ci>
<cn cellml:units="millivolt">106.157</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>pi_infinity</ci>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="millivolt">28.6</cn>
</apply>
<cn cellml:units="millivolt">17.1</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>piy</ci>
</apply>
<apply>
<divide/>
<apply>
<minus/>
<ci>pi_infinity</ci>
<ci>piy</ci>
</apply>
<ci>tau_pi</ci>
</apply>
</apply>
</math>
</component>
<component name="i_Ks">
<variable cmeta:id="i_Ks.i_Ks" name="i_Ks" public_interface="out" units="nanoA"/>
<variable cmeta:id="i_Ks.g_Ks" name="g_Ks" public_interface="out" units="microS"/>
<variable cmeta:id="i_Ks.g_Ks_" initial_value="0.00065" name="g_Ks_" units="microS"/>
<variable cmeta:id="i_Ks.E_Ks" name="E_Ks" units="millivolt"/>
<variable cmeta:id="i_Ks.time" name="time" private_interface="out" public_interface="in" units="second"/>
<variable cmeta:id="i_Ks.V" name="V" private_interface="out" public_interface="in" units="millivolt"/>
<variable cmeta:id="i_Ks.Ki" name="Ki" public_interface="in" units="millimolar"/>
<variable cmeta:id="i_Ks.Ko" name="Ko" public_interface="in" units="millimolar"/>
<variable cmeta:id="i_Ks.Nai" name="Nai" public_interface="in" units="millimolar"/>
<variable cmeta:id="i_Ks.Nao" name="Nao" public_interface="in" units="millimolar"/>
<variable cmeta:id="i_Ks.RTONF" name="RTONF" public_interface="in" units="millivolt"/>
<variable cmeta:id="i_Ks.n" name="n" private_interface="in" units="dimensionless"/>
<variable cmeta:id="i_Ks.Iso_1_uM" name="Iso_1_uM" private_interface="out" public_interface="in" units="dimensionless"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>g_Ks</ci>
<piecewise>
<piece>
<apply>
<times/>
<cn cellml:units="dimensionless">1.2</cn>
<ci>g_Ks_</ci>
</apply>
<apply>
<gt/>
<ci>Iso_1_uM</ci>
<cn cellml:units="dimensionless">0</cn>
</apply>
</piece>
<otherwise>
<ci>g_Ks_</ci>
</otherwise>
</piecewise>
</apply>
<apply>
<eq/>
<ci>E_Ks</ci>
<apply>
<times/>
<ci>RTONF</ci>
<apply>
<ln/>
<apply>
<divide/>
<apply>
<plus/>
<ci>Ko</ci>
<apply>
<times/>
<cn cellml:units="dimensionless">0.12</cn>
<ci>Nao</ci>
</apply>
</apply>
<apply>
<plus/>
<ci>Ki</ci>
<apply>
<times/>
<cn cellml:units="dimensionless">0.12</cn>
<ci>Nai</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>i_Ks</ci>
<apply>
<times/>
<ci>g_Ks</ci>
<apply>
<minus/>
<ci>V</ci>
<ci>E_Ks</ci>
</apply>
<apply>
<power/>
<ci>n</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
</apply>
</math>
</component>
<component name="i_Ks_n_gate">
<variable cmeta:id="i_Ks_n_gate.n" initial_value="0.1162" name="n" public_interface="out" units="dimensionless"/>
<variable cmeta:id="i_Ks_n_gate.n_infinity" name="n_infinity" units="dimensionless"/>
<variable cmeta:id="i_Ks_n_gate.tau_n" name="tau_n" units="second"/>
<variable cmeta:id="i_Ks_n_gate.V" name="V" public_interface="in" units="millivolt"/>
<variable cmeta:id="i_Ks_n_gate.time" name="time" public_interface="in" units="second"/>
<variable cmeta:id="i_Ks_n_gate.Iso_1_uM" name="Iso_1_uM" public_interface="in" units="dimensionless"/>
<variable cmeta:id="i_Ks_n_gate.Iso_shift" name="Iso_shift" units="millivolt"/>
<variable cmeta:id="i_Ks_n_gate.alpha_n" name="alpha_n" units="per_second"/>
<variable cmeta:id="i_Ks_n_gate.beta_n" name="beta_n" units="per_second"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>Iso_shift</ci>
<piecewise>
<piece>
<apply>
<minus/>
<cn cellml:units="millivolt">14</cn>
</apply>
<apply>
<gt/>
<ci>Iso_1_uM</ci>
<cn cellml:units="dimensionless">0</cn>
</apply>
</piece>
<otherwise>
<cn cellml:units="millivolt">0</cn>
</otherwise>
</piecewise>
</apply>
<apply>
<eq/>
<ci>n_infinity</ci>
<apply>
<root/>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<minus/>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="millivolt">0.6383</cn>
</apply>
<ci>Iso_shift</ci>
</apply>
</apply>
<cn cellml:units="millivolt">10.7071</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>tau_n</ci>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<plus/>
<ci>alpha_n</ci>
<ci>beta_n</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>alpha_n</ci>
<apply>
<divide/>
<cn cellml:units="per_second">28</cn>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<minus/>
<apply>
<minus/>
<ci>V</ci>
<cn cellml:units="millivolt">40</cn>
</apply>
<ci>Iso_shift</ci>
</apply>
</apply>
<cn cellml:units="millivolt">3</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>beta_n</ci>
<apply>
<times/>
<cn cellml:units="per_second">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<minus/>
<apply>
<minus/>
<ci>V</ci>
<ci>Iso_shift</ci>
</apply>
<cn cellml:units="millivolt">5</cn>
</apply>
</apply>
<cn cellml:units="millivolt">25</cn>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>n</ci>
</apply>
<apply>
<divide/>
<apply>
<minus/>
<ci>n_infinity</ci>
<ci>n</ci>
</apply>
<ci>tau_n</ci>
</apply>
</apply>
</math>
</component>
<component name="i_KACh">
<variable cmeta:id="i_KACh.ACh" name="ACh" private_interface="out" public_interface="in" units="millimolar"/>
<variable cmeta:id="i_KACh.i_KACh" name="i_KACh" public_interface="out" units="nanoA"/>
<variable cmeta:id="i_KACh.g_KACh" initial_value="0.00345" name="g_KACh" units="microS"/>
<variable cmeta:id="i_KACh.E_K" name="E_K" public_interface="in" units="millivolt"/>
<variable cmeta:id="i_KACh.time" name="time" private_interface="out" public_interface="in" units="second"/>
<variable cmeta:id="i_KACh.V" name="V" private_interface="out" public_interface="in" units="millivolt"/>
<variable cmeta:id="i_KACh.a" name="a" private_interface="in" units="dimensionless"/>
<variable cmeta:id="i_KACh.ACh_on" initial_value="1" name="ACh_on" units="dimensionless"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>i_KACh</ci>
<piecewise>
<piece>
<apply>
<times/>
<ci>ACh_on</ci>
<ci>g_KACh</ci>
<apply>
<minus/>
<ci>V</ci>
<ci>E_K</ci>
</apply>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="millivolt">20</cn>
</apply>
<cn cellml:units="millivolt">20</cn>
</apply>
</apply>
</apply>
<ci>a</ci>
</apply>
<apply>
<gt/>
<ci>ACh</ci>
<cn cellml:units="millimolar">0</cn>
</apply>
</piece>
<otherwise>
<cn cellml:units="nanoA">0</cn>
</otherwise>
</piecewise>
</apply>
</math>
</component>
<component name="i_KACh_a_gate">
<variable cmeta:id="i_KACh_a_gate.a" initial_value="0.00277" name="a" public_interface="out" units="dimensionless"/>
<variable cmeta:id="i_KACh_a_gate.alpha_a" name="alpha_a" units="per_second"/>
<variable cmeta:id="i_KACh_a_gate.beta_a" name="beta_a" units="per_second"/>
<variable cmeta:id="i_KACh_a_gate.a_infinity" name="a_infinity" units="dimensionless"/>
<variable cmeta:id="i_KACh_a_gate.tau_a" name="tau_a" units="second"/>
<variable cmeta:id="i_KACh_a_gate.ACh" name="ACh" public_interface="in" units="millimolar"/>
<variable cmeta:id="i_KACh_a_gate.V" name="V" public_interface="in" units="millivolt"/>
<variable cmeta:id="i_KACh_a_gate.time" name="time" public_interface="in" units="second"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>alpha_a</ci>
<apply>
<plus/>
<apply>
<divide/>
<apply>
<minus/>
<cn cellml:units="per_second">3.5988</cn>
<cn cellml:units="per_second">0.025641</cn>
</apply>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<divide/>
<cn cellml:units="dimensionless">0.0000012155</cn>
<apply>
<power/>
<apply>
<times/>
<cn cellml:units="per_millimolar">1</cn>
<ci>ACh</ci>
</apply>
<cn cellml:units="dimensionless">1.6951</cn>
</apply>
</apply>
</apply>
</apply>
<cn cellml:units="per_second">0.025641</cn>
</apply>
</apply>
<apply>
<eq/>
<ci>beta_a</ci>
<apply>
<times/>
<cn cellml:units="per_second">10</cn>
<apply>
<exp/>
<apply>
<times/>
<cn cellml:units="per_millivolt">0.0133</cn>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="millivolt">40</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>a_infinity</ci>
<apply>
<divide/>
<ci>alpha_a</ci>
<apply>
<plus/>
<ci>alpha_a</ci>
<ci>beta_a</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>tau_a</ci>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<plus/>
<ci>alpha_a</ci>
<ci>beta_a</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>a</ci>
</apply>
<apply>
<divide/>
<apply>
<minus/>
<ci>a_infinity</ci>
<ci>a</ci>
</apply>
<ci>tau_a</ci>
</apply>
</apply>
</math>
</component>
<group>
<relationship_ref relationship="containment"/>
<component_ref component="Membrane">
<component_ref component="Nai_concentration"/>
<component_ref component="Cell_parameters"/>
<component_ref component="Voltage_clamp"/>
<component_ref component="i_Kur">
<component_ref component="i_Kur_sKur_gate"/>
<component_ref component="i_Kur_rKur_gate"/>
</component_ref>
<component_ref component="Ca_buffering"/>
<component_ref component="Ca_dynamics"/>
<component_ref component="Rate_modulation_experiments"/>
<component_ref component="i_KACh">
<component_ref component="i_KACh_a_gate"/>
</component_ref>
<component_ref component="i_NaK"/>
<component_ref component="i_CaL">
<component_ref component="i_CaL_fCa_gate"/>
<component_ref component="i_CaL_dL_gate"/>
<component_ref component="i_CaL_fL_gate"/>
</component_ref>
<component_ref component="i_NaCa"/>
<component_ref component="i_Kr">
<component_ref component="i_Kr_pi_gate"/>
<component_ref component="i_Kr_pa_gate"/>
</component_ref>
<component_ref component="i_CaT">
<component_ref component="i_CaT_fT_gate"/>
<component_ref component="i_CaT_dT_gate"/>
</component_ref>
<component_ref component="i_Ks">
<component_ref component="i_Ks_n_gate"/>
</component_ref>
<component_ref component="i_to">
<component_ref component="i_to_q_gate"/>
<component_ref component="i_to_r_gate"/>
</component_ref>
<component_ref component="i_Na">
<component_ref component="i_Na_m_gate"/>
<component_ref component="i_Na_h_gate"/>
</component_ref>
<component_ref component="Ionic_values"/>
<component_ref component="Ca_intracellular_fluxes"/>
<component_ref component="i_f">
<component_ref component="i_f_y_gate"/>
</component_ref>
<component_ref component="Ca_SR_release"/>
</component_ref>
</group>
<group>
<relationship_ref relationship="encapsulation"/>
<component_ref component="i_CaL">
<component_ref component="i_CaL_fCa_gate"/>
<component_ref component="i_CaL_dL_gate"/>
<component_ref component="i_CaL_fL_gate"/>
</component_ref>
<component_ref component="i_Kr">
<component_ref component="i_Kr_pi_gate"/>
<component_ref component="i_Kr_pa_gate"/>
</component_ref>
<component_ref component="i_CaT">
<component_ref component="i_CaT_fT_gate"/>
<component_ref component="i_CaT_dT_gate"/>
</component_ref>
<component_ref component="i_Kur">
<component_ref component="i_Kur_sKur_gate"/>
<component_ref component="i_Kur_rKur_gate"/>
</component_ref>
<component_ref component="i_KACh">
<component_ref component="i_KACh_a_gate"/>
</component_ref>
<component_ref component="Membrane">
<component_ref component="Voltage_clamp"/>
</component_ref>
<component_ref component="i_Ks">
<component_ref component="i_Ks_n_gate"/>
</component_ref>
<component_ref component="i_Na">
<component_ref component="i_Na_m_gate"/>
<component_ref component="i_Na_h_gate"/>
</component_ref>
<component_ref component="i_to">
<component_ref component="i_to_q_gate"/>
<component_ref component="i_to_r_gate"/>
</component_ref>
<component_ref component="i_f">
<component_ref component="i_f_y_gate"/>
</component_ref>
</group>
<connection>
<map_components component_1="environment" component_2="i_Ks"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection>
<map_components component_1="Membrane" component_2="Voltage_clamp"/>
<map_variables variable_1="time" variable_2="time"/>
<map_variables variable_1="V_clamp" variable_2="V_clamp"/>
</connection>
<connection>
<map_components component_1="i_f" component_2="i_f_y_gate"/>
<map_variables variable_1="time" variable_2="time"/>
<map_variables variable_1="V" variable_2="V"/>
<map_variables variable_1="ACh" variable_2="ACh"/>
<map_variables variable_1="Iso_1_uM" variable_2="Iso_1_uM"/>
<map_variables variable_1="y" variable_2="y"/>
</connection>
<connection>
<map_components component_1="i_to_r_gate" component_2="i_to"/>
<map_variables variable_1="r" variable_2="r"/>
</connection>
<connection>
<map_components component_1="Nai_concentration" component_2="i_NaCa"/>
<map_variables variable_1="Nai" variable_2="Nai"/>
<map_variables variable_1="i_NaCa" variable_2="i_NaCa"/>
</connection>
<connection>
<map_components component_1="Ionic_values" component_2="i_f"/>
<map_variables variable_1="E_Na" variable_2="E_Na"/>
<map_variables variable_1="E_K" variable_2="E_K"/>
<map_variables variable_1="Ko" variable_2="Ko"/>
</connection>
<connection>
<map_components component_1="Ca_intracellular_fluxes" component_2="Ca_dynamics"/>
<map_variables variable_1="j_tr" variable_2="j_tr"/>
<map_variables variable_1="j_up" variable_2="j_up"/>
<map_variables variable_1="j_Ca_dif" variable_2="j_Ca_dif"/>
</connection>
<connection>
<map_components component_1="i_NaK" component_2="Membrane"/>
<map_variables variable_1="i_NaK" variable_2="i_NaK"/>
</connection>
<connection>
<map_components component_1="environment" component_2="i_to"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection>
<map_components component_1="Ionic_values" component_2="i_Kur"/>
<map_variables variable_1="E_K" variable_2="E_K"/>
</connection>
<connection>
<map_components component_1="Nai_concentration" component_2="i_NaK"/>
<map_variables variable_1="Nai" variable_2="Nai"/>
<map_variables variable_1="i_NaK" variable_2="i_NaK"/>
</connection>
<connection>
<map_components component_1="i_Na_h_gate" component_2="i_Na"/>
<map_variables variable_1="h" variable_2="h"/>
</connection>
<connection>
<map_components component_1="Cell_parameters" component_2="Nai_concentration"/>
<map_variables variable_1="V_sub" variable_2="V_sub"/>
<map_variables variable_1="V_i" variable_2="V_i"/>
</connection>
<connection>
<map_components component_1="i_Ks" component_2="Membrane"/>
<map_variables variable_1="i_Ks" variable_2="i_Ks"/>
</connection>
<connection>
<map_components component_1="Ionic_values" component_2="i_to"/>
<map_variables variable_1="E_K" variable_2="E_K"/>
</connection>
<connection>
<map_components component_1="environment" component_2="i_CaT"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection>
<map_components component_1="i_CaT" component_2="Membrane"/>
<map_variables variable_1="i_CaT" variable_2="i_CaT"/>
</connection>
<connection>
<map_components component_1="Ionic_values" component_2="i_Ks"/>
<map_variables variable_1="Nao" variable_2="Nao"/>
<map_variables variable_1="Ko" variable_2="Ko"/>
<map_variables variable_1="Ki" variable_2="Ki"/>
</connection>
<connection>
<map_components component_1="Nai_concentration" component_2="i_Ks"/>
<map_variables variable_1="Nai" variable_2="Nai"/>
</connection>
<connection>
<map_components component_1="i_Kur" component_2="i_Kur_rKur_gate"/>
<map_variables variable_1="time" variable_2="time"/>
<map_variables variable_1="V" variable_2="V"/>
<map_variables variable_1="r_Kur" variable_2="r_Kur"/>
</connection>
<connection>
<map_components component_1="i_CaL" component_2="Membrane"/>
<map_variables variable_1="i_CaL" variable_2="i_CaL"/>
</connection>
<connection>
<map_components component_1="Membrane" component_2="i_CaT"/>
<map_variables variable_1="RTONF" variable_2="RTONF"/>
<map_variables variable_1="V" variable_2="V"/>
</connection>
<connection>
<map_components component_1="Rate_modulation_experiments" component_2="i_Ks"/>
<map_variables variable_1="Iso_1_uM" variable_2="Iso_1_uM"/>
</connection>
<connection>
<map_components component_1="i_Ks" component_2="i_Ks_n_gate"/>
<map_variables variable_1="time" variable_2="time"/>
<map_variables variable_1="V" variable_2="V"/>
<map_variables variable_1="Iso_1_uM" variable_2="Iso_1_uM"/>
<map_variables variable_1="n" variable_2="n"/>
</connection>
<connection>
<map_components component_1="Membrane" component_2="i_Na"/>
<map_variables variable_1="RTONF" variable_2="RTONF"/>
<map_variables variable_1="V" variable_2="V"/>
<map_variables variable_1="i_Na" variable_2="i_Na"/>
</connection>
<connection>
<map_components component_1="i_f" component_2="Membrane"/>
<map_variables variable_1="i_f" variable_2="i_f"/>
</connection>
<connection>
<map_components component_1="Membrane" component_2="i_Kr"/>
<map_variables variable_1="RTONF" variable_2="RTONF"/>
<map_variables variable_1="V" variable_2="V"/>
<map_variables variable_1="i_Kr" variable_2="i_Kr"/>
</connection>
<connection>
<map_components component_1="Ca_dynamics" component_2="Ca_intracellular_fluxes"/>
<map_variables variable_1="Cai" variable_2="Cai"/>
<map_variables variable_1="Ca_jsr" variable_2="Ca_jsr"/>
<map_variables variable_1="Ca_sub" variable_2="Ca_sub"/>
<map_variables variable_1="Ca_nsr" variable_2="Ca_nsr"/>
</connection>
<connection>
<map_components component_1="Membrane" component_2="Nai_concentration"/>
<map_variables variable_1="F" variable_2="F"/>
</connection>
<connection>
<map_components component_1="i_Na" component_2="i_Na_m_gate"/>
<map_variables variable_1="V" variable_2="V"/>
<map_variables variable_1="time" variable_2="time"/>
<map_variables variable_1="m" variable_2="m"/>
</connection>
<connection>
<map_components component_1="environment" component_2="Membrane"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection>
<map_components component_1="Ca_dynamics" component_2="i_CaL"/>
<map_variables variable_1="Ca_sub" variable_2="Ca_sub"/>
<map_variables variable_1="i_siCa" variable_2="i_siCa"/>
</connection>
<connection>
<map_components component_1="i_KACh" component_2="Membrane"/>
<map_variables variable_1="i_KACh" variable_2="i_KACh"/>
</connection>
<connection>
<map_components component_1="i_CaT_fT_gate" component_2="i_CaT"/>
<map_variables variable_1="fT" variable_2="fT"/>
</connection>
<connection>
<map_components component_1="i_CaT" component_2="i_CaT_fT_gate"/>
<map_variables variable_1="V" variable_2="V"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection>
<map_components component_1="Ionic_values" component_2="i_NaCa"/>
<map_variables variable_1="Cao" variable_2="Cao"/>
<map_variables variable_1="Nao" variable_2="Nao"/>
</connection>
<connection>
<map_components component_1="i_Kur_sKur_gate" component_2="i_Kur"/>
<map_variables variable_1="s_Kur" variable_2="s_Kur"/>
</connection>
<connection>
<map_components component_1="i_CaL" component_2="Nai_concentration"/>
<map_variables variable_1="i_siNa" variable_2="i_siNa"/>
</connection>
<connection>
<map_components component_1="environment" component_2="Nai_concentration"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection>
<map_components component_1="Rate_modulation_experiments" component_2="i_KACh"/>
<map_variables variable_1="ACh" variable_2="ACh"/>
</connection>
<connection>
<map_components component_1="i_Na" component_2="i_Na_h_gate"/>
<map_variables variable_1="V" variable_2="V"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection>
<map_components component_1="Ionic_values" component_2="i_Na"/>
<map_variables variable_1="Nao" variable_2="Nao"/>
<map_variables variable_1="Ko" variable_2="Ko"/>
<map_variables variable_1="Ki" variable_2="Ki"/>
</connection>
<connection>
<map_components component_1="Ca_dynamics" component_2="Ca_buffering"/>
<map_variables variable_1="Cai" variable_2="Cai"/>
<map_variables variable_1="Ca_jsr" variable_2="Ca_jsr"/>
<map_variables variable_1="Ca_sub" variable_2="Ca_sub"/>
</connection>
<connection>
<map_components component_1="Ca_dynamics" component_2="i_NaCa"/>
<map_variables variable_1="Ca_sub" variable_2="Ca_sub"/>
<map_variables variable_1="i_NaCa" variable_2="i_NaCa"/>
</connection>
<connection>
<map_components component_1="environment" component_2="Ca_dynamics"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection>
<map_components component_1="Membrane" component_2="i_NaK"/>
<map_variables variable_1="V" variable_2="V"/>
</connection>
<connection>
<map_components component_1="Nai_concentration" component_2="i_CaL"/>
<map_variables variable_1="Nai" variable_2="Nai"/>
</connection>
<connection>
<map_components component_1="i_Kr_pa_gate" component_2="i_Kr"/>
<map_variables variable_1="paF" variable_2="paF"/>
<map_variables variable_1="paS" variable_2="paS"/>
</connection>
<connection>
<map_components component_1="environment" component_2="i_KACh"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection>
<map_components component_1="Rate_modulation_experiments" component_2="Ca_intracellular_fluxes"/>
<map_variables variable_1="ACh" variable_2="ACh"/>
<map_variables variable_1="Iso_1_uM" variable_2="Iso_1_uM"/>
</connection>
<connection>
<map_components component_1="i_Na" component_2="Nai_concentration"/>
<map_variables variable_1="i_Na" variable_2="i_Na"/>
</connection>
<connection>
<map_components component_1="environment" component_2="i_CaL"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection>
<map_components component_1="i_CaL" component_2="i_CaL_fCa_gate"/>
<map_variables variable_1="V" variable_2="V"/>
<map_variables variable_1="Ca_sub" variable_2="Ca_sub"/>
<map_variables variable_1="time" variable_2="time"/>
<map_variables variable_1="fCa" variable_2="fCa"/>
</connection>
<connection>
<map_components component_1="Cell_parameters" component_2="Ca_dynamics"/>
<map_variables variable_1="V_sub" variable_2="V_sub"/>
<map_variables variable_1="V_jsr" variable_2="V_jsr"/>
<map_variables variable_1="V_nsr" variable_2="V_nsr"/>
<map_variables variable_1="V_i" variable_2="V_i"/>
</connection>
<connection>
<map_components component_1="i_Kur" component_2="i_Kur_sKur_gate"/>
<map_variables variable_1="time" variable_2="time"/>
<map_variables variable_1="V" variable_2="V"/>
</connection>
<connection>
<map_components component_1="Ca_dynamics" component_2="i_CaT"/>
<map_variables variable_1="Ca_sub" variable_2="Ca_sub"/>
<map_variables variable_1="i_CaT" variable_2="i_CaT"/>
</connection>
<connection>
<map_components component_1="i_Kr" component_2="i_Kr_pi_gate"/>
<map_variables variable_1="V" variable_2="V"/>
<map_variables variable_1="time" variable_2="time"/>
<map_variables variable_1="piy" variable_2="piy"/>
</connection>
<connection>
<map_components component_1="i_to_q_gate" component_2="i_to"/>
<map_variables variable_1="q" variable_2="q"/>
</connection>
<connection>
<map_components component_1="i_CaL" component_2="i_CaL_fL_gate"/>
<map_variables variable_1="V" variable_2="V"/>
<map_variables variable_1="time" variable_2="time"/>
<map_variables variable_1="fL" variable_2="fL"/>
</connection>
<connection>
<map_components component_1="Rate_modulation_experiments" component_2="i_f"/>
<map_variables variable_1="ACh" variable_2="ACh"/>
<map_variables variable_1="Iso_1_uM" variable_2="Iso_1_uM"/>
</connection>
<connection>
<map_components component_1="environment" component_2="i_f"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection>
<map_components component_1="Ca_dynamics" component_2="Ionic_values"/>
<map_variables variable_1="Ca_sub" variable_2="Ca_sub"/>
</connection>
<connection>
<map_components component_1="Ca_SR_release" component_2="Ca_dynamics"/>
<map_variables variable_1="j_SRCarel" variable_2="j_SRCarel"/>
<map_variables variable_1="Ca_jsr" variable_2="Ca_jsr"/>
<map_variables variable_1="Ca_sub" variable_2="Ca_sub"/>
</connection>
<connection>
<map_components component_1="i_Kr" component_2="i_Kr_pa_gate"/>
<map_variables variable_1="V" variable_2="V"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection>
<map_components component_1="i_KACh_a_gate" component_2="i_KACh"/>
<map_variables variable_1="a" variable_2="a"/>
</connection>
<connection>
<map_components component_1="i_f" component_2="Nai_concentration"/>
<map_variables variable_1="i_fNa" variable_2="i_fNa"/>
</connection>
<connection>
<map_components component_1="environment" component_2="i_Na"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection>
<map_components component_1="Rate_modulation_experiments" component_2="i_NaK"/>
<map_variables variable_1="Iso_1_uM" variable_2="Iso_1_uM"/>
</connection>
<connection>
<map_components component_1="Rate_modulation_experiments" component_2="i_CaL"/>
<map_variables variable_1="ACh" variable_2="ACh"/>
<map_variables variable_1="Iso_1_uM" variable_2="Iso_1_uM"/>
</connection>
<connection>
<map_components component_1="Membrane" component_2="i_Ks"/>
<map_variables variable_1="RTONF" variable_2="RTONF"/>
<map_variables variable_1="V" variable_2="V"/>
</connection>
<connection>
<map_components component_1="Ionic_values" component_2="i_Kr"/>
<map_variables variable_1="E_K" variable_2="E_K"/>
<map_variables variable_1="Ko" variable_2="Ko"/>
</connection>
<connection>
<map_components component_1="Membrane" component_2="Ionic_values"/>
<map_variables variable_1="RTONF" variable_2="RTONF"/>
</connection>
<connection>
<map_components component_1="Membrane" component_2="i_Kur"/>
<map_variables variable_1="V" variable_2="V"/>
<map_variables variable_1="i_Kur" variable_2="i_Kur"/>
</connection>
<connection>
<map_components component_1="environment" component_2="Ca_SR_release"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection>
<map_components component_1="Ionic_values" component_2="i_CaT"/>
<map_variables variable_1="Cao" variable_2="Cao"/>
</connection>
<connection>
<map_components component_1="i_to" component_2="i_to_q_gate"/>
<map_variables variable_1="V" variable_2="V"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection>
<map_components component_1="Ionic_values" component_2="i_CaL"/>
<map_variables variable_1="Cao" variable_2="Cao"/>
<map_variables variable_1="Nao" variable_2="Nao"/>
<map_variables variable_1="Ko" variable_2="Ko"/>
<map_variables variable_1="Ki" variable_2="Ki"/>
</connection>
<connection>
<map_components component_1="i_CaL_dL_gate" component_2="i_CaL"/>
<map_variables variable_1="dL" variable_2="dL"/>
</connection>
<connection>
<map_components component_1="Membrane" component_2="i_CaL"/>
<map_variables variable_1="RTONF" variable_2="RTONF"/>
<map_variables variable_1="V" variable_2="V"/>
</connection>
<connection>
<map_components component_1="i_to" component_2="i_to_r_gate"/>
<map_variables variable_1="V" variable_2="V"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection>
<map_components component_1="environment" component_2="i_NaK"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection>
<map_components component_1="Membrane" component_2="i_KACh"/>
<map_variables variable_1="V" variable_2="V"/>
</connection>
<connection>
<map_components component_1="Ionic_values" component_2="i_NaK"/>
<map_variables variable_1="E_Na" variable_2="E_Na"/>
<map_variables variable_1="Ko" variable_2="Ko"/>
</connection>
<connection>
<map_components component_1="Membrane" component_2="i_f"/>
<map_variables variable_1="RTONF" variable_2="RTONF"/>
<map_variables variable_1="V" variable_2="V"/>
</connection>
<connection>
<map_components component_1="Nai_concentration" component_2="Ionic_values"/>
<map_variables variable_1="Nai" variable_2="Nai"/>
</connection>
<connection>
<map_components component_1="i_NaCa" component_2="Membrane"/>
<map_variables variable_1="i_NaCa" variable_2="i_NaCa"/>
</connection>
<connection>
<map_components component_1="i_KACh" component_2="i_KACh_a_gate"/>
<map_variables variable_1="ACh" variable_2="ACh"/>
<map_variables variable_1="time" variable_2="time"/>
<map_variables variable_1="V" variable_2="V"/>
</connection>
<connection>
<map_components component_1="environment" component_2="i_NaCa"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection>
<map_components component_1="environment" component_2="i_Kr"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection>
<map_components component_1="Ca_buffering" component_2="Ca_dynamics"/>
<map_variables variable_1="delta_fTMC" variable_2="delta_fTMC"/>
<map_variables variable_1="delta_fCMi" variable_2="delta_fCMi"/>
<map_variables variable_1="delta_fCMs" variable_2="delta_fCMs"/>
<map_variables variable_1="CQ_tot" variable_2="CQ_tot"/>
<map_variables variable_1="TC_tot" variable_2="TC_tot"/>
<map_variables variable_1="CM_tot" variable_2="CM_tot"/>
<map_variables variable_1="TMC_tot" variable_2="TMC_tot"/>
<map_variables variable_1="delta_fCQ" variable_2="delta_fCQ"/>
<map_variables variable_1="delta_fTC" variable_2="delta_fTC"/>
</connection>
<connection>
<map_components component_1="i_CaT" component_2="i_CaT_dT_gate"/>
<map_variables variable_1="V" variable_2="V"/>
<map_variables variable_1="time" variable_2="time"/>
<map_variables variable_1="dT" variable_2="dT"/>
</connection>
<connection>
<map_components component_1="Membrane" component_2="i_NaCa"/>
<map_variables variable_1="RTONF" variable_2="RTONF"/>
<map_variables variable_1="V" variable_2="V"/>
</connection>
<connection>
<map_components component_1="i_CaL" component_2="i_CaL_dL_gate"/>
<map_variables variable_1="V" variable_2="V"/>
<map_variables variable_1="Iso_1_uM" variable_2="Iso_1_uM"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection>
<map_components component_1="Membrane" component_2="i_to"/>
<map_variables variable_1="V" variable_2="V"/>
<map_variables variable_1="i_to" variable_2="i_to"/>
</connection>
<connection>
<map_components component_1="Ionic_values" component_2="i_KACh"/>
<map_variables variable_1="E_K" variable_2="E_K"/>
</connection>
<connection>
<map_components component_1="environment" component_2="i_Kur"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection>
<map_components component_1="environment" component_2="Ca_buffering"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection>
<map_components component_1="Membrane" component_2="Ca_dynamics"/>
<map_variables variable_1="F" variable_2="F"/>
</connection>
<connection>
<map_components component_1="Nai_concentration" component_2="i_Na"/>
<map_variables variable_1="Nai" variable_2="Nai"/>
</connection>
<connection>
<map_components component_1="Membrane" component_2="Brain_stem"/>
<map_variables variable_1="time" variable_2="t"/>
<map_variables variable_1="g" variable_2="g"/>
</connection>
</model>
