Location: Cardiovascular model with the feedback loops @ ed61dccc717c / models / CellMLV2 / Gee_whiz_flat.cellml

Author:
Weiwei Ai <wai484@aucklanduni.ac.nz>
Date:
2024-07-17 10:48:12+12:00
Desc:
Working version from Alan (removed the connection duplications)
Permanent Source URI:
https://models.fieldml.org/workspace/b6e/rawfile/ed61dccc717c423c88ae5090bb3cb5759e57372f/models/CellMLV2/Gee_whiz_flat.cellml

<?xml version="1.0" encoding="UTF-8"?>
<model xmlns="http://www.cellml.org/cellml/2.0#" name="Gee_whiz">
  <units name="mmHg">
    <unit multiplier="133.32" units="pascal"/>
  </units>
  <units name="millimolar">
    <unit prefix="milli" units="mole"/>
    <unit exponent="-1" units="litre"/>
  </units>
  <units name="per_mmHg">
    <unit exponent="-1" units="mmHg"/>
  </units>
  <units name="millisecond">
    <unit prefix="milli" units="second"/>
  </units>
  <units name="per_millimolar">
    <unit exponent="-1" units="millimolar"/>
  </units>
  <units name="per_s_per_mM">
    <unit exponent="-1" units="second"/>
    <unit units="per_millimolar"/>
  </units>
  <units name="per_s">
    <unit exponent="-1" units="second"/>
  </units>
  <units name="mL">
    <unit prefix="milli" units="litre"/>
  </units>
  <units name="mmHg_per_s">
    <unit units="mmHg"/>
    <unit exponent="-1" units="second"/>
  </units>
  <units name="mL_per_s">
    <unit units="mL"/>
    <unit exponent="-1" units="second"/>
  </units>
  <units name="mL_per_mmHg">
    <unit units="mL"/>
    <unit exponent="-1" units="mmHg"/>
  </units>
  <units name="mmHg_s_per_mL">
    <unit units="mmHg"/>
    <unit exponent="-1" units="mL"/>
    <unit units="second"/>
  </units>
  <units name="mmHg_s2_per_mL">
    <unit units="mmHg"/>
    <unit exponent="-1" units="mL"/>
    <unit exponent="2" units="second"/>
  </units>
  <units name="mL_per_s_1">
    <unit units="mL"/>
    <unit exponent="-1" units="second"/>
  </units>
  <units name="s_per_mL">
    <unit exponent="-1" units="mL_per_s_1"/>
  </units>
  <units name="per_mL">
    <unit exponent="-1" units="mL"/>
  </units>
  <units name="mmHg_per_mL">
    <unit units="mmHg"/>
    <unit exponent="-1" units="mL"/>
  </units>
  <units name="mmHg_per_s_1">
    <unit units="mmHg"/>
    <unit exponent="-1" units="second"/>
  </units>
  <units name="mL_per_s_2">
    <unit units="mL"/>
    <unit exponent="-1" units="second"/>
  </units>
  <units name="mL_per_mmHg_1">
    <unit units="mL"/>
    <unit exponent="-1" units="mmHg"/>
  </units>
  <units name="mmHg_s_per_mL_1">
    <unit units="mmHg"/>
    <unit exponent="-1" units="mL"/>
    <unit units="second"/>
  </units>
  <units name="mmHg_s2_per_mL_1">
    <unit units="mmHg"/>
    <unit exponent="-1" units="mL"/>
    <unit exponent="2" units="second"/>
  </units>
  <units name="s_per_mL_1">
    <unit exponent="-1" units="mL_per_s_2"/>
  </units>
  <units name="per_mL_1">
    <unit exponent="-1" units="mL"/>
  </units>
  <units name="mmHg_per_mL_1">
    <unit units="mmHg"/>
    <unit exponent="-1" units="mL"/>
  </units>
  <units name="s2">
    <unit exponent="2" units="second"/>
  </units>
  <units name="millivolt">
    <unit prefix="milli" units="volt"/>
  </units>
  <units name="joule_per_mole_kelvin">
    <unit units="joule"/>
    <unit exponent="-1" units="mole"/>
    <unit exponent="-1" units="kelvin"/>
  </units>
  <units name="coulomb_per_millimole">
    <unit units="coulomb"/>
    <unit exponent="-1" prefix="milli" units="mole"/>
  </units>
  <units name="picoF">
    <unit prefix="pico" units="farad"/>
  </units>
  <units name="picoA_per_picoF">
    <unit prefix="pico" units="ampere"/>
    <unit exponent="-1" prefix="pico" units="farad"/>
  </units>
  <units name="nanoS_per_picoF">
    <unit prefix="nano" units="siemens"/>
    <unit exponent="-1" prefix="pico" units="farad"/>
  </units>
  <units name="per_millisecond">
    <unit exponent="-1" prefix="milli" units="second"/>
  </units>
  <units name="per_millivolt">
    <unit exponent="-1" prefix="milli" units="volt"/>
  </units>
  <units name="litre_per_farad_second">
    <unit units="litre"/>
    <unit exponent="-1" units="farad"/>
    <unit exponent="-1" units="second"/>
  </units>
  <units name="millivolt_square">
    <unit exponent="2" prefix="milli" units="volt"/>
  </units>
  <units name="millimolar_per_millisecond">
    <unit units="millimolar"/>
    <unit exponent="-1" prefix="milli" units="second"/>
  </units>
  <component name="cvs">
    <variable name="t" units="second" interface="public"/>
    <variable name="phi" units="dimensionless" interface="public"/>
    <variable name="q_ra" units="mL" initial_value="4"/>
    <variable name="q_rv" units="mL" initial_value="5"/>
    <variable name="q_la" units="mL" initial_value="4"/>
    <variable name="q_lv" units="mL" initial_value="10"/>
    <variable name="q_pa" units="mL" initial_value="4"/>
    <variable name="q_pc" units="mL" initial_value="5"/>
    <variable name="q_pv" units="mL" initial_value="4"/>
    <variable name="q_sa" units="mL" initial_value="5"/>
    <variable name="q_sc" units="mL" initial_value="10"/>
    <variable name="q_sv" units="mL"/>
    <variable name="q_tot" units="mL" initial_value="5800"/>
    <variable name="q_tot_u" units="mL" initial_value="3500"/>
    <variable name="u_ra" units="mmHg"/>
    <variable name="u_rv" units="mmHg" interface="public"/>
    <variable name="u_la" units="mmHg"/>
    <variable name="u_lv" units="mmHg" interface="public"/>
    <variable name="u_pa" units="mmHg"/>
    <variable name="u_pc" units="mmHg"/>
    <variable name="u_pv" units="mmHg" interface="public"/>
    <variable name="u_sa" units="mmHg" interface="public"/>
    <variable name="du_sa" units="mmHg_per_s" interface="public"/>
    <variable name="u_sc" units="mmHg"/>
    <variable name="u_sv" units="mmHg"/>
    <variable name="v_TV" units="mL_per_s_1"/>
    <variable name="v_PV" units="mL_per_s_1"/>
    <variable name="v_MV" units="mL_per_s_1"/>
    <variable name="v_AV" units="mL_per_s_1"/>
    <variable name="v_pa" units="mL_per_s_1" initial_value="5"/>
    <variable name="v_pv" units="mL_per_s_1"/>
    <variable name="v_la" units="mL_per_s_1"/>
    <variable name="v_sa" units="mL_per_s_1" initial_value="5"/>
    <variable name="v_sc" units="mL_per_s_1"/>
    <variable name="v_sv" units="mL_per_s_1"/>
    <variable name="C_ra" units="mL_per_mmHg" initial_value="31.25"/>
    <variable name="C_pa" units="mL_per_mmHg" initial_value="0.76"/>
    <variable name="C_pc" units="mL_per_mmHg" initial_value="5.8"/>
    <variable name="C_pv" units="mL_per_mmHg" initial_value="25.37"/>
    <variable name="C_la" units="mL_per_mmHg" initial_value="19.23"/>
    <variable name="C_sa" units="mL_per_mmHg" initial_value="0.28"/>
    <variable name="C_sc" units="mL_per_mmHg" initial_value="3.72"/>
    <variable name="C_sv" units="mL_per_mmHg" initial_value="111.11"/>
    <variable name="R_TV" units="mmHg_s_per_mL" initial_value="2.5e-3"/>
    <variable name="R_PV" units="mmHg_s_per_mL"/>
    <variable name="R_MV" units="mmHg_s_per_mL" initial_value="2.5e-3"/>
    <variable name="R_AV" units="mmHg_s_per_mL"/>
    <variable name="R_pa" units="mmHg_s_per_mL" initial_value="0.023"/>
    <variable name="R_pv" units="mmHg_s_per_mL" initial_value="0.0894"/>
    <variable name="R_la" units="mmHg_s_per_mL" initial_value="0.0056"/>
    <variable name="R_sv" units="mmHg_s_per_mL" initial_value="0.0054"/>
    <variable name="R_sc" units="mmHg_s_per_mL" initial_value="0.8895"/>
    <variable name="R_sa" units="mmHg_s_per_mL" initial_value="0.06"/>
    <variable name="L_pa" units="mmHg_s2_per_mL" initial_value="1.8e-4"/>
    <variable name="L_sa" units="mmHg_s2_per_mL" initial_value="2.2e-4"/>
    <variable name="k_PV_R" units="s_per_mL" initial_value="1.4e-3"/>
    <variable name="k_AV_R" units="s_per_mL" initial_value="3.75e-4"/>
    <variable name="k_lv_E" units="per_mL" initial_value="0.014"/>
    <variable name="k_rv_E" units="per_mL" initial_value="0.011"/>
    <variable name="q_lv_u" units="mL" initial_value="16.77"/>
    <variable name="q_rv_u" units="mL" initial_value="40.8"/>
    <variable name="u_lv_0" units="mmHg" initial_value="1.5"/>
    <variable name="u_rv_0" units="mmHg" initial_value="1.5"/>
    <variable name="E_lv_max" units="mmHg_per_mL" initial_value="1.283"/>
    <variable name="E_rv_max" units="mmHg_per_mL" initial_value="0.757"/>
    <variable name="u_rv_max" units="mmHg"/>
    <variable name="u_lv_max" units="mmHg"/>
    <math xmlns="http://www.w3.org/1998/Math/MathML" xmlns:cellml="http://www.cellml.org/cellml/2.0#">
      <apply>
        <eq/>
        <ci>u_lv_max</ci>
        <apply>
          <plus/>
          <apply>
            <times/>
            <ci>phi</ci>
            <ci>E_lv_max</ci>
            <apply>
              <minus/>
              <ci>q_lv</ci>
              <ci>q_lv_u</ci>
            </apply>
          </apply>
          <apply>
            <times/>
            <apply>
              <minus/>
              <cn cellml:units="dimensionless">1</cn>
              <ci>phi</ci>
            </apply>
            <ci>u_lv_0</ci>
            <apply>
              <minus/>
              <apply>
                <exp/>
                <apply>
                  <times/>
                  <ci>k_lv_E</ci>
                  <ci>q_lv</ci>
                </apply>
              </apply>
              <cn cellml:units="dimensionless">1</cn>
            </apply>
          </apply>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>u_rv_max</ci>
        <apply>
          <plus/>
          <apply>
            <times/>
            <ci>phi</ci>
            <ci>E_rv_max</ci>
            <apply>
              <minus/>
              <ci>q_rv</ci>
              <ci>q_rv_u</ci>
            </apply>
          </apply>
          <apply>
            <times/>
            <apply>
              <minus/>
              <cn cellml:units="dimensionless">1</cn>
              <ci>phi</ci>
            </apply>
            <ci>u_rv_0</ci>
            <apply>
              <minus/>
              <apply>
                <exp/>
                <apply>
                  <times/>
                  <ci>k_rv_E</ci>
                  <ci>q_rv</ci>
                </apply>
              </apply>
              <cn cellml:units="dimensionless">1</cn>
            </apply>
          </apply>
        </apply>
      </apply>
      <apply>
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci>t</ci>
          </bvar>
          <ci>q_ra</ci>
        </apply>
        <apply>
          <minus/>
          <ci>v_sv</ci>
          <ci>v_TV</ci>
        </apply>
      </apply>
      <apply>
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci>t</ci>
          </bvar>
          <ci>q_rv</ci>
        </apply>
        <apply>
          <minus/>
          <ci>v_TV</ci>
          <ci>v_PV</ci>
        </apply>
      </apply>
      <apply>
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci>t</ci>
          </bvar>
          <ci>q_pa</ci>
        </apply>
        <apply>
          <minus/>
          <ci>v_PV</ci>
          <ci>v_pa</ci>
        </apply>
      </apply>
      <apply>
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci>t</ci>
          </bvar>
          <ci>q_pc</ci>
        </apply>
        <apply>
          <minus/>
          <ci>v_pa</ci>
          <ci>v_pv</ci>
        </apply>
      </apply>
      <apply>
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci>t</ci>
          </bvar>
          <ci>q_pv</ci>
        </apply>
        <apply>
          <minus/>
          <ci>v_pv</ci>
          <ci>v_la</ci>
        </apply>
      </apply>
      <apply>
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci>t</ci>
          </bvar>
          <ci>q_la</ci>
        </apply>
        <apply>
          <minus/>
          <ci>v_la</ci>
          <ci>v_MV</ci>
        </apply>
      </apply>
      <apply>
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci>t</ci>
          </bvar>
          <ci>q_lv</ci>
        </apply>
        <apply>
          <minus/>
          <ci>v_MV</ci>
          <ci>v_AV</ci>
        </apply>
      </apply>
      <apply>
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci>t</ci>
          </bvar>
          <ci>q_sa</ci>
        </apply>
        <apply>
          <minus/>
          <ci>v_AV</ci>
          <ci>v_sa</ci>
        </apply>
      </apply>
      <apply>
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci>t</ci>
          </bvar>
          <ci>q_sc</ci>
        </apply>
        <apply>
          <minus/>
          <ci>v_sa</ci>
          <ci>v_sc</ci>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>q_sv</ci>
        <apply>
          <minus/>
          <apply>
            <minus/>
            <ci>q_tot</ci>
            <ci>q_tot_u</ci>
          </apply>
          <apply>
            <plus/>
            <ci>q_ra</ci>
            <ci>q_rv</ci>
            <ci>q_pa</ci>
            <ci>q_pc</ci>
            <ci>q_pv</ci>
            <ci>q_la</ci>
            <ci>q_lv</ci>
            <ci>q_sa</ci>
            <ci>q_sc</ci>
          </apply>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>u_ra</ci>
        <apply>
          <divide/>
          <ci>q_ra</ci>
          <ci>C_ra</ci>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>u_rv</ci>
        <apply>
          <minus/>
          <ci>u_rv_max</ci>
          <apply>
            <times/>
            <ci>R_PV</ci>
            <ci>v_PV</ci>
          </apply>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>u_pa</ci>
        <apply>
          <divide/>
          <ci>q_pa</ci>
          <ci>C_pa</ci>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>u_pc</ci>
        <apply>
          <divide/>
          <ci>q_pc</ci>
          <ci>C_pc</ci>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>u_pv</ci>
        <apply>
          <divide/>
          <ci>q_pv</ci>
          <ci>C_pv</ci>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>u_la</ci>
        <apply>
          <divide/>
          <ci>q_la</ci>
          <ci>C_la</ci>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>u_lv</ci>
        <apply>
          <minus/>
          <ci>u_lv_max</ci>
          <apply>
            <times/>
            <ci>R_AV</ci>
            <ci>v_AV</ci>
          </apply>
        </apply>
      </apply>
      <!--       ode(u_sa1,t) = du_sa;-->
      <apply>
        <eq/>
        <ci>u_sa</ci>
        <apply>
          <divide/>
          <ci>q_sa</ci>
          <ci>C_sa</ci>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>u_sc</ci>
        <apply>
          <divide/>
          <ci>q_sc</ci>
          <ci>C_sc</ci>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>u_sv</ci>
        <apply>
          <divide/>
          <ci>q_sv</ci>
          <ci>C_sv</ci>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>du_sa</ci>
        <apply>
          <divide/>
          <apply>
            <minus/>
            <ci>v_AV</ci>
            <ci>v_sa</ci>
          </apply>
          <ci>C_sa</ci>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>v_TV</ci>
        <piecewise>
          <piece>
            <cn cellml:units="mL_per_s_1">0</cn>
            <apply>
              <lt/>
              <ci>u_ra</ci>
              <ci>u_rv</ci>
            </apply>
          </piece>
          <otherwise>
            <apply>
              <divide/>
              <apply>
                <minus/>
                <ci>u_ra</ci>
                <ci>u_rv</ci>
              </apply>
              <ci>R_TV</ci>
            </apply>
          </otherwise>
        </piecewise>
      </apply>
      <apply>
        <eq/>
        <ci>v_PV</ci>
        <piecewise>
          <piece>
            <cn cellml:units="mL_per_s_1">0</cn>
            <apply>
              <lt/>
              <ci>u_rv_max</ci>
              <ci>u_pa</ci>
            </apply>
          </piece>
          <otherwise>
            <apply>
              <divide/>
              <apply>
                <minus/>
                <ci>u_rv_max</ci>
                <ci>u_pa</ci>
              </apply>
              <ci>R_PV</ci>
            </apply>
          </otherwise>
        </piecewise>
      </apply>
      <apply>
        <eq/>
        <ci>R_PV</ci>
        <apply>
          <times/>
          <ci>k_PV_R</ci>
          <ci>u_rv_max</ci>
        </apply>
      </apply>
      <apply>
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci>t</ci>
          </bvar>
          <ci>v_pa</ci>
        </apply>
        <apply>
          <divide/>
          <apply>
            <minus/>
            <apply>
              <minus/>
              <ci>u_pa</ci>
              <ci>u_pc</ci>
            </apply>
            <apply>
              <times/>
              <ci>R_pa</ci>
              <ci>v_pa</ci>
            </apply>
          </apply>
          <ci>L_pa</ci>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>v_pv</ci>
        <apply>
          <divide/>
          <apply>
            <minus/>
            <ci>u_pc</ci>
            <ci>u_pv</ci>
          </apply>
          <ci>R_pv</ci>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>v_la</ci>
        <apply>
          <divide/>
          <apply>
            <minus/>
            <ci>u_pv</ci>
            <ci>u_la</ci>
          </apply>
          <ci>R_la</ci>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>v_MV</ci>
        <piecewise>
          <piece>
            <cn cellml:units="mL_per_s_1">0</cn>
            <apply>
              <lt/>
              <ci>u_la</ci>
              <ci>u_lv</ci>
            </apply>
          </piece>
          <otherwise>
            <apply>
              <divide/>
              <apply>
                <minus/>
                <ci>u_la</ci>
                <ci>u_lv</ci>
              </apply>
              <ci>R_MV</ci>
            </apply>
          </otherwise>
        </piecewise>
      </apply>
      <apply>
        <eq/>
        <ci>v_AV</ci>
        <piecewise>
          <piece>
            <cn cellml:units="mL_per_s_1">0</cn>
            <apply>
              <lt/>
              <ci>u_lv_max</ci>
              <ci>u_sa</ci>
            </apply>
          </piece>
          <otherwise>
            <apply>
              <divide/>
              <apply>
                <minus/>
                <ci>u_lv_max</ci>
                <ci>u_sa</ci>
              </apply>
              <ci>R_AV</ci>
            </apply>
          </otherwise>
        </piecewise>
      </apply>
      <apply>
        <eq/>
        <ci>R_AV</ci>
        <apply>
          <times/>
          <ci>k_AV_R</ci>
          <ci>u_lv_max</ci>
        </apply>
      </apply>
      <apply>
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci>t</ci>
          </bvar>
          <ci>v_sa</ci>
        </apply>
        <apply>
          <divide/>
          <apply>
            <minus/>
            <apply>
              <minus/>
              <ci>u_sa</ci>
              <ci>u_sc</ci>
            </apply>
            <apply>
              <times/>
              <ci>R_sa</ci>
              <ci>v_sa</ci>
            </apply>
          </apply>
          <ci>L_sa</ci>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>v_sc</ci>
        <apply>
          <divide/>
          <apply>
            <minus/>
            <ci>u_sc</ci>
            <ci>u_sv</ci>
          </apply>
          <ci>R_sc</ci>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>v_sv</ci>
        <apply>
          <divide/>
          <apply>
            <minus/>
            <ci>u_sv</ci>
            <ci>u_ra</ci>
          </apply>
          <ci>R_sv</ci>
        </apply>
      </apply>
    </math>
  </component>
  <component name="cvs0">
    <variable name="t" units="second" interface="public"/>
    <variable name="phi" units="dimensionless" interface="public"/>
    <variable name="q_ra" units="mL" initial_value="4"/>
    <variable name="q_rv" units="mL" initial_value="5"/>
    <variable name="q_la" units="mL" initial_value="4"/>
    <variable name="q_lv" units="mL" initial_value="10"/>
    <variable name="q_pa" units="mL" initial_value="4"/>
    <variable name="q_pc" units="mL" initial_value="5"/>
    <variable name="q_pv" units="mL" initial_value="4"/>
    <variable name="q_sa" units="mL" initial_value="5"/>
    <variable name="q_sc" units="mL" initial_value="10"/>
    <variable name="q_sv" units="mL"/>
    <variable name="q_tot" units="mL" initial_value="5800"/>
    <variable name="q_tot_u" units="mL" initial_value="3500"/>
    <variable name="u_ra" units="mmHg"/>
    <variable name="u_rv" units="mmHg" interface="public"/>
    <variable name="u_la" units="mmHg"/>
    <variable name="u_lv" units="mmHg" interface="public"/>
    <variable name="u_pa" units="mmHg"/>
    <variable name="u_pc" units="mmHg"/>
    <variable name="u_pv" units="mmHg" interface="public"/>
    <variable name="u_sa" units="mmHg" interface="public"/>
    <variable name="du_sa" units="mmHg_per_s_1" interface="public"/>
    <variable name="u_sc" units="mmHg"/>
    <variable name="u_sv" units="mmHg"/>
    <variable name="v_TV" units="mL_per_s_2"/>
    <variable name="v_PV" units="mL_per_s_2"/>
    <variable name="v_MV" units="mL_per_s_2"/>
    <variable name="v_AV" units="mL_per_s_2"/>
    <variable name="v_pa" units="mL_per_s_2" initial_value="5"/>
    <variable name="v_pv" units="mL_per_s_2"/>
    <variable name="v_la" units="mL_per_s_2"/>
    <variable name="v_sa" units="mL_per_s_2" initial_value="5"/>
    <variable name="v_sc" units="mL_per_s_2"/>
    <variable name="v_sv" units="mL_per_s_2"/>
    <variable name="C_ra" units="mL_per_mmHg_1" initial_value="31.25"/>
    <variable name="C_pa" units="mL_per_mmHg_1" initial_value="0.76"/>
    <variable name="C_pc" units="mL_per_mmHg_1" initial_value="5.8"/>
    <variable name="C_pv" units="mL_per_mmHg_1" initial_value="25.37"/>
    <variable name="C_la" units="mL_per_mmHg_1" initial_value="19.23"/>
    <variable name="C_sa" units="mL_per_mmHg_1" initial_value="0.28"/>
    <variable name="C_sc" units="mL_per_mmHg_1" initial_value="3.72"/>
    <variable name="C_sv" units="mL_per_mmHg_1" initial_value="111.11"/>
    <variable name="R_TV" units="mmHg_s_per_mL_1" initial_value="2.5e-3"/>
    <variable name="R_PV" units="mmHg_s_per_mL_1"/>
    <variable name="R_MV" units="mmHg_s_per_mL_1" initial_value="2.5e-3"/>
    <variable name="R_AV" units="mmHg_s_per_mL_1"/>
    <variable name="R_pa" units="mmHg_s_per_mL_1" initial_value="0.023"/>
    <variable name="R_pv" units="mmHg_s_per_mL_1" initial_value="0.0894"/>
    <variable name="R_la" units="mmHg_s_per_mL_1" initial_value="0.0056"/>
    <variable name="R_sv" units="mmHg_s_per_mL_1" initial_value="0.0054"/>
    <variable name="R_sc" units="mmHg_s_per_mL_1" initial_value="0.8895"/>
    <variable name="R_sa" units="mmHg_s_per_mL_1" initial_value="0.06"/>
    <variable name="L_pa" units="mmHg_s2_per_mL_1" initial_value="1.8e-4"/>
    <variable name="L_sa" units="mmHg_s2_per_mL_1" initial_value="2.2e-4"/>
    <variable name="k_PV_R" units="s_per_mL_1" initial_value="1.4e-3"/>
    <variable name="k_AV_R" units="s_per_mL_1" initial_value="3.75e-4"/>
    <variable name="k_lv_E" units="per_mL_1" initial_value="0.014"/>
    <variable name="k_rv_E" units="per_mL_1" initial_value="0.011"/>
    <variable name="q_lv_u" units="mL" initial_value="16.77"/>
    <variable name="q_rv_u" units="mL" initial_value="40.8"/>
    <variable name="u_lv_0" units="mmHg" initial_value="1.5"/>
    <variable name="u_rv_0" units="mmHg" initial_value="1.5"/>
    <variable name="E_lv_max" units="mmHg_per_mL_1" initial_value="1.283"/>
    <variable name="E_rv_max" units="mmHg_per_mL_1" initial_value="0.757"/>
    <variable name="u_rv_max" units="mmHg"/>
    <variable name="u_lv_max" units="mmHg"/>
    <math xmlns="http://www.w3.org/1998/Math/MathML" xmlns:cellml="http://www.cellml.org/cellml/2.0#">
      <apply>
        <eq/>
        <ci>u_lv_max</ci>
        <apply>
          <plus/>
          <apply>
            <times/>
            <ci>phi</ci>
            <ci>E_lv_max</ci>
            <apply>
              <minus/>
              <ci>q_lv</ci>
              <ci>q_lv_u</ci>
            </apply>
          </apply>
          <apply>
            <times/>
            <apply>
              <minus/>
              <cn cellml:units="dimensionless">1</cn>
              <ci>phi</ci>
            </apply>
            <ci>u_lv_0</ci>
            <apply>
              <minus/>
              <apply>
                <exp/>
                <apply>
                  <times/>
                  <ci>k_lv_E</ci>
                  <ci>q_lv</ci>
                </apply>
              </apply>
              <cn cellml:units="dimensionless">1</cn>
            </apply>
          </apply>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>u_rv_max</ci>
        <apply>
          <plus/>
          <apply>
            <times/>
            <ci>phi</ci>
            <ci>E_rv_max</ci>
            <apply>
              <minus/>
              <ci>q_rv</ci>
              <ci>q_rv_u</ci>
            </apply>
          </apply>
          <apply>
            <times/>
            <apply>
              <minus/>
              <cn cellml:units="dimensionless">1</cn>
              <ci>phi</ci>
            </apply>
            <ci>u_rv_0</ci>
            <apply>
              <minus/>
              <apply>
                <exp/>
                <apply>
                  <times/>
                  <ci>k_rv_E</ci>
                  <ci>q_rv</ci>
                </apply>
              </apply>
              <cn cellml:units="dimensionless">1</cn>
            </apply>
          </apply>
        </apply>
      </apply>
      <apply>
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci>t</ci>
          </bvar>
          <ci>q_ra</ci>
        </apply>
        <apply>
          <minus/>
          <ci>v_sv</ci>
          <ci>v_TV</ci>
        </apply>
      </apply>
      <apply>
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci>t</ci>
          </bvar>
          <ci>q_rv</ci>
        </apply>
        <apply>
          <minus/>
          <ci>v_TV</ci>
          <ci>v_PV</ci>
        </apply>
      </apply>
      <apply>
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci>t</ci>
          </bvar>
          <ci>q_pa</ci>
        </apply>
        <apply>
          <minus/>
          <ci>v_PV</ci>
          <ci>v_pa</ci>
        </apply>
      </apply>
      <apply>
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci>t</ci>
          </bvar>
          <ci>q_pc</ci>
        </apply>
        <apply>
          <minus/>
          <ci>v_pa</ci>
          <ci>v_pv</ci>
        </apply>
      </apply>
      <apply>
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci>t</ci>
          </bvar>
          <ci>q_pv</ci>
        </apply>
        <apply>
          <minus/>
          <ci>v_pv</ci>
          <ci>v_la</ci>
        </apply>
      </apply>
      <apply>
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci>t</ci>
          </bvar>
          <ci>q_la</ci>
        </apply>
        <apply>
          <minus/>
          <ci>v_la</ci>
          <ci>v_MV</ci>
        </apply>
      </apply>
      <apply>
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci>t</ci>
          </bvar>
          <ci>q_lv</ci>
        </apply>
        <apply>
          <minus/>
          <ci>v_MV</ci>
          <ci>v_AV</ci>
        </apply>
      </apply>
      <apply>
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci>t</ci>
          </bvar>
          <ci>q_sa</ci>
        </apply>
        <apply>
          <minus/>
          <ci>v_AV</ci>
          <ci>v_sa</ci>
        </apply>
      </apply>
      <apply>
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci>t</ci>
          </bvar>
          <ci>q_sc</ci>
        </apply>
        <apply>
          <minus/>
          <ci>v_sa</ci>
          <ci>v_sc</ci>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>q_sv</ci>
        <apply>
          <minus/>
          <apply>
            <minus/>
            <ci>q_tot</ci>
            <ci>q_tot_u</ci>
          </apply>
          <apply>
            <plus/>
            <ci>q_ra</ci>
            <ci>q_rv</ci>
            <ci>q_pa</ci>
            <ci>q_pc</ci>
            <ci>q_pv</ci>
            <ci>q_la</ci>
            <ci>q_lv</ci>
            <ci>q_sa</ci>
            <ci>q_sc</ci>
          </apply>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>u_ra</ci>
        <apply>
          <divide/>
          <ci>q_ra</ci>
          <ci>C_ra</ci>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>u_rv</ci>
        <apply>
          <minus/>
          <ci>u_rv_max</ci>
          <apply>
            <times/>
            <ci>R_PV</ci>
            <ci>v_PV</ci>
          </apply>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>u_pa</ci>
        <apply>
          <divide/>
          <ci>q_pa</ci>
          <ci>C_pa</ci>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>u_pc</ci>
        <apply>
          <divide/>
          <ci>q_pc</ci>
          <ci>C_pc</ci>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>u_pv</ci>
        <apply>
          <divide/>
          <ci>q_pv</ci>
          <ci>C_pv</ci>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>u_la</ci>
        <apply>
          <divide/>
          <ci>q_la</ci>
          <ci>C_la</ci>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>u_lv</ci>
        <apply>
          <minus/>
          <ci>u_lv_max</ci>
          <apply>
            <times/>
            <ci>R_AV</ci>
            <ci>v_AV</ci>
          </apply>
        </apply>
      </apply>
      <!--       ode(u_sa1,t) = du_sa;-->
      <apply>
        <eq/>
        <ci>u_sa</ci>
        <apply>
          <divide/>
          <ci>q_sa</ci>
          <ci>C_sa</ci>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>u_sc</ci>
        <apply>
          <divide/>
          <ci>q_sc</ci>
          <ci>C_sc</ci>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>u_sv</ci>
        <apply>
          <divide/>
          <ci>q_sv</ci>
          <ci>C_sv</ci>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>du_sa</ci>
        <apply>
          <divide/>
          <apply>
            <minus/>
            <ci>v_AV</ci>
            <ci>v_sa</ci>
          </apply>
          <ci>C_sa</ci>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>v_TV</ci>
        <piecewise>
          <piece>
            <cn cellml:units="mL_per_s_2">0</cn>
            <apply>
              <lt/>
              <ci>u_ra</ci>
              <ci>u_rv</ci>
            </apply>
          </piece>
          <otherwise>
            <apply>
              <divide/>
              <apply>
                <minus/>
                <ci>u_ra</ci>
                <ci>u_rv</ci>
              </apply>
              <ci>R_TV</ci>
            </apply>
          </otherwise>
        </piecewise>
      </apply>
      <apply>
        <eq/>
        <ci>v_PV</ci>
        <piecewise>
          <piece>
            <cn cellml:units="mL_per_s_2">0</cn>
            <apply>
              <lt/>
              <ci>u_rv_max</ci>
              <ci>u_pa</ci>
            </apply>
          </piece>
          <otherwise>
            <apply>
              <divide/>
              <apply>
                <minus/>
                <ci>u_rv_max</ci>
                <ci>u_pa</ci>
              </apply>
              <ci>R_PV</ci>
            </apply>
          </otherwise>
        </piecewise>
      </apply>
      <apply>
        <eq/>
        <ci>R_PV</ci>
        <apply>
          <times/>
          <ci>k_PV_R</ci>
          <ci>u_rv_max</ci>
        </apply>
      </apply>
      <apply>
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci>t</ci>
          </bvar>
          <ci>v_pa</ci>
        </apply>
        <apply>
          <divide/>
          <apply>
            <minus/>
            <apply>
              <minus/>
              <ci>u_pa</ci>
              <ci>u_pc</ci>
            </apply>
            <apply>
              <times/>
              <ci>R_pa</ci>
              <ci>v_pa</ci>
            </apply>
          </apply>
          <ci>L_pa</ci>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>v_pv</ci>
        <apply>
          <divide/>
          <apply>
            <minus/>
            <ci>u_pc</ci>
            <ci>u_pv</ci>
          </apply>
          <ci>R_pv</ci>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>v_la</ci>
        <apply>
          <divide/>
          <apply>
            <minus/>
            <ci>u_pv</ci>
            <ci>u_la</ci>
          </apply>
          <ci>R_la</ci>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>v_MV</ci>
        <piecewise>
          <piece>
            <cn cellml:units="mL_per_s_2">0</cn>
            <apply>
              <lt/>
              <ci>u_la</ci>
              <ci>u_lv</ci>
            </apply>
          </piece>
          <otherwise>
            <apply>
              <divide/>
              <apply>
                <minus/>
                <ci>u_la</ci>
                <ci>u_lv</ci>
              </apply>
              <ci>R_MV</ci>
            </apply>
          </otherwise>
        </piecewise>
      </apply>
      <apply>
        <eq/>
        <ci>v_AV</ci>
        <piecewise>
          <piece>
            <cn cellml:units="mL_per_s_2">0</cn>
            <apply>
              <lt/>
              <ci>u_lv_max</ci>
              <ci>u_sa</ci>
            </apply>
          </piece>
          <otherwise>
            <apply>
              <divide/>
              <apply>
                <minus/>
                <ci>u_lv_max</ci>
                <ci>u_sa</ci>
              </apply>
              <ci>R_AV</ci>
            </apply>
          </otherwise>
        </piecewise>
      </apply>
      <apply>
        <eq/>
        <ci>R_AV</ci>
        <apply>
          <times/>
          <ci>k_AV_R</ci>
          <ci>u_lv_max</ci>
        </apply>
      </apply>
      <apply>
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci>t</ci>
          </bvar>
          <ci>v_sa</ci>
        </apply>
        <apply>
          <divide/>
          <apply>
            <minus/>
            <apply>
              <minus/>
              <ci>u_sa</ci>
              <ci>u_sc</ci>
            </apply>
            <apply>
              <times/>
              <ci>R_sa</ci>
              <ci>v_sa</ci>
            </apply>
          </apply>
          <ci>L_sa</ci>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>v_sc</ci>
        <apply>
          <divide/>
          <apply>
            <minus/>
            <ci>u_sc</ci>
            <ci>u_sv</ci>
          </apply>
          <ci>R_sc</ci>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>v_sv</ci>
        <apply>
          <divide/>
          <apply>
            <minus/>
            <ci>u_sv</ci>
            <ci>u_ra</ci>
          </apply>
          <ci>R_sv</ci>
        </apply>
      </apply>
    </math>
  </component>
  <component name="activation">
    <variable name="t" units="second" interface="public"/>
    <variable name="T" units="second" interface="public"/>
    <variable name="T_0" units="second" interface="public"/>
    <variable name="T_sys_0" units="second" initial_value="0.4"/>
    <variable name="k_sys" units="s2" initial_value="0.075"/>
    <variable name="T_sys" units="second"/>
    <variable name="u" units="dimensionless" interface="public"/>
    <variable name="phi" units="dimensionless" interface="public"/>
    <math xmlns="http://www.w3.org/1998/Math/MathML" xmlns:cellml="http://www.cellml.org/cellml/2.0#">
      <apply>
        <eq/>
        <ci>u</ci>
        <apply>
          <minus/>
          <apply>
            <divide/>
            <ci>t</ci>
            <ci>T_0</ci>
          </apply>
          <apply>
            <floor/>
            <apply>
              <divide/>
              <ci>t</ci>
              <ci>T_0</ci>
            </apply>
          </apply>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>T</ci>
        <apply>
          <minus/>
          <ci>t</ci>
          <apply>
            <times/>
            <ci>T_0</ci>
            <apply>
              <floor/>
              <apply>
                <divide/>
                <ci>t</ci>
                <ci>T_0</ci>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>T_sys</ci>
        <piecewise>
          <piece>
            <apply>
              <minus/>
              <ci>T_sys_0</ci>
              <apply>
                <divide/>
                <ci>k_sys</ci>
                <ci>T</ci>
              </apply>
            </apply>
            <apply>
              <gt/>
              <ci>T</ci>
              <cn cellml:units="second">0</cn>
            </apply>
          </piece>
          <otherwise>
            <apply>
              <minus/>
              <ci>T_sys_0</ci>
              <apply>
                <divide/>
                <ci>k_sys</ci>
                <ci>T_0</ci>
              </apply>
            </apply>
          </otherwise>
        </piecewise>
      </apply>
      <apply>
        <eq/>
        <ci>phi</ci>
        <piecewise>
          <piece>
            <apply>
              <power/>
              <apply>
                <sin/>
                <apply>
                  <divide/>
                  <apply>
                    <times/>
                    <pi/>
                    <ci>u</ci>
                    <ci>T</ci>
                  </apply>
                  <ci>T_sys_0</ci>
                </apply>
              </apply>
              <cn cellml:units="dimensionless">2</cn>
            </apply>
            <apply>
              <and/>
              <apply>
                <geq/>
                <ci>u</ci>
                <cn cellml:units="dimensionless">0</cn>
              </apply>
              <apply>
                <leq/>
                <ci>u</ci>
                <apply>
                  <divide/>
                  <ci>T_sys_0</ci>
                  <ci>T</ci>
                </apply>
              </apply>
            </apply>
          </piece>
          <otherwise>
            <cn cellml:units="dimensionless">0</cn>
          </otherwise>
        </piecewise>
      </apply>
    </math>
  </component>
  <component name="environment">
    <variable name="time" units="millisecond" interface="public" id="environment_time"/>
  </component>
  <component name="membrane" id="id_00001">
    <variable name="V" units="millivolt" initial_value="-86.2" interface="public" id="membrane_V"/>
    <variable name="R" units="joule_per_mole_kelvin" initial_value="8.314" interface="public"/>
    <variable name="T" units="kelvin" initial_value="310" interface="public"/>
    <variable name="F" units="coulomb_per_millimole" initial_value="96.485" interface="public"/>
    <variable name="Cm" units="picoF" initial_value="185" interface="public"/>
    <variable name="V_c" units="mL" initial_value="16404" interface="public"/>
    <variable name="time" units="millisecond" interface="public"/>
    <variable name="i_K1" units="picoA_per_picoF" interface="public"/>
    <variable name="i_to" units="picoA_per_picoF" interface="public"/>
    <variable name="i_Kr" units="picoA_per_picoF" interface="public"/>
    <variable name="i_Ks" units="picoA_per_picoF" interface="public"/>
    <variable name="i_CaL" units="picoA_per_picoF" interface="public"/>
    <variable name="i_NaK" units="picoA_per_picoF" interface="public"/>
    <variable name="i_Na" units="picoA_per_picoF" interface="public"/>
    <variable name="i_b_Na" units="picoA_per_picoF" interface="public"/>
    <variable name="i_NaCa" units="picoA_per_picoF" interface="public"/>
    <variable name="i_b_Ca" units="picoA_per_picoF" interface="public"/>
    <variable name="i_p_K" units="picoA_per_picoF" interface="public"/>
    <variable name="i_p_Ca" units="picoA_per_picoF" interface="public"/>
    <variable name="i_Stim" units="picoA_per_picoF" interface="public"/>
    <variable name="stim_start" units="millisecond" initial_value="10"/>
    <variable name="stim_period" units="millisecond" initial_value="1000" interface="public"/>
    <variable name="stim_duration" units="millisecond" initial_value="1"/>
    <variable name="stim_amplitude" units="picoA_per_picoF" initial_value="-52"/>
    <math xmlns="http://www.w3.org/1998/Math/MathML" xmlns:cellml="http://www.cellml.org/cellml/2.0#">
      <apply>
        <eq/>
        <ci>i_Stim</ci>
        <piecewise>
          <piece>
            <ci>stim_amplitude</ci>
            <apply>
              <and/>
              <apply>
                <geq/>
                <apply>
                  <minus/>
                  <ci>time</ci>
                  <apply>
                    <times/>
                    <apply>
                      <floor/>
                      <apply>
                        <divide/>
                        <ci>time</ci>
                        <ci>stim_period</ci>
                      </apply>
                    </apply>
                    <ci>stim_period</ci>
                  </apply>
                </apply>
                <ci>stim_start</ci>
              </apply>
              <apply>
                <leq/>
                <apply>
                  <minus/>
                  <ci>time</ci>
                  <apply>
                    <times/>
                    <apply>
                      <floor/>
                      <apply>
                        <divide/>
                        <ci>time</ci>
                        <ci>stim_period</ci>
                      </apply>
                    </apply>
                    <ci>stim_period</ci>
                  </apply>
                </apply>
                <apply>
                  <plus/>
                  <ci>stim_start</ci>
                  <ci>stim_duration</ci>
                </apply>
              </apply>
            </apply>
          </piece>
          <otherwise>
            <cn cellml:units="picoA_per_picoF">0</cn>
          </otherwise>
        </piecewise>
      </apply>
      <apply>
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci>time</ci>
          </bvar>
          <ci>V</ci>
        </apply>
        <apply>
          <minus/>
          <apply>
            <plus/>
            <ci>i_K1</ci>
            <ci>i_to</ci>
            <ci>i_Kr</ci>
            <ci>i_Ks</ci>
            <ci>i_CaL</ci>
            <ci>i_NaK</ci>
            <ci>i_Na</ci>
            <ci>i_b_Na</ci>
            <ci>i_NaCa</ci>
            <ci>i_b_Ca</ci>
            <ci>i_p_K</ci>
            <ci>i_p_Ca</ci>
            <ci>i_Stim</ci>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  <component name="fast_sodium_current" id="id_00009">
    <variable name="i_Na" units="picoA_per_picoF" interface="public" id="fast_sodium_current_i_Na"/>
    <variable name="g_Na" units="nanoS_per_picoF" initial_value="14.838"/>
    <variable name="E_Na" units="millivolt" interface="public"/>
    <variable name="time" units="millisecond" interface="public_and_private"/>
    <variable name="V" units="millivolt" interface="public_and_private"/>
    <variable name="m" units="dimensionless" interface="private"/>
    <variable name="h" units="dimensionless" interface="private"/>
    <variable name="j" units="dimensionless" interface="private"/>
    <math xmlns="http://www.w3.org/1998/Math/MathML" xmlns:cellml="http://www.cellml.org/cellml/2.0#">
      <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>
          <ci>j</ci>
          <apply>
            <minus/>
            <ci>V</ci>
            <ci>E_Na</ci>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  <component name="fast_sodium_current_m_gate" id="id_00010">
    <variable name="m" units="dimensionless" initial_value="0" interface="public"/>
    <variable name="V" units="millivolt" interface="public"/>
    <variable name="time" units="millisecond" interface="public"/>
    <variable name="m_inf" units="dimensionless"/>
    <variable name="alpha_m" units="dimensionless"/>
    <variable name="beta_m" units="dimensionless"/>
    <variable name="tau_m" units="millisecond"/>
    <math xmlns="http://www.w3.org/1998/Math/MathML" xmlns:cellml="http://www.cellml.org/cellml/2.0#">
      <apply>
        <eq/>
        <ci>m_inf</ci>
        <apply>
          <divide/>
          <cn cellml:units="dimensionless">1</cn>
          <apply>
            <power/>
            <apply>
              <plus/>
              <cn cellml:units="dimensionless">1</cn>
              <apply>
                <exp/>
                <apply>
                  <divide/>
                  <apply>
                    <minus/>
                    <apply>
                      <minus/>
                      <cn cellml:units="millivolt">56.86</cn>
                    </apply>
                    <ci>V</ci>
                  </apply>
                  <cn cellml:units="millivolt">9.03</cn>
                </apply>
              </apply>
            </apply>
            <cn cellml:units="dimensionless">2</cn>
          </apply>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>alpha_m</ci>
        <apply>
          <divide/>
          <cn cellml:units="dimensionless">1</cn>
          <apply>
            <plus/>
            <cn cellml:units="dimensionless">1</cn>
            <apply>
              <exp/>
              <apply>
                <divide/>
                <apply>
                  <minus/>
                  <apply>
                    <minus/>
                    <cn cellml:units="millivolt">60</cn>
                  </apply>
                  <ci>V</ci>
                </apply>
                <cn cellml:units="millivolt">5</cn>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>beta_m</ci>
        <apply>
          <plus/>
          <apply>
            <divide/>
            <cn cellml:units="dimensionless">0.1</cn>
            <apply>
              <plus/>
              <cn cellml:units="dimensionless">1</cn>
              <apply>
                <exp/>
                <apply>
                  <divide/>
                  <apply>
                    <plus/>
                    <ci>V</ci>
                    <cn cellml:units="millivolt">35</cn>
                  </apply>
                  <cn cellml:units="millivolt">5</cn>
                </apply>
              </apply>
            </apply>
          </apply>
          <apply>
            <divide/>
            <cn cellml:units="dimensionless">0.1</cn>
            <apply>
              <plus/>
              <cn cellml:units="dimensionless">1</cn>
              <apply>
                <exp/>
                <apply>
                  <divide/>
                  <apply>
                    <minus/>
                    <ci>V</ci>
                    <cn cellml:units="millivolt">50</cn>
                  </apply>
                  <cn cellml:units="millivolt">200</cn>
                </apply>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>tau_m</ci>
        <apply>
          <times/>
          <cn cellml:units="millisecond">1</cn>
          <ci>alpha_m</ci>
          <ci>beta_m</ci>
        </apply>
      </apply>
      <apply>
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci>time</ci>
          </bvar>
          <ci>m</ci>
        </apply>
        <apply>
          <divide/>
          <apply>
            <minus/>
            <ci>m_inf</ci>
            <ci>m</ci>
          </apply>
          <ci>tau_m</ci>
        </apply>
      </apply>
    </math>
  </component>
  <component name="fast_sodium_current_h_gate" id="id_00011">
    <variable name="h" units="dimensionless" initial_value="0.75" interface="public"/>
    <variable name="V" units="millivolt" interface="public"/>
    <variable name="time" units="millisecond" interface="public"/>
    <variable name="h_inf" units="dimensionless"/>
    <variable name="alpha_h" units="per_millisecond"/>
    <variable name="beta_h" units="per_millisecond"/>
    <variable name="tau_h" units="millisecond"/>
    <math xmlns="http://www.w3.org/1998/Math/MathML" xmlns:cellml="http://www.cellml.org/cellml/2.0#">
      <apply>
        <eq/>
        <ci>h_inf</ci>
        <apply>
          <divide/>
          <cn cellml:units="dimensionless">1</cn>
          <apply>
            <power/>
            <apply>
              <plus/>
              <cn cellml:units="dimensionless">1</cn>
              <apply>
                <exp/>
                <apply>
                  <divide/>
                  <apply>
                    <plus/>
                    <ci>V</ci>
                    <cn cellml:units="millivolt">71.55</cn>
                  </apply>
                  <cn cellml:units="millivolt">7.43</cn>
                </apply>
              </apply>
            </apply>
            <cn cellml:units="dimensionless">2</cn>
          </apply>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>alpha_h</ci>
        <piecewise>
          <piece>
            <apply>
              <times/>
              <cn cellml:units="per_millisecond">0.057</cn>
              <apply>
                <exp/>
                <apply>
                  <divide/>
                  <apply>
                    <minus/>
                    <apply>
                      <plus/>
                      <ci>V</ci>
                      <cn cellml:units="millivolt">80</cn>
                    </apply>
                  </apply>
                  <cn cellml:units="millivolt">6.8</cn>
                </apply>
              </apply>
            </apply>
            <apply>
              <lt/>
              <ci>V</ci>
              <apply>
                <minus/>
                <cn cellml:units="millivolt">40</cn>
              </apply>
            </apply>
          </piece>
          <otherwise>
            <cn cellml:units="per_millisecond">0</cn>
          </otherwise>
        </piecewise>
      </apply>
      <apply>
        <eq/>
        <ci>beta_h</ci>
        <piecewise>
          <piece>
            <apply>
              <plus/>
              <apply>
                <times/>
                <cn cellml:units="per_millisecond">2.7</cn>
                <apply>
                  <exp/>
                  <apply>
                    <times/>
                    <cn cellml:units="per_millivolt">0.079</cn>
                    <ci>V</ci>
                  </apply>
                </apply>
              </apply>
              <apply>
                <times/>
                <cn cellml:units="per_millisecond">310000</cn>
                <apply>
                  <exp/>
                  <apply>
                    <times/>
                    <cn cellml:units="per_millivolt">0.3485</cn>
                    <ci>V</ci>
                  </apply>
                </apply>
              </apply>
            </apply>
            <apply>
              <lt/>
              <ci>V</ci>
              <apply>
                <minus/>
                <cn cellml:units="millivolt">40</cn>
              </apply>
            </apply>
          </piece>
          <otherwise>
            <apply>
              <divide/>
              <cn cellml:units="per_millisecond">0.77</cn>
              <apply>
                <times/>
                <cn cellml:units="dimensionless">0.13</cn>
                <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                    <exp/>
                    <apply>
                      <divide/>
                      <apply>
                        <plus/>
                        <ci>V</ci>
                        <cn cellml:units="millivolt">10.66</cn>
                      </apply>
                      <apply>
                        <minus/>
                        <cn cellml:units="millivolt">11.1</cn>
                      </apply>
                    </apply>
                  </apply>
                </apply>
              </apply>
            </apply>
          </otherwise>
        </piecewise>
      </apply>
      <apply>
        <eq/>
        <ci>tau_h</ci>
        <apply>
          <divide/>
          <cn cellml:units="dimensionless">1</cn>
          <apply>
            <plus/>
            <ci>alpha_h</ci>
            <ci>beta_h</ci>
          </apply>
        </apply>
      </apply>
      <apply>
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci>time</ci>
          </bvar>
          <ci>h</ci>
        </apply>
        <apply>
          <divide/>
          <apply>
            <minus/>
            <ci>h_inf</ci>
            <ci>h</ci>
          </apply>
          <ci>tau_h</ci>
        </apply>
      </apply>
    </math>
  </component>
  <component name="fast_sodium_current_j_gate" id="id_00012">
    <variable name="j" units="dimensionless" initial_value="0.75" interface="public"/>
    <variable name="V" units="millivolt" interface="public"/>
    <variable name="time" units="millisecond" interface="public"/>
    <variable name="j_inf" units="dimensionless"/>
    <variable name="alpha_j" units="per_millisecond"/>
    <variable name="beta_j" units="per_millisecond"/>
    <variable name="tau_j" units="millisecond"/>
    <math xmlns="http://www.w3.org/1998/Math/MathML" xmlns:cellml="http://www.cellml.org/cellml/2.0#">
      <apply>
        <eq/>
        <ci>j_inf</ci>
        <apply>
          <divide/>
          <cn cellml:units="dimensionless">1</cn>
          <apply>
            <power/>
            <apply>
              <plus/>
              <cn cellml:units="dimensionless">1</cn>
              <apply>
                <exp/>
                <apply>
                  <divide/>
                  <apply>
                    <plus/>
                    <ci>V</ci>
                    <cn cellml:units="millivolt">71.55</cn>
                  </apply>
                  <cn cellml:units="millivolt">7.43</cn>
                </apply>
              </apply>
            </apply>
            <cn cellml:units="dimensionless">2</cn>
          </apply>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>alpha_j</ci>
        <piecewise>
          <piece>
            <apply>
              <divide/>
              <apply>
                <divide/>
                <apply>
                  <times/>
                  <apply>
                    <minus/>
                    <apply>
                      <times/>
                      <apply>
                        <minus/>
                        <cn cellml:units="per_millisecond">25428</cn>
                      </apply>
                      <apply>
                        <exp/>
                        <apply>
                          <times/>
                          <cn cellml:units="per_millivolt">0.2444</cn>
                          <ci>V</ci>
                        </apply>
                      </apply>
                    </apply>
                    <apply>
                      <times/>
                      <cn type="e-notation" cellml:units="per_millisecond">6.948<sep/>-6</cn>
                      <apply>
                        <exp/>
                        <apply>
                          <times/>
                          <apply>
                            <minus/>
                            <cn cellml:units="per_millivolt">0.04391</cn>
                          </apply>
                          <ci>V</ci>
                        </apply>
                      </apply>
                    </apply>
                  </apply>
                  <apply>
                    <plus/>
                    <ci>V</ci>
                    <cn cellml:units="millivolt">37.78</cn>
                  </apply>
                </apply>
                <cn cellml:units="millivolt">1</cn>
              </apply>
              <apply>
                <plus/>
                <cn cellml:units="dimensionless">1</cn>
                <apply>
                  <exp/>
                  <apply>
                    <times/>
                    <cn cellml:units="per_millivolt">0.311</cn>
                    <apply>
                      <plus/>
                      <ci>V</ci>
                      <cn cellml:units="millivolt">79.23</cn>
                    </apply>
                  </apply>
                </apply>
              </apply>
            </apply>
            <apply>
              <lt/>
              <ci>V</ci>
              <apply>
                <minus/>
                <cn cellml:units="millivolt">40</cn>
              </apply>
            </apply>
          </piece>
          <otherwise>
            <cn cellml:units="per_millisecond">0</cn>
          </otherwise>
        </piecewise>
      </apply>
      <apply>
        <eq/>
        <ci>beta_j</ci>
        <piecewise>
          <piece>
            <apply>
              <divide/>
              <apply>
                <times/>
                <cn cellml:units="per_millisecond">0.02424</cn>
                <apply>
                  <exp/>
                  <apply>
                    <times/>
                    <apply>
                      <minus/>
                      <cn cellml:units="per_millivolt">0.01052</cn>
                    </apply>
                    <ci>V</ci>
                  </apply>
                </apply>
              </apply>
              <apply>
                <plus/>
                <cn cellml:units="dimensionless">1</cn>
                <apply>
                  <exp/>
                  <apply>
                    <times/>
                    <apply>
                      <minus/>
                      <cn cellml:units="per_millivolt">0.1378</cn>
                    </apply>
                    <apply>
                      <plus/>
                      <ci>V</ci>
                      <cn cellml:units="millivolt">40.14</cn>
                    </apply>
                  </apply>
                </apply>
              </apply>
            </apply>
            <apply>
              <lt/>
              <ci>V</ci>
              <apply>
                <minus/>
                <cn cellml:units="millivolt">40</cn>
              </apply>
            </apply>
          </piece>
          <otherwise>
            <apply>
              <divide/>
              <apply>
                <times/>
                <cn cellml:units="per_millisecond">0.6</cn>
                <apply>
                  <exp/>
                  <apply>
                    <times/>
                    <cn cellml:units="per_millivolt">0.057</cn>
                    <ci>V</ci>
                  </apply>
                </apply>
              </apply>
              <apply>
                <plus/>
                <cn cellml:units="dimensionless">1</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">32</cn>
                    </apply>
                  </apply>
                </apply>
              </apply>
            </apply>
          </otherwise>
        </piecewise>
      </apply>
      <apply>
        <eq/>
        <ci>tau_j</ci>
        <apply>
          <divide/>
          <cn cellml:units="dimensionless">1</cn>
          <apply>
            <plus/>
            <ci>alpha_j</ci>
            <ci>beta_j</ci>
          </apply>
        </apply>
      </apply>
      <apply>
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci>time</ci>
          </bvar>
          <ci>j</ci>
        </apply>
        <apply>
          <divide/>
          <apply>
            <minus/>
            <ci>j_inf</ci>
            <ci>j</ci>
          </apply>
          <ci>tau_j</ci>
        </apply>
      </apply>
    </math>
  </component>
  <component name="L_type_Ca_current" id="id_00014">
    <variable name="i_CaL" units="picoA_per_picoF" interface="public" id="L_type_Ca_current_i_CaL"/>
    <variable name="g_CaL" units="litre_per_farad_second" initial_value="0.175"/>
    <variable name="Ca_o" units="millimolar" interface="public"/>
    <variable name="Ca_i" units="millimolar" interface="public_and_private"/>
    <variable name="R" units="joule_per_mole_kelvin" interface="public"/>
    <variable name="F" units="coulomb_per_millimole" interface="public"/>
    <variable name="T" units="kelvin" interface="public"/>
    <variable name="time" units="millisecond" interface="public_and_private"/>
    <variable name="V" units="millivolt" interface="public_and_private"/>
    <variable name="d" units="dimensionless" interface="public_and_private"/>
    <variable name="f" units="dimensionless" interface="private"/>
    <variable name="fCa" units="dimensionless" interface="private"/>
    <math xmlns="http://www.w3.org/1998/Math/MathML" xmlns:cellml="http://www.cellml.org/cellml/2.0#">
      <apply>
        <eq/>
        <ci>i_CaL</ci>
        <apply>
          <divide/>
          <apply>
            <times/>
            <apply>
              <divide/>
              <apply>
                <times/>
                <ci>g_CaL</ci>
                <ci>d</ci>
                <ci>f</ci>
                <ci>fCa</ci>
                <cn cellml:units="dimensionless">4</cn>
                <ci>V</ci>
                <apply>
                  <power/>
                  <ci>F</ci>
                  <cn cellml:units="dimensionless">2</cn>
                </apply>
              </apply>
              <apply>
                <times/>
                <ci>R</ci>
                <ci>T</ci>
              </apply>
            </apply>
            <apply>
              <minus/>
              <apply>
                <times/>
                <ci>Ca_i</ci>
                <apply>
                  <exp/>
                  <apply>
                    <divide/>
                    <apply>
                      <times/>
                      <cn cellml:units="dimensionless">2</cn>
                      <ci>V</ci>
                      <ci>F</ci>
                    </apply>
                    <apply>
                      <times/>
                      <ci>R</ci>
                      <ci>T</ci>
                    </apply>
                  </apply>
                </apply>
              </apply>
              <apply>
                <times/>
                <cn cellml:units="dimensionless">0.341</cn>
                <ci>Ca_o</ci>
              </apply>
            </apply>
          </apply>
          <apply>
            <minus/>
            <apply>
              <exp/>
              <apply>
                <divide/>
                <apply>
                  <times/>
                  <cn cellml:units="dimensionless">2</cn>
                  <ci>V</ci>
                  <ci>F</ci>
                </apply>
                <apply>
                  <times/>
                  <ci>R</ci>
                  <ci>T</ci>
                </apply>
              </apply>
            </apply>
            <cn cellml:units="dimensionless">1</cn>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  <component name="L_type_Ca_current_d_gate" id="id_00015">
    <variable name="d" units="dimensionless" initial_value="0" interface="public"/>
    <variable name="V" units="millivolt" interface="public"/>
    <variable name="time" units="millisecond" interface="public"/>
    <variable name="d_inf" units="dimensionless"/>
    <variable name="alpha_d" units="dimensionless"/>
    <variable name="beta_d" units="dimensionless"/>
    <variable name="gamma_d" units="millisecond"/>
    <variable name="tau_d" units="millisecond"/>
    <math xmlns="http://www.w3.org/1998/Math/MathML" xmlns:cellml="http://www.cellml.org/cellml/2.0#">
      <apply>
        <eq/>
        <ci>d_inf</ci>
        <apply>
          <divide/>
          <cn cellml:units="dimensionless">1</cn>
          <apply>
            <plus/>
            <cn cellml:units="dimensionless">1</cn>
            <apply>
              <exp/>
              <apply>
                <divide/>
                <apply>
                  <minus/>
                  <apply>
                    <minus/>
                    <cn cellml:units="millivolt">5</cn>
                  </apply>
                  <ci>V</ci>
                </apply>
                <cn cellml:units="millivolt">7.5</cn>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>alpha_d</ci>
        <apply>
          <plus/>
          <apply>
            <divide/>
            <cn cellml:units="dimensionless">1.4</cn>
            <apply>
              <plus/>
              <cn cellml:units="dimensionless">1</cn>
              <apply>
                <exp/>
                <apply>
                  <divide/>
                  <apply>
                    <minus/>
                    <apply>
                      <minus/>
                      <cn cellml:units="millivolt">35</cn>
                    </apply>
                    <ci>V</ci>
                  </apply>
                  <cn cellml:units="millivolt">13</cn>
                </apply>
              </apply>
            </apply>
          </apply>
          <cn cellml:units="dimensionless">0.25</cn>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>beta_d</ci>
        <apply>
          <divide/>
          <cn cellml:units="dimensionless">1.4</cn>
          <apply>
            <plus/>
            <cn cellml:units="dimensionless">1</cn>
            <apply>
              <exp/>
              <apply>
                <divide/>
                <apply>
                  <plus/>
                  <ci>V</ci>
                  <cn cellml:units="millivolt">5</cn>
                </apply>
                <cn cellml:units="millivolt">5</cn>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>gamma_d</ci>
        <apply>
          <divide/>
          <cn cellml:units="millisecond">1</cn>
          <apply>
            <plus/>
            <cn cellml:units="dimensionless">1</cn>
            <apply>
              <exp/>
              <apply>
                <divide/>
                <apply>
                  <minus/>
                  <cn cellml:units="millivolt">50</cn>
                  <ci>V</ci>
                </apply>
                <cn cellml:units="millivolt">20</cn>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>tau_d</ci>
        <apply>
          <plus/>
          <apply>
            <times/>
            <cn cellml:units="millisecond">1</cn>
            <ci>alpha_d</ci>
            <ci>beta_d</ci>
          </apply>
          <ci>gamma_d</ci>
        </apply>
      </apply>
      <apply>
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci>time</ci>
          </bvar>
          <ci>d</ci>
        </apply>
        <apply>
          <divide/>
          <apply>
            <minus/>
            <ci>d_inf</ci>
            <ci>d</ci>
          </apply>
          <ci>tau_d</ci>
        </apply>
      </apply>
    </math>
  </component>
  <component name="L_type_Ca_current_f_gate" id="id_00016">
    <variable name="f" units="dimensionless" initial_value="1" interface="public"/>
    <variable name="V" units="millivolt" interface="public"/>
    <variable name="time" units="millisecond" interface="public"/>
    <variable name="f_inf" units="dimensionless"/>
    <variable name="tau_f" units="millisecond"/>
    <math xmlns="http://www.w3.org/1998/Math/MathML" xmlns:cellml="http://www.cellml.org/cellml/2.0#">
      <apply>
        <eq/>
        <ci>f_inf</ci>
        <apply>
          <divide/>
          <cn cellml:units="dimensionless">1</cn>
          <apply>
            <plus/>
            <cn cellml:units="dimensionless">1</cn>
            <apply>
              <exp/>
              <apply>
                <divide/>
                <apply>
                  <plus/>
                  <ci>V</ci>
                  <cn cellml:units="millivolt">20</cn>
                </apply>
                <cn cellml:units="millivolt">7</cn>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>tau_f</ci>
        <apply>
          <plus/>
          <apply>
            <times/>
            <cn cellml:units="millisecond">1125</cn>
            <apply>
              <exp/>
              <apply>
                <divide/>
                <apply>
                  <minus/>
                  <apply>
                    <power/>
                    <apply>
                      <plus/>
                      <ci>V</ci>
                      <cn cellml:units="millivolt">27</cn>
                    </apply>
                    <cn cellml:units="dimensionless">2</cn>
                  </apply>
                </apply>
                <cn cellml:units="millivolt_square">240</cn>
              </apply>
            </apply>
          </apply>
          <cn cellml:units="millisecond">80</cn>
          <apply>
            <divide/>
            <cn cellml:units="millisecond">165</cn>
            <apply>
              <plus/>
              <cn cellml:units="dimensionless">1</cn>
              <apply>
                <exp/>
                <apply>
                  <divide/>
                  <apply>
                    <minus/>
                    <cn cellml:units="millivolt">25</cn>
                    <ci>V</ci>
                  </apply>
                  <cn cellml:units="millivolt">10</cn>
                </apply>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      <apply>
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci>time</ci>
          </bvar>
          <ci>f</ci>
        </apply>
        <apply>
          <divide/>
          <apply>
            <minus/>
            <ci>f_inf</ci>
            <ci>f</ci>
          </apply>
          <ci>tau_f</ci>
        </apply>
      </apply>
    </math>
  </component>
  <component name="L_type_Ca_current_fCa_gate" id="id_00017">
    <variable name="fCa" units="dimensionless" initial_value="1" interface="public"/>
    <variable name="Ca_i" units="millimolar" interface="public"/>
    <variable name="V" units="millivolt" interface="public"/>
    <variable name="time" units="millisecond" interface="public"/>
    <variable name="alpha_fCa" units="dimensionless"/>
    <variable name="beta_fCa" units="dimensionless"/>
    <variable name="gama_fCa" units="dimensionless"/>
    <variable name="fCa_inf" units="dimensionless"/>
    <variable name="tau_fCa" units="millisecond"/>
    <variable name="d_fCa" units="per_millisecond"/>
    <math xmlns="http://www.w3.org/1998/Math/MathML" xmlns:cellml="http://www.cellml.org/cellml/2.0#">
      <apply>
        <eq/>
        <ci>alpha_fCa</ci>
        <apply>
          <divide/>
          <cn cellml:units="dimensionless">1</cn>
          <apply>
            <plus/>
            <cn cellml:units="dimensionless">1</cn>
            <apply>
              <power/>
              <apply>
                <divide/>
                <ci>Ca_i</ci>
                <cn cellml:units="millimolar">0.000325</cn>
              </apply>
              <cn cellml:units="dimensionless">8</cn>
            </apply>
          </apply>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>beta_fCa</ci>
        <apply>
          <divide/>
          <cn cellml:units="dimensionless">0.1</cn>
          <apply>
            <plus/>
            <cn cellml:units="dimensionless">1</cn>
            <apply>
              <exp/>
              <apply>
                <divide/>
                <apply>
                  <minus/>
                  <ci>Ca_i</ci>
                  <cn cellml:units="millimolar">0.0005</cn>
                </apply>
                <cn cellml:units="millimolar">0.0001</cn>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>gama_fCa</ci>
        <apply>
          <divide/>
          <cn cellml:units="dimensionless">0.2</cn>
          <apply>
            <plus/>
            <cn cellml:units="dimensionless">1</cn>
            <apply>
              <exp/>
              <apply>
                <divide/>
                <apply>
                  <minus/>
                  <ci>Ca_i</ci>
                  <cn cellml:units="millimolar">0.00075</cn>
                </apply>
                <cn cellml:units="millimolar">0.0008</cn>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>fCa_inf</ci>
        <apply>
          <divide/>
          <apply>
            <plus/>
            <ci>alpha_fCa</ci>
            <ci>beta_fCa</ci>
            <ci>gama_fCa</ci>
            <cn cellml:units="dimensionless">0.23</cn>
          </apply>
          <cn cellml:units="dimensionless">1.46</cn>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>tau_fCa</ci>
        <cn cellml:units="millisecond">2</cn>
      </apply>
      <apply>
        <eq/>
        <ci>d_fCa</ci>
        <apply>
          <divide/>
          <apply>
            <minus/>
            <ci>fCa_inf</ci>
            <ci>fCa</ci>
          </apply>
          <ci>tau_fCa</ci>
        </apply>
      </apply>
      <apply>
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci>time</ci>
          </bvar>
          <ci>fCa</ci>
        </apply>
        <piecewise>
          <piece>
            <cn cellml:units="per_millisecond">0</cn>
            <apply>
              <and/>
              <apply>
                <gt/>
                <ci>fCa_inf</ci>
                <ci>fCa</ci>
              </apply>
              <apply>
                <gt/>
                <ci>V</ci>
                <apply>
                  <minus/>
                  <cn cellml:units="millivolt">60</cn>
                </apply>
              </apply>
            </apply>
          </piece>
          <otherwise>
            <ci>d_fCa</ci>
          </otherwise>
        </piecewise>
      </apply>
    </math>
  </component>
  <component name="rapid_time_dependent_potassium_current" id="id_00004">
    <variable name="i_Kr" units="picoA_per_picoF" interface="public" id="rapid_time_dependent_potassium_current_i_Kr"/>
    <variable name="g_Kr" units="nanoS_per_picoF" initial_value="0.096"/>
    <variable name="E_K" units="millivolt" interface="public"/>
    <variable name="K_o" units="millimolar" interface="public"/>
    <variable name="time" units="millisecond" interface="public_and_private"/>
    <variable name="V" units="millivolt" interface="public_and_private"/>
    <variable name="Xr1" units="dimensionless" interface="private"/>
    <variable name="Xr2" units="dimensionless" interface="private"/>
    <math xmlns="http://www.w3.org/1998/Math/MathML" xmlns:cellml="http://www.cellml.org/cellml/2.0#">
      <apply>
        <eq/>
        <ci>i_Kr</ci>
        <apply>
          <times/>
          <ci>g_Kr</ci>
          <apply>
            <root/>
            <apply>
              <divide/>
              <ci>K_o</ci>
              <cn cellml:units="millimolar">5.4</cn>
            </apply>
          </apply>
          <ci>Xr1</ci>
          <ci>Xr2</ci>
          <apply>
            <minus/>
            <ci>V</ci>
            <ci>E_K</ci>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  <component name="rapid_time_dependent_potassium_current_Xr1_gate" id="id_00005">
    <variable name="Xr1" units="dimensionless" initial_value="0" interface="public"/>
    <variable name="V" units="millivolt" interface="public"/>
    <variable name="time" units="millisecond" interface="public"/>
    <variable name="xr1_inf" units="dimensionless"/>
    <variable name="alpha_xr1" units="dimensionless"/>
    <variable name="beta_xr1" units="dimensionless"/>
    <variable name="tau_xr1" units="millisecond"/>
    <math xmlns="http://www.w3.org/1998/Math/MathML" xmlns:cellml="http://www.cellml.org/cellml/2.0#">
      <apply>
        <eq/>
        <ci>xr1_inf</ci>
        <apply>
          <divide/>
          <cn cellml:units="dimensionless">1</cn>
          <apply>
            <plus/>
            <cn cellml:units="dimensionless">1</cn>
            <apply>
              <exp/>
              <apply>
                <divide/>
                <apply>
                  <minus/>
                  <apply>
                    <minus/>
                    <cn cellml:units="millivolt">26</cn>
                  </apply>
                  <ci>V</ci>
                </apply>
                <cn cellml:units="millivolt">7</cn>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>alpha_xr1</ci>
        <apply>
          <divide/>
          <cn cellml:units="dimensionless">450</cn>
          <apply>
            <plus/>
            <cn cellml:units="dimensionless">1</cn>
            <apply>
              <exp/>
              <apply>
                <divide/>
                <apply>
                  <minus/>
                  <apply>
                    <minus/>
                    <cn cellml:units="millivolt">45</cn>
                  </apply>
                  <ci>V</ci>
                </apply>
                <cn cellml:units="millivolt">10</cn>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>beta_xr1</ci>
        <apply>
          <divide/>
          <cn cellml:units="dimensionless">6</cn>
          <apply>
            <plus/>
            <cn cellml:units="dimensionless">1</cn>
            <apply>
              <exp/>
              <apply>
                <divide/>
                <apply>
                  <plus/>
                  <ci>V</ci>
                  <cn cellml:units="millivolt">30</cn>
                </apply>
                <cn cellml:units="millivolt">11.5</cn>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>tau_xr1</ci>
        <apply>
          <times/>
          <cn cellml:units="millisecond">1</cn>
          <ci>alpha_xr1</ci>
          <ci>beta_xr1</ci>
        </apply>
      </apply>
      <apply>
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci>time</ci>
          </bvar>
          <ci>Xr1</ci>
        </apply>
        <apply>
          <divide/>
          <apply>
            <minus/>
            <ci>xr1_inf</ci>
            <ci>Xr1</ci>
          </apply>
          <ci>tau_xr1</ci>
        </apply>
      </apply>
    </math>
  </component>
  <component name="rapid_time_dependent_potassium_current_Xr2_gate" id="id_00006">
    <variable name="Xr2" units="dimensionless" initial_value="1" interface="public"/>
    <variable name="V" units="millivolt" interface="public"/>
    <variable name="time" units="millisecond" interface="public"/>
    <variable name="xr2_inf" units="dimensionless"/>
    <variable name="alpha_xr2" units="dimensionless"/>
    <variable name="beta_xr2" units="dimensionless"/>
    <variable name="tau_xr2" units="millisecond"/>
    <math xmlns="http://www.w3.org/1998/Math/MathML" xmlns:cellml="http://www.cellml.org/cellml/2.0#">
      <apply>
        <eq/>
        <ci>xr2_inf</ci>
        <apply>
          <divide/>
          <cn cellml:units="dimensionless">1</cn>
          <apply>
            <plus/>
            <cn cellml:units="dimensionless">1</cn>
            <apply>
              <exp/>
              <apply>
                <divide/>
                <apply>
                  <plus/>
                  <ci>V</ci>
                  <cn cellml:units="millivolt">88</cn>
                </apply>
                <cn cellml:units="millivolt">24</cn>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>alpha_xr2</ci>
        <apply>
          <divide/>
          <cn cellml:units="dimensionless">3</cn>
          <apply>
            <plus/>
            <cn cellml:units="dimensionless">1</cn>
            <apply>
              <exp/>
              <apply>
                <divide/>
                <apply>
                  <minus/>
                  <apply>
                    <minus/>
                    <cn cellml:units="millivolt">60</cn>
                  </apply>
                  <ci>V</ci>
                </apply>
                <cn cellml:units="millivolt">20</cn>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>beta_xr2</ci>
        <apply>
          <divide/>
          <cn cellml:units="dimensionless">1.12</cn>
          <apply>
            <plus/>
            <cn cellml:units="dimensionless">1</cn>
            <apply>
              <exp/>
              <apply>
                <divide/>
                <apply>
                  <minus/>
                  <ci>V</ci>
                  <cn cellml:units="millivolt">60</cn>
                </apply>
                <cn cellml:units="millivolt">20</cn>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>tau_xr2</ci>
        <apply>
          <times/>
          <cn cellml:units="millisecond">1</cn>
          <ci>alpha_xr2</ci>
          <ci>beta_xr2</ci>
        </apply>
      </apply>
      <apply>
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci>time</ci>
          </bvar>
          <ci>Xr2</ci>
        </apply>
        <apply>
          <divide/>
          <apply>
            <minus/>
            <ci>xr2_inf</ci>
            <ci>Xr2</ci>
          </apply>
          <ci>tau_xr2</ci>
        </apply>
      </apply>
    </math>
  </component>
  <component name="slow_time_dependent_potassium_current" id="id_00007">
    <variable name="i_Ks" units="picoA_per_picoF" interface="public" id="slow_time_dependent_potassium_current_i_Ks"/>
    <variable name="g_Ks" units="nanoS_per_picoF" initial_value="0.062"/>
    <variable name="E_Ks" units="millivolt" interface="public"/>
    <variable name="time" units="millisecond" interface="public_and_private"/>
    <variable name="V" units="millivolt" interface="public_and_private"/>
    <variable name="Xs" units="dimensionless" interface="private"/>
    <math xmlns="http://www.w3.org/1998/Math/MathML" xmlns:cellml="http://www.cellml.org/cellml/2.0#">
      <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_time_dependent_potassium_current_Xs_gate" id="id_00008">
    <variable name="Xs" units="dimensionless" initial_value="0" interface="public"/>
    <variable name="V" units="millivolt" interface="public"/>
    <variable name="time" units="millisecond" interface="public"/>
    <variable name="xs_inf" units="dimensionless"/>
    <variable name="alpha_xs" units="dimensionless"/>
    <variable name="beta_xs" units="dimensionless"/>
    <variable name="tau_xs" units="millisecond"/>
    <math xmlns="http://www.w3.org/1998/Math/MathML" xmlns:cellml="http://www.cellml.org/cellml/2.0#">
      <apply>
        <eq/>
        <ci>xs_inf</ci>
        <apply>
          <divide/>
          <cn cellml:units="dimensionless">1</cn>
          <apply>
            <plus/>
            <cn cellml:units="dimensionless">1</cn>
            <apply>
              <exp/>
              <apply>
                <divide/>
                <apply>
                  <minus/>
                  <apply>
                    <minus/>
                    <cn cellml:units="millivolt">5</cn>
                  </apply>
                  <ci>V</ci>
                </apply>
                <cn cellml:units="millivolt">14</cn>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>alpha_xs</ci>
        <apply>
          <divide/>
          <cn cellml:units="dimensionless">1100</cn>
          <apply>
            <root/>
            <apply>
              <plus/>
              <cn cellml:units="dimensionless">1</cn>
              <apply>
                <exp/>
                <apply>
                  <divide/>
                  <apply>
                    <minus/>
                    <apply>
                      <minus/>
                      <cn cellml:units="millivolt">10</cn>
                    </apply>
                    <ci>V</ci>
                  </apply>
                  <cn cellml:units="millivolt">6</cn>
                </apply>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>beta_xs</ci>
        <apply>
          <divide/>
          <cn cellml:units="dimensionless">1</cn>
          <apply>
            <plus/>
            <cn cellml:units="dimensionless">1</cn>
            <apply>
              <exp/>
              <apply>
                <divide/>
                <apply>
                  <minus/>
                  <ci>V</ci>
                  <cn cellml:units="millivolt">60</cn>
                </apply>
                <cn cellml:units="millivolt">20</cn>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>tau_xs</ci>
        <apply>
          <times/>
          <cn cellml:units="millisecond">1</cn>
          <ci>alpha_xs</ci>
          <ci>beta_xs</ci>
        </apply>
      </apply>
      <apply>
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci>time</ci>
          </bvar>
          <ci>Xs</ci>
        </apply>
        <apply>
          <divide/>
          <apply>
            <minus/>
            <ci>xs_inf</ci>
            <ci>Xs</ci>
          </apply>
          <ci>tau_xs</ci>
        </apply>
      </apply>
    </math>
  </component>
  <component name="sodium_background_current" id="id_00013">
    <variable name="i_b_Na" units="picoA_per_picoF" interface="public" id="sodium_background_current_i_b_Na"/>
    <variable name="g_bna" units="nanoS_per_picoF" initial_value="0.00029"/>
    <variable name="E_Na" units="millivolt" interface="public"/>
    <variable name="time" units="millisecond" interface="public"/>
    <variable name="V" units="millivolt" interface="public"/>
    <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="calcium_background_current" id="id_00018">
    <variable name="i_b_Ca" units="picoA_per_picoF" interface="public" id="calcium_background_current_i_b_Ca"/>
    <variable name="g_bca" units="nanoS_per_picoF" initial_value="0.000592"/>
    <variable name="E_Ca" units="millivolt" interface="public"/>
    <variable name="time" units="millisecond" interface="public"/>
    <variable name="V" units="millivolt" interface="public"/>
    <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" id="id_00019">
    <variable name="i_to" units="picoA_per_picoF" interface="public" id="transient_outward_current_i_to"/>
    <variable name="g_to" units="nanoS_per_picoF" initial_value="0.294"/>
    <variable name="E_K" units="millivolt" interface="public"/>
    <variable name="V" units="millivolt" interface="public_and_private"/>
    <variable name="time" units="millisecond" interface="public_and_private"/>
    <variable name="s" units="dimensionless" interface="private"/>
    <variable name="r" units="dimensionless" interface="private"/>
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply>
        <eq/>
        <ci>i_to</ci>
        <apply>
          <times/>
          <ci>g_to</ci>
          <ci>r</ci>
          <ci>s</ci>
          <apply>
            <minus/>
            <ci>V</ci>
            <ci>E_K</ci>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  <component name="transient_outward_current_s_gate" id="id_00020">
    <variable name="s" units="dimensionless" initial_value="1" interface="public"/>
    <variable name="V" units="millivolt" interface="public"/>
    <variable name="time" units="millisecond" interface="public"/>
    <variable name="s_inf" units="dimensionless"/>
    <variable name="tau_s" units="millisecond"/>
    <math xmlns="http://www.w3.org/1998/Math/MathML" xmlns:cellml="http://www.cellml.org/cellml/2.0#">
      <apply>
        <eq/>
        <ci>s_inf</ci>
        <apply>
          <divide/>
          <cn cellml:units="dimensionless">1</cn>
          <apply>
            <plus/>
            <cn cellml:units="dimensionless">1</cn>
            <apply>
              <exp/>
              <apply>
                <divide/>
                <apply>
                  <plus/>
                  <ci>V</ci>
                  <cn cellml:units="millivolt">20</cn>
                </apply>
                <cn cellml:units="millivolt">5</cn>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>tau_s</ci>
        <apply>
          <plus/>
          <apply>
            <times/>
            <cn cellml:units="millisecond">85</cn>
            <apply>
              <exp/>
              <apply>
                <divide/>
                <apply>
                  <minus/>
                  <apply>
                    <power/>
                    <apply>
                      <plus/>
                      <ci>V</ci>
                      <cn cellml:units="millivolt">45</cn>
                    </apply>
                    <cn cellml:units="dimensionless">2</cn>
                  </apply>
                </apply>
                <cn cellml:units="millivolt_square">320</cn>
              </apply>
            </apply>
          </apply>
          <apply>
            <divide/>
            <cn cellml:units="millisecond">5</cn>
            <apply>
              <plus/>
              <cn cellml:units="dimensionless">1</cn>
              <apply>
                <exp/>
                <apply>
                  <divide/>
                  <apply>
                    <minus/>
                    <ci>V</ci>
                    <cn cellml:units="millivolt">20</cn>
                  </apply>
                  <cn cellml:units="millivolt">5</cn>
                </apply>
              </apply>
            </apply>
          </apply>
          <cn cellml:units="millisecond">3</cn>
        </apply>
      </apply>
      <apply>
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci>time</ci>
          </bvar>
          <ci>s</ci>
        </apply>
        <apply>
          <divide/>
          <apply>
            <minus/>
            <ci>s_inf</ci>
            <ci>s</ci>
          </apply>
          <ci>tau_s</ci>
        </apply>
      </apply>
    </math>
  </component>
  <component name="transient_outward_current_r_gate" id="id_00021">
    <variable name="r" units="dimensionless" initial_value="0" interface="public"/>
    <variable name="V" units="millivolt" interface="public"/>
    <variable name="time" units="millisecond" interface="public"/>
    <variable name="r_inf" units="dimensionless"/>
    <variable name="tau_r" units="millisecond"/>
    <math xmlns="http://www.w3.org/1998/Math/MathML" xmlns:cellml="http://www.cellml.org/cellml/2.0#">
      <apply>
        <eq/>
        <ci>r_inf</ci>
        <apply>
          <divide/>
          <cn cellml:units="dimensionless">1</cn>
          <apply>
            <plus/>
            <cn cellml:units="dimensionless">1</cn>
            <apply>
              <exp/>
              <apply>
                <divide/>
                <apply>
                  <minus/>
                  <cn cellml:units="millivolt">20</cn>
                  <ci>V</ci>
                </apply>
                <cn cellml:units="millivolt">6</cn>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>tau_r</ci>
        <apply>
          <plus/>
          <apply>
            <times/>
            <cn cellml:units="millisecond">9.5</cn>
            <apply>
              <exp/>
              <apply>
                <divide/>
                <apply>
                  <minus/>
                  <apply>
                    <power/>
                    <apply>
                      <plus/>
                      <ci>V</ci>
                      <cn cellml:units="millivolt">40</cn>
                    </apply>
                    <cn cellml:units="dimensionless">2</cn>
                  </apply>
                </apply>
                <cn cellml:units="millivolt_square">1800</cn>
              </apply>
            </apply>
          </apply>
          <cn cellml:units="millisecond">0.8</cn>
        </apply>
      </apply>
      <apply>
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci>time</ci>
          </bvar>
          <ci>r</ci>
        </apply>
        <apply>
          <divide/>
          <apply>
            <minus/>
            <ci>r_inf</ci>
            <ci>r</ci>
          </apply>
          <ci>tau_r</ci>
        </apply>
      </apply>
    </math>
  </component>
  <component name="sodium_potassium_pump_current" id="id_00022">
    <variable name="i_NaK" units="picoA_per_picoF" interface="public" id="sodium_potassium_pump_current_i_NaK"/>
    <variable name="P_NaK" units="picoA_per_picoF" initial_value="1.362"/>
    <variable name="K_mk" units="millimolar" initial_value="1"/>
    <variable name="K_mNa" units="millimolar" initial_value="40"/>
    <variable name="K_o" units="millimolar" interface="public"/>
    <variable name="Na_i" units="millimolar" interface="public"/>
    <variable name="time" units="millisecond" interface="public"/>
    <variable name="R" units="joule_per_mole_kelvin" interface="public"/>
    <variable name="F" units="coulomb_per_millimole" interface="public"/>
    <variable name="T" units="kelvin" interface="public"/>
    <variable name="V" units="millivolt" interface="public_and_private"/>
    <math xmlns="http://www.w3.org/1998/Math/MathML" xmlns:cellml="http://www.cellml.org/cellml/2.0#">
      <apply>
        <eq/>
        <ci>i_NaK</ci>
        <apply>
          <divide/>
          <apply>
            <divide/>
            <apply>
              <times/>
              <apply>
                <divide/>
                <apply>
                  <times/>
                  <ci>P_NaK</ci>
                  <ci>K_o</ci>
                </apply>
                <apply>
                  <plus/>
                  <ci>K_o</ci>
                  <ci>K_mk</ci>
                </apply>
              </apply>
              <ci>Na_i</ci>
            </apply>
            <apply>
              <plus/>
              <ci>Na_i</ci>
              <ci>K_mNa</ci>
            </apply>
          </apply>
          <apply>
            <plus/>
            <cn cellml:units="dimensionless">1</cn>
            <apply>
              <times/>
              <cn cellml:units="dimensionless">0.1245</cn>
              <apply>
                <exp/>
                <apply>
                  <divide/>
                  <apply>
                    <times/>
                    <apply>
                      <minus/>
                      <cn cellml:units="dimensionless">0.1</cn>
                    </apply>
                    <ci>V</ci>
                    <ci>F</ci>
                  </apply>
                  <apply>
                    <times/>
                    <ci>R</ci>
                    <ci>T</ci>
                  </apply>
                </apply>
              </apply>
            </apply>
            <apply>
              <times/>
              <cn cellml:units="dimensionless">0.0353</cn>
              <apply>
                <exp/>
                <apply>
                  <divide/>
                  <apply>
                    <times/>
                    <apply>
                      <minus/>
                      <ci>V</ci>
                    </apply>
                    <ci>F</ci>
                  </apply>
                  <apply>
                    <times/>
                    <ci>R</ci>
                    <ci>T</ci>
                  </apply>
                </apply>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  <component name="sodium_calcium_exchanger_current" id="id_00023">
    <variable name="i_NaCa" units="picoA_per_picoF" interface="public" id="sodium_calcium_exchanger_current_i_NaCa"/>
    <variable name="K_NaCa" units="picoA_per_picoF" initial_value="1000"/>
    <variable name="K_sat" units="dimensionless" initial_value="0.1"/>
    <variable name="alpha" units="dimensionless" initial_value="2.5"/>
    <variable name="gamma" units="dimensionless" initial_value="0.35"/>
    <variable name="Km_Ca" units="millimolar" initial_value="1.38"/>
    <variable name="Km_Nai" units="millimolar" initial_value="87.5"/>
    <variable name="Na_i" units="millimolar" interface="public"/>
    <variable name="Na_o" units="millimolar" interface="public"/>
    <variable name="Ca_i" units="millimolar" interface="public"/>
    <variable name="Ca_o" units="millimolar" interface="public"/>
    <variable name="R" units="joule_per_mole_kelvin" interface="public"/>
    <variable name="F" units="coulomb_per_millimole" interface="public"/>
    <variable name="T" units="kelvin" interface="public"/>
    <variable name="V" units="millivolt" interface="public"/>
    <variable name="time" units="millisecond" interface="public"/>
    <math xmlns="http://www.w3.org/1998/Math/MathML" xmlns:cellml="http://www.cellml.org/cellml/2.0#">
      <apply>
        <eq/>
        <ci>i_NaCa</ci>
        <apply>
          <divide/>
          <apply>
            <times/>
            <ci>K_NaCa</ci>
            <apply>
              <minus/>
              <apply>
                <times/>
                <apply>
                  <exp/>
                  <apply>
                    <divide/>
                    <apply>
                      <times/>
                      <ci>gamma</ci>
                      <ci>V</ci>
                      <ci>F</ci>
                    </apply>
                    <apply>
                      <times/>
                      <ci>R</ci>
                      <ci>T</ci>
                    </apply>
                  </apply>
                </apply>
                <apply>
                  <power/>
                  <ci>Na_i</ci>
                  <cn cellml:units="dimensionless">3</cn>
                </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>
                      <ci>V</ci>
                      <ci>F</ci>
                    </apply>
                    <apply>
                      <times/>
                      <ci>R</ci>
                      <ci>T</ci>
                    </apply>
                  </apply>
                </apply>
                <apply>
                  <power/>
                  <ci>Na_o</ci>
                  <cn cellml:units="dimensionless">3</cn>
                </apply>
                <ci>Ca_i</ci>
                <ci>alpha</ci>
              </apply>
            </apply>
          </apply>
          <apply>
            <times/>
            <apply>
              <plus/>
              <apply>
                <power/>
                <ci>Km_Nai</ci>
                <cn cellml:units="dimensionless">3</cn>
              </apply>
              <apply>
                <power/>
                <ci>Na_o</ci>
                <cn cellml:units="dimensionless">3</cn>
              </apply>
            </apply>
            <apply>
              <plus/>
              <ci>Km_Ca</ci>
              <ci>Ca_o</ci>
            </apply>
            <apply>
              <plus/>
              <cn cellml:units="dimensionless">1</cn>
              <apply>
                <times/>
                <ci>K_sat</ci>
                <apply>
                  <exp/>
                  <apply>
                    <divide/>
                    <apply>
                      <times/>
                      <apply>
                        <minus/>
                        <ci>gamma</ci>
                        <cn cellml:units="dimensionless">1</cn>
                      </apply>
                      <ci>V</ci>
                      <ci>F</ci>
                    </apply>
                    <apply>
                      <times/>
                      <ci>R</ci>
                      <ci>T</ci>
                    </apply>
                  </apply>
                </apply>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  <component name="calcium_pump_current" id="id_00027">
    <variable name="i_p_Ca" units="picoA_per_picoF" interface="public" id="calcium_pump_current_i_p_Ca"/>
    <variable name="g_pCa" units="picoA_per_picoF" initial_value="0.825"/>
    <variable name="K_pCa" units="millimolar" initial_value="0.0005"/>
    <variable name="Ca_i" units="millimolar" interface="public"/>
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply>
        <eq/>
        <ci>i_p_Ca</ci>
        <apply>
          <divide/>
          <apply>
            <times/>
            <ci>g_pCa</ci>
            <ci>Ca_i</ci>
          </apply>
          <apply>
            <plus/>
            <ci>Ca_i</ci>
            <ci>K_pCa</ci>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  <component name="potassium_pump_current" id="id_00069">
    <variable name="i_p_K" units="picoA_per_picoF" interface="public" id="potassium_pump_current_i_p_K"/>
    <variable name="g_pK" units="nanoS_per_picoF" initial_value="0.0146"/>
    <variable name="E_K" units="millivolt" interface="public"/>
    <variable name="V" units="millivolt" interface="public"/>
    <math xmlns="http://www.w3.org/1998/Math/MathML" xmlns:cellml="http://www.cellml.org/cellml/2.0#">
      <apply>
        <eq/>
        <ci>i_p_K</ci>
        <apply>
          <divide/>
          <apply>
            <times/>
            <ci>g_pK</ci>
            <apply>
              <minus/>
              <ci>V</ci>
              <ci>E_K</ci>
            </apply>
          </apply>
          <apply>
            <plus/>
            <cn cellml:units="dimensionless">1</cn>
            <apply>
              <exp/>
              <apply>
                <divide/>
                <apply>
                  <minus/>
                  <cn cellml:units="millivolt">25</cn>
                  <ci>V</ci>
                </apply>
                <cn cellml:units="millivolt">5.98</cn>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  <component name="calcium_dynamics" id="id_00024">
    <variable name="Ca_i" units="millimolar" initial_value="0.0002" interface="public"/>
    <variable name="Ca_o" units="millimolar" initial_value="2" interface="public"/>
    <variable name="Ca_SR" units="millimolar" initial_value="0.2"/>
    <variable name="i_rel" units="millimolar_per_millisecond" id="calcium_dynamics_i_rel"/>
    <variable name="i_up" units="millimolar_per_millisecond" id="calcium_dynamics_i_up"/>
    <variable name="i_leak" units="millimolar_per_millisecond" id="calcium_dynamics_i_leak"/>
    <variable name="g" units="dimensionless" initial_value="1"/>
    <variable name="tau_g" units="millisecond" initial_value="2"/>
    <variable name="g_inf" units="dimensionless"/>
    <variable name="a_rel" units="millimolar_per_millisecond" initial_value="0.016464"/>
    <variable name="b_rel" units="millimolar" initial_value="0.25"/>
    <variable name="c_rel" units="millimolar_per_millisecond" initial_value="0.008232"/>
    <variable name="K_up" units="millimolar" initial_value="0.00025"/>
    <variable name="V_leak" units="per_millisecond" initial_value="8e-5"/>
    <variable name="Vmax_up" units="millimolar_per_millisecond" initial_value="0.000425"/>
    <variable name="ddt_Ca_i_total" units="millimolar_per_millisecond"/>
    <variable name="ddt_Ca_sr_total" units="millimolar_per_millisecond"/>
    <variable name="f_JCa_i_free" units="dimensionless"/>
    <variable name="f_JCa_sr_free" units="dimensionless"/>
    <variable name="Buf_c" units="millimolar" initial_value="0.15"/>
    <variable name="K_buf_c" units="millimolar" initial_value="0.001"/>
    <variable name="Buf_sr" units="millimolar" initial_value="10"/>
    <variable name="K_buf_sr" units="millimolar" initial_value="0.3"/>
    <variable name="V_sr" units="mL" initial_value="1094"/>
    <variable name="V" units="millivolt" interface="public"/>
    <variable name="V_c" units="mL" interface="public"/>
    <variable name="time" units="millisecond" interface="public"/>
    <variable name="d" units="dimensionless" interface="public"/>
    <variable name="F" units="coulomb_per_millimole" interface="public"/>
    <variable name="Cm" units="picoF" interface="public"/>
    <variable name="i_CaL" units="picoA_per_picoF" interface="public"/>
    <variable name="i_NaCa" units="picoA_per_picoF" interface="public"/>
    <variable name="i_p_Ca" units="picoA_per_picoF" interface="public"/>
    <variable name="i_b_Ca" units="picoA_per_picoF" interface="public"/>
    <variable name="d_g" units="per_millisecond"/>
    <math xmlns="http://www.w3.org/1998/Math/MathML" xmlns:cellml="http://www.cellml.org/cellml/2.0#">
      <apply>
        <eq/>
        <ci>i_rel</ci>
        <apply>
          <times/>
          <apply>
            <plus/>
            <apply>
              <divide/>
              <apply>
                <times/>
                <ci>a_rel</ci>
                <apply>
                  <power/>
                  <ci>Ca_SR</ci>
                  <cn cellml:units="dimensionless">2</cn>
                </apply>
              </apply>
              <apply>
                <plus/>
                <apply>
                  <power/>
                  <ci>b_rel</ci>
                  <cn cellml:units="dimensionless">2</cn>
                </apply>
                <apply>
                  <power/>
                  <ci>Ca_SR</ci>
                  <cn cellml:units="dimensionless">2</cn>
                </apply>
              </apply>
            </apply>
            <ci>c_rel</ci>
          </apply>
          <ci>d</ci>
          <ci>g</ci>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>i_up</ci>
        <apply>
          <divide/>
          <ci>Vmax_up</ci>
          <apply>
            <plus/>
            <cn cellml:units="dimensionless">1</cn>
            <apply>
              <divide/>
              <apply>
                <power/>
                <ci>K_up</ci>
                <cn cellml:units="dimensionless">2</cn>
              </apply>
              <apply>
                <power/>
                <ci>Ca_i</ci>
                <cn cellml:units="dimensionless">2</cn>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>i_leak</ci>
        <apply>
          <times/>
          <ci>V_leak</ci>
          <apply>
            <minus/>
            <ci>Ca_SR</ci>
            <ci>Ca_i</ci>
          </apply>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>g_inf</ci>
        <piecewise>
          <piece>
            <apply>
              <divide/>
              <cn cellml:units="dimensionless">1</cn>
              <apply>
                <plus/>
                <cn cellml:units="dimensionless">1</cn>
                <apply>
                  <power/>
                  <apply>
                    <divide/>
                    <ci>Ca_i</ci>
                    <cn cellml:units="millimolar">0.00035</cn>
                  </apply>
                  <cn cellml:units="dimensionless">6</cn>
                </apply>
              </apply>
            </apply>
            <apply>
              <lt/>
              <ci>Ca_i</ci>
              <cn cellml:units="millimolar">0.00035</cn>
            </apply>
          </piece>
          <otherwise>
            <apply>
              <divide/>
              <cn cellml:units="dimensionless">1</cn>
              <apply>
                <plus/>
                <cn cellml:units="dimensionless">1</cn>
                <apply>
                  <power/>
                  <apply>
                    <divide/>
                    <ci>Ca_i</ci>
                    <cn cellml:units="millimolar">0.00035</cn>
                  </apply>
                  <cn cellml:units="dimensionless">16</cn>
                </apply>
              </apply>
            </apply>
          </otherwise>
        </piecewise>
      </apply>
      <apply>
        <eq/>
        <ci>d_g</ci>
        <apply>
          <divide/>
          <apply>
            <minus/>
            <ci>g_inf</ci>
            <ci>g</ci>
          </apply>
          <ci>tau_g</ci>
        </apply>
      </apply>
      <apply>
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci>time</ci>
          </bvar>
          <ci>g</ci>
        </apply>
        <piecewise>
          <piece>
            <cn cellml:units="per_millisecond">0</cn>
            <apply>
              <and/>
              <apply>
                <gt/>
                <ci>g_inf</ci>
                <ci>g</ci>
              </apply>
              <apply>
                <gt/>
                <ci>V</ci>
                <apply>
                  <minus/>
                  <cn cellml:units="millivolt">60</cn>
                </apply>
              </apply>
            </apply>
          </piece>
          <otherwise>
            <ci>d_g</ci>
          </otherwise>
        </piecewise>
      </apply>
      <apply>
        <eq/>
        <ci>ddt_Ca_i_total</ci>
        <apply>
          <plus/>
          <apply>
            <minus/>
            <apply>
              <plus/>
              <apply>
                <times/>
                <apply>
                  <divide/>
                  <apply>
                    <minus/>
                    <apply>
                      <minus/>
                      <apply>
                        <plus/>
                        <ci>i_CaL</ci>
                        <ci>i_b_Ca</ci>
                        <ci>i_p_Ca</ci>
                      </apply>
                      <apply>
                        <times/>
                        <cn cellml:units="dimensionless">2</cn>
                        <ci>i_NaCa</ci>
                      </apply>
                    </apply>
                  </apply>
                  <apply>
                    <times/>
                    <cn cellml:units="dimensionless">2</cn>
                    <ci>V_c</ci>
                    <ci>F</ci>
                  </apply>
                </apply>
                <ci>Cm</ci>
              </apply>
              <ci>i_leak</ci>
            </apply>
            <ci>i_up</ci>
          </apply>
          <ci>i_rel</ci>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>ddt_Ca_sr_total</ci>
        <apply>
          <times/>
          <apply>
            <divide/>
            <ci>V_c</ci>
            <ci>V_sr</ci>
          </apply>
          <apply>
            <minus/>
            <ci>i_up</ci>
            <apply>
              <plus/>
              <ci>i_rel</ci>
              <ci>i_leak</ci>
            </apply>
          </apply>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>f_JCa_i_free</ci>
        <apply>
          <divide/>
          <cn cellml:units="dimensionless">1</cn>
          <apply>
            <plus/>
            <cn cellml:units="dimensionless">1</cn>
            <apply>
              <divide/>
              <apply>
                <times/>
                <ci>Buf_c</ci>
                <ci>K_buf_c</ci>
              </apply>
              <apply>
                <power/>
                <apply>
                  <plus/>
                  <ci>Ca_i</ci>
                  <ci>K_buf_c</ci>
                </apply>
                <cn cellml:units="dimensionless">2</cn>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>f_JCa_sr_free</ci>
        <apply>
          <divide/>
          <cn cellml:units="dimensionless">1</cn>
          <apply>
            <plus/>
            <cn cellml:units="dimensionless">1</cn>
            <apply>
              <divide/>
              <apply>
                <times/>
                <ci>Buf_sr</ci>
                <ci>K_buf_sr</ci>
              </apply>
              <apply>
                <power/>
                <apply>
                  <plus/>
                  <ci>Ca_SR</ci>
                  <ci>K_buf_sr</ci>
                </apply>
                <cn cellml:units="dimensionless">2</cn>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      <apply>
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci>time</ci>
          </bvar>
          <ci>Ca_i</ci>
        </apply>
        <apply>
          <times/>
          <ci>ddt_Ca_i_total</ci>
          <ci>f_JCa_i_free</ci>
        </apply>
      </apply>
      <apply>
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci>time</ci>
          </bvar>
          <ci>Ca_SR</ci>
        </apply>
        <apply>
          <times/>
          <ci>ddt_Ca_sr_total</ci>
          <ci>f_JCa_sr_free</ci>
        </apply>
      </apply>
    </math>
  </component>
  <component name="sodium_dynamics" id="id_00025">
    <variable name="Na_i" units="millimolar" initial_value="11.6" interface="public"/>
    <variable name="Na_o" units="millimolar" initial_value="140" interface="public"/>
    <variable name="time" units="millisecond" interface="public"/>
    <variable name="F" units="coulomb_per_millimole" interface="public"/>
    <variable name="Cm" units="picoF" interface="public"/>
    <variable name="V_c" units="mL" interface="public"/>
    <variable name="i_Na" units="picoA_per_picoF" interface="public"/>
    <variable name="i_NaCa" units="picoA_per_picoF" interface="public"/>
    <variable name="i_NaK" units="picoA_per_picoF" interface="public"/>
    <variable name="i_b_Na" units="picoA_per_picoF" interface="public"/>
    <math xmlns="http://www.w3.org/1998/Math/MathML" xmlns:cellml="http://www.cellml.org/cellml/2.0#">
      <apply>
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci>time</ci>
          </bvar>
          <ci>Na_i</ci>
        </apply>
        <apply>
          <divide/>
          <apply>
            <times/>
            <apply>
              <minus/>
              <apply>
                <plus/>
                <ci>i_Na</ci>
                <ci>i_b_Na</ci>
                <apply>
                  <times/>
                  <cn cellml:units="dimensionless">3</cn>
                  <ci>i_NaK</ci>
                </apply>
                <apply>
                  <times/>
                  <cn cellml:units="dimensionless">3</cn>
                  <ci>i_NaCa</ci>
                </apply>
              </apply>
            </apply>
            <ci>Cm</ci>
          </apply>
          <apply>
            <times/>
            <ci>V_c</ci>
            <ci>F</ci>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  <component name="potassium_dynamics" id="id_00026">
    <variable name="K_i" units="millimolar" initial_value="138.3" interface="public"/>
    <variable name="K_o" units="millimolar" initial_value="5.4" interface="public"/>
    <variable name="time" units="millisecond" interface="public"/>
    <variable name="F" units="coulomb_per_millimole" interface="public"/>
    <variable name="Cm" units="picoF" interface="public"/>
    <variable name="V_c" units="mL" interface="public"/>
    <variable name="i_K1" units="picoA_per_picoF" interface="public"/>
    <variable name="i_to" units="picoA_per_picoF" interface="public"/>
    <variable name="i_NaK" units="picoA_per_picoF" interface="public"/>
    <variable name="i_Kr" units="picoA_per_picoF" interface="public"/>
    <variable name="i_Ks" units="picoA_per_picoF" interface="public"/>
    <variable name="i_p_K" units="picoA_per_picoF" interface="public"/>
    <variable name="i_Stim" units="picoA_per_picoF" interface="public"/>
    <math xmlns="http://www.w3.org/1998/Math/MathML" xmlns:cellml="http://www.cellml.org/cellml/2.0#">
      <apply>
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci>time</ci>
          </bvar>
          <ci>K_i</ci>
        </apply>
        <apply>
          <divide/>
          <apply>
            <times/>
            <apply>
              <minus/>
              <apply>
                <minus/>
                <apply>
                  <plus/>
                  <ci>i_K1</ci>
                  <ci>i_to</ci>
                  <ci>i_Kr</ci>
                  <ci>i_Ks</ci>
                  <ci>i_p_K</ci>
                  <ci>i_Stim</ci>
                </apply>
                <apply>
                  <times/>
                  <cn cellml:units="dimensionless">2</cn>
                  <ci>i_NaK</ci>
                </apply>
              </apply>
            </apply>
            <ci>Cm</ci>
          </apply>
          <apply>
            <times/>
            <ci>V_c</ci>
            <ci>F</ci>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  <component name="reversal_potentials" id="id_00002">
    <variable name="E_Na" units="millivolt" interface="public"/>
    <variable name="E_K" units="millivolt" interface="public"/>
    <variable name="E_Ks" units="millivolt" interface="public"/>
    <variable name="E_Ca" units="millivolt" interface="public"/>
    <variable name="P_kna" units="dimensionless" initial_value="0.03"/>
    <variable name="K_o" units="millimolar" interface="public"/>
    <variable name="Na_o" units="millimolar" interface="public"/>
    <variable name="K_i" units="millimolar" interface="public"/>
    <variable name="Na_i" units="millimolar" interface="public"/>
    <variable name="Ca_o" units="millimolar" interface="public"/>
    <variable name="Ca_i" units="millimolar" interface="public"/>
    <variable name="R" units="joule_per_mole_kelvin" interface="public"/>
    <variable name="F" units="coulomb_per_millimole" interface="public"/>
    <variable name="T" units="kelvin" interface="public"/>
    <variable name="time" units="millisecond" interface="public"/>
    <math xmlns="http://www.w3.org/1998/Math/MathML" xmlns:cellml="http://www.cellml.org/cellml/2.0#">
      <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>
    </math>
  </component>
  <component name="inward_rectifier_potassium_current" id="id_00003">
    <variable name="i_K1" units="picoA_per_picoF" interface="public" id="inward_rectifier_potassium_current_i_K1"/>
    <variable name="g_K1" units="nanoS_per_picoF" initial_value="5.405"/>
    <variable name="V" units="millivolt" interface="public"/>
    <variable name="K_o" units="millimolar" interface="public"/>
    <variable name="E_K" units="millivolt" interface="public"/>
    <variable name="time" units="millisecond" interface="public"/>
    <variable name="xK1_inf" units="dimensionless"/>
    <variable name="alpha_K1" units="per_millisecond"/>
    <variable name="beta_K1" units="per_millisecond"/>
    <math xmlns="http://www.w3.org/1998/Math/MathML" xmlns:cellml="http://www.cellml.org/cellml/2.0#">
      <apply>
        <eq/>
        <ci>alpha_K1</ci>
        <apply>
          <divide/>
          <cn cellml:units="per_millisecond">0.1</cn>
          <apply>
            <plus/>
            <cn cellml:units="dimensionless">1</cn>
            <apply>
              <exp/>
              <apply>
                <times/>
                <cn cellml:units="per_millivolt">0.06</cn>
                <apply>
                  <minus/>
                  <apply>
                    <minus/>
                    <ci>V</ci>
                    <ci>E_K</ci>
                  </apply>
                  <cn cellml:units="millivolt">200</cn>
                </apply>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>beta_K1</ci>
        <apply>
          <divide/>
          <apply>
            <plus/>
            <apply>
              <times/>
              <cn cellml:units="per_millisecond">3</cn>
              <apply>
                <exp/>
                <apply>
                  <times/>
                  <cn cellml:units="per_millivolt">0.0002</cn>
                  <apply>
                    <plus/>
                    <apply>
                      <minus/>
                      <ci>V</ci>
                      <ci>E_K</ci>
                    </apply>
                    <cn cellml:units="millivolt">100</cn>
                  </apply>
                </apply>
              </apply>
            </apply>
            <apply>
              <times/>
              <cn cellml:units="per_millisecond">1</cn>
              <apply>
                <exp/>
                <apply>
                  <times/>
                  <cn cellml:units="per_millivolt">0.1</cn>
                  <apply>
                    <minus/>
                    <apply>
                      <minus/>
                      <ci>V</ci>
                      <ci>E_K</ci>
                    </apply>
                    <cn cellml:units="millivolt">10</cn>
                  </apply>
                </apply>
              </apply>
            </apply>
          </apply>
          <apply>
            <plus/>
            <cn cellml:units="dimensionless">1</cn>
            <apply>
              <exp/>
              <apply>
                <times/>
                <apply>
                  <minus/>
                  <cn cellml:units="per_millivolt">0.5</cn>
                </apply>
                <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>E_K</ci>
                </apply>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>xK1_inf</ci>
        <apply>
          <divide/>
          <ci>alpha_K1</ci>
          <apply>
            <plus/>
            <ci>alpha_K1</ci>
            <ci>beta_K1</ci>
          </apply>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>i_K1</ci>
        <apply>
          <times/>
          <ci>g_K1</ci>
          <ci>xK1_inf</ci>
          <apply>
            <root/>
            <apply>
              <divide/>
              <ci>K_o</ci>
              <cn cellml:units="millimolar">5.4</cn>
            </apply>
          </apply>
          <apply>
            <minus/>
            <ci>V</ci>
            <ci>E_K</ci>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  <component name="main">
    <variable name="phi_Cai" units="dimensionless" interface="public"/>
    <variable name="time" units="millisecond" interface="public"/>
    <variable name="phi" units="dimensionless" interface="public"/>
    <variable name="u_pv" units="mmHg" interface="public"/>
    <variable name="u_sa" units="mmHg" interface="public"/>
    <variable name="u_pv0" units="mmHg" interface="public"/>
    <variable name="u_sa0" units="mmHg" interface="public"/>
    <variable name="u_rv" units="mmHg" interface="public"/>
    <variable name="u_lv" units="mmHg" interface="public"/>
    <variable name="u_rv0" units="mmHg" interface="public"/>
    <variable name="u_lv0" units="mmHg" interface="public"/>
  </component>
  <component name="TT04">
    <variable name="Ca_i" units="millimolar" interface="public"/>
    <variable name="K_Cai" units="per_mmHg" initial_value="0.2"/>
    <variable name="phi_Cai" units="dimensionless" interface="public"/>
    <variable name="time" units="second" interface="public"/>
    <variable name="T_ref" units="mmHg" initial_value="125.0"/>
    <variable name="rho_0" units="per_s_per_mM" initial_value="100.0e3"/>
    <variable name="rho_1" units="per_s" initial_value="163.0"/>
    <variable name="Ca_b_max" units="millimolar" initial_value="2.26e-3"/>
    <variable name="pC_50_ref" units="dimensionless" initial_value="6"/>
    <variable name="n_ref" units="dimensionless" initial_value="4.25"/>
    <variable name="alpha_0" units="per_s" initial_value="2.0"/>
    <variable name="beta_0" units="dimensionless" initial_value="1.45"/>
    <variable name="beta_1" units="dimensionless" initial_value="1.95"/>
    <variable name="beta_2" units="dimensionless" initial_value="0.31"/>
    <variable name="lambda" units="dimensionless" initial_value="1"/>
    <variable name="C_50" units="millimolar"/>
    <variable name="pC_50" units="dimensionless"/>
    <variable name="n" units="dimensionless"/>
    <variable name="T_o" units="mmHg"/>
    <variable name="Ca_b" units="millimolar" initial_value="0"/>
    <variable name="z" units="dimensionless" initial_value="0"/>
    <math xmlns="http://www.w3.org/1998/Math/MathML" xmlns:cellml="http://www.cellml.org/cellml/2.0#">
      <apply>
        <eq/>
        <ci>phi_Cai</ci>
        <apply>
          <times/>
          <apply>
            <minus/>
            <ci>T_o</ci>
            <cn cellml:units="mmHg">1.2</cn>
          </apply>
          <ci>K_Cai</ci>
        </apply>
      </apply>
    </math>
    <math xmlns="http://www.w3.org/1998/Math/MathML" xmlns:cellml="http://www.cellml.org/cellml/2.0#">
      <apply>
        <eq/>
        <ci>pC_50</ci>
        <apply>
          <times/>
          <ci>pC_50_ref</ci>
          <apply>
            <plus/>
            <cn cellml:units="dimensionless">1.0</cn>
            <apply>
              <times/>
              <ci>beta_2</ci>
              <apply>
                <minus/>
                <ci>lambda</ci>
                <cn cellml:units="dimensionless">1.0</cn>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>C_50</ci>
        <apply>
          <times/>
          <cn type="e-notation" cellml:units="millimolar">1.0<sep/>-3</cn>
          <apply>
            <power/>
            <cn cellml:units="dimensionless">10</cn>
            <apply>
              <minus/>
              <cn cellml:units="dimensionless">6</cn>
              <ci>pC_50</ci>
            </apply>
          </apply>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>n</ci>
        <apply>
          <times/>
          <ci>n_ref</ci>
          <apply>
            <plus/>
            <cn cellml:units="dimensionless">1.0</cn>
            <apply>
              <times/>
              <ci>beta_1</ci>
              <apply>
                <minus/>
                <ci>lambda</ci>
                <cn cellml:units="dimensionless">1.0</cn>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      <apply>
        <eq/>
        <ci>T_o</ci>
        <apply>
          <times/>
          <ci>T_ref</ci>
          <apply>
            <plus/>
            <cn cellml:units="dimensionless">1.0</cn>
            <apply>
              <times/>
              <ci>beta_0</ci>
              <apply>
                <minus/>
                <ci>lambda</ci>
                <cn cellml:units="dimensionless">1.0</cn>
              </apply>
            </apply>
          </apply>
          <ci>z</ci>
        </apply>
      </apply>
      <apply>
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci>time</ci>
          </bvar>
          <ci>Ca_b</ci>
        </apply>
        <apply>
          <minus/>
          <apply>
            <times/>
            <ci>rho_0</ci>
            <ci>Ca_i</ci>
            <apply>
              <minus/>
              <ci>Ca_b_max</ci>
              <ci>Ca_b</ci>
            </apply>
          </apply>
          <apply>
            <times/>
            <ci>rho_1</ci>
            <ci>Ca_b</ci>
          </apply>
        </apply>
      </apply>
      <apply>
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci>time</ci>
          </bvar>
          <ci>z</ci>
        </apply>
        <apply>
          <times/>
          <ci>alpha_0</ci>
          <apply>
            <minus/>
            <apply>
              <times/>
              <apply>
                <power/>
                <apply>
                  <divide/>
                  <ci>Ca_b</ci>
                  <ci>C_50</ci>
                </apply>
                <ci>n</ci>
              </apply>
              <apply>
                <minus/>
                <cn cellml:units="dimensionless">1.0</cn>
                <ci>z</ci>
              </apply>
            </apply>
            <ci>z</ci>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  <connection component_1="cvs" component_2="main">
    <map_variables variable_1="t" variable_2="time"/>
    <map_variables variable_1="u_rv" variable_2="u_rv"/>
    <map_variables variable_1="u_lv" variable_2="u_lv"/>
    <map_variables variable_1="u_pv" variable_2="u_pv"/>
    <map_variables variable_1="u_sa" variable_2="u_sa"/>
  </connection>
  <connection component_1="cvs" component_2="TT04">
    <map_variables variable_1="phi" variable_2="phi_Cai"/>
  </connection>
  <connection component_1="cvs0" component_2="main">
    <map_variables variable_1="t" variable_2="time"/>
    <map_variables variable_1="u_rv" variable_2="u_rv0"/>
    <map_variables variable_1="u_lv" variable_2="u_lv0"/>
    <map_variables variable_1="u_pv" variable_2="u_pv0"/>
    <map_variables variable_1="u_sa" variable_2="u_sa0"/>
  </connection>
  <connection component_1="cvs0" component_2="activation">
    <map_variables variable_1="phi" variable_2="phi"/>
  </connection>
  <connection component_1="activation" component_2="main">
    <map_variables variable_1="t" variable_2="time"/>
    <map_variables variable_1="phi" variable_2="phi"/>
  </connection>
  <connection component_1="activation" component_2="membrane">
    <map_variables variable_1="T_0" variable_2="stim_period"/>
  </connection>
  <connection component_1="environment" component_2="main">
    <map_variables variable_1="time" variable_2="time"/>
  </connection>
  <connection component_1="environment" component_2="membrane">
    <map_variables variable_1="time" variable_2="time"/>
  </connection>
  <connection component_1="environment" component_2="fast_sodium_current">
    <map_variables variable_1="time" variable_2="time"/>
  </connection>
  <connection component_1="environment" component_2="L_type_Ca_current">
    <map_variables variable_1="time" variable_2="time"/>
  </connection>
  <connection component_1="environment" component_2="rapid_time_dependent_potassium_current">
    <map_variables variable_1="time" variable_2="time"/>
  </connection>
  <connection component_1="environment" component_2="slow_time_dependent_potassium_current">
    <map_variables variable_1="time" variable_2="time"/>
  </connection>
  <connection component_1="environment" component_2="sodium_background_current">
    <map_variables variable_1="time" variable_2="time"/>
  </connection>
  <connection component_1="environment" component_2="calcium_background_current">
    <map_variables variable_1="time" variable_2="time"/>
  </connection>
  <connection component_1="environment" component_2="transient_outward_current">
    <map_variables variable_1="time" variable_2="time"/>
  </connection>
  <connection component_1="environment" component_2="sodium_potassium_pump_current">
    <map_variables variable_1="time" variable_2="time"/>
  </connection>
  <connection component_1="environment" component_2="sodium_calcium_exchanger_current">
    <map_variables variable_1="time" variable_2="time"/>
  </connection>
  <connection component_1="environment" component_2="calcium_dynamics">
    <map_variables variable_1="time" variable_2="time"/>
  </connection>
  <connection component_1="environment" component_2="sodium_dynamics">
    <map_variables variable_1="time" variable_2="time"/>
  </connection>
  <connection component_1="environment" component_2="potassium_dynamics">
    <map_variables variable_1="time" variable_2="time"/>
  </connection>
  <connection component_1="environment" component_2="reversal_potentials">
    <map_variables variable_1="time" variable_2="time"/>
  </connection>
  <connection component_1="environment" component_2="inward_rectifier_potassium_current">
    <map_variables variable_1="time" variable_2="time"/>
  </connection>
  <connection 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 component_1="membrane" component_2="L_type_Ca_current">
    <map_variables variable_1="V" variable_2="V"/>
    <map_variables variable_1="R" variable_2="R"/>
    <map_variables variable_1="T" variable_2="T"/>
    <map_variables variable_1="F" variable_2="F"/>
    <map_variables variable_1="i_CaL" variable_2="i_CaL"/>
  </connection>
  <connection component_1="membrane" component_2="rapid_time_dependent_potassium_current">
    <map_variables variable_1="V" variable_2="V"/>
    <map_variables variable_1="i_Kr" variable_2="i_Kr"/>
  </connection>
  <connection component_1="membrane" component_2="slow_time_dependent_potassium_current">
    <map_variables variable_1="V" variable_2="V"/>
    <map_variables variable_1="i_Ks" variable_2="i_Ks"/>
  </connection>
  <connection 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 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 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 component_1="membrane" component_2="sodium_potassium_pump_current">
    <map_variables variable_1="V" variable_2="V"/>
    <map_variables variable_1="R" variable_2="R"/>
    <map_variables variable_1="T" variable_2="T"/>
    <map_variables variable_1="F" variable_2="F"/>
    <map_variables variable_1="i_NaK" variable_2="i_NaK"/>
  </connection>
  <connection component_1="membrane" component_2="sodium_calcium_exchanger_current">
    <map_variables variable_1="V" variable_2="V"/>
    <map_variables variable_1="R" variable_2="R"/>
    <map_variables variable_1="T" variable_2="T"/>
    <map_variables variable_1="F" variable_2="F"/>
    <map_variables variable_1="i_NaCa" variable_2="i_NaCa"/>
  </connection>
  <connection component_1="membrane" component_2="potassium_pump_current">
    <map_variables variable_1="V" variable_2="V"/>
    <map_variables variable_1="i_p_K" variable_2="i_p_K"/>
  </connection>
  <connection component_1="membrane" component_2="calcium_dynamics">
    <map_variables variable_1="V" variable_2="V"/>
    <map_variables variable_1="F" variable_2="F"/>
    <map_variables variable_1="Cm" variable_2="Cm"/>
    <map_variables variable_1="V_c" variable_2="V_c"/>
  </connection>
  <connection component_1="membrane" component_2="inward_rectifier_potassium_current">
    <map_variables variable_1="V" variable_2="V"/>
    <map_variables variable_1="i_K1" variable_2="i_K1"/>
  </connection>
  <connection 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 component_1="membrane" component_2="sodium_dynamics">
    <map_variables variable_1="F" variable_2="F"/>
    <map_variables variable_1="Cm" variable_2="Cm"/>
    <map_variables variable_1="V_c" variable_2="V_c"/>
  </connection>
  <connection component_1="membrane" component_2="potassium_dynamics">
    <map_variables variable_1="F" variable_2="F"/>
    <map_variables variable_1="Cm" variable_2="Cm"/>
    <map_variables variable_1="V_c" variable_2="V_c"/>
    <map_variables variable_1="i_Stim" variable_2="i_Stim"/>
  </connection>
  <connection component_1="membrane" component_2="calcium_pump_current">
    <map_variables variable_1="i_p_Ca" variable_2="i_p_Ca"/>
  </connection>
  <connection component_1="fast_sodium_current" component_2="sodium_dynamics">
    <map_variables variable_1="i_Na" variable_2="i_Na"/>
  </connection>
  <connection component_1="fast_sodium_current" component_2="reversal_potentials">
    <map_variables variable_1="E_Na" variable_2="E_Na"/>
  </connection>
  <connection component_1="fast_sodium_current" component_2="fast_sodium_current_m_gate">
    <map_variables variable_1="time" variable_2="time"/>
    <map_variables variable_1="V" variable_2="V"/>
    <map_variables variable_1="m" variable_2="m"/>
  </connection>
  <connection component_1="fast_sodium_current" component_2="fast_sodium_current_h_gate">
    <map_variables variable_1="time" variable_2="time"/>
    <map_variables variable_1="V" variable_2="V"/>
    <map_variables variable_1="h" variable_2="h"/>
  </connection>
  <connection component_1="fast_sodium_current" component_2="fast_sodium_current_j_gate">
    <map_variables variable_1="time" variable_2="time"/>
    <map_variables variable_1="V" variable_2="V"/>
    <map_variables variable_1="j" variable_2="j"/>
  </connection>
  <connection component_1="L_type_Ca_current" component_2="calcium_dynamics">
    <map_variables variable_1="i_CaL" variable_2="i_CaL"/>
    <map_variables variable_1="Ca_o" variable_2="Ca_o"/>
    <map_variables variable_1="Ca_i" variable_2="Ca_i"/>
    <map_variables variable_1="d" variable_2="d"/>
  </connection>
  <connection component_1="L_type_Ca_current" component_2="L_type_Ca_current_fCa_gate">
    <map_variables variable_1="Ca_i" variable_2="Ca_i"/>
    <map_variables variable_1="time" variable_2="time"/>
    <map_variables variable_1="V" variable_2="V"/>
    <map_variables variable_1="fCa" variable_2="fCa"/>
  </connection>
  <connection component_1="L_type_Ca_current" component_2="L_type_Ca_current_d_gate">
    <map_variables variable_1="time" variable_2="time"/>
    <map_variables variable_1="V" variable_2="V"/>
    <map_variables variable_1="d" variable_2="d"/>
  </connection>
  <connection component_1="L_type_Ca_current" component_2="L_type_Ca_current_f_gate">
    <map_variables variable_1="time" variable_2="time"/>
    <map_variables variable_1="V" variable_2="V"/>
    <map_variables variable_1="f" variable_2="f"/>
  </connection>
  <connection component_1="rapid_time_dependent_potassium_current" component_2="potassium_dynamics">
    <map_variables variable_1="i_Kr" variable_2="i_Kr"/>
    <map_variables variable_1="K_o" variable_2="K_o"/>
  </connection>
  <connection component_1="rapid_time_dependent_potassium_current" component_2="reversal_potentials">
    <map_variables variable_1="E_K" variable_2="E_K"/>
  </connection>
  <connection component_1="rapid_time_dependent_potassium_current" component_2="rapid_time_dependent_potassium_current_Xr1_gate">
    <map_variables variable_1="time" variable_2="time"/>
    <map_variables variable_1="V" variable_2="V"/>
    <map_variables variable_1="Xr1" variable_2="Xr1"/>
  </connection>
  <connection component_1="rapid_time_dependent_potassium_current" component_2="rapid_time_dependent_potassium_current_Xr2_gate">
    <map_variables variable_1="time" variable_2="time"/>
    <map_variables variable_1="V" variable_2="V"/>
    <map_variables variable_1="Xr2" variable_2="Xr2"/>
  </connection>
  <connection component_1="slow_time_dependent_potassium_current" component_2="potassium_dynamics">
    <map_variables variable_1="i_Ks" variable_2="i_Ks"/>
  </connection>
  <connection component_1="slow_time_dependent_potassium_current" component_2="reversal_potentials">
    <map_variables variable_1="E_Ks" variable_2="E_Ks"/>
  </connection>
  <connection component_1="slow_time_dependent_potassium_current" component_2="slow_time_dependent_potassium_current_Xs_gate">
    <map_variables variable_1="time" variable_2="time"/>
    <map_variables variable_1="V" variable_2="V"/>
    <map_variables variable_1="Xs" variable_2="Xs"/>
  </connection>
  <connection component_1="sodium_background_current" component_2="sodium_dynamics">
    <map_variables variable_1="i_b_Na" variable_2="i_b_Na"/>
  </connection>
  <connection component_1="sodium_background_current" component_2="reversal_potentials">
    <map_variables variable_1="E_Na" variable_2="E_Na"/>
  </connection>
  <connection component_1="calcium_background_current" component_2="calcium_dynamics">
    <map_variables variable_1="i_b_Ca" variable_2="i_b_Ca"/>
  </connection>
  <connection component_1="calcium_background_current" component_2="reversal_potentials">
    <map_variables variable_1="E_Ca" variable_2="E_Ca"/>
  </connection>
  <connection component_1="transient_outward_current" component_2="potassium_dynamics">
    <map_variables variable_1="i_to" variable_2="i_to"/>
  </connection>
  <connection component_1="transient_outward_current" component_2="reversal_potentials">
    <map_variables variable_1="E_K" variable_2="E_K"/>
  </connection>
  <connection component_1="transient_outward_current" component_2="transient_outward_current_s_gate">
    <map_variables variable_1="V" variable_2="V"/>
    <map_variables variable_1="time" variable_2="time"/>
    <map_variables variable_1="s" variable_2="s"/>
  </connection>
  <connection component_1="transient_outward_current" component_2="transient_outward_current_r_gate">
    <map_variables variable_1="V" variable_2="V"/>
    <map_variables variable_1="time" variable_2="time"/>
    <map_variables variable_1="r" variable_2="r"/>
  </connection>
  <connection component_1="sodium_potassium_pump_current" component_2="sodium_dynamics">
    <map_variables variable_1="i_NaK" variable_2="i_NaK"/>
    <map_variables variable_1="Na_i" variable_2="Na_i"/>
  </connection>
  <connection component_1="sodium_potassium_pump_current" component_2="potassium_dynamics">
    <map_variables variable_1="i_NaK" variable_2="i_NaK"/>
    <map_variables variable_1="K_o" variable_2="K_o"/>
  </connection>
  <connection component_1="sodium_calcium_exchanger_current" component_2="calcium_dynamics">
    <map_variables variable_1="i_NaCa" variable_2="i_NaCa"/>
    <map_variables variable_1="Ca_i" variable_2="Ca_i"/>
    <map_variables variable_1="Ca_o" variable_2="Ca_o"/>
  </connection>
  <connection component_1="sodium_calcium_exchanger_current" component_2="sodium_dynamics">
    <map_variables variable_1="i_NaCa" variable_2="i_NaCa"/>
    <map_variables variable_1="Na_i" variable_2="Na_i"/>
    <map_variables variable_1="Na_o" variable_2="Na_o"/>
  </connection>
  <connection component_1="calcium_pump_current" component_2="calcium_dynamics">
    <map_variables variable_1="i_p_Ca" variable_2="i_p_Ca"/>
    <map_variables variable_1="Ca_i" variable_2="Ca_i"/>
  </connection>
  <connection component_1="potassium_pump_current" component_2="potassium_dynamics">
    <map_variables variable_1="i_p_K" variable_2="i_p_K"/>
  </connection>
  <connection component_1="potassium_pump_current" component_2="reversal_potentials">
    <map_variables variable_1="E_K" variable_2="E_K"/>
  </connection>
  <connection component_1="calcium_dynamics" component_2="TT04">
    <map_variables variable_1="Ca_i" variable_2="Ca_i"/>
  </connection>
  <connection component_1="calcium_dynamics" component_2="reversal_potentials">
    <map_variables variable_1="Ca_i" variable_2="Ca_i"/>
    <map_variables variable_1="Ca_o" variable_2="Ca_o"/>
  </connection>
  <connection component_1="sodium_dynamics" component_2="reversal_potentials">
    <map_variables variable_1="Na_i" variable_2="Na_i"/>
    <map_variables variable_1="Na_o" variable_2="Na_o"/>
  </connection>
  <connection component_1="potassium_dynamics" component_2="reversal_potentials">
    <map_variables variable_1="K_i" variable_2="K_i"/>
    <map_variables variable_1="K_o" variable_2="K_o"/>
  </connection>
  <connection component_1="potassium_dynamics" component_2="inward_rectifier_potassium_current">
    <map_variables variable_1="K_o" variable_2="K_o"/>
    <map_variables variable_1="i_K1" variable_2="i_K1"/>
  </connection>
  <connection component_1="reversal_potentials" component_2="inward_rectifier_potassium_current">
    <map_variables variable_1="E_K" variable_2="E_K"/>
  </connection>
  <connection component_1="main" component_2="TT04">
    <map_variables variable_1="phi_Cai" variable_2="phi_Cai"/>
    <map_variables variable_1="time" variable_2="time"/>
  </connection>
  <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="fast_sodium_current_j_gate"/>
    </component_ref>
    <component_ref component="L_type_Ca_current">
      <component_ref component="L_type_Ca_current_d_gate"/>
      <component_ref component="L_type_Ca_current_f_gate"/>
      <component_ref component="L_type_Ca_current_fCa_gate"/>
    </component_ref>
    <component_ref component="rapid_time_dependent_potassium_current">
      <component_ref component="rapid_time_dependent_potassium_current_Xr1_gate"/>
      <component_ref component="rapid_time_dependent_potassium_current_Xr2_gate"/>
    </component_ref>
    <component_ref component="slow_time_dependent_potassium_current">
      <component_ref component="slow_time_dependent_potassium_current_Xs_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>
  </encapsulation>
</model>