<?xml version="1.0" encoding="utf-8"?>
<!--
This CellML file was generated on 10.02.2012 at 13:45:07 using:
COR (0.9.31.1409)
Copyright 2002-2012 Dr Alan Garny
http://cor.physiol.ox.ac.uk/ - cor@physiol.ox.ac.uk
CellML 1.0 was used to generate this model
http://www.cellml.org/
-->
<model xmlns:dc="http://purl.org/dc/elements/1.1/" xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#" xmlns:bqs="http://www.cellml.org/bqs/1.0#" xmlns:dcterms="http://purl.org/dc/terms/" xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" xmlns:ns7="http://www.cellml.org/metadata/simulation/1.0#" xmlns:xlink="http://www.w3.org/1999/xlink" name="ekaterinburg_rat" cmeta:id="ekaterinburg_rat" xmlns="http://www.cellml.org/cellml/1.0#" xmlns:cellml="http://www.cellml.org/cellml/1.0#" xmlns:cmeta="http://www.cellml.org/metadata/1.0#">
<units name="per_millimolar">
<unit units="millimolar" exponent="-1"/>
</units>
<units name="per_second2">
<unit units="second" exponent="-2"/>
</units>
<units name="micrometre_per_second">
<unit units="micrometre"/>
<unit units="per_second"/>
</units>
<units name="micrometre_per_second2">
<unit units="micrometre"/>
<unit units="per_second2"/>
</units>
<units name="second_per_micrometre">
<unit units="second"/>
<unit units="per_micrometre"/>
</units>
<units name="millinewton">
<unit units="newton" prefix="milli"/>
</units>
<units name="millinewton_per_second">
<unit units="millinewton"/>
<unit units="per_second"/>
</units>
<units name="millinewton_second_per_micrometre">
<unit units="millinewton"/>
<unit units="second_per_micrometre"/>
</units>
<units name="micrometre">
<unit units="metre" prefix="micro"/>
</units>
<units name="per_micrometre">
<unit units="micrometre" exponent="-1"/>
</units>
<units name="per_second">
<unit units="second" exponent="-1"/>
</units>
<units name="millivolt">
<unit units="volt" prefix="milli"/>
</units>
<units name="per_millivolt">
<unit units="volt" prefix="milli" exponent="-1"/>
</units>
<units name="per_millivolt_second">
<unit units="millivolt" exponent="-1"/>
<unit units="second" exponent="-1"/>
</units>
<units name="microS">
<unit units="siemens" prefix="micro"/>
</units>
<units name="microF">
<unit units="farad" prefix="micro"/>
</units>
<units name="nanoA">
<unit units="ampere" prefix="nano"/>
</units>
<units name="mA_nA">
<unit units="ampere" prefix="milli"/>
<unit units="nanoA" exponent="-1"/>
</units>
<units name="nanoA_per_millimolar">
<unit units="ampere" prefix="nano"/>
<unit units="millimolar" exponent="-1"/>
</units>
<units name="millimolar">
<unit units="mole" prefix="milli"/>
<unit units="litre" exponent="-1"/>
</units>
<units name="millimolar4">
<unit units="millimolar" exponent="4"/>
</units>
<units name="millimolar_per_second">
<unit units="millimolar"/>
<unit units="second" exponent="-1"/>
</units>
<units name="per_millimolar_second">
<unit units="millimolar" exponent="-1"/>
<unit units="second" exponent="-1"/>
</units>
<units name="joule_per_kilomole_kelvin">
<unit units="joule"/>
<unit units="mole" prefix="kilo" exponent="-1"/>
<unit units="kelvin" exponent="-1"/>
</units>
<units name="coulomb_per_mole">
<unit units="coulomb"/>
<unit units="mole" exponent="-1"/>
</units>
<units name="micrometre3">
<unit units="metre" prefix="micro" exponent="3"/>
</units>
<units name="litre_micrometre3">
<unit units="litre"/>
<unit units="micrometre3" exponent="-1"/>
</units>
<component name="environment">
<variable name="time" units="second" public_interface="out"/>
</component>
<component name="parameters">
<variable name="isotonic" units="dimensionless" initial_value="0" public_interface="out"/>
<variable name="alpha_1" units="per_micrometre" initial_value="14.6" public_interface="out"/>
<variable name="beta_1" units="millinewton" initial_value="0.84" public_interface="out"/>
<variable name="alpha_2" units="per_micrometre" initial_value="14.6" public_interface="out"/>
<variable name="beta_2" units="millinewton" initial_value="0.0018" public_interface="out"/>
<variable name="alpha_3" units="per_micrometre" initial_value="48" public_interface="out"/>
<variable name="beta_3" units="millinewton" initial_value="0.015" public_interface="out"/>
<variable name="lambda" units="millinewton" initial_value="30" public_interface="out"/>
<variable name="mu" units="dimensionless" initial_value="3" public_interface="out"/>
<variable name="k_mu" units="dimensionless" initial_value="0.6" public_interface="out"/>
<variable name="kappa" units="dimensionless" initial_value="0.705" public_interface="out"/>
<variable name="kappa_0" units="dimensionless" initial_value="3" public_interface="out"/>
<variable name="m_0" units="dimensionless" initial_value="0.9" public_interface="out"/>
<variable name="v_max" units="micrometre_per_second" initial_value="5.5" public_interface="out"/>
<variable name="a" units="dimensionless" initial_value="0.25" public_interface="out"/>
<variable name="d_h" units="dimensionless" initial_value="0.5" public_interface="out"/>
<variable name="alpha_P" units="dimensionless" initial_value="4" public_interface="out"/>
</component>
<component name="isotonic">
<variable name="isotonic" units="dimensionless" public_interface="in"/>
<variable name="l" units="micrometre" public_interface="in"/>
<variable name="F_muscle" units="millinewton" public_interface="in"/>
<variable name="F_afterload" units="millinewton" initial_value="2"/>
<variable name="isotonic_mode" units="dimensionless" public_interface="out"/>
<variable name="l_0" units="micrometre" public_interface="in"/>
<variable name="stim_period" units="second" public_interface="in"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>isotonic_mode</ci>
<piecewise>
<piece>
<cn cellml:units="dimensionless">1</cn>
<apply>
<and/>
<apply>
<eq/>
<ci>isotonic</ci>
<cn cellml:units="dimensionless">1</cn>
</apply>
<apply>
<gt/>
<ci>F_muscle</ci>
<ci>F_afterload</ci>
</apply>
<apply>
<leq/>
<ci>l</ci>
<apply>
<times/>
<ci>l_0</ci>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<cn type="e-notation" cellml:units="dimensionless">1 <sep/>
-4</cn>
</apply>
</apply>
</apply>
</apply>
</piece>
<otherwise>
<cn cellml:units="dimensionless">0</cn>
</otherwise>
</piecewise>
</apply>
</math>
</component>
<component name="parameters_izakov_et_al_1991">
<variable name="S_0" units="micrometre" initial_value="1.14" public_interface="out"/>
<variable name="q_v" units="per_second" public_interface="out"/>
<variable name="q_1" units="per_second" initial_value="17.3"/>
<variable name="q_2" units="per_second" initial_value="259"/>
<variable name="q_3" units="per_second" initial_value="17.3"/>
<variable name="q_4" units="per_second" initial_value="15" public_interface="out"/>
<variable name="x_st" units="dimensionless" initial_value="0.964285"/>
<variable name="v_st" units="micrometre_per_second" public_interface="out"/>
<variable name="v_1" units="micrometre_per_second" public_interface="out"/>
<variable name="alpha_G" units="dimensionless" initial_value="1" public_interface="out"/>
<variable name="k_A" units="per_millimolar" initial_value="40" public_interface="out"/>
<variable name="v" units="micrometre_per_second" public_interface="in"/>
<variable name="v_max" units="micrometre_per_second" public_interface="in"/>
<variable name="alpha_Q" units="dimensionless" initial_value="10"/>
<variable name="beta_Q" units="dimensionless" initial_value="5"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>v_st</ci>
<apply>
<times/>
<ci>x_st</ci>
<ci>v_max</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>q_v</ci>
<piecewise>
<piece>
<apply>
<minus/>
<ci>q_1</ci>
<apply>
<divide/>
<apply>
<times/>
<ci>q_2</ci>
<ci>v</ci>
</apply>
<ci>v_max</ci>
</apply>
</apply>
<apply>
<leq/>
<ci>v</ci>
<cn cellml:units="micrometre_per_second">0</cn>
</apply>
</piece>
<piece>
<apply>
<plus/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<ci>q_4</ci>
<ci>q_3</ci>
</apply>
<ci>v</ci>
</apply>
<ci>v_st</ci>
</apply>
<ci>q_3</ci>
</apply>
<apply>
<and/>
<apply>
<leq/>
<ci>v</ci>
<ci>v_st</ci>
</apply>
<apply>
<lt/>
<cn cellml:units="micrometre_per_second">0</cn>
<ci>v</ci>
</apply>
</apply>
</piece>
<otherwise>
<apply>
<divide/>
<ci>q_4</ci>
<apply>
<power/>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<divide/>
<apply>
<times/>
<ci>beta_Q</ci>
<apply>
<minus/>
<ci>v</ci>
<ci>v_st</ci>
</apply>
</apply>
<ci>v_max</ci>
</apply>
</apply>
<ci>alpha_Q</ci>
</apply>
</apply>
</otherwise>
</piecewise>
</apply>
<apply>
<eq/>
<ci>v_1</ci>
<apply>
<divide/>
<ci>v_max</ci>
<cn cellml:units="dimensionless">10</cn>
</apply>
</apply>
</math>
</component>
<component name="force">
<variable name="F_CE" units="millinewton"/>
<variable name="F_muscle" units="millinewton" public_interface="out"/>
<variable name="F_XSE" units="millinewton"/>
<variable name="F_SE" units="millinewton"/>
<variable name="F_PE" units="millinewton"/>
<variable name="N" units="dimensionless" public_interface="in"/>
<variable name="k_P_vis" units="millinewton_second_per_micrometre" public_interface="in"/>
<variable name="k_S_vis" units="millinewton_second_per_micrometre" public_interface="in"/>
<variable name="w" units="micrometre_per_second" public_interface="in"/>
<variable name="v" units="micrometre_per_second" public_interface="in"/>
<variable name="l_1" units="micrometre" public_interface="in"/>
<variable name="l_2" units="micrometre" public_interface="in"/>
<variable name="l_3" units="micrometre" public_interface="in"/>
<variable name="p_v" units="dimensionless" public_interface="in"/>
<variable name="alpha_1" units="per_micrometre" public_interface="in"/>
<variable name="beta_1" units="millinewton" public_interface="in"/>
<variable name="alpha_2" units="per_micrometre" public_interface="in"/>
<variable name="beta_2" units="millinewton" public_interface="in"/>
<variable name="alpha_3" units="per_micrometre" public_interface="in"/>
<variable name="beta_3" units="millinewton" public_interface="in"/>
<variable name="lambda" units="millinewton" public_interface="in"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>F_CE</ci>
<apply>
<times/>
<ci>lambda</ci>
<ci>p_v</ci>
<ci>N</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>F_muscle</ci>
<ci>F_XSE</ci>
</apply>
<apply>
<eq/>
<ci>F_SE</ci>
<apply>
<times/>
<ci>beta_1</ci>
<apply>
<minus/>
<apply>
<exp/>
<apply>
<times/>
<ci>alpha_1</ci>
<apply>
<minus/>
<ci>l_2</ci>
<ci>l_1</ci>
</apply>
</apply>
</apply>
<cn cellml:units="dimensionless">1</cn>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>F_PE</ci>
<apply>
<times/>
<ci>beta_2</ci>
<apply>
<minus/>
<apply>
<exp/>
<apply>
<times/>
<ci>alpha_2</ci>
<ci>l_2</ci>
</apply>
</apply>
<cn cellml:units="dimensionless">1</cn>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>F_XSE</ci>
<apply>
<times/>
<ci>beta_3</ci>
<apply>
<minus/>
<apply>
<exp/>
<apply>
<times/>
<ci>alpha_3</ci>
<ci>l_3</ci>
</apply>
</apply>
<cn cellml:units="dimensionless">1</cn>
</apply>
</apply>
</apply>
</math>
</component>
<component name="crossbridge_kinetics">
<variable name="N" units="dimensionless" initial_value="7.455e-8" public_interface="out"/>
<variable name="K_kappa" units="per_second" public_interface="out"/>
<variable name="M_A" units="dimensionless"/>
<variable name="n_1" units="dimensionless"/>
<variable name="L_oz" units="dimensionless" public_interface="out"/>
<variable name="k_p_v" units="per_second"/>
<variable name="k_m_v" units="per_second"/>
<variable name="time" units="second" public_interface="in"/>
<variable name="A" units="millimolar" public_interface="in"/>
<variable name="l_1" units="micrometre" public_interface="in"/>
<variable name="A_tot" units="millimolar" public_interface="in"/>
<variable name="mu" units="dimensionless" public_interface="in"/>
<variable name="k_mu" units="dimensionless" public_interface="in"/>
<variable name="kappa" units="dimensionless" public_interface="in"/>
<variable name="kappa_0" units="dimensionless" public_interface="in"/>
<variable name="m_0" units="dimensionless" public_interface="in"/>
<variable name="S_0" units="micrometre" public_interface="in"/>
<variable name="q_v" units="per_second" public_interface="in"/>
<variable name="v" units="micrometre_per_second" public_interface="in"/>
<variable name="v_st" units="micrometre_per_second" public_interface="in"/>
<variable name="G_star" units="dimensionless" public_interface="in"/>
<variable name="P_star" units="dimensionless" public_interface="in"/>
<variable name="q_4" units="per_second" public_interface="in"/>
<variable name="g_1" units="per_micrometre" initial_value="0.6"/>
<variable name="g_2" units="dimensionless" initial_value="0.52"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>M_A</ci>
<apply>
<divide/>
<apply>
<times/>
<apply>
<power/>
<apply>
<divide/>
<ci>A</ci>
<ci>A_tot</ci>
</apply>
<ci>mu</ci>
</apply>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<power/>
<ci>k_mu</ci>
<ci>mu</ci>
</apply>
</apply>
</apply>
<apply>
<plus/>
<apply>
<power/>
<apply>
<divide/>
<ci>A</ci>
<ci>A_tot</ci>
</apply>
<ci>mu</ci>
</apply>
<apply>
<power/>
<ci>k_mu</ci>
<ci>mu</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>n_1</ci>
<piecewise>
<piece>
<cn cellml:units="dimensionless">0</cn>
<apply>
<lt/>
<apply>
<plus/>
<apply>
<times/>
<ci>g_1</ci>
<ci>l_1</ci>
</apply>
<ci>g_2</ci>
</apply>
<cn cellml:units="dimensionless">0</cn>
</apply>
</piece>
<piece>
<apply>
<plus/>
<apply>
<times/>
<ci>g_1</ci>
<ci>l_1</ci>
</apply>
<ci>g_2</ci>
</apply>
<apply>
<lt/>
<apply>
<plus/>
<apply>
<times/>
<ci>g_1</ci>
<ci>l_1</ci>
</apply>
<ci>g_2</ci>
</apply>
<cn cellml:units="dimensionless">1</cn>
</apply>
</piece>
<otherwise>
<cn cellml:units="dimensionless">1</cn>
</otherwise>
</piecewise>
</apply>
<apply>
<eq/>
<ci>L_oz</ci>
<piecewise>
<piece>
<apply>
<divide/>
<apply>
<plus/>
<ci>l_1</ci>
<ci>S_0</ci>
</apply>
<apply>
<plus/>
<cn cellml:units="micrometre">0.46</cn>
<ci>S_0</ci>
</apply>
</apply>
<apply>
<gt/>
<ci>l_1</ci>
<cn cellml:units="micrometre">0.55</cn>
</apply>
</piece>
<otherwise>
<apply>
<times/>
<apply>
<plus/>
<ci>S_0</ci>
<cn cellml:units="micrometre">0.55</cn>
</apply>
<cn cellml:units="per_micrometre">1</cn>
</apply>
</otherwise>
</piecewise>
</apply>
<apply>
<eq/>
<ci>k_p_v</ci>
<apply>
<times/>
<ci>kappa</ci>
<ci>kappa_0</ci>
<ci>q_v</ci>
<ci>m_0</ci>
<ci>G_star</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>k_m_v</ci>
<apply>
<times/>
<ci>kappa_0</ci>
<ci>q_v</ci>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<times/>
<ci>kappa</ci>
<ci>m_0</ci>
<ci>G_star</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>K_kappa</ci>
<apply>
<minus/>
<apply>
<times/>
<ci>k_p_v</ci>
<ci>M_A</ci>
<ci>n_1</ci>
<ci>L_oz</ci>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<ci>N</ci>
</apply>
</apply>
<apply>
<times/>
<ci>k_m_v</ci>
<ci>N</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>N</ci>
</apply>
<ci>K_kappa</ci>
</apply>
</math>
</component>
<component name="length">
<variable name="dl_1_dt" units="micrometre_per_second" public_interface="out"/>
<variable name="l_1" units="micrometre" initial_value="0.436333342969918" public_interface="out"/>
<variable name="dl_2_dt" units="micrometre_per_second" public_interface="out"/>
<variable name="l_2" units="micrometre" initial_value="0.436333525334166" public_interface="out"/>
<variable name="dl_3_dt" units="micrometre_per_second"/>
<variable name="l_3" units="micrometre" initial_value="0.088805830771694" public_interface="out"/>
<variable name="time" units="second" public_interface="in"/>
<variable name="isotonic_mode" units="dimensionless" public_interface="in"/>
<variable name="v" units="micrometre_per_second" public_interface="in"/>
<variable name="w" units="micrometre_per_second" public_interface="in"/>
<variable name="l" units="micrometre" public_interface="out"/>
<variable name="k_S_vis" units="millinewton_second_per_micrometre" public_interface="in"/>
<variable name="phi_chi_2" units="micrometre_per_second" public_interface="in"/>
<variable name="stim_start" units="second" public_interface="in"/>
<variable name="stim_end" units="second" public_interface="in"/>
<variable name="stim_period" units="second" public_interface="in"/>
<variable name="stim_duration" units="second" public_interface="in"/>
<variable name="stim_amplitude" units="nanoA" public_interface="in"/>
<variable name="l_0" units="micrometre" initial_value="0.525139356105856" public_interface="out"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>l</ci>
<apply>
<plus/>
<ci>l_2</ci>
<ci>l_3</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>dl_1_dt</ci>
<ci>v</ci>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>l_1</ci>
</apply>
<ci>dl_1_dt</ci>
</apply>
<apply>
<eq/>
<ci>dl_2_dt</ci>
<piecewise>
<piece>
<ci>phi_chi_2</ci>
<apply>
<eq/>
<ci>k_S_vis</ci>
<cn cellml:units="millinewton_second_per_micrometre">0</cn>
</apply>
</piece>
<otherwise>
<ci>w</ci>
</otherwise>
</piecewise>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>l_2</ci>
</apply>
<ci>dl_2_dt</ci>
</apply>
<apply>
<eq/>
<ci>dl_3_dt</ci>
<piecewise>
<piece>
<cn cellml:units="micrometre_per_second">0</cn>
<apply>
<eq/>
<ci>isotonic_mode</ci>
<cn cellml:units="dimensionless">1</cn>
</apply>
</piece>
<piece>
<apply>
<minus/>
<ci>phi_chi_2</ci>
</apply>
<apply>
<and/>
<apply>
<eq/>
<ci>isotonic_mode</ci>
<cn cellml:units="dimensionless">0</cn>
</apply>
<apply>
<eq/>
<ci>k_S_vis</ci>
<cn cellml:units="millinewton_second_per_micrometre">0</cn>
</apply>
</apply>
</piece>
<piece>
<apply>
<minus/>
<ci>w</ci>
</apply>
<apply>
<and/>
<apply>
<eq/>
<ci>isotonic_mode</ci>
<cn cellml:units="dimensionless">0</cn>
</apply>
<apply>
<neq/>
<ci>k_S_vis</ci>
<cn cellml:units="millinewton_second_per_micrometre">0</cn>
</apply>
</apply>
</piece>
</piecewise>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>l_3</ci>
</apply>
<piecewise>
<piece>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<minus/>
<apply>
<plus/>
<ci>l_2</ci>
<ci>l_3</ci>
</apply>
<ci>l_0</ci>
</apply>
</apply>
<ci>stim_duration</ci>
</apply>
<apply>
<and/>
<apply>
<geq/>
<ci>time</ci>
<ci>stim_start</ci>
</apply>
<apply>
<leq/>
<ci>time</ci>
<ci>stim_end</ci>
</apply>
<apply>
<leq/>
<apply>
<minus/>
<apply>
<minus/>
<ci>time</ci>
<ci>stim_start</ci>
</apply>
<apply>
<times/>
<apply>
<floor/>
<apply>
<divide/>
<apply>
<minus/>
<ci>time</ci>
<ci>stim_start</ci>
</apply>
<ci>stim_period</ci>
</apply>
</apply>
<ci>stim_period</ci>
</apply>
</apply>
<ci>stim_duration</ci>
</apply>
</apply>
</piece>
<otherwise>
<ci>dl_3_dt</ci>
</otherwise>
</piecewise>
</apply>
</math>
</component>
<component name="CE_velocity">
<variable name="v" units="micrometre_per_second" initial_value="0" public_interface="out"/>
<variable name="phi_chi" units="micrometre_per_second2" public_interface="out"/>
<variable name="phi_chi_2" units="micrometre_per_second" public_interface="out"/>
<variable name="k_P_vis" units="millinewton_second_per_micrometre" public_interface="out"/>
<variable name="k_S_vis" units="millinewton_second_per_micrometre" public_interface="in"/>
<variable name="time" units="second" public_interface="in"/>
<variable name="w" units="micrometre_per_second" public_interface="in"/>
<variable name="K_kappa" units="per_second" public_interface="in"/>
<variable name="N" units="dimensionless" public_interface="in"/>
<variable name="l_1" units="micrometre" public_interface="in"/>
<variable name="dl_1_dt" units="micrometre_per_second" public_interface="in"/>
<variable name="l_2" units="micrometre" public_interface="in"/>
<variable name="l_3" units="micrometre" public_interface="in"/>
<variable name="alpha_1" units="per_micrometre" public_interface="in"/>
<variable name="beta_1" units="millinewton" public_interface="in"/>
<variable name="alpha_2" units="per_micrometre" public_interface="in"/>
<variable name="beta_2" units="millinewton" public_interface="in"/>
<variable name="alpha_3" units="per_micrometre" public_interface="in"/>
<variable name="beta_3" units="millinewton" public_interface="in"/>
<variable name="lambda" units="millinewton" public_interface="in"/>
<variable name="isotonic_mode" units="dimensionless" public_interface="in"/>
<variable name="p_v" units="dimensionless" public_interface="in"/>
<variable name="p_prime_v" units="second_per_micrometre" public_interface="in"/>
<variable name="alpha_P_lengthening" units="per_micrometre" initial_value="16"/>
<variable name="beta_P_lengthening" units="millinewton_second_per_micrometre" initial_value="0.0015"/>
<variable name="alpha_P_shortening" units="per_micrometre" initial_value="16"/>
<variable name="beta_P_shortening" units="millinewton_second_per_micrometre" initial_value="0.0015"/>
<variable name="alp_p" units="per_micrometre"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>alp_p</ci>
<piecewise>
<piece>
<ci>alpha_P_lengthening</ci>
<apply>
<leq/>
<ci>v</ci>
<cn cellml:units="micrometre_per_second">0</cn>
</apply>
</piece>
<otherwise>
<ci>alpha_P_shortening</ci>
</otherwise>
</piecewise>
</apply>
<apply>
<eq/>
<ci>k_P_vis</ci>
<piecewise>
<piece>
<apply>
<times/>
<ci>beta_P_lengthening</ci>
<apply>
<exp/>
<apply>
<times/>
<ci>alpha_P_lengthening</ci>
<ci>l_1</ci>
</apply>
</apply>
</apply>
<apply>
<leq/>
<ci>v</ci>
<cn cellml:units="micrometre_per_second">0</cn>
</apply>
</piece>
<otherwise>
<apply>
<times/>
<ci>beta_P_shortening</ci>
<apply>
<exp/>
<apply>
<times/>
<ci>alpha_P_shortening</ci>
<ci>l_1</ci>
</apply>
</apply>
</apply>
</otherwise>
</piecewise>
</apply>
<apply>
<eq/>
<ci>phi_chi</ci>
<piecewise>
<piece>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<plus/>
<apply>
<times/>
<ci>lambda</ci>
<ci>K_kappa</ci>
<ci>p_v</ci>
</apply>
<apply>
<times/>
<ci>alp_p</ci>
<ci>k_P_vis</ci>
<apply>
<power/>
<ci>v</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
<apply>
<times/>
<ci>alpha_2</ci>
<ci>beta_2</ci>
<apply>
<exp/>
<apply>
<times/>
<ci>alpha_2</ci>
<ci>l_2</ci>
</apply>
</apply>
<ci>w</ci>
</apply>
</apply>
</apply>
<apply>
<plus/>
<apply>
<times/>
<ci>lambda</ci>
<ci>N</ci>
<ci>p_prime_v</ci>
</apply>
<ci>k_P_vis</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>isotonic_mode</ci>
<cn cellml:units="dimensionless">1</cn>
</apply>
</piece>
<otherwise>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<plus/>
<apply>
<times/>
<ci>lambda</ci>
<ci>K_kappa</ci>
<ci>p_v</ci>
</apply>
<apply>
<times/>
<ci>alp_p</ci>
<ci>k_P_vis</ci>
<apply>
<power/>
<ci>v</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
<apply>
<times/>
<apply>
<plus/>
<apply>
<times/>
<ci>alpha_2</ci>
<ci>beta_2</ci>
<apply>
<exp/>
<apply>
<times/>
<ci>alpha_2</ci>
<ci>l_2</ci>
</apply>
</apply>
</apply>
<apply>
<times/>
<ci>alpha_3</ci>
<ci>beta_3</ci>
<apply>
<exp/>
<apply>
<times/>
<ci>alpha_3</ci>
<ci>l_3</ci>
</apply>
</apply>
</apply>
</apply>
<ci>w</ci>
</apply>
</apply>
</apply>
<apply>
<plus/>
<apply>
<times/>
<ci>lambda</ci>
<ci>N</ci>
<ci>p_prime_v</ci>
</apply>
<ci>k_P_vis</ci>
</apply>
</apply>
</otherwise>
</piecewise>
</apply>
<apply>
<eq/>
<ci>phi_chi_2</ci>
<piecewise>
<piece>
<apply>
<divide/>
<apply>
<times/>
<ci>alpha_1</ci>
<ci>beta_1</ci>
<apply>
<exp/>
<apply>
<times/>
<ci>alpha_1</ci>
<apply>
<minus/>
<ci>l_2</ci>
<ci>l_1</ci>
</apply>
</apply>
</apply>
<ci>v</ci>
</apply>
<apply>
<plus/>
<apply>
<times/>
<ci>alpha_1</ci>
<ci>beta_1</ci>
<apply>
<exp/>
<apply>
<times/>
<ci>alpha_1</ci>
<apply>
<minus/>
<ci>l_2</ci>
<ci>l_1</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<times/>
<ci>alpha_2</ci>
<ci>beta_2</ci>
<apply>
<exp/>
<apply>
<times/>
<ci>alpha_2</ci>
<ci>l_2</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>isotonic_mode</ci>
<cn cellml:units="dimensionless">1</cn>
</apply>
</piece>
<otherwise>
<apply>
<divide/>
<apply>
<times/>
<ci>alpha_1</ci>
<ci>beta_1</ci>
<apply>
<exp/>
<apply>
<times/>
<ci>alpha_1</ci>
<apply>
<minus/>
<ci>l_2</ci>
<ci>l_1</ci>
</apply>
</apply>
</apply>
<ci>v</ci>
</apply>
<apply>
<plus/>
<apply>
<times/>
<ci>alpha_1</ci>
<ci>beta_1</ci>
<apply>
<exp/>
<apply>
<times/>
<ci>alpha_1</ci>
<apply>
<minus/>
<ci>l_2</ci>
<ci>l_1</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<times/>
<ci>alpha_2</ci>
<ci>beta_2</ci>
<apply>
<exp/>
<apply>
<times/>
<ci>alpha_2</ci>
<ci>l_2</ci>
</apply>
</apply>
</apply>
<apply>
<times/>
<ci>alpha_3</ci>
<ci>beta_3</ci>
<apply>
<exp/>
<apply>
<times/>
<ci>alpha_3</ci>
<ci>l_3</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</otherwise>
</piecewise>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>v</ci>
</apply>
<piecewise>
<piece>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<times/>
<ci>alpha_1</ci>
<ci>beta_1</ci>
<apply>
<exp/>
<apply>
<times/>
<ci>alpha_1</ci>
<apply>
<minus/>
<ci>l_2</ci>
<ci>l_1</ci>
</apply>
</apply>
</apply>
<apply>
<minus/>
<ci>phi_chi_2</ci>
<ci>v</ci>
</apply>
</apply>
<apply>
<plus/>
<apply>
<times/>
<ci>lambda</ci>
<ci>K_kappa</ci>
<ci>p_v</ci>
</apply>
<apply>
<times/>
<ci>alp_p</ci>
<ci>k_P_vis</ci>
<apply>
<power/>
<ci>v</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
</apply>
</apply>
<apply>
<plus/>
<apply>
<times/>
<ci>lambda</ci>
<ci>N</ci>
<ci>p_prime_v</ci>
</apply>
<ci>k_P_vis</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>k_S_vis</ci>
<cn cellml:units="millinewton_second_per_micrometre">0</cn>
</apply>
</piece>
<otherwise>
<ci>phi_chi</ci>
</otherwise>
</piecewise>
</apply>
</math>
</component>
<component name="PE_velocity">
<variable name="w" units="micrometre_per_second" initial_value="0" public_interface="out"/>
<variable name="k_S_vis" units="millinewton_second_per_micrometre" public_interface="out"/>
<variable name="time" units="second" public_interface="in"/>
<variable name="v" units="micrometre_per_second" public_interface="in"/>
<variable name="phi_chi" units="micrometre_per_second2" public_interface="in"/>
<variable name="l_1" units="micrometre" public_interface="in"/>
<variable name="dl_1_dt" units="micrometre_per_second" public_interface="in"/>
<variable name="l_2" units="micrometre" public_interface="in"/>
<variable name="dl_2_dt" units="micrometre_per_second" public_interface="in"/>
<variable name="l_3" units="micrometre" public_interface="in"/>
<variable name="alpha_1" units="per_micrometre" public_interface="in"/>
<variable name="beta_1" units="millinewton" public_interface="in"/>
<variable name="alpha_2" units="per_micrometre" public_interface="in"/>
<variable name="beta_2" units="millinewton" public_interface="in"/>
<variable name="alpha_3" units="per_micrometre" public_interface="in"/>
<variable name="beta_3" units="millinewton" public_interface="in"/>
<variable name="lambda" units="millinewton" public_interface="in"/>
<variable name="isotonic_mode" units="dimensionless" public_interface="in"/>
<variable name="alpha_S_lengthening" units="per_micrometre" initial_value="46"/>
<variable name="beta_S_lengthening" units="millinewton_second_per_micrometre" initial_value="0"/>
<variable name="alpha_S_shortening" units="per_micrometre" initial_value="39"/>
<variable name="beta_S_shortening" units="millinewton_second_per_micrometre" initial_value="0"/>
<variable name="alp_s" units="per_micrometre"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>alp_s</ci>
<piecewise>
<piece>
<ci>alpha_S_lengthening</ci>
<apply>
<leq/>
<ci>w</ci>
<ci>v</ci>
</apply>
</piece>
<otherwise>
<ci>alpha_S_shortening</ci>
</otherwise>
</piecewise>
</apply>
<apply>
<eq/>
<ci>k_S_vis</ci>
<piecewise>
<piece>
<apply>
<times/>
<ci>beta_S_lengthening</ci>
<apply>
<exp/>
<apply>
<times/>
<ci>alpha_S_lengthening</ci>
<apply>
<minus/>
<ci>l_2</ci>
<ci>l_1</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<leq/>
<ci>w</ci>
<ci>v</ci>
</apply>
</piece>
<otherwise>
<apply>
<times/>
<ci>beta_S_shortening</ci>
<apply>
<exp/>
<apply>
<times/>
<ci>alpha_S_shortening</ci>
<apply>
<minus/>
<ci>l_2</ci>
<ci>l_1</ci>
</apply>
</apply>
</apply>
</apply>
</otherwise>
</piecewise>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>w</ci>
</apply>
<piecewise>
<piece>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<minus/>
<apply>
<times/>
<ci>k_S_vis</ci>
<apply>
<minus/>
<ci>phi_chi</ci>
<apply>
<times/>
<ci>alp_s</ci>
<apply>
<power/>
<apply>
<minus/>
<ci>w</ci>
<ci>v</ci>
</apply>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
</apply>
</apply>
<apply>
<times/>
<ci>alpha_1</ci>
<ci>beta_1</ci>
<apply>
<exp/>
<apply>
<times/>
<ci>alpha_1</ci>
<apply>
<minus/>
<ci>l_2</ci>
<ci>l_1</ci>
</apply>
</apply>
</apply>
<apply>
<minus/>
<ci>w</ci>
<ci>v</ci>
</apply>
</apply>
</apply>
<apply>
<times/>
<ci>alpha_2</ci>
<ci>beta_2</ci>
<apply>
<exp/>
<apply>
<times/>
<ci>alpha_2</ci>
<ci>l_2</ci>
</apply>
</apply>
<ci>w</ci>
</apply>
</apply>
<ci>k_S_vis</ci>
</apply>
<apply>
<and/>
<apply>
<eq/>
<ci>isotonic_mode</ci>
<cn cellml:units="dimensionless">1</cn>
</apply>
<apply>
<neq/>
<ci>k_S_vis</ci>
<cn cellml:units="millinewton_second_per_micrometre">0</cn>
</apply>
</apply>
</piece>
<piece>
<apply>
<minus/>
<apply>
<minus/>
<ci>phi_chi</ci>
<apply>
<times/>
<ci>alp_s</ci>
<apply>
<power/>
<apply>
<minus/>
<ci>w</ci>
<ci>v</ci>
</apply>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
</apply>
<apply>
<divide/>
<apply>
<plus/>
<apply>
<times/>
<ci>alpha_1</ci>
<ci>beta_1</ci>
<apply>
<exp/>
<apply>
<times/>
<ci>alpha_1</ci>
<apply>
<minus/>
<ci>l_2</ci>
<ci>l_1</ci>
</apply>
</apply>
</apply>
<apply>
<minus/>
<ci>w</ci>
<ci>v</ci>
</apply>
</apply>
<apply>
<times/>
<apply>
<plus/>
<apply>
<times/>
<ci>alpha_2</ci>
<ci>beta_2</ci>
<apply>
<exp/>
<apply>
<times/>
<ci>alpha_2</ci>
<ci>l_2</ci>
</apply>
</apply>
</apply>
<apply>
<times/>
<ci>alpha_3</ci>
<ci>beta_3</ci>
<apply>
<exp/>
<apply>
<times/>
<ci>alpha_3</ci>
<ci>l_3</ci>
</apply>
</apply>
</apply>
</apply>
<ci>w</ci>
</apply>
</apply>
<ci>k_S_vis</ci>
</apply>
</apply>
<apply>
<and/>
<apply>
<eq/>
<ci>isotonic_mode</ci>
<cn cellml:units="dimensionless">0</cn>
</apply>
<apply>
<neq/>
<ci>k_S_vis</ci>
<cn cellml:units="millinewton_second_per_micrometre">0</cn>
</apply>
</apply>
</piece>
<piece>
<cn cellml:units="micrometre_per_second2">0</cn>
<apply>
<eq/>
<ci>k_S_vis</ci>
<cn cellml:units="millinewton_second_per_micrometre">0</cn>
</apply>
</piece>
</piecewise>
</apply>
</math>
</component>
<component name="average_crossbridge_force">
<variable name="p_v" units="dimensionless" public_interface="out"/>
<variable name="p_prime_v" units="second_per_micrometre" public_interface="out"/>
<variable name="P_star" units="dimensionless" public_interface="out"/>
<variable name="G_star" units="dimensionless" public_interface="out"/>
<variable name="gamma" units="dimensionless"/>
<variable name="case_1" units="second_per_micrometre"/>
<variable name="case_2" units="second_per_micrometre"/>
<variable name="case_3" units="second_per_micrometre"/>
<variable name="case_4" units="second_per_micrometre"/>
<variable name="v" units="micrometre_per_second" public_interface="in"/>
<variable name="phi_chi" units="micrometre_per_second2" public_interface="in"/>
<variable name="v_1" units="micrometre_per_second" public_interface="in"/>
<variable name="v_max" units="micrometre_per_second" public_interface="in"/>
<variable name="a" units="dimensionless" public_interface="in"/>
<variable name="d_h" units="dimensionless" public_interface="in"/>
<variable name="alpha_G" units="dimensionless" public_interface="in"/>
<variable name="alpha_P" units="dimensionless" public_interface="in"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>gamma</ci>
<apply>
<divide/>
<apply>
<times/>
<ci>a</ci>
<ci>d_h</ci>
<apply>
<power/>
<apply>
<divide/>
<ci>v_1</ci>
<ci>v_max</ci>
</apply>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
<apply>
<minus/>
<apply>
<times/>
<cn cellml:units="dimensionless">3</cn>
<ci>a</ci>
<ci>d_h</ci>
</apply>
<apply>
<divide/>
<apply>
<times/>
<apply>
<plus/>
<ci>a</ci>
<cn cellml:units="dimensionless">1</cn>
</apply>
<ci>v_1</ci>
</apply>
<ci>v_max</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>P_star</ci>
<piecewise>
<piece>
<apply>
<divide/>
<apply>
<times/>
<ci>a</ci>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<divide/>
<ci>v</ci>
<ci>v_max</ci>
</apply>
</apply>
</apply>
<apply>
<minus/>
<ci>a</ci>
<apply>
<divide/>
<ci>v</ci>
<ci>v_max</ci>
</apply>
</apply>
</apply>
<apply>
<leq/>
<ci>v</ci>
<cn cellml:units="micrometre_per_second">0</cn>
</apply>
</piece>
<otherwise>
<apply>
<minus/>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<ci>d_h</ci>
</apply>
<apply>
<divide/>
<apply>
<times/>
<apply>
<power/>
<ci>d_h</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
<ci>a</ci>
</apply>
<apply>
<plus/>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<ci>a</ci>
<ci>d_h</ci>
</apply>
<ci>gamma</ci>
</apply>
<apply>
<power/>
<apply>
<divide/>
<ci>v</ci>
<ci>v_max</ci>
</apply>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
<apply>
<divide/>
<apply>
<times/>
<apply>
<plus/>
<ci>a</ci>
<cn cellml:units="dimensionless">1</cn>
</apply>
<ci>v</ci>
</apply>
<ci>v_max</ci>
</apply>
<apply>
<times/>
<ci>a</ci>
<ci>d_h</ci>
</apply>
</apply>
</apply>
</apply>
</otherwise>
</piecewise>
</apply>
<apply>
<eq/>
<ci>G_star</ci>
<piecewise>
<piece>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<divide/>
<apply>
<times/>
<cn cellml:units="dimensionless">0.6</cn>
<ci>v</ci>
</apply>
<ci>v_max</ci>
</apply>
</apply>
<apply>
<and/>
<apply>
<leq/>
<apply>
<minus/>
<ci>v_max</ci>
</apply>
<ci>v</ci>
</apply>
<apply>
<leq/>
<ci>v</ci>
<cn cellml:units="micrometre_per_second">0</cn>
</apply>
</apply>
</piece>
<piece>
<apply>
<divide/>
<ci>P_star</ci>
<apply>
<plus/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="dimensionless">0.4</cn>
<ci>a</ci>
</apply>
<cn cellml:units="dimensionless">1</cn>
</apply>
<ci>v</ci>
</apply>
<apply>
<times/>
<ci>a</ci>
<ci>v_max</ci>
</apply>
</apply>
<cn cellml:units="dimensionless">1</cn>
</apply>
</apply>
<apply>
<and/>
<apply>
<lt/>
<cn cellml:units="micrometre_per_second">0</cn>
<ci>v</ci>
</apply>
<apply>
<leq/>
<ci>v</ci>
<ci>v_1</ci>
</apply>
</apply>
</piece>
<otherwise>
<apply>
<divide/>
<apply>
<times/>
<ci>P_star</ci>
<apply>
<exp/>
<apply>
<times/>
<apply>
<minus/>
<ci>alpha_G</ci>
</apply>
<apply>
<power/>
<apply>
<divide/>
<apply>
<minus/>
<ci>v</ci>
<ci>v_1</ci>
</apply>
<ci>v_max</ci>
</apply>
<ci>alpha_P</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<plus/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="dimensionless">0.4</cn>
<ci>a</ci>
</apply>
<cn cellml:units="dimensionless">1</cn>
</apply>
<ci>v</ci>
</apply>
<apply>
<times/>
<ci>a</ci>
<ci>v_max</ci>
</apply>
</apply>
<cn cellml:units="dimensionless">1</cn>
</apply>
</apply>
</otherwise>
</piecewise>
</apply>
<apply>
<eq/>
<ci>case_1</ci>
<apply>
<divide/>
<apply>
<times/>
<ci>a</ci>
<apply>
<plus/>
<cn cellml:units="dimensionless">0.4</cn>
<apply>
<times/>
<cn cellml:units="dimensionless">0.4</cn>
<ci>a</ci>
</apply>
</apply>
</apply>
<apply>
<times/>
<ci>v_max</ci>
<apply>
<power/>
<apply>
<times/>
<apply>
<plus/>
<ci>a</ci>
<cn cellml:units="dimensionless">1</cn>
</apply>
<cn cellml:units="dimensionless">0.4</cn>
</apply>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>case_2</ci>
<apply>
<divide/>
<apply>
<times/>
<ci>a</ci>
<cn cellml:units="dimensionless">1</cn>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<times/>
<cn cellml:units="dimensionless">0.4</cn>
<ci>a</ci>
</apply>
<apply>
<divide/>
<apply>
<times/>
<cn cellml:units="dimensionless">1.2</cn>
<ci>v</ci>
</apply>
<ci>v_max</ci>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">0.6</cn>
<apply>
<power/>
<apply>
<divide/>
<ci>v</ci>
<ci>v_max</ci>
</apply>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
</apply>
</apply>
<apply>
<times/>
<ci>v_max</ci>
<apply>
<power/>
<apply>
<times/>
<apply>
<minus/>
<ci>a</ci>
<apply>
<divide/>
<ci>v</ci>
<ci>v_max</ci>
</apply>
</apply>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<divide/>
<apply>
<times/>
<cn cellml:units="dimensionless">0.6</cn>
<ci>v</ci>
</apply>
<ci>v_max</ci>
</apply>
</apply>
</apply>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>case_3</ci>
<apply>
<divide/>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="dimensionless">0.4</cn>
<ci>a</ci>
</apply>
<cn cellml:units="dimensionless">1</cn>
</apply>
<apply>
<times/>
<ci>a</ci>
<ci>v_max</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>case_4</ci>
<apply>
<times/>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<ci>v_max</ci>
</apply>
<apply>
<exp/>
<apply>
<times/>
<apply>
<minus/>
<ci>alpha_G</ci>
</apply>
<apply>
<power/>
<apply>
<minus/>
<apply>
<divide/>
<ci>v</ci>
<ci>v_max</ci>
</apply>
<apply>
<divide/>
<ci>v_1</ci>
<ci>v_max</ci>
</apply>
</apply>
<ci>alpha_P</ci>
</apply>
</apply>
</apply>
<apply>
<plus/>
<apply>
<divide/>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="dimensionless">0.4</cn>
<ci>a</ci>
</apply>
<cn cellml:units="dimensionless">1</cn>
</apply>
<ci>a</ci>
</apply>
<apply>
<times/>
<ci>alpha_G</ci>
<ci>alpha_P</ci>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<divide/>
<apply>
<times/>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="dimensionless">0.4</cn>
<ci>a</ci>
</apply>
<cn cellml:units="dimensionless">1</cn>
</apply>
<ci>v</ci>
</apply>
<apply>
<times/>
<ci>a</ci>
<ci>v_max</ci>
</apply>
</apply>
</apply>
<apply>
<power/>
<apply>
<minus/>
<apply>
<divide/>
<ci>v</ci>
<ci>v_max</ci>
</apply>
<apply>
<divide/>
<ci>v_1</ci>
<ci>v_max</ci>
</apply>
</apply>
<apply>
<minus/>
<ci>alpha_P</ci>
<cn cellml:units="dimensionless">1</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>p_v</ci>
<apply>
<divide/>
<ci>P_star</ci>
<ci>G_star</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>p_prime_v</ci>
<piecewise>
<piece>
<ci>case_1</ci>
<apply>
<leq/>
<ci>v</ci>
<apply>
<minus/>
<ci>v_max</ci>
</apply>
</apply>
</piece>
<piece>
<ci>case_2</ci>
<apply>
<and/>
<apply>
<lt/>
<apply>
<minus/>
<ci>v_max</ci>
</apply>
<ci>v</ci>
</apply>
<apply>
<leq/>
<ci>v</ci>
<cn cellml:units="micrometre_per_second">0</cn>
</apply>
</apply>
</piece>
<piece>
<ci>case_3</ci>
<apply>
<and/>
<apply>
<lt/>
<cn cellml:units="micrometre_per_second">0</cn>
<ci>v</ci>
</apply>
<apply>
<leq/>
<ci>v</ci>
<ci>v_1</ci>
</apply>
</apply>
</piece>
<otherwise>
<ci>case_4</ci>
</otherwise>
</piecewise>
</apply>
</math>
</component>
<component name="membrane">
<variable name="V" units="millivolt" initial_value="-93.658148" public_interface="out"/>
<variable name="R" units="joule_per_kilomole_kelvin" initial_value="8314.472" public_interface="out"/>
<variable name="T" units="kelvin" initial_value="310" public_interface="out"/>
<variable name="F" units="coulomb_per_mole" initial_value="96485.3415" public_interface="out"/>
<variable name="Cm" units="microF" initial_value="9.5e-5"/>
<variable name="time" units="second" public_interface="in"/>
<variable name="i_K1" units="nanoA" public_interface="in"/>
<variable name="i_to" units="nanoA" public_interface="in"/>
<variable name="i_Kr" units="nanoA" public_interface="in"/>
<variable name="i_Ks" units="nanoA" public_interface="in"/>
<variable name="i_Ca_L_K_cyt" units="nanoA" public_interface="in"/>
<variable name="i_Ca_L_K_ds" units="nanoA" public_interface="in"/>
<variable name="i_NaK" units="nanoA" public_interface="in"/>
<variable name="i_Na" units="nanoA" public_interface="in"/>
<variable name="i_b_Na" units="nanoA" public_interface="in"/>
<variable name="i_p_Na" units="nanoA" public_interface="in"/>
<variable name="i_Ca_L_Na_cyt" units="nanoA" public_interface="in"/>
<variable name="i_Ca_L_Na_ds" units="nanoA" public_interface="in"/>
<variable name="i_NaCa_cyt" units="nanoA" public_interface="in"/>
<variable name="i_NaCa_ds" units="nanoA" public_interface="in"/>
<variable name="i_Ca_L_Ca_cyt" units="nanoA" public_interface="in"/>
<variable name="i_Ca_L_Ca_ds" units="nanoA" public_interface="in"/>
<variable name="i_b_Ca" units="nanoA" public_interface="in"/>
<variable name="i_Stim" units="nanoA"/>
<variable name="stim_start" units="second" initial_value="0.06" public_interface="out"/>
<variable name="stim_end" units="second" initial_value="10000" public_interface="out"/>
<variable name="stim_period" units="second" initial_value="1" public_interface="out"/>
<variable name="stim_duration" units="second" initial_value="0.0025" public_interface="out"/>
<variable name="stim_amplitude" units="nanoA" initial_value="-3" public_interface="out"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>i_Stim</ci>
<piecewise>
<piece>
<ci>stim_amplitude</ci>
<apply>
<and/>
<apply>
<geq/>
<ci>time</ci>
<ci>stim_start</ci>
</apply>
<apply>
<leq/>
<ci>time</ci>
<ci>stim_end</ci>
</apply>
<apply>
<leq/>
<apply>
<minus/>
<apply>
<minus/>
<ci>time</ci>
<ci>stim_start</ci>
</apply>
<apply>
<times/>
<apply>
<floor/>
<apply>
<divide/>
<apply>
<minus/>
<ci>time</ci>
<ci>stim_start</ci>
</apply>
<ci>stim_period</ci>
</apply>
</apply>
<ci>stim_period</ci>
</apply>
</apply>
<ci>stim_duration</ci>
</apply>
</apply>
</piece>
<otherwise>
<cn cellml:units="nanoA">0</cn>
</otherwise>
</piecewise>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>V</ci>
</apply>
<apply>
<times/>
<apply>
<divide/>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
</apply>
<ci>Cm</ci>
</apply>
<apply>
<plus/>
<ci>i_Stim</ci>
<ci>i_K1</ci>
<ci>i_to</ci>
<ci>i_Kr</ci>
<ci>i_Ks</ci>
<ci>i_NaK</ci>
<ci>i_Na</ci>
<ci>i_b_Na</ci>
<ci>i_p_Na</ci>
<ci>i_Ca_L_Na_cyt</ci>
<ci>i_Ca_L_Na_ds</ci>
<ci>i_NaCa_cyt</ci>
<ci>i_NaCa_ds</ci>
<ci>i_Ca_L_Ca_cyt</ci>
<ci>i_Ca_L_Ca_ds</ci>
<ci>i_Ca_L_K_cyt</ci>
<ci>i_Ca_L_K_ds</ci>
<ci>i_b_Ca</ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="reversal_potentials">
<variable name="E_Na" units="millivolt" public_interface="out"/>
<variable name="E_K" units="millivolt" public_interface="out"/>
<variable name="E_Ks" units="millivolt" public_interface="out"/>
<variable name="E_Ca" units="millivolt" public_interface="out"/>
<variable name="E_mh" units="millivolt" public_interface="out"/>
<variable name="P_kna" units="dimensionless" initial_value="0.03"/>
<variable name="K_o" units="millimolar" public_interface="in"/>
<variable name="Na_o" units="millimolar" public_interface="in"/>
<variable name="K_i" units="millimolar" public_interface="in"/>
<variable name="Na_i" units="millimolar" public_interface="in"/>
<variable name="Ca_o" units="millimolar" public_interface="in"/>
<variable name="Ca_i" units="millimolar" public_interface="in"/>
<variable name="R" units="joule_per_kilomole_kelvin" public_interface="in"/>
<variable name="F" units="coulomb_per_mole" public_interface="in"/>
<variable name="T" units="kelvin" public_interface="in"/>
<variable name="time" units="second" public_interface="in"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>E_Na</ci>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
<ci>F</ci>
</apply>
<apply>
<ln/>
<apply>
<divide/>
<ci>Na_o</ci>
<ci>Na_i</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>E_K</ci>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
<ci>F</ci>
</apply>
<apply>
<ln/>
<apply>
<divide/>
<ci>K_o</ci>
<ci>K_i</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>E_Ks</ci>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
<ci>F</ci>
</apply>
<apply>
<ln/>
<apply>
<divide/>
<apply>
<plus/>
<ci>K_o</ci>
<apply>
<times/>
<ci>P_kna</ci>
<ci>Na_o</ci>
</apply>
</apply>
<apply>
<plus/>
<ci>K_i</ci>
<apply>
<times/>
<ci>P_kna</ci>
<ci>Na_i</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>E_Ca</ci>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<cn cellml:units="dimensionless">0.5</cn>
<ci>R</ci>
<ci>T</ci>
</apply>
<ci>F</ci>
</apply>
<apply>
<ln/>
<apply>
<divide/>
<ci>Ca_o</ci>
<ci>Ca_i</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>E_mh</ci>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
<ci>F</ci>
</apply>
<apply>
<ln/>
<apply>
<divide/>
<apply>
<plus/>
<ci>Na_o</ci>
<apply>
<times/>
<cn cellml:units="dimensionless">0.12</cn>
<ci>K_o</ci>
</apply>
</apply>
<apply>
<plus/>
<ci>Na_i</ci>
<apply>
<times/>
<cn cellml:units="dimensionless">0.12</cn>
<ci>K_i</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<component name="time_independent_potassium_current">
<variable name="i_K1" units="nanoA" public_interface="out"/>
<variable name="K_mk1" units="millimolar" initial_value="10" public_interface="out"/>
<variable name="g_K1" units="microS" initial_value="0.5"/>
<variable name="K_o" units="millimolar" public_interface="in"/>
<variable name="V" units="millivolt" public_interface="in"/>
<variable name="E_K" units="millivolt" public_interface="in"/>
<variable name="R" units="joule_per_kilomole_kelvin" public_interface="in"/>
<variable name="F" units="coulomb_per_mole" public_interface="in"/>
<variable name="T" units="kelvin" public_interface="in"/>
<variable name="time" units="second" public_interface="in"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>i_K1</ci>
<apply>
<divide/>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<ci>g_K1</ci>
<ci>K_o</ci>
</apply>
<apply>
<plus/>
<ci>K_o</ci>
<ci>K_mk1</ci>
</apply>
</apply>
<apply>
<minus/>
<ci>V</ci>
<ci>E_K</ci>
</apply>
</apply>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<apply>
<minus/>
<ci>V</ci>
<ci>E_K</ci>
</apply>
<cn cellml:units="millivolt">10</cn>
</apply>
<ci>F</ci>
<cn cellml:units="dimensionless">1.25</cn>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<component name="rapid_delayed_rectifier_potassium_current">
<variable name="i_Kr" units="nanoA" public_interface="out"/>
<variable name="g_Kr1" units="microS" initial_value="0.0021"/>
<variable name="g_Kr2" units="microS" initial_value="0.0013"/>
<variable name="E_K" units="millivolt" public_interface="in"/>
<variable name="time" units="second" public_interface="in" private_interface="out"/>
<variable name="V" units="millivolt" public_interface="in" private_interface="out"/>
<variable name="xr1" units="dimensionless" private_interface="in"/>
<variable name="xr2" units="dimensionless" private_interface="in"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>i_Kr</ci>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<plus/>
<apply>
<times/>
<ci>g_Kr1</ci>
<ci>xr1</ci>
</apply>
<apply>
<times/>
<ci>g_Kr2</ci>
<ci>xr2</ci>
</apply>
</apply>
<cn cellml:units="dimensionless">1</cn>
</apply>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="millivolt">9</cn>
</apply>
<cn cellml:units="millivolt">22.4</cn>
</apply>
</apply>
</apply>
</apply>
<apply>
<minus/>
<ci>V</ci>
<ci>E_K</ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="rapid_delayed_rectifier_potassium_current_xr1_gate">
<variable name="xr1" units="dimensionless" initial_value="8.88859784542779e-6" public_interface="out"/>
<variable name="V" units="millivolt" public_interface="in"/>
<variable name="alpha_xr1" units="per_second"/>
<variable name="beta_xr1" units="per_second"/>
<variable name="time" units="second" public_interface="in"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>alpha_xr1</ci>
<apply>
<divide/>
<cn cellml:units="per_second">50</cn>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<minus/>
<ci>V</ci>
<cn cellml:units="millivolt">5</cn>
</apply>
</apply>
<cn cellml:units="millivolt">9</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>beta_xr1</ci>
<apply>
<times/>
<cn cellml:units="per_second">0.05</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<minus/>
<ci>V</ci>
<cn cellml:units="millivolt">20</cn>
</apply>
</apply>
<cn cellml:units="millivolt">15</cn>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>xr1</ci>
</apply>
<apply>
<minus/>
<apply>
<times/>
<ci>alpha_xr1</ci>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<ci>xr1</ci>
</apply>
</apply>
<apply>
<times/>
<ci>beta_xr1</ci>
<ci>xr1</ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="rapid_delayed_rectifier_potassium_current_xr2_gate">
<variable name="xr2" units="dimensionless" initial_value="1.53745791069154e-7" public_interface="out"/>
<variable name="V" units="millivolt" public_interface="in"/>
<variable name="alpha_xr2" units="per_second"/>
<variable name="beta_xr2" units="per_second"/>
<variable name="time" units="second" public_interface="in"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>alpha_xr2</ci>
<apply>
<divide/>
<cn cellml:units="per_second">50</cn>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<minus/>
<ci>V</ci>
<cn cellml:units="millivolt">5</cn>
</apply>
</apply>
<cn cellml:units="millivolt">9</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>beta_xr2</ci>
<apply>
<times/>
<cn cellml:units="per_second">0.4</cn>
<apply>
<exp/>
<apply>
<minus/>
<apply>
<power/>
<apply>
<divide/>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="millivolt">30</cn>
</apply>
<cn cellml:units="millivolt">30</cn>
</apply>
<cn cellml:units="dimensionless">3</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>xr2</ci>
</apply>
<apply>
<minus/>
<apply>
<times/>
<ci>alpha_xr2</ci>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<ci>xr2</ci>
</apply>
</apply>
<apply>
<times/>
<ci>beta_xr2</ci>
<ci>xr2</ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="slow_delayed_rectifier_potassium_current">
<variable name="i_Ks" units="nanoA" public_interface="out"/>
<variable name="g_Ks" units="microS" initial_value="0.0026"/>
<variable name="E_Ks" units="millivolt" public_interface="in"/>
<variable name="time" units="second" public_interface="in" private_interface="out"/>
<variable name="V" units="millivolt" public_interface="in" private_interface="out"/>
<variable name="xs" units="dimensionless" private_interface="in"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>i_Ks</ci>
<apply>
<times/>
<ci>g_Ks</ci>
<apply>
<power/>
<ci>xs</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
<apply>
<minus/>
<ci>V</ci>
<ci>E_Ks</ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="slow_delayed_rectifier_potassium_current_xs_gate">
<variable name="xs" units="dimensionless" initial_value="0.001" public_interface="out"/>
<variable name="V" units="millivolt" public_interface="in"/>
<variable name="alpha_xs" units="per_second"/>
<variable name="beta_xs" units="per_second"/>
<variable name="time" units="second" public_interface="in"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>alpha_xs</ci>
<apply>
<divide/>
<cn cellml:units="per_second">14</cn>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<minus/>
<ci>V</ci>
<cn cellml:units="millivolt">40</cn>
</apply>
</apply>
<cn cellml:units="millivolt">9</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>beta_xs</ci>
<apply>
<times/>
<cn cellml:units="per_second">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<ci>V</ci>
</apply>
<cn cellml:units="millivolt">45</cn>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>xs</ci>
</apply>
<apply>
<minus/>
<apply>
<times/>
<ci>alpha_xs</ci>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<ci>xs</ci>
</apply>
</apply>
<apply>
<times/>
<ci>beta_xs</ci>
<ci>xs</ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="sodium_activated_potassium_current">
<variable name="i_KNa" units="nanoA" public_interface="out"/>
<variable name="g_K_Na" units="microS" initial_value="0"/>
<variable name="K_kna" units="millimolar" initial_value="20"/>
<variable name="Na_i" units="millimolar" public_interface="in"/>
<variable name="E_K" units="millivolt" public_interface="in"/>
<variable name="V" units="millivolt" public_interface="in"/>
<variable name="time" units="second" public_interface="in"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>i_KNa</ci>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<ci>g_K_Na</ci>
<ci>Na_i</ci>
</apply>
<apply>
<plus/>
<ci>Na_i</ci>
<ci>K_kna</ci>
</apply>
</apply>
<apply>
<minus/>
<ci>V</ci>
<ci>E_K</ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="fast_sodium_current">
<variable name="i_Na" units="nanoA" public_interface="out"/>
<variable name="g_Na" units="microS" initial_value="2.5"/>
<variable name="E_mh" units="millivolt" public_interface="in"/>
<variable name="time" units="second" public_interface="in" private_interface="out"/>
<variable name="V" units="millivolt" public_interface="in" private_interface="out"/>
<variable name="m" units="dimensionless" private_interface="in"/>
<variable name="h" units="dimensionless" private_interface="in"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>i_Na</ci>
<apply>
<times/>
<ci>g_Na</ci>
<apply>
<power/>
<ci>m</ci>
<cn cellml:units="dimensionless">3</cn>
</apply>
<ci>h</ci>
<apply>
<minus/>
<ci>V</ci>
<ci>E_mh</ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="fast_sodium_current_m_gate">
<variable name="m" units="dimensionless" initial_value="0.0015" public_interface="out"/>
<variable name="alpha_m" units="per_second"/>
<variable name="beta_m" units="per_second"/>
<variable name="V" units="millivolt" public_interface="in"/>
<variable name="time" units="second" public_interface="in"/>
<variable name="delta_m" units="millivolt" initial_value="1e-5"/>
<variable name="E0_m" units="millivolt"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>E0_m</ci>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="millivolt">41</cn>
</apply>
</apply>
<apply>
<eq/>
<ci>alpha_m</ci>
<piecewise>
<piece>
<cn cellml:units="per_second">2000</cn>
<apply>
<lt/>
<apply>
<abs/>
<ci>E0_m</ci>
</apply>
<ci>delta_m</ci>
</apply>
</piece>
<otherwise>
<apply>
<divide/>
<apply>
<times/>
<cn cellml:units="per_millivolt_second">200</cn>
<ci>E0_m</ci>
</apply>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="per_millivolt">0.1</cn>
</apply>
<ci>E0_m</ci>
</apply>
</apply>
</apply>
</apply>
</otherwise>
</piecewise>
</apply>
<apply>
<eq/>
<ci>beta_m</ci>
<apply>
<times/>
<cn cellml:units="per_second">8000</cn>
<apply>
<exp/>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="per_millivolt">0.056</cn>
</apply>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="millivolt">66</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>m</ci>
</apply>
<apply>
<minus/>
<apply>
<times/>
<ci>alpha_m</ci>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<ci>m</ci>
</apply>
</apply>
<apply>
<times/>
<ci>beta_m</ci>
<ci>m</ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="fast_sodium_current_h_gate">
<variable name="h" units="dimensionless" initial_value="0.995" public_interface="out"/>
<variable name="alpha_h" units="per_second"/>
<variable name="beta_h" units="per_second"/>
<variable name="V" units="millivolt" public_interface="in"/>
<variable name="time" units="second" public_interface="in"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>alpha_h</ci>
<apply>
<times/>
<cn cellml:units="per_second">20</cn>
<apply>
<exp/>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="per_millivolt">0.125</cn>
</apply>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="millivolt">75</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>beta_h</ci>
<apply>
<divide/>
<cn cellml:units="per_second">2000</cn>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<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>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>h</ci>
</apply>
<apply>
<minus/>
<apply>
<times/>
<ci>alpha_h</ci>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<ci>h</ci>
</apply>
</apply>
<apply>
<times/>
<ci>beta_h</ci>
<ci>h</ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="persistent_sodium_current">
<variable name="i_p_Na" units="nanoA" public_interface="out"/>
<variable name="g_pna" units="microS" initial_value="0.004"/>
<variable name="E_Na" units="millivolt" public_interface="in"/>
<variable name="time" units="second" public_interface="in"/>
<variable name="V" units="millivolt" public_interface="in"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>i_p_Na</ci>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<ci>g_pna</ci>
<cn cellml:units="dimensionless">1</cn>
</apply>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="millivolt">52</cn>
</apply>
</apply>
<cn cellml:units="millivolt">8</cn>
</apply>
</apply>
</apply>
</apply>
<apply>
<minus/>
<ci>V</ci>
<ci>E_Na</ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="sodium_background_current">
<variable name="i_b_Na" units="nanoA" public_interface="out"/>
<variable name="g_bna" units="microS" initial_value="0.0006"/>
<variable name="Na_i" units="millimolar" public_interface="in"/>
<variable name="E_Na" units="millivolt" public_interface="in"/>
<variable name="time" units="second" public_interface="in"/>
<variable name="V" units="millivolt" public_interface="in"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>i_b_Na</ci>
<apply>
<times/>
<ci>g_bna</ci>
<apply>
<minus/>
<ci>V</ci>
<ci>E_Na</ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="L_type_Ca_channel">
<variable name="i_Ca_L_Ca_cyt" units="nanoA" public_interface="out"/>
<variable name="i_Ca_L_K_cyt" units="nanoA" public_interface="out"/>
<variable name="i_Ca_L_Na_cyt" units="nanoA" public_interface="out"/>
<variable name="i_Ca_L_Ca_ds" units="nanoA" public_interface="out"/>
<variable name="i_Ca_L_K_ds" units="nanoA" public_interface="out"/>
<variable name="i_Ca_L_Na_ds" units="nanoA" public_interface="out"/>
<variable name="i_Ca_L" units="nanoA" public_interface="out"/>
<variable name="P_Ca_L" units="nanoA_per_millimolar" initial_value="0.1"/>
<variable name="P_CaK" units="dimensionless" initial_value="0.002"/>
<variable name="P_CaNa" units="dimensionless" initial_value="0.01"/>
<variable name="Ca_o" units="millimolar" public_interface="in"/>
<variable name="Ca_i" units="millimolar" public_interface="in" private_interface="out"/>
<variable name="Ca_ds" units="millimolar" public_interface="in" private_interface="out"/>
<variable name="K_o" units="millimolar" public_interface="in"/>
<variable name="K_i" units="millimolar" public_interface="in"/>
<variable name="Na_o" units="millimolar" public_interface="in"/>
<variable name="Na_i" units="millimolar" public_interface="in"/>
<variable name="R" units="joule_per_kilomole_kelvin" public_interface="in"/>
<variable name="F" units="coulomb_per_mole" public_interface="in"/>
<variable name="T" units="kelvin" public_interface="in"/>
<variable name="time" units="second" public_interface="in" private_interface="out"/>
<variable name="V" units="millivolt" public_interface="in" private_interface="out"/>
<variable name="d" units="dimensionless" private_interface="in"/>
<variable name="f" units="dimensionless" private_interface="in"/>
<variable name="f2" units="dimensionless" private_interface="in"/>
<variable name="f2ds" units="dimensionless" private_interface="in"/>
<variable name="Km_f2" units="millimolar" initial_value="100000" private_interface="out"/>
<variable name="Km_f2ds" units="millimolar" initial_value="0.001" private_interface="out"/>
<variable name="R_decay" units="per_second" initial_value="20" public_interface="out" private_interface="out"/>
<variable name="FrICa" units="dimensionless" initial_value="1"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>i_Ca_L_Ca_cyt</ci>
<apply>
<times/>
<apply>
<divide/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<ci>FrICa</ci>
</apply>
<cn cellml:units="dimensionless">4</cn>
<ci>P_Ca_L</ci>
<ci>d</ci>
<ci>f</ci>
<ci>f2</ci>
<apply>
<minus/>
<ci>V</ci>
<cn cellml:units="millivolt">50</cn>
</apply>
<ci>F</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<apply>
<minus/>
<ci>V</ci>
<cn cellml:units="millivolt">50</cn>
</apply>
</apply>
<ci>F</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<minus/>
<apply>
<times/>
<ci>Ca_i</ci>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<cn cellml:units="millivolt">100</cn>
<ci>F</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<times/>
<ci>Ca_o</ci>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<apply>
<minus/>
<ci>V</ci>
<cn cellml:units="millivolt">50</cn>
</apply>
</apply>
<ci>F</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>i_Ca_L_K_cyt</ci>
<apply>
<times/>
<apply>
<divide/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<ci>FrICa</ci>
</apply>
<ci>P_CaK</ci>
<ci>P_Ca_L</ci>
<ci>d</ci>
<ci>f</ci>
<ci>f2</ci>
<apply>
<minus/>
<ci>V</ci>
<cn cellml:units="millivolt">50</cn>
</apply>
<ci>F</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<apply>
<minus/>
<ci>V</ci>
<cn cellml:units="millivolt">50</cn>
</apply>
</apply>
<ci>F</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<minus/>
<apply>
<times/>
<ci>K_i</ci>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<cn cellml:units="millivolt">50</cn>
<ci>F</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<times/>
<ci>K_o</ci>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<apply>
<minus/>
<ci>V</ci>
<cn cellml:units="millivolt">50</cn>
</apply>
</apply>
<ci>F</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>i_Ca_L_Na_cyt</ci>
<apply>
<times/>
<apply>
<divide/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<ci>FrICa</ci>
</apply>
<ci>P_CaNa</ci>
<ci>P_Ca_L</ci>
<ci>d</ci>
<ci>f</ci>
<ci>f2</ci>
<apply>
<minus/>
<ci>V</ci>
<cn cellml:units="millivolt">50</cn>
</apply>
<ci>F</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<apply>
<minus/>
<ci>V</ci>
<cn cellml:units="millivolt">50</cn>
</apply>
</apply>
<ci>F</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<minus/>
<apply>
<times/>
<ci>Na_i</ci>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<cn cellml:units="millivolt">50</cn>
<ci>F</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<times/>
<ci>Na_o</ci>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<apply>
<minus/>
<ci>V</ci>
<cn cellml:units="millivolt">50</cn>
</apply>
</apply>
<ci>F</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>i_Ca_L_Ca_ds</ci>
<apply>
<times/>
<apply>
<divide/>
<apply>
<divide/>
<apply>
<times/>
<ci>FrICa</ci>
<cn cellml:units="dimensionless">4</cn>
<ci>P_Ca_L</ci>
<ci>d</ci>
<ci>f</ci>
<ci>f2ds</ci>
<apply>
<minus/>
<ci>V</ci>
<cn cellml:units="millivolt">50</cn>
</apply>
<ci>F</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<apply>
<minus/>
<ci>V</ci>
<cn cellml:units="millivolt">50</cn>
</apply>
</apply>
<ci>F</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<minus/>
<apply>
<times/>
<ci>Ca_i</ci>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<cn cellml:units="millivolt">100</cn>
<ci>F</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<times/>
<ci>Ca_o</ci>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<apply>
<minus/>
<ci>V</ci>
<cn cellml:units="millivolt">50</cn>
</apply>
</apply>
<ci>F</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>i_Ca_L_K_ds</ci>
<apply>
<times/>
<apply>
<divide/>
<apply>
<divide/>
<apply>
<times/>
<ci>FrICa</ci>
<ci>P_CaK</ci>
<ci>P_Ca_L</ci>
<ci>d</ci>
<ci>f</ci>
<ci>f2ds</ci>
<apply>
<minus/>
<ci>V</ci>
<cn cellml:units="millivolt">50</cn>
</apply>
<ci>F</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<apply>
<minus/>
<ci>V</ci>
<cn cellml:units="millivolt">50</cn>
</apply>
</apply>
<ci>F</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<minus/>
<apply>
<times/>
<ci>K_i</ci>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<cn cellml:units="millivolt">50</cn>
<ci>F</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<times/>
<ci>K_o</ci>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<apply>
<minus/>
<ci>V</ci>
<cn cellml:units="millivolt">50</cn>
</apply>
</apply>
<ci>F</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>i_Ca_L_Na_ds</ci>
<apply>
<times/>
<apply>
<divide/>
<apply>
<divide/>
<apply>
<times/>
<ci>FrICa</ci>
<ci>P_CaNa</ci>
<ci>P_Ca_L</ci>
<ci>d</ci>
<ci>f</ci>
<ci>f2ds</ci>
<apply>
<minus/>
<ci>V</ci>
<cn cellml:units="millivolt">50</cn>
</apply>
<ci>F</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<apply>
<minus/>
<ci>V</ci>
<cn cellml:units="millivolt">50</cn>
</apply>
</apply>
<ci>F</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<minus/>
<apply>
<times/>
<ci>Na_i</ci>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<cn cellml:units="millivolt">50</cn>
<ci>F</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<times/>
<ci>Na_o</ci>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<apply>
<minus/>
<ci>V</ci>
<cn cellml:units="millivolt">50</cn>
</apply>
</apply>
<ci>F</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>i_Ca_L</ci>
<apply>
<plus/>
<ci>i_Ca_L_Ca_cyt</ci>
<ci>i_Ca_L_K_cyt</ci>
<ci>i_Ca_L_Na_cyt</ci>
<ci>i_Ca_L_Ca_ds</ci>
<ci>i_Ca_L_K_ds</ci>
<ci>i_Ca_L_Na_ds</ci>
</apply>
</apply>
</math>
</component>
<component name="L_type_Ca_channel_d_gate">
<variable name="d" units="dimensionless" initial_value="0" public_interface="out"/>
<variable name="alpha_d" units="per_second"/>
<variable name="beta_d" units="per_second"/>
<variable name="V" units="millivolt" public_interface="in"/>
<variable name="time" units="second" public_interface="in"/>
<variable name="E0_d" units="millivolt"/>
<variable name="speed_d" units="dimensionless" initial_value="3"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>E0_d</ci>
<apply>
<minus/>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="millivolt">24</cn>
</apply>
<cn cellml:units="millivolt">5</cn>
</apply>
</apply>
<apply>
<eq/>
<ci>alpha_d</ci>
<piecewise>
<piece>
<cn cellml:units="per_second">120</cn>
<apply>
<lt/>
<apply>
<abs/>
<ci>E0_d</ci>
</apply>
<cn cellml:units="millivolt">0.0001</cn>
</apply>
</piece>
<otherwise>
<apply>
<divide/>
<apply>
<times/>
<cn cellml:units="per_millivolt_second">30</cn>
<ci>E0_d</ci>
</apply>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<ci>E0_d</ci>
</apply>
<cn cellml:units="millivolt">4</cn>
</apply>
</apply>
</apply>
</apply>
</otherwise>
</piecewise>
</apply>
<apply>
<eq/>
<ci>beta_d</ci>
<piecewise>
<piece>
<cn cellml:units="per_second">120</cn>
<apply>
<lt/>
<apply>
<abs/>
<ci>E0_d</ci>
</apply>
<cn cellml:units="millivolt">0.0001</cn>
</apply>
</piece>
<otherwise>
<apply>
<divide/>
<apply>
<times/>
<cn cellml:units="per_millivolt_second">12</cn>
<ci>E0_d</ci>
</apply>
<apply>
<minus/>
<apply>
<exp/>
<apply>
<divide/>
<ci>E0_d</ci>
<cn cellml:units="millivolt">10</cn>
</apply>
</apply>
<cn cellml:units="dimensionless">1</cn>
</apply>
</apply>
</otherwise>
</piecewise>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>d</ci>
</apply>
<apply>
<times/>
<ci>speed_d</ci>
<apply>
<minus/>
<apply>
<times/>
<ci>alpha_d</ci>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<ci>d</ci>
</apply>
</apply>
<apply>
<times/>
<ci>beta_d</ci>
<ci>d</ci>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<component name="L_type_Ca_channel_f_gate">
<variable name="f" units="dimensionless" initial_value="1" public_interface="out"/>
<variable name="alpha_f" units="per_second"/>
<variable name="beta_f" units="per_second"/>
<variable name="V" units="millivolt" public_interface="in"/>
<variable name="time" units="second" public_interface="in"/>
<variable name="speed_f" units="dimensionless" initial_value="0.3"/>
<variable name="delta_f" units="millivolt" initial_value="0.0001"/>
<variable name="E0_f" units="millivolt"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>E0_f</ci>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="millivolt">34</cn>
</apply>
</apply>
<apply>
<eq/>
<ci>alpha_f</ci>
<piecewise>
<piece>
<cn cellml:units="per_second">25</cn>
<apply>
<lt/>
<apply>
<abs/>
<ci>E0_f</ci>
</apply>
<ci>delta_f</ci>
</apply>
</piece>
<otherwise>
<apply>
<divide/>
<apply>
<times/>
<cn cellml:units="per_millivolt_second">6.25</cn>
<ci>E0_f</ci>
</apply>
<apply>
<minus/>
<apply>
<exp/>
<apply>
<divide/>
<ci>E0_f</ci>
<cn cellml:units="millivolt">4</cn>
</apply>
</apply>
<cn cellml:units="dimensionless">1</cn>
</apply>
</apply>
</otherwise>
</piecewise>
</apply>
<apply>
<eq/>
<ci>beta_f</ci>
<apply>
<divide/>
<cn cellml:units="per_second">12</cn>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
</apply>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="millivolt">34</cn>
</apply>
</apply>
<cn cellml:units="millivolt">4</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>f</ci>
</apply>
<apply>
<times/>
<ci>speed_f</ci>
<apply>
<minus/>
<apply>
<times/>
<ci>alpha_f</ci>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<ci>f</ci>
</apply>
</apply>
<apply>
<times/>
<ci>beta_f</ci>
<ci>f</ci>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<component name="L_type_Ca_channel_f2_gate">
<variable name="f2" units="dimensionless" initial_value="1" public_interface="out"/>
<variable name="Km_f2" units="millimolar" public_interface="in"/>
<variable name="Ca_i" units="millimolar" public_interface="in"/>
<variable name="V" units="millivolt" public_interface="in"/>
<variable name="time" units="second" public_interface="in"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>f2</ci>
</apply>
<apply>
<minus/>
<cn cellml:units="per_second">1</cn>
<apply>
<times/>
<cn cellml:units="per_second">1</cn>
<apply>
<plus/>
<apply>
<divide/>
<ci>Ca_i</ci>
<apply>
<plus/>
<ci>Km_f2</ci>
<ci>Ca_i</ci>
</apply>
</apply>
<ci>f2</ci>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<component name="L_type_Ca_channel_f2ds_gate">
<variable name="f2ds" units="dimensionless" initial_value="0.997" public_interface="out"/>
<variable name="Km_f2ds" units="millimolar" public_interface="in"/>
<variable name="R_decay" units="per_second" public_interface="in"/>
<variable name="Ca_ds" units="millimolar" public_interface="in"/>
<variable name="V" units="millivolt" public_interface="in"/>
<variable name="time" units="second" public_interface="in"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>f2ds</ci>
</apply>
<apply>
<times/>
<ci>R_decay</ci>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<plus/>
<apply>
<divide/>
<ci>Ca_ds</ci>
<apply>
<plus/>
<ci>Km_f2ds</ci>
<ci>Ca_ds</ci>
</apply>
</apply>
<ci>f2ds</ci>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<component name="calcium_background_current">
<variable name="i_b_Ca" units="nanoA" public_interface="out"/>
<variable name="g_bca" units="microS" initial_value="0.00025"/>
<variable name="E_Ca" units="millivolt" public_interface="in"/>
<variable name="time" units="second" public_interface="in"/>
<variable name="V" units="millivolt" public_interface="in"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>i_b_Ca</ci>
<apply>
<times/>
<ci>g_bca</ci>
<apply>
<minus/>
<ci>V</ci>
<ci>E_Ca</ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="transient_outward_current">
<variable name="i_to" units="nanoA" public_interface="out"/>
<variable name="g_to" units="microS" initial_value="0.006"/>
<variable name="g_tos" units="dimensionless" initial_value="0"/>
<variable name="E_K" units="millivolt" public_interface="in"/>
<variable name="V" units="millivolt" public_interface="in" private_interface="out"/>
<variable name="time" units="second" public_interface="in" private_interface="out"/>
<variable name="s" units="dimensionless" private_interface="in"/>
<variable name="r" units="dimensionless" private_interface="in"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>i_to</ci>
<apply>
<times/>
<ci>g_to</ci>
<apply>
<plus/>
<ci>g_tos</ci>
<apply>
<times/>
<ci>s</ci>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<ci>g_tos</ci>
</apply>
</apply>
</apply>
<ci>r</ci>
<apply>
<minus/>
<ci>V</ci>
<ci>E_K</ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="transient_outward_current_s_gate">
<variable name="s" units="dimensionless" initial_value="0.997044616031121" public_interface="out"/>
<variable name="alpha_s" units="per_second"/>
<variable name="beta_s" units="per_second"/>
<variable name="V" units="millivolt" public_interface="in"/>
<variable name="time" units="second" public_interface="in"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>alpha_s</ci>
<apply>
<times/>
<cn cellml:units="per_second">0.033</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<ci>V</ci>
</apply>
<cn cellml:units="millivolt">17</cn>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>beta_s</ci>
<apply>
<divide/>
<cn cellml:units="per_second">33</cn>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</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">10</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>s</ci>
</apply>
<apply>
<minus/>
<apply>
<times/>
<ci>alpha_s</ci>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<ci>s</ci>
</apply>
</apply>
<apply>
<times/>
<ci>beta_s</ci>
<ci>s</ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="transient_outward_current_r_gate">
<variable name="r" units="dimensionless" initial_value="1.63117173173398e-8" public_interface="out"/>
<variable name="V" units="millivolt" public_interface="in"/>
<variable name="time" units="second" public_interface="in"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>r</ci>
</apply>
<apply>
<times/>
<cn cellml:units="per_second">333</cn>
<apply>
<minus/>
<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">4</cn>
</apply>
</apply>
<cn cellml:units="millivolt">5</cn>
</apply>
</apply>
</apply>
</apply>
<ci>r</ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="sodium_potassium_pump">
<variable name="i_NaK" units="nanoA" public_interface="out"/>
<variable name="i_NaK_max" units="nanoA" initial_value="0.7"/>
<variable name="K_mK" units="millimolar" initial_value="1"/>
<variable name="K_mNa" units="millimolar" initial_value="24.2"/>
<variable name="K_o" units="millimolar" public_interface="in"/>
<variable name="Na_i" units="millimolar" public_interface="in"/>
<variable name="time" units="second" public_interface="in"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>i_NaK</ci>
<apply>
<divide/>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<ci>i_NaK_max</ci>
<ci>K_o</ci>
</apply>
<apply>
<plus/>
<ci>K_mK</ci>
<ci>K_o</ci>
</apply>
</apply>
<ci>Na_i</ci>
</apply>
<apply>
<plus/>
<ci>K_mNa</ci>
<ci>Na_i</ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="sodium_calcium_exchanger">
<variable name="i_NaCa" units="nanoA" public_interface="out"/>
<variable name="i_NaCa_cyt" units="nanoA" public_interface="out"/>
<variable name="i_NaCa_ds" units="nanoA" public_interface="out"/>
<variable name="k_NaCa" units="nanoA" initial_value="0.0005"/>
<variable name="n_NaCa" units="dimensionless" initial_value="3" public_interface="out"/>
<variable name="d_NaCa" units="dimensionless" initial_value="0"/>
<variable name="gamma" units="dimensionless" initial_value="0.5"/>
<variable name="FRiNaCa" units="dimensionless" initial_value="0.001"/>
<variable name="Na_i" units="millimolar" public_interface="in"/>
<variable name="Na_o" units="millimolar" public_interface="in"/>
<variable name="Ca_i" units="millimolar" public_interface="in"/>
<variable name="Ca_ds" units="millimolar" public_interface="in"/>
<variable name="Ca_o" units="millimolar" public_interface="in"/>
<variable name="R" units="joule_per_kilomole_kelvin" public_interface="in"/>
<variable name="F" units="coulomb_per_mole" public_interface="in"/>
<variable name="T" units="kelvin" public_interface="in"/>
<variable name="V" units="millivolt" public_interface="in"/>
<variable name="time" units="second" public_interface="in"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>i_NaCa_cyt</ci>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<ci>FRiNaCa</ci>
</apply>
<ci>k_NaCa</ci>
<apply>
<minus/>
<apply>
<times/>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<ci>gamma</ci>
<apply>
<minus/>
<ci>n_NaCa</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
<ci>V</ci>
<ci>F</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
<apply>
<power/>
<ci>Na_i</ci>
<ci>n_NaCa</ci>
</apply>
<ci>Ca_o</ci>
</apply>
<apply>
<times/>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<ci>gamma</ci>
<cn cellml:units="dimensionless">1</cn>
</apply>
<apply>
<minus/>
<ci>n_NaCa</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
<ci>V</ci>
<ci>F</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
<apply>
<power/>
<ci>Na_o</ci>
<ci>n_NaCa</ci>
</apply>
<ci>Ca_i</ci>
</apply>
</apply>
</apply>
<apply>
<times/>
<apply>
<plus/>
<cn cellml:units="millimolar4">1</cn>
<apply>
<times/>
<ci>d_NaCa</ci>
<apply>
<plus/>
<apply>
<times/>
<ci>Ca_i</ci>
<apply>
<power/>
<ci>Na_o</ci>
<ci>n_NaCa</ci>
</apply>
</apply>
<apply>
<times/>
<ci>Ca_o</ci>
<apply>
<power/>
<ci>Na_i</ci>
<ci>n_NaCa</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<divide/>
<ci>Ca_i</ci>
<cn cellml:units="millimolar">0.0069</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>i_NaCa_ds</ci>
<apply>
<divide/>
<apply>
<times/>
<ci>FRiNaCa</ci>
<ci>k_NaCa</ci>
<apply>
<minus/>
<apply>
<times/>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<ci>gamma</ci>
<apply>
<minus/>
<ci>n_NaCa</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
<ci>V</ci>
<ci>F</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
<apply>
<power/>
<ci>Na_i</ci>
<ci>n_NaCa</ci>
</apply>
<ci>Ca_o</ci>
</apply>
<apply>
<times/>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<ci>gamma</ci>
<cn cellml:units="dimensionless">1</cn>
</apply>
<apply>
<minus/>
<ci>n_NaCa</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
<ci>V</ci>
<ci>F</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
<apply>
<power/>
<ci>Na_o</ci>
<ci>n_NaCa</ci>
</apply>
<ci>Ca_ds</ci>
</apply>
</apply>
</apply>
<apply>
<times/>
<apply>
<plus/>
<cn cellml:units="millimolar4">1</cn>
<apply>
<times/>
<ci>d_NaCa</ci>
<apply>
<plus/>
<apply>
<times/>
<ci>Ca_ds</ci>
<apply>
<power/>
<ci>Na_o</ci>
<ci>n_NaCa</ci>
</apply>
</apply>
<apply>
<times/>
<ci>Ca_o</ci>
<apply>
<power/>
<ci>Na_i</ci>
<ci>n_NaCa</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<divide/>
<ci>Ca_ds</ci>
<cn cellml:units="millimolar">0.0069</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>i_NaCa</ci>
<apply>
<plus/>
<ci>i_NaCa_cyt</ci>
<ci>i_NaCa_ds</ci>
</apply>
</apply>
</math>
</component>
<component name="sarcoplasmic_reticulum_calcium_pump">
<variable name="i_up" units="millimolar_per_second" public_interface="out"/>
<variable name="K_1" units="dimensionless"/>
<variable name="K_2" units="millimolar"/>
<variable name="K_cyca" units="millimolar" initial_value="0.00015"/>
<variable name="K_xcs" units="dimensionless" initial_value="0.4"/>
<variable name="K_srca" units="millimolar" initial_value="0.5"/>
<variable name="alpha_up" units="millimolar_per_second" initial_value="1"/>
<variable name="beta_up" units="millimolar_per_second" initial_value="0.03"/>
<variable name="Ca_i" units="millimolar" public_interface="in"/>
<variable name="Ca_up" units="millimolar" public_interface="in"/>
<variable name="time" units="second" public_interface="in"/>
<variable name="flag_ingib" units="dimensionless" initial_value="0"/>
<variable name="K_inh" units="millimolar" initial_value="4"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>K_1</ci>
<apply>
<divide/>
<apply>
<times/>
<ci>K_cyca</ci>
<ci>K_xcs</ci>
</apply>
<ci>K_srca</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>K_2</ci>
<apply>
<plus/>
<ci>Ca_i</ci>
<apply>
<times/>
<ci>Ca_up</ci>
<ci>K_1</ci>
</apply>
<apply>
<times/>
<ci>K_cyca</ci>
<ci>K_xcs</ci>
</apply>
<ci>K_cyca</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>i_up</ci>
<piecewise>
<piece>
<apply>
<minus/>
<apply>
<times/>
<apply>
<divide/>
<ci>Ca_i</ci>
<ci>K_2</ci>
</apply>
<ci>alpha_up</ci>
</apply>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<ci>Ca_up</ci>
<ci>K_1</ci>
</apply>
<ci>K_2</ci>
</apply>
<ci>beta_up</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>flag_ingib</ci>
<cn cellml:units="dimensionless">0</cn>
</apply>
</piece>
<otherwise>
<apply>
<minus/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<divide/>
<ci>Ca_i</ci>
<ci>K_2</ci>
</apply>
<ci>alpha_up</ci>
</apply>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<divide/>
<ci>Ca_up</ci>
<ci>K_inh</ci>
</apply>
</apply>
</apply>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<ci>Ca_up</ci>
<ci>K_1</ci>
</apply>
<ci>K_2</ci>
</apply>
<ci>beta_up</ci>
</apply>
</apply>
</otherwise>
</piecewise>
</apply>
</math>
</component>
<component name="calcium_translocation">
<variable name="i_trans" units="millimolar_per_second" public_interface="out"/>
<variable name="Ca_rel" units="millimolar" public_interface="in"/>
<variable name="Ca_up" units="millimolar" public_interface="in"/>
<variable name="time" units="second" public_interface="in"/>
<variable name="a_tr" units="per_second" initial_value="15"/>
<variable name="alpha_CaS" units="per_millimolar_second" initial_value="50000"/>
<variable name="beta_CaS" units="per_second" initial_value="32500"/>
<variable name="CaS_tot" units="millimolar" initial_value="40" public_interface="out"/>
<variable name="beta" units="millimolar" public_interface="out"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>i_trans</ci>
<apply>
<times/>
<ci>a_tr</ci>
<apply>
<minus/>
<ci>Ca_up</ci>
<ci>Ca_rel</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>beta</ci>
<apply>
<divide/>
<ci>beta_CaS</ci>
<ci>alpha_CaS</ci>
</apply>
</apply>
</math>
</component>
<component name="calcium_release">
<variable name="i_rel" units="millimolar_per_second" public_interface="out"/>
<variable name="VoltDep" units="dimensionless"/>
<variable name="RegBindSite" units="dimensionless"/>
<variable name="CaiReg" units="dimensionless"/>
<variable name="CadsReg" units="dimensionless"/>
<variable name="ActRate" units="per_second"/>
<variable name="InactRate" units="per_second"/>
<variable name="SRLeak" units="per_second" initial_value="0.05"/>
<variable name="K_m_rel" units="per_second" initial_value="10000"/>
<variable name="K_m_Ca_cyt" units="millimolar" initial_value="0.0005"/>
<variable name="K_m_Ca_ds" units="millimolar" initial_value="0.01"/>
<variable name="PrecFrac" units="dimensionless"/>
<variable name="ActFrac" units="dimensionless" initial_value="0.001914"/>
<variable name="ProdFrac" units="dimensionless" initial_value="0.2854569"/>
<variable name="SpeedRel" units="dimensionless"/>
<variable name="time" units="second" public_interface="in"/>
<variable name="V" units="millivolt" public_interface="in"/>
<variable name="Ca_i" units="millimolar" public_interface="in"/>
<variable name="Ca_ds" units="millimolar" public_interface="in"/>
<variable name="Ca_rel" units="millimolar" public_interface="in"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>VoltDep</ci>
<apply>
<exp/>
<apply>
<times/>
<cn cellml:units="per_millivolt">0.08</cn>
<apply>
<minus/>
<ci>V</ci>
<cn cellml:units="millivolt">40</cn>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>CaiReg</ci>
<apply>
<divide/>
<ci>Ca_i</ci>
<apply>
<plus/>
<ci>Ca_i</ci>
<ci>K_m_Ca_cyt</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>CadsReg</ci>
<apply>
<divide/>
<ci>Ca_ds</ci>
<apply>
<plus/>
<ci>Ca_ds</ci>
<ci>K_m_Ca_ds</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>RegBindSite</ci>
<apply>
<plus/>
<ci>CaiReg</ci>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<ci>CaiReg</ci>
</apply>
<ci>CadsReg</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>ActRate</ci>
<apply>
<times/>
<cn cellml:units="per_second">500</cn>
<apply>
<power/>
<ci>RegBindSite</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>InactRate</ci>
<apply>
<plus/>
<cn cellml:units="per_second">60</cn>
<apply>
<times/>
<cn cellml:units="per_second">500</cn>
<apply>
<power/>
<ci>RegBindSite</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>SpeedRel</ci>
<piecewise>
<piece>
<cn cellml:units="dimensionless">5</cn>
<apply>
<lt/>
<ci>V</ci>
<apply>
<minus/>
<cn cellml:units="millivolt">50</cn>
</apply>
</apply>
</piece>
<otherwise>
<cn cellml:units="dimensionless">1</cn>
</otherwise>
</piecewise>
</apply>
<apply>
<eq/>
<ci>PrecFrac</ci>
<apply>
<minus/>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<ci>ActFrac</ci>
</apply>
<ci>ProdFrac</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>ActFrac</ci>
</apply>
<apply>
<minus/>
<apply>
<times/>
<ci>PrecFrac</ci>
<ci>SpeedRel</ci>
<ci>ActRate</ci>
</apply>
<apply>
<times/>
<ci>ActFrac</ci>
<ci>SpeedRel</ci>
<ci>InactRate</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>ProdFrac</ci>
</apply>
<apply>
<minus/>
<apply>
<times/>
<ci>ActFrac</ci>
<ci>SpeedRel</ci>
<ci>InactRate</ci>
</apply>
<apply>
<times/>
<ci>SpeedRel</ci>
<cn cellml:units="per_second">0.6</cn>
<ci>ProdFrac</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>i_rel</ci>
<apply>
<times/>
<apply>
<plus/>
<apply>
<times/>
<apply>
<power/>
<apply>
<divide/>
<ci>ActFrac</ci>
<apply>
<plus/>
<ci>ActFrac</ci>
<cn cellml:units="dimensionless">0.25</cn>
</apply>
</apply>
<cn cellml:units="dimensionless">2</cn>
</apply>
<ci>K_m_rel</ci>
</apply>
<ci>SRLeak</ci>
</apply>
<ci>Ca_rel</ci>
</apply>
</apply>
</math>
</component>
<component name="extracellular_sodium_concentration">
<variable name="Na_o" units="millimolar" initial_value="140" public_interface="out"/>
</component>
<component name="intracellular_sodium_concentration">
<variable name="Na_i" units="millimolar" initial_value="5.18787513289509" public_interface="out"/>
<variable name="V_i" units="micrometre3" public_interface="in"/>
<variable name="time" units="second" public_interface="in"/>
<variable name="F" units="coulomb_per_mole" public_interface="in"/>
<variable name="i_Na" units="nanoA" public_interface="in"/>
<variable name="i_b_Na" units="nanoA" public_interface="in"/>
<variable name="i_p_Na" units="nanoA" public_interface="in"/>
<variable name="i_Ca_L_Na_cyt" units="nanoA" public_interface="in"/>
<variable name="i_Ca_L_Na_ds" units="nanoA" public_interface="in"/>
<variable name="i_NaK" units="nanoA" public_interface="in"/>
<variable name="i_NaCa_cyt" units="nanoA" public_interface="in"/>
<variable name="n_NaK" units="dimensionless" initial_value="1.5" public_interface="out"/>
<variable name="n_NaCa" units="dimensionless" public_interface="in"/>
<variable name="i_NaCa" units="nanoA" public_interface="in"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>Na_i</ci>
</apply>
<apply>
<times/>
<apply>
<divide/>
<apply>
<minus/>
<cn cellml:units="mA_nA">1</cn>
</apply>
<apply>
<times/>
<cn cellml:units="litre_micrometre3">1</cn>
<ci>V_i</ci>
<ci>F</ci>
</apply>
</apply>
<apply>
<plus/>
<ci>i_Na</ci>
<ci>i_p_Na</ci>
<ci>i_b_Na</ci>
<ci>i_Ca_L_Na_cyt</ci>
<ci>i_Ca_L_Na_ds</ci>
<apply>
<times/>
<apply>
<divide/>
<ci>n_NaK</ci>
<apply>
<minus/>
<ci>n_NaK</ci>
<cn cellml:units="dimensionless">1</cn>
</apply>
</apply>
<ci>i_NaK</ci>
</apply>
<apply>
<times/>
<apply>
<divide/>
<ci>n_NaCa</ci>
<apply>
<minus/>
<ci>n_NaCa</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
<ci>i_NaCa</ci>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<component name="extracellular_calcium_concentration">
<variable name="Ca_o" units="millimolar" initial_value="2" public_interface="out"/>
</component>
<component name="extracellular_potassium_concentration">
<variable name="K_o" units="millimolar" initial_value="3.988" public_interface="out"/>
<variable name="K_b" units="millimolar" initial_value="4"/>
<variable name="i_Kr" units="nanoA" public_interface="in"/>
<variable name="i_Ks" units="nanoA" public_interface="in"/>
<variable name="i_K1" units="nanoA" public_interface="in"/>
<variable name="i_to" units="nanoA" public_interface="in"/>
<variable name="i_NaK" units="nanoA" public_interface="in"/>
<variable name="i_Ca_L_K_cyt" units="nanoA" public_interface="in"/>
<variable name="i_Ca_L_K_ds" units="nanoA" public_interface="in"/>
<variable name="time" units="second" public_interface="in"/>
<variable name="n_NaK" units="dimensionless" public_interface="in"/>
<variable name="F" units="coulomb_per_mole" public_interface="in"/>
<variable name="V_e" units="micrometre3" public_interface="in"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>K_o</ci>
</apply>
<apply>
<minus/>
<apply>
<divide/>
<apply>
<times/>
<cn cellml:units="mA_nA">1</cn>
<apply>
<plus/>
<apply>
<minus/>
<apply>
<plus/>
<ci>i_Kr</ci>
<ci>i_Ks</ci>
<ci>i_K1</ci>
<ci>i_to</ci>
</apply>
<apply>
<times/>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<minus/>
<ci>n_NaK</ci>
<cn cellml:units="dimensionless">1</cn>
</apply>
</apply>
<ci>i_NaK</ci>
</apply>
</apply>
<ci>i_Ca_L_K_cyt</ci>
<ci>i_Ca_L_K_ds</ci>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="litre_micrometre3">1</cn>
<ci>V_e</ci>
<ci>F</ci>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="per_second">0.7</cn>
<apply>
<minus/>
<ci>K_o</ci>
<ci>K_b</ci>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<component name="intracellular_potassium_concentration">
<variable name="K_i" units="millimolar" initial_value="139.054" public_interface="out"/>
<variable name="V_i" units="micrometre3" public_interface="in"/>
<variable name="i_K1" units="nanoA" public_interface="in"/>
<variable name="i_Kr" units="nanoA" public_interface="in"/>
<variable name="i_Ks" units="nanoA" public_interface="in"/>
<variable name="i_Ca_L_K_cyt" units="nanoA" public_interface="in"/>
<variable name="i_Ca_L_K_ds" units="nanoA" public_interface="in"/>
<variable name="i_to" units="nanoA" public_interface="in"/>
<variable name="i_NaK" units="nanoA" public_interface="in"/>
<variable name="time" units="second" public_interface="in"/>
<variable name="F" units="coulomb_per_mole" public_interface="in"/>
<variable name="n_NaK" units="dimensionless" public_interface="in"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>K_i</ci>
</apply>
<apply>
<times/>
<apply>
<divide/>
<apply>
<minus/>
<cn cellml:units="mA_nA">1</cn>
</apply>
<apply>
<times/>
<cn cellml:units="litre_micrometre3">1</cn>
<ci>V_i</ci>
<ci>F</ci>
</apply>
</apply>
<apply>
<minus/>
<apply>
<plus/>
<ci>i_K1</ci>
<ci>i_Kr</ci>
<ci>i_Ks</ci>
<ci>i_Ca_L_K_cyt</ci>
<ci>i_Ca_L_K_ds</ci>
<ci>i_to</ci>
</apply>
<apply>
<times/>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<minus/>
<ci>n_NaK</ci>
<cn cellml:units="dimensionless">1</cn>
</apply>
</apply>
<ci>i_NaK</ci>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<component name="intracellular_calcium_concentration">
<variable name="Ca_i" units="millimolar" initial_value="6.15e-6" public_interface="out"/>
<variable name="Ca_ds" units="millimolar" initial_value="2.55e-6" public_interface="out"/>
<variable name="Ca_up" units="millimolar" initial_value="0.994579" public_interface="out"/>
<variable name="Ca_rel" units="millimolar" initial_value="0.989665" public_interface="out"/>
<variable name="radius" units="micrometre" initial_value="12"/>
<variable name="length" units="micrometre" initial_value="74"/>
<variable name="V_i" units="micrometre3" public_interface="out"/>
<variable name="V_e" units="micrometre3" public_interface="out"/>
<variable name="V_Cell" units="micrometre3"/>
<variable name="V_i_ratio" units="dimensionless" initial_value="0.49"/>
<variable name="V_ds_ratio" units="dimensionless" initial_value="0.1"/>
<variable name="V_rel_ratio" units="dimensionless" initial_value="0.003"/>
<variable name="V_e_ratio" units="dimensionless" initial_value="0.4"/>
<variable name="V_up_ratio" units="dimensionless" initial_value="0.03"/>
<variable name="Kdecay" units="per_second" initial_value="10"/>
<variable name="i_up" units="millimolar_per_second" public_interface="in"/>
<variable name="i_trans" units="millimolar_per_second" public_interface="in"/>
<variable name="i_rel" units="millimolar_per_second" public_interface="in"/>
<variable name="i_NaCa_cyt" units="nanoA" public_interface="in"/>
<variable name="i_NaCa_ds" units="nanoA" public_interface="in"/>
<variable name="i_Ca_L_Ca_cyt" units="nanoA" public_interface="in"/>
<variable name="i_Ca_L_Ca_ds" units="nanoA" public_interface="in"/>
<variable name="i_b_Ca" units="nanoA" public_interface="in"/>
<variable name="F" units="coulomb_per_mole" public_interface="in"/>
<variable name="time" units="second" public_interface="in"/>
<variable name="A" units="millimolar" initial_value="0.00015" public_interface="out"/>
<variable name="N_A" units="dimensionless"/>
<variable name="pi_N_A" units="dimensionless"/>
<variable name="pi_min" units="dimensionless" initial_value="0.03"/>
<variable name="N" units="dimensionless" public_interface="in"/>
<variable name="a_on" units="per_millimolar_second" initial_value="70000"/>
<variable name="a_off" units="per_second" initial_value="200"/>
<variable name="k_A" units="per_millimolar" public_interface="in"/>
<variable name="L_oz" units="dimensionless" public_interface="in"/>
<variable name="A_tot" units="millimolar" initial_value="0.07" public_interface="out"/>
<variable name="B_1" units="millimolar" initial_value="0" public_interface="out"/>
<variable name="B_1_tot" units="millimolar" initial_value="0.08"/>
<variable name="b_1_on" units="per_millimolar_second" initial_value="100000"/>
<variable name="b_1_off" units="per_second" initial_value="182"/>
<variable name="B_2" units="millimolar" initial_value="0" public_interface="out"/>
<variable name="B_2_tot" units="millimolar" initial_value="0.1"/>
<variable name="b_2_on" units="per_millimolar_second" initial_value="1000"/>
<variable name="b_2_off" units="per_second" initial_value="3"/>
<variable name="CaS_tot" units="millimolar" public_interface="in"/>
<variable name="beta" units="millimolar" public_interface="in"/>
<variable name="n_NaCa" units="dimensionless" public_interface="in"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>N_A</ci>
<apply>
<divide/>
<apply>
<times/>
<ci>A_tot</ci>
<ci>N</ci>
</apply>
<apply>
<times/>
<ci>L_oz</ci>
<ci>A</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>pi_N_A</ci>
<piecewise>
<piece>
<ci>pi_min</ci>
<apply>
<leq/>
<ci>N_A</ci>
<cn cellml:units="dimensionless">0</cn>
</apply>
</piece>
<piece>
<apply>
<power/>
<ci>pi_min</ci>
<ci>N_A</ci>
</apply>
<apply>
<leq/>
<ci>N_A</ci>
<cn cellml:units="dimensionless">1</cn>
</apply>
</piece>
<otherwise>
<cn cellml:units="dimensionless">1</cn>
</otherwise>
</piecewise>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>A</ci>
</apply>
<apply>
<minus/>
<apply>
<times/>
<ci>a_on</ci>
<apply>
<minus/>
<ci>A_tot</ci>
<ci>A</ci>
</apply>
<ci>Ca_i</ci>
</apply>
<apply>
<times/>
<ci>a_off</ci>
<apply>
<exp/>
<apply>
<times/>
<apply>
<minus/>
<ci>k_A</ci>
</apply>
<ci>A</ci>
</apply>
</apply>
<ci>pi_N_A</ci>
<ci>A</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>B_1</ci>
</apply>
<apply>
<minus/>
<apply>
<times/>
<ci>b_1_on</ci>
<apply>
<minus/>
<ci>B_1_tot</ci>
<ci>B_1</ci>
</apply>
<ci>Ca_i</ci>
</apply>
<apply>
<times/>
<ci>b_1_off</ci>
<ci>B_1</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>B_2</ci>
</apply>
<apply>
<minus/>
<apply>
<times/>
<ci>b_2_on</ci>
<apply>
<minus/>
<ci>B_2_tot</ci>
<ci>B_2</ci>
</apply>
<ci>Ca_i</ci>
</apply>
<apply>
<times/>
<ci>b_2_off</ci>
<ci>B_2</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>V_Cell</ci>
<apply>
<divide/>
<apply>
<times/>
<cn cellml:units="dimensionless">3.141592654</cn>
<apply>
<power/>
<apply>
<divide/>
<ci>radius</ci>
<cn cellml:units="dimensionless">1000</cn>
</apply>
<cn cellml:units="dimensionless">2</cn>
</apply>
<ci>length</ci>
</apply>
<cn cellml:units="dimensionless">1000</cn>
</apply>
</apply>
<apply>
<eq/>
<ci>V_i</ci>
<apply>
<times/>
<ci>V_Cell</ci>
<ci>V_i_ratio</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>V_e</ci>
<apply>
<times/>
<ci>V_Cell</ci>
<ci>V_e_ratio</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>Ca_i</ci>
</apply>
<apply>
<minus/>
<apply>
<minus/>
<apply>
<minus/>
<apply>
<minus/>
<apply>
<plus/>
<apply>
<times/>
<apply>
<divide/>
<apply>
<minus/>
<cn cellml:units="mA_nA">1</cn>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">2</cn>
<cn cellml:units="litre_micrometre3">1</cn>
<ci>V_i</ci>
<ci>F</ci>
</apply>
</apply>
<apply>
<minus/>
<apply>
<plus/>
<ci>i_Ca_L_Ca_cyt</ci>
<ci>i_b_Ca</ci>
</apply>
<apply>
<times/>
<apply>
<divide/>
<cn cellml:units="dimensionless">2</cn>
<apply>
<minus/>
<ci>n_NaCa</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
<ci>i_NaCa_cyt</ci>
</apply>
</apply>
</apply>
<apply>
<times/>
<ci>Ca_ds</ci>
<ci>V_ds_ratio</ci>
<ci>Kdecay</ci>
</apply>
<apply>
<divide/>
<apply>
<times/>
<ci>i_rel</ci>
<ci>V_rel_ratio</ci>
</apply>
<ci>V_i_ratio</ci>
</apply>
</apply>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>A</ci>
</apply>
</apply>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>B_1</ci>
</apply>
</apply>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>B_2</ci>
</apply>
</apply>
<ci>i_up</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>Ca_ds</ci>
</apply>
<apply>
<minus/>
<apply>
<divide/>
<apply>
<plus/>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="mA_nA">1</cn>
</apply>
<ci>i_Ca_L_Ca_ds</ci>
</apply>
<apply>
<divide/>
<apply>
<times/>
<cn cellml:units="mA_nA">2</cn>
<ci>i_NaCa_ds</ci>
</apply>
<apply>
<minus/>
<ci>n_NaCa</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">2</cn>
<cn cellml:units="litre_micrometre3">1</cn>
<ci>V_ds_ratio</ci>
<ci>V_i</ci>
<ci>F</ci>
</apply>
</apply>
<apply>
<times/>
<ci>Ca_ds</ci>
<ci>Kdecay</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>Ca_up</ci>
</apply>
<apply>
<minus/>
<apply>
<times/>
<apply>
<divide/>
<ci>V_i_ratio</ci>
<ci>V_up_ratio</ci>
</apply>
<ci>i_up</ci>
</apply>
<ci>i_trans</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>Ca_rel</ci>
</apply>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<times/>
<apply>
<divide/>
<ci>V_up_ratio</ci>
<ci>V_rel_ratio</ci>
</apply>
<ci>i_trans</ci>
</apply>
<ci>i_rel</ci>
</apply>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<divide/>
<apply>
<times/>
<ci>beta</ci>
<ci>CaS_tot</ci>
</apply>
<apply>
<power/>
<apply>
<plus/>
<ci>Ca_rel</ci>
<ci>beta</ci>
</apply>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<component name="fibroblast">
<variable name="g_fibro_junct" units="microS" initial_value="2.9e-4"/>
<variable name="g_fibro" units="microS" initial_value="2e-4"/>
<variable name="c_fibro" units="microF" initial_value="1e-5"/>
<variable name="g_fibro_stretch" units="microS" initial_value="0"/>
<variable name="E_fibro_stretch" units="millivolt" initial_value="0"/>
<variable name="V_fibro" units="millivolt" initial_value="-20"/>
<variable name="i_fibro" units="nanoA"/>
<variable name="i_fibro_junct" units="nanoA"/>
<variable name="V" units="millivolt" public_interface="in"/>
<variable name="time" units="second" public_interface="in"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>i_fibro</ci>
<apply>
<plus/>
<apply>
<times/>
<ci>g_fibro</ci>
<apply>
<plus/>
<ci>V_fibro</ci>
<cn cellml:units="millivolt">20</cn>
</apply>
</apply>
<apply>
<times/>
<ci>g_fibro_stretch</ci>
<apply>
<minus/>
<ci>V_fibro</ci>
<ci>E_fibro_stretch</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>i_fibro_junct</ci>
<apply>
<times/>
<apply>
<minus/>
<ci>g_fibro_junct</ci>
</apply>
<apply>
<minus/>
<ci>V</ci>
<ci>V_fibro</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>V_fibro</ci>
</apply>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<plus/>
<ci>i_fibro</ci>
<ci>i_fibro_junct</ci>
</apply>
</apply>
<ci>c_fibro</ci>
</apply>
</apply>
</math>
</component>
<group>
<relationship_ref relationship="containment"/>
<component_ref component="membrane">
<component_ref component="fast_sodium_current">
<component_ref component="fast_sodium_current_m_gate"/>
<component_ref component="fast_sodium_current_h_gate"/>
</component_ref>
<component_ref component="L_type_Ca_channel">
<component_ref component="L_type_Ca_channel_d_gate"/>
<component_ref component="L_type_Ca_channel_f_gate"/>
<component_ref component="L_type_Ca_channel_f2_gate"/>
<component_ref component="L_type_Ca_channel_f2ds_gate"/>
</component_ref>
<component_ref component="rapid_delayed_rectifier_potassium_current">
<component_ref component="rapid_delayed_rectifier_potassium_current_xr1_gate"/>
<component_ref component="rapid_delayed_rectifier_potassium_current_xr2_gate"/>
</component_ref>
<component_ref component="slow_delayed_rectifier_potassium_current">
<component_ref component="slow_delayed_rectifier_potassium_current_xs_gate"/>
</component_ref>
<component_ref component="time_independent_potassium_current"/>
<component_ref component="sodium_activated_potassium_current"/>
<component_ref component="persistent_sodium_current"/>
<component_ref component="sodium_background_current"/>
<component_ref component="calcium_background_current"/>
<component_ref component="transient_outward_current">
<component_ref component="transient_outward_current_s_gate"/>
<component_ref component="transient_outward_current_r_gate"/>
</component_ref>
<component_ref component="sodium_potassium_pump"/>
<component_ref component="sodium_calcium_exchanger"/>
<component_ref component="sarcoplasmic_reticulum_calcium_pump"/>
<component_ref component="calcium_translocation"/>
<component_ref component="calcium_release"/>
<component_ref component="extracellular_sodium_concentration"/>
<component_ref component="intracellular_sodium_concentration"/>
<component_ref component="extracellular_potassium_concentration"/>
<component_ref component="intracellular_potassium_concentration"/>
<component_ref component="extracellular_calcium_concentration"/>
<component_ref component="intracellular_calcium_concentration"/>
<component_ref component="reversal_potentials"/>
<component_ref component="fibroblast"/>
</component_ref>
</group>
<group>
<relationship_ref relationship="encapsulation"/>
<component_ref component="fast_sodium_current">
<component_ref component="fast_sodium_current_m_gate"/>
<component_ref component="fast_sodium_current_h_gate"/>
</component_ref>
<component_ref component="L_type_Ca_channel">
<component_ref component="L_type_Ca_channel_d_gate"/>
<component_ref component="L_type_Ca_channel_f_gate"/>
<component_ref component="L_type_Ca_channel_f2_gate"/>
<component_ref component="L_type_Ca_channel_f2ds_gate"/>
</component_ref>
<component_ref component="transient_outward_current">
<component_ref component="transient_outward_current_s_gate"/>
<component_ref component="transient_outward_current_r_gate"/>
</component_ref>
<component_ref component="rapid_delayed_rectifier_potassium_current">
<component_ref component="rapid_delayed_rectifier_potassium_current_xr1_gate"/>
<component_ref component="rapid_delayed_rectifier_potassium_current_xr2_gate"/>
</component_ref>
<component_ref component="slow_delayed_rectifier_potassium_current">
<component_ref component="slow_delayed_rectifier_potassium_current_xs_gate"/>
</component_ref>
</group>
<connection>
<map_components component_1="membrane" component_2="environment"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection>
<map_components component_1="reversal_potentials" component_2="environment"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection>
<map_components component_1="fast_sodium_current" component_2="environment"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection>
<map_components component_1="L_type_Ca_channel" component_2="environment"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection>
<map_components component_1="rapid_delayed_rectifier_potassium_current" component_2="environment"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection>
<map_components component_1="slow_delayed_rectifier_potassium_current" component_2="environment"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection>
<map_components component_1="sodium_activated_potassium_current" component_2="environment"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection>
<map_components component_1="time_independent_potassium_current" component_2="environment"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection>
<map_components component_1="persistent_sodium_current" component_2="environment"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection>
<map_components component_1="sodium_background_current" component_2="environment"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection>
<map_components component_1="calcium_background_current" component_2="environment"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection>
<map_components component_1="transient_outward_current" component_2="environment"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection>
<map_components component_1="sodium_potassium_pump" component_2="environment"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection>
<map_components component_1="sodium_calcium_exchanger" component_2="environment"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection>
<map_components component_1="intracellular_sodium_concentration" component_2="environment"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection>
<map_components component_1="intracellular_potassium_concentration" component_2="environment"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection>
<map_components component_1="intracellular_calcium_concentration" component_2="environment"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection>
<map_components component_1="sarcoplasmic_reticulum_calcium_pump" component_2="environment"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection>
<map_components component_1="calcium_translocation" component_2="environment"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection>
<map_components component_1="calcium_release" component_2="environment"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection>
<map_components component_1="fibroblast" component_2="environment"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection>
<map_components component_1="calcium_release" component_2="membrane"/>
<map_variables variable_1="V" variable_2="V"/>
</connection>
<connection>
<map_components component_1="membrane" component_2="reversal_potentials"/>
<map_variables variable_1="R" variable_2="R"/>
<map_variables variable_1="T" variable_2="T"/>
<map_variables variable_1="F" variable_2="F"/>
</connection>
<connection>
<map_components component_1="membrane" component_2="fast_sodium_current"/>
<map_variables variable_1="V" variable_2="V"/>
<map_variables variable_1="i_Na" variable_2="i_Na"/>
</connection>
<connection>
<map_components component_1="membrane" component_2="L_type_Ca_channel"/>
<map_variables variable_1="V" variable_2="V"/>
<map_variables variable_1="i_Ca_L_Ca_cyt" variable_2="i_Ca_L_Ca_cyt"/>
<map_variables variable_1="i_Ca_L_Ca_ds" variable_2="i_Ca_L_Ca_ds"/>
<map_variables variable_1="i_Ca_L_Na_cyt" variable_2="i_Ca_L_Na_cyt"/>
<map_variables variable_1="i_Ca_L_Na_ds" variable_2="i_Ca_L_Na_ds"/>
<map_variables variable_1="i_Ca_L_K_cyt" variable_2="i_Ca_L_K_cyt"/>
<map_variables variable_1="i_Ca_L_K_ds" variable_2="i_Ca_L_K_ds"/>
<map_variables variable_1="R" variable_2="R"/>
<map_variables variable_1="T" variable_2="T"/>
<map_variables variable_1="F" variable_2="F"/>
</connection>
<connection>
<map_components component_1="membrane" component_2="rapid_delayed_rectifier_potassium_current"/>
<map_variables variable_1="V" variable_2="V"/>
<map_variables variable_1="i_Kr" variable_2="i_Kr"/>
</connection>
<connection>
<map_components component_1="membrane" component_2="slow_delayed_rectifier_potassium_current"/>
<map_variables variable_1="V" variable_2="V"/>
<map_variables variable_1="i_Ks" variable_2="i_Ks"/>
</connection>
<connection>
<map_components component_1="membrane" component_2="transient_outward_current"/>
<map_variables variable_1="V" variable_2="V"/>
<map_variables variable_1="i_to" variable_2="i_to"/>
</connection>
<connection>
<map_components component_1="membrane" component_2="time_independent_potassium_current"/>
<map_variables variable_1="V" variable_2="V"/>
<map_variables variable_1="i_K1" variable_2="i_K1"/>
<map_variables variable_1="R" variable_2="R"/>
<map_variables variable_1="T" variable_2="T"/>
<map_variables variable_1="F" variable_2="F"/>
</connection>
<connection>
<map_components component_1="membrane" component_2="sodium_activated_potassium_current"/>
<map_variables variable_1="V" variable_2="V"/>
</connection>
<connection>
<map_components component_1="membrane" component_2="persistent_sodium_current"/>
<map_variables variable_1="V" variable_2="V"/>
<map_variables variable_1="i_p_Na" variable_2="i_p_Na"/>
</connection>
<connection>
<map_components component_1="membrane" component_2="sodium_background_current"/>
<map_variables variable_1="V" variable_2="V"/>
<map_variables variable_1="i_b_Na" variable_2="i_b_Na"/>
</connection>
<connection>
<map_components component_1="membrane" component_2="calcium_background_current"/>
<map_variables variable_1="V" variable_2="V"/>
<map_variables variable_1="i_b_Ca" variable_2="i_b_Ca"/>
</connection>
<connection>
<map_components component_1="membrane" component_2="sodium_potassium_pump"/>
<map_variables variable_1="i_NaK" variable_2="i_NaK"/>
</connection>
<connection>
<map_components component_1="membrane" component_2="sodium_calcium_exchanger"/>
<map_variables variable_1="V" variable_2="V"/>
<map_variables variable_1="i_NaCa_cyt" variable_2="i_NaCa_cyt"/>
<map_variables variable_1="i_NaCa_ds" variable_2="i_NaCa_ds"/>
<map_variables variable_1="R" variable_2="R"/>
<map_variables variable_1="T" variable_2="T"/>
<map_variables variable_1="F" variable_2="F"/>
</connection>
<connection>
<map_components component_1="membrane" component_2="fibroblast"/>
<map_variables variable_1="V" variable_2="V"/>
</connection>
<connection>
<map_components component_1="membrane" component_2="intracellular_sodium_concentration"/>
<map_variables variable_1="F" variable_2="F"/>
</connection>
<connection>
<map_components component_1="membrane" component_2="intracellular_potassium_concentration"/>
<map_variables variable_1="F" variable_2="F"/>
</connection>
<connection>
<map_components component_1="membrane" component_2="intracellular_calcium_concentration"/>
<map_variables variable_1="F" variable_2="F"/>
</connection>
<connection>
<map_components component_1="reversal_potentials" component_2="sodium_background_current"/>
<map_variables variable_1="E_Na" variable_2="E_Na"/>
</connection>
<connection>
<map_components component_1="reversal_potentials" component_2="persistent_sodium_current"/>
<map_variables variable_1="E_Na" variable_2="E_Na"/>
</connection>
<connection>
<map_components component_1="reversal_potentials" component_2="time_independent_potassium_current"/>
<map_variables variable_1="E_K" variable_2="E_K"/>
</connection>
<connection>
<map_components component_1="reversal_potentials" component_2="rapid_delayed_rectifier_potassium_current"/>
<map_variables variable_1="E_K" variable_2="E_K"/>
</connection>
<connection>
<map_components component_1="reversal_potentials" component_2="slow_delayed_rectifier_potassium_current"/>
<map_variables variable_1="E_Ks" variable_2="E_Ks"/>
</connection>
<connection>
<map_components component_1="reversal_potentials" component_2="sodium_activated_potassium_current"/>
<map_variables variable_1="E_K" variable_2="E_K"/>
</connection>
<connection>
<map_components component_1="reversal_potentials" component_2="transient_outward_current"/>
<map_variables variable_1="E_K" variable_2="E_K"/>
</connection>
<connection>
<map_components component_1="reversal_potentials" component_2="calcium_background_current"/>
<map_variables variable_1="E_Ca" variable_2="E_Ca"/>
</connection>
<connection>
<map_components component_1="reversal_potentials" component_2="fast_sodium_current"/>
<map_variables variable_1="E_mh" variable_2="E_mh"/>
</connection>
<connection>
<map_components component_1="reversal_potentials" component_2="extracellular_potassium_concentration"/>
<map_variables variable_1="K_o" variable_2="K_o"/>
</connection>
<connection>
<map_components component_1="reversal_potentials" component_2="extracellular_sodium_concentration"/>
<map_variables variable_1="Na_o" variable_2="Na_o"/>
</connection>
<connection>
<map_components component_1="reversal_potentials" component_2="intracellular_potassium_concentration"/>
<map_variables variable_1="K_i" variable_2="K_i"/>
</connection>
<connection>
<map_components component_1="reversal_potentials" component_2="intracellular_sodium_concentration"/>
<map_variables variable_1="Na_i" variable_2="Na_i"/>
</connection>
<connection>
<map_components component_1="reversal_potentials" component_2="extracellular_calcium_concentration"/>
<map_variables variable_1="Ca_o" variable_2="Ca_o"/>
</connection>
<connection>
<map_components component_1="reversal_potentials" component_2="intracellular_calcium_concentration"/>
<map_variables variable_1="Ca_i" variable_2="Ca_i"/>
</connection>
<connection>
<map_components component_1="time_independent_potassium_current" component_2="intracellular_potassium_concentration"/>
<map_variables variable_1="i_K1" variable_2="i_K1"/>
</connection>
<connection>
<map_components component_1="rapid_delayed_rectifier_potassium_current" component_2="intracellular_potassium_concentration"/>
<map_variables variable_1="i_Kr" variable_2="i_Kr"/>
</connection>
<connection>
<map_components component_1="slow_delayed_rectifier_potassium_current" component_2="intracellular_potassium_concentration"/>
<map_variables variable_1="i_Ks" variable_2="i_Ks"/>
</connection>
<connection>
<map_components component_1="fast_sodium_current" component_2="intracellular_sodium_concentration"/>
<map_variables variable_1="i_Na" variable_2="i_Na"/>
</connection>
<connection>
<map_components component_1="sodium_background_current" component_2="intracellular_sodium_concentration"/>
<map_variables variable_1="i_b_Na" variable_2="i_b_Na"/>
<map_variables variable_1="Na_i" variable_2="Na_i"/>
</connection>
<connection>
<map_components component_1="persistent_sodium_current" component_2="intracellular_sodium_concentration"/>
<map_variables variable_1="i_p_Na" variable_2="i_p_Na"/>
</connection>
<connection>
<map_components component_1="L_type_Ca_channel" component_2="intracellular_sodium_concentration"/>
<map_variables variable_1="i_Ca_L_Na_cyt" variable_2="i_Ca_L_Na_cyt"/>
<map_variables variable_1="i_Ca_L_Na_ds" variable_2="i_Ca_L_Na_ds"/>
<map_variables variable_1="Na_i" variable_2="Na_i"/>
</connection>
<connection>
<map_components component_1="L_type_Ca_channel" component_2="intracellular_calcium_concentration"/>
<map_variables variable_1="i_Ca_L_Ca_cyt" variable_2="i_Ca_L_Ca_cyt"/>
<map_variables variable_1="i_Ca_L_Ca_ds" variable_2="i_Ca_L_Ca_ds"/>
<map_variables variable_1="Ca_i" variable_2="Ca_i"/>
<map_variables variable_1="Ca_ds" variable_2="Ca_ds"/>
</connection>
<connection>
<map_components component_1="L_type_Ca_channel" component_2="intracellular_potassium_concentration"/>
<map_variables variable_1="i_Ca_L_K_cyt" variable_2="i_Ca_L_K_cyt"/>
<map_variables variable_1="i_Ca_L_K_ds" variable_2="i_Ca_L_K_ds"/>
<map_variables variable_1="K_i" variable_2="K_i"/>
</connection>
<connection>
<map_components component_1="L_type_Ca_channel" component_2="extracellular_sodium_concentration"/>
<map_variables variable_1="Na_o" variable_2="Na_o"/>
</connection>
<connection>
<map_components component_1="L_type_Ca_channel" component_2="extracellular_calcium_concentration"/>
<map_variables variable_1="Ca_o" variable_2="Ca_o"/>
</connection>
<connection>
<map_components component_1="calcium_background_current" component_2="intracellular_calcium_concentration"/>
<map_variables variable_1="i_b_Ca" variable_2="i_b_Ca"/>
</connection>
<connection>
<map_components component_1="sodium_activated_potassium_current" component_2="intracellular_sodium_concentration"/>
<map_variables variable_1="Na_i" variable_2="Na_i"/>
</connection>
<connection>
<map_components component_1="sodium_potassium_pump" component_2="intracellular_sodium_concentration"/>
<map_variables variable_1="i_NaK" variable_2="i_NaK"/>
<map_variables variable_1="Na_i" variable_2="Na_i"/>
</connection>
<connection>
<map_components component_1="sodium_potassium_pump" component_2="intracellular_potassium_concentration"/>
<map_variables variable_1="i_NaK" variable_2="i_NaK"/>
</connection>
<connection>
<map_components component_1="transient_outward_current" component_2="intracellular_potassium_concentration"/>
<map_variables variable_1="i_to" variable_2="i_to"/>
</connection>
<connection>
<map_components component_1="sodium_calcium_exchanger" component_2="intracellular_calcium_concentration"/>
<map_variables variable_1="i_NaCa_cyt" variable_2="i_NaCa_cyt"/>
<map_variables variable_1="i_NaCa_ds" variable_2="i_NaCa_ds"/>
<map_variables variable_1="Ca_i" variable_2="Ca_i"/>
<map_variables variable_1="Ca_ds" variable_2="Ca_ds"/>
<map_variables variable_1="n_NaCa" variable_2="n_NaCa"/>
</connection>
<connection>
<map_components component_1="sodium_calcium_exchanger" component_2="intracellular_sodium_concentration"/>
<map_variables variable_1="i_NaCa_cyt" variable_2="i_NaCa_cyt"/>
<map_variables variable_1="i_NaCa" variable_2="i_NaCa"/>
<map_variables variable_1="n_NaCa" variable_2="n_NaCa"/>
<map_variables variable_1="Na_i" variable_2="Na_i"/>
</connection>
<connection>
<map_components component_1="sodium_calcium_exchanger" component_2="extracellular_sodium_concentration"/>
<map_variables variable_1="Na_o" variable_2="Na_o"/>
</connection>
<connection>
<map_components component_1="sodium_calcium_exchanger" component_2="extracellular_calcium_concentration"/>
<map_variables variable_1="Ca_o" variable_2="Ca_o"/>
</connection>
<connection>
<map_components component_1="sarcoplasmic_reticulum_calcium_pump" component_2="intracellular_calcium_concentration"/>
<map_variables variable_1="i_up" variable_2="i_up"/>
<map_variables variable_1="Ca_i" variable_2="Ca_i"/>
<map_variables variable_1="Ca_up" variable_2="Ca_up"/>
</connection>
<connection>
<map_components component_1="calcium_translocation" component_2="intracellular_calcium_concentration"/>
<map_variables variable_1="i_trans" variable_2="i_trans"/>
<map_variables variable_1="Ca_rel" variable_2="Ca_rel"/>
<map_variables variable_1="Ca_up" variable_2="Ca_up"/>
<map_variables variable_1="CaS_tot" variable_2="CaS_tot"/>
<map_variables variable_1="beta" variable_2="beta"/>
</connection>
<connection>
<map_components component_1="calcium_release" component_2="intracellular_calcium_concentration"/>
<map_variables variable_1="i_rel" variable_2="i_rel"/>
<map_variables variable_1="Ca_rel" variable_2="Ca_rel"/>
<map_variables variable_1="Ca_i" variable_2="Ca_i"/>
<map_variables variable_1="Ca_ds" variable_2="Ca_ds"/>
</connection>
<connection>
<map_components component_1="rapid_delayed_rectifier_potassium_current" component_2="rapid_delayed_rectifier_potassium_current_xr1_gate"/>
<map_variables variable_1="xr1" variable_2="xr1"/>
<map_variables variable_1="time" variable_2="time"/>
<map_variables variable_1="V" variable_2="V"/>
</connection>
<connection>
<map_components component_1="rapid_delayed_rectifier_potassium_current" component_2="rapid_delayed_rectifier_potassium_current_xr2_gate"/>
<map_variables variable_1="xr2" variable_2="xr2"/>
<map_variables variable_1="time" variable_2="time"/>
<map_variables variable_1="V" variable_2="V"/>
</connection>
<connection>
<map_components component_1="slow_delayed_rectifier_potassium_current" component_2="slow_delayed_rectifier_potassium_current_xs_gate"/>
<map_variables variable_1="xs" variable_2="xs"/>
<map_variables variable_1="time" variable_2="time"/>
<map_variables variable_1="V" variable_2="V"/>
</connection>
<connection>
<map_components component_1="fast_sodium_current" component_2="fast_sodium_current_m_gate"/>
<map_variables variable_1="m" variable_2="m"/>
<map_variables variable_1="time" variable_2="time"/>
<map_variables variable_1="V" variable_2="V"/>
</connection>
<connection>
<map_components component_1="fast_sodium_current" component_2="fast_sodium_current_h_gate"/>
<map_variables variable_1="h" variable_2="h"/>
<map_variables variable_1="time" variable_2="time"/>
<map_variables variable_1="V" variable_2="V"/>
</connection>
<connection>
<map_components component_1="L_type_Ca_channel" component_2="L_type_Ca_channel_d_gate"/>
<map_variables variable_1="d" variable_2="d"/>
<map_variables variable_1="time" variable_2="time"/>
<map_variables variable_1="V" variable_2="V"/>
</connection>
<connection>
<map_components component_1="L_type_Ca_channel" component_2="L_type_Ca_channel_f_gate"/>
<map_variables variable_1="f" variable_2="f"/>
<map_variables variable_1="time" variable_2="time"/>
<map_variables variable_1="V" variable_2="V"/>
</connection>
<connection>
<map_components component_1="L_type_Ca_channel" component_2="L_type_Ca_channel_f2_gate"/>
<map_variables variable_1="f2" variable_2="f2"/>
<map_variables variable_1="time" variable_2="time"/>
<map_variables variable_1="V" variable_2="V"/>
<map_variables variable_1="Ca_i" variable_2="Ca_i"/>
<map_variables variable_1="Km_f2" variable_2="Km_f2"/>
</connection>
<connection>
<map_components component_1="L_type_Ca_channel" component_2="L_type_Ca_channel_f2ds_gate"/>
<map_variables variable_1="f2ds" variable_2="f2ds"/>
<map_variables variable_1="time" variable_2="time"/>
<map_variables variable_1="V" variable_2="V"/>
<map_variables variable_1="Ca_ds" variable_2="Ca_ds"/>
<map_variables variable_1="Km_f2ds" variable_2="Km_f2ds"/>
<map_variables variable_1="R_decay" variable_2="R_decay"/>
</connection>
<connection>
<map_components component_1="transient_outward_current" component_2="transient_outward_current_s_gate"/>
<map_variables variable_1="s" variable_2="s"/>
<map_variables variable_1="time" variable_2="time"/>
<map_variables variable_1="V" variable_2="V"/>
</connection>
<connection>
<map_components component_1="transient_outward_current" component_2="transient_outward_current_r_gate"/>
<map_variables variable_1="r" variable_2="r"/>
<map_variables variable_1="time" variable_2="time"/>
<map_variables variable_1="V" variable_2="V"/>
</connection>
<connection>
<map_components component_1="intracellular_sodium_concentration" component_2="intracellular_calcium_concentration"/>
<map_variables variable_1="V_i" variable_2="V_i"/>
</connection>
<connection>
<map_components component_1="intracellular_potassium_concentration" component_2="intracellular_calcium_concentration"/>
<map_variables variable_1="V_i" variable_2="V_i"/>
</connection>
<connection>
<map_components component_1="parameters" component_2="isotonic"/>
<map_variables variable_1="isotonic" variable_2="isotonic"/>
</connection>
<connection>
<map_components component_1="force" component_2="isotonic"/>
<map_variables variable_1="F_muscle" variable_2="F_muscle"/>
</connection>
<connection>
<map_components component_1="length" component_2="isotonic"/>
<map_variables variable_1="isotonic_mode" variable_2="isotonic_mode"/>
<map_variables variable_1="l" variable_2="l"/>
<map_variables variable_1="l_0" variable_2="l_0"/>
</connection>
<connection>
<map_components component_1="PE_velocity" component_2="isotonic"/>
<map_variables variable_1="isotonic_mode" variable_2="isotonic_mode"/>
</connection>
<connection>
<map_components component_1="CE_velocity" component_2="isotonic"/>
<map_variables variable_1="isotonic_mode" variable_2="isotonic_mode"/>
</connection>
<connection>
<map_components component_1="membrane" component_2="isotonic"/>
<map_variables variable_1="stim_period" variable_2="stim_period"/>
</connection>
<connection>
<map_components component_1="environment" component_2="crossbridge_kinetics"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection>
<map_components component_1="environment" component_2="length"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection>
<map_components component_1="environment" component_2="CE_velocity"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection>
<map_components component_1="environment" component_2="PE_velocity"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection>
<map_components component_1="parameters" component_2="parameters_izakov_et_al_1991"/>
<map_variables variable_1="v_max" variable_2="v_max"/>
</connection>
<connection>
<map_components component_1="force" component_2="parameters"/>
<map_variables variable_1="alpha_1" variable_2="alpha_1"/>
<map_variables variable_1="beta_1" variable_2="beta_1"/>
<map_variables variable_1="alpha_2" variable_2="alpha_2"/>
<map_variables variable_1="beta_2" variable_2="beta_2"/>
<map_variables variable_1="alpha_3" variable_2="alpha_3"/>
<map_variables variable_1="beta_3" variable_2="beta_3"/>
<map_variables variable_1="lambda" variable_2="lambda"/>
</connection>
<connection>
<map_components component_1="force" component_2="crossbridge_kinetics"/>
<map_variables variable_1="N" variable_2="N"/>
</connection>
<connection>
<map_components component_1="force" component_2="length"/>
<map_variables variable_1="l_1" variable_2="l_1"/>
<map_variables variable_1="l_2" variable_2="l_2"/>
<map_variables variable_1="l_3" variable_2="l_3"/>
</connection>
<connection>
<map_components component_1="force" component_2="CE_velocity"/>
<map_variables variable_1="v" variable_2="v"/>
<map_variables variable_1="k_P_vis" variable_2="k_P_vis"/>
</connection>
<connection>
<map_components component_1="force" component_2="PE_velocity"/>
<map_variables variable_1="w" variable_2="w"/>
<map_variables variable_1="k_S_vis" variable_2="k_S_vis"/>
</connection>
<connection>
<map_components component_1="force" component_2="average_crossbridge_force"/>
<map_variables variable_1="p_v" variable_2="p_v"/>
</connection>
<connection>
<map_components component_1="crossbridge_kinetics" component_2="parameters"/>
<map_variables variable_1="mu" variable_2="mu"/>
<map_variables variable_1="k_mu" variable_2="k_mu"/>
<map_variables variable_1="kappa" variable_2="kappa"/>
<map_variables variable_1="kappa_0" variable_2="kappa_0"/>
<map_variables variable_1="m_0" variable_2="m_0"/>
</connection>
<connection>
<map_components component_1="crossbridge_kinetics" component_2="intracellular_calcium_concentration"/>
<map_variables variable_1="A" variable_2="A"/>
<map_variables variable_1="A_tot" variable_2="A_tot"/>
<map_variables variable_1="N" variable_2="N"/>
<map_variables variable_1="L_oz" variable_2="L_oz"/>
</connection>
<connection>
<map_components component_1="crossbridge_kinetics" component_2="average_crossbridge_force"/>
<map_variables variable_1="G_star" variable_2="G_star"/>
<map_variables variable_1="P_star" variable_2="P_star"/>
</connection>
<connection>
<map_components component_1="crossbridge_kinetics" component_2="length"/>
<map_variables variable_1="l_1" variable_2="l_1"/>
</connection>
<connection>
<map_components component_1="crossbridge_kinetics" component_2="CE_velocity"/>
<map_variables variable_1="K_kappa" variable_2="K_kappa"/>
<map_variables variable_1="v" variable_2="v"/>
<map_variables variable_1="N" variable_2="N"/>
</connection>
<connection>
<map_components component_1="crossbridge_kinetics" component_2="parameters_izakov_et_al_1991"/>
<map_variables variable_1="S_0" variable_2="S_0"/>
<map_variables variable_1="q_v" variable_2="q_v"/>
<map_variables variable_1="v_st" variable_2="v_st"/>
<map_variables variable_1="q_4" variable_2="q_4"/>
</connection>
<connection>
<map_components component_1="length" component_2="CE_velocity"/>
<map_variables variable_1="v" variable_2="v"/>
<map_variables variable_1="phi_chi_2" variable_2="phi_chi_2"/>
<map_variables variable_1="dl_1_dt" variable_2="dl_1_dt"/>
<map_variables variable_1="l_1" variable_2="l_1"/>
<map_variables variable_1="l_2" variable_2="l_2"/>
<map_variables variable_1="l_3" variable_2="l_3"/>
</connection>
<connection>
<map_components component_1="length" component_2="PE_velocity"/>
<map_variables variable_1="w" variable_2="w"/>
<map_variables variable_1="k_S_vis" variable_2="k_S_vis"/>
<map_variables variable_1="dl_1_dt" variable_2="dl_1_dt"/>
<map_variables variable_1="l_1" variable_2="l_1"/>
<map_variables variable_1="dl_2_dt" variable_2="dl_2_dt"/>
<map_variables variable_1="l_2" variable_2="l_2"/>
<map_variables variable_1="l_3" variable_2="l_3"/>
</connection>
<connection>
<map_components component_1="length" component_2="membrane"/>
<map_variables variable_1="stim_start" variable_2="stim_start"/>
<map_variables variable_1="stim_end" variable_2="stim_end"/>
<map_variables variable_1="stim_period" variable_2="stim_period"/>
<map_variables variable_1="stim_duration" variable_2="stim_duration"/>
<map_variables variable_1="stim_amplitude" variable_2="stim_amplitude"/>
</connection>
<connection>
<map_components component_1="CE_velocity" component_2="parameters"/>
<map_variables variable_1="alpha_1" variable_2="alpha_1"/>
<map_variables variable_1="beta_1" variable_2="beta_1"/>
<map_variables variable_1="alpha_2" variable_2="alpha_2"/>
<map_variables variable_1="beta_2" variable_2="beta_2"/>
<map_variables variable_1="alpha_3" variable_2="alpha_3"/>
<map_variables variable_1="beta_3" variable_2="beta_3"/>
<map_variables variable_1="lambda" variable_2="lambda"/>
</connection>
<connection>
<map_components component_1="CE_velocity" component_2="parameters_izakov_et_al_1991"/>
<map_variables variable_1="v" variable_2="v"/>
</connection>
<connection>
<map_components component_1="CE_velocity" component_2="PE_velocity"/>
<map_variables variable_1="v" variable_2="v"/>
<map_variables variable_1="w" variable_2="w"/>
<map_variables variable_1="k_S_vis" variable_2="k_S_vis"/>
<map_variables variable_1="phi_chi" variable_2="phi_chi"/>
</connection>
<connection>
<map_components component_1="PE_velocity" component_2="parameters"/>
<map_variables variable_1="alpha_1" variable_2="alpha_1"/>
<map_variables variable_1="beta_1" variable_2="beta_1"/>
<map_variables variable_1="alpha_2" variable_2="alpha_2"/>
<map_variables variable_1="beta_2" variable_2="beta_2"/>
<map_variables variable_1="alpha_3" variable_2="alpha_3"/>
<map_variables variable_1="beta_3" variable_2="beta_3"/>
<map_variables variable_1="lambda" variable_2="lambda"/>
</connection>
<connection>
<map_components component_1="average_crossbridge_force" component_2="parameters"/>
<map_variables variable_1="v_max" variable_2="v_max"/>
<map_variables variable_1="a" variable_2="a"/>
<map_variables variable_1="d_h" variable_2="d_h"/>
<map_variables variable_1="alpha_P" variable_2="alpha_P"/>
</connection>
<connection>
<map_components component_1="average_crossbridge_force" component_2="parameters_izakov_et_al_1991"/>
<map_variables variable_1="v_1" variable_2="v_1"/>
<map_variables variable_1="alpha_G" variable_2="alpha_G"/>
</connection>
<connection>
<map_components component_1="average_crossbridge_force" component_2="CE_velocity"/>
<map_variables variable_1="v" variable_2="v"/>
<map_variables variable_1="phi_chi" variable_2="phi_chi"/>
<map_variables variable_1="p_v" variable_2="p_v"/>
<map_variables variable_1="p_prime_v" variable_2="p_prime_v"/>
</connection>
<connection>
<map_components component_1="intracellular_calcium_concentration" component_2="parameters_izakov_et_al_1991"/>
<map_variables variable_1="k_A" variable_2="k_A"/>
</connection>
<connection>
<map_components component_1="extracellular_potassium_concentration" component_2="rapid_delayed_rectifier_potassium_current"/>
<map_variables variable_1="i_Kr" variable_2="i_Kr"/>
</connection>
<connection>
<map_components component_1="extracellular_potassium_concentration" component_2="slow_delayed_rectifier_potassium_current"/>
<map_variables variable_1="i_Ks" variable_2="i_Ks"/>
</connection>
<connection>
<map_components component_1="extracellular_potassium_concentration" component_2="time_independent_potassium_current"/>
<map_variables variable_1="i_K1" variable_2="i_K1"/>
<map_variables variable_1="K_o" variable_2="K_o"/>
</connection>
<connection>
<map_components component_1="extracellular_potassium_concentration" component_2="transient_outward_current"/>
<map_variables variable_1="i_to" variable_2="i_to"/>
</connection>
<connection>
<map_components component_1="extracellular_potassium_concentration" component_2="environment"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection>
<map_components component_1="extracellular_potassium_concentration" component_2="sodium_potassium_pump"/>
<map_variables variable_1="i_NaK" variable_2="i_NaK"/>
<map_variables variable_1="K_o" variable_2="K_o"/>
</connection>
<connection>
<map_components component_1="extracellular_potassium_concentration" component_2="L_type_Ca_channel"/>
<map_variables variable_1="i_Ca_L_K_ds" variable_2="i_Ca_L_K_ds"/>
<map_variables variable_1="i_Ca_L_K_cyt" variable_2="i_Ca_L_K_cyt"/>
<map_variables variable_1="K_o" variable_2="K_o"/>
</connection>
<connection>
<map_components component_1="extracellular_potassium_concentration" component_2="membrane"/>
<map_variables variable_1="F" variable_2="F"/>
</connection>
<connection>
<map_components component_1="extracellular_potassium_concentration" component_2="intracellular_sodium_concentration"/>
<map_variables variable_1="n_NaK" variable_2="n_NaK"/>
</connection>
<connection>
<map_components component_1="intracellular_potassium_concentration" component_2="intracellular_sodium_concentration"/>
<map_variables variable_1="n_NaK" variable_2="n_NaK"/>
</connection>
<connection>
<map_components component_1="extracellular_potassium_concentration" component_2="intracellular_calcium_concentration"/>
<map_variables variable_1="V_e" variable_2="V_e"/>
</connection>
<rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#" xmlns:bqs="http://www.cellml.org/bqs/1.0#" xmlns:dc="http://purl.org/dc/elements/1.1/" xmlns:dcterms="http://purl.org/dc/terms/" xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">
<foo xmlns="http://www.example.org/foo/bar/baz" />
<!--
The following RDF block contains metadata that applies to this document
as a whole, as indicated by the empty about attribute on the
<rdf:Description> element.
-->
<rdf:Description rdf:about="">
<!--
The Model Builder Metadata. The Dublin Core "creator" element is used
to indicate the person who translated the model into CellML.
-->
<dc:creator rdf:parseType="Resource">
<vCard:N rdf:parseType="Resource">
<vCard:Family>Sulman</vCard:Family>
<vCard:Given>Tatiana</vCard:Given>
</vCard:N>
<vCard:ORG rdf:parseType="Resource">
<vCard:Orgname>Ural Division of the Russian Academy of Sciences.</vCard:Orgname>
<vCard:Orgunit>Mathematical Physilogy Lab, Institute of Immunology and Physiology</vCard:Orgunit>
</vCard:ORG>
</dc:creator>
<!--
The Creation Date metadata. This is the date on which the model
was translated into CellML.
<dcterms:created rdf:parseType="Resource">
<dcterms:W3CDTF>2009-09-28</dcterms:W3CDTF>
</dcterms:created>
-->
<!--
The Modification History metadata. This lists the changes that have been
made to the document, who made the changes, and when they were made.
-->
</rdf:Description>
<!--
The following metadata refers to the model itself, as indicated by the
reference to the ID "ekaterinburg_rat", which is
declared on the <model> element.
-->
<rdf:Description rdf:about="#ekaterinburg_rat">
<!-- A human readable name for the model. -->
<dc:title>
Mathematical modeling of mechanically modulated rhythm disturbances in homogeneous and heterogeneous myocardium with attenuated activity of na+-k+ pump
</dc:title>
<!-- Keyword(s) -->
<bqs:reference rdf:parseType="Resource">
<dc:subject rdf:parseType="Resource">
<bqs:subject_type>keyword</bqs:subject_type>
<rdf:value>
<rdf:Bag>
<rdf:li>cardiac</rdf:li>
<rdf:li>electrophysiology</rdf:li>
</rdf:Bag>
</rdf:value>
</dc:subject>
</bqs:reference>
<bqs:reference rdf:parseType="Resource">
<bqs:Pubmed_id>18259823</bqs:Pubmed_id>
<bqs:JournalArticle rdf:parseType="Resource">
<dc:creator>
<rdf:Seq>
<rdf:li rdf:parseType="Resource">
<bqs:Person rdf:parseType="Resource">
<vCard:N rdf:parseType="Resource">
<vCard:Family>Sulman</vCard:Family>
<vCard:Given>Tatiana</vCard:Given>
</vCard:N>
</bqs:Person>
</rdf:li>
<rdf:li rdf:parseType="Resource">
<bqs:Person rdf:parseType="Resource">
<vCard:N rdf:parseType="Resource">
<vCard:Family>Katnelson</vCard:Family>
<vCard:Given>Leonid</vCard:Given>
<vCard:Other>B</vCard:Other>
</vCard:N>
</bqs:Person>
</rdf:li>
<rdf:li rdf:parseType="Resource">
<bqs:Person rdf:parseType="Resource">
<vCard:N rdf:parseType="Resource">
<vCard:Family>Solovyova</vCard:Family>
<vCard:Given>Olga</vCard:Given>
</vCard:N>
</bqs:Person>
</rdf:li>
<rdf:li rdf:parseType="Resource">
<bqs:Person rdf:parseType="Resource">
<vCard:N rdf:parseType="Resource">
<vCard:Family>Markhasin</vCard:Family>
<vCard:Given>Vladimir</vCard:Given>
<vCard:Other>S</vCard:Other>
</vCard:N>
</bqs:Person>
</rdf:li>
</rdf:Seq>
</dc:creator>
<dc:title>
Mathematical modeling of mechanically modulated rhythm disturbances in homogeneous and heterogeneous myocardium with attenuated activity of na+-k+ pump
</dc:title>
<bqs:Journal rdf:parseType="Resource">
<dc:title>Bulletin of Mathematical Biology</dc:title>
</bqs:Journal>
<bqs:volume>70</bqs:volume>
<bqs:first_page>910</bqs:first_page>
<bqs:last_page>949</bqs:last_page>
</bqs:JournalArticle>
</bqs:reference>
</rdf:Description>
</rdf:RDF>
</model>
