Rendering of the source text

<?xml version="1.0" encoding="utf-8"?>
<!--
This CellML file was generated on 10.02.2012 at 13:45:07 using:

COR (0.9.31.1409)
Copyright 2002-2012 Dr Alan Garny
http://cor.physiol.ox.ac.uk/ - cor@physiol.ox.ac.uk

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

      <dcterms:created rdf:parseType="Resource">
        <dcterms:W3CDTF>2009-09-28</dcterms:W3CDTF>
      </dcterms:created>
      -->
      <!--
        The Modification History metadata. This lists the changes that have been
        made to the document, who made the changes, and when they were made.
      -->
     
    </rdf:Description>
    
    <!--
      The following metadata refers to the model itself, as indicated by the
      reference to the ID "ekaterinburg_rat", which is    
      declared on the <model> element.
    -->
    <rdf:Description rdf:about="#ekaterinburg_rat">
      <!-- A human readable name for the model. -->
      <dc:title>
        Mathematical modeling of mechanically modulated rhythm disturbances in homogeneous and heterogeneous myocardium with attenuated activity of na+-k+ pump
      </dc:title>
      
      <!--  Keyword(s) -->
      <bqs:reference rdf:parseType="Resource">
        <dc:subject rdf:parseType="Resource">
          <bqs:subject_type>keyword</bqs:subject_type>
          <rdf:value>
            <rdf:Bag>
              <rdf:li>cardiac</rdf:li>
              <rdf:li>electrophysiology</rdf:li>
            </rdf:Bag>
          </rdf:value>
        </dc:subject>
      </bqs:reference>

      <bqs:reference rdf:parseType="Resource">
        <bqs:Pubmed_id>18259823</bqs:Pubmed_id>
        <bqs:JournalArticle rdf:parseType="Resource">
          <dc:creator>
            <rdf:Seq>
              <rdf:li rdf:parseType="Resource">
                <bqs:Person rdf:parseType="Resource">
                  <vCard:N rdf:parseType="Resource">
                    <vCard:Family>Sulman</vCard:Family>
                    <vCard:Given>Tatiana</vCard:Given>
                  </vCard:N>
                </bqs:Person>
              </rdf:li>
              <rdf:li rdf:parseType="Resource">
                <bqs:Person rdf:parseType="Resource">
                  <vCard:N rdf:parseType="Resource">
                    <vCard:Family>Katnelson</vCard:Family>
                    <vCard:Given>Leonid</vCard:Given>
                    <vCard:Other>B</vCard:Other>
                  </vCard:N>
                </bqs:Person>
              </rdf:li>
			  <rdf:li rdf:parseType="Resource">
                <bqs:Person rdf:parseType="Resource">
                  <vCard:N rdf:parseType="Resource">
                    <vCard:Family>Solovyova</vCard:Family>
                    <vCard:Given>Olga</vCard:Given>
                  </vCard:N>
                </bqs:Person>
              </rdf:li>
			  <rdf:li rdf:parseType="Resource">
                <bqs:Person rdf:parseType="Resource">
                  <vCard:N rdf:parseType="Resource">
                    <vCard:Family>Markhasin</vCard:Family>
                    <vCard:Given>Vladimir</vCard:Given>
                    <vCard:Other>S</vCard:Other>
				  </vCard:N>
                </bqs:Person>
              </rdf:li>
            </rdf:Seq>
          </dc:creator>
          <dc:title>
            Mathematical modeling of mechanically modulated rhythm disturbances in homogeneous and heterogeneous myocardium with attenuated activity of na+-k+ pump
          </dc:title>
          <bqs:Journal rdf:parseType="Resource">
            <dc:title>Bulletin of Mathematical Biology</dc:title>
          </bqs:Journal>
          <bqs:volume>70</bqs:volume>
          <bqs:first_page>910</bqs:first_page>
          <bqs:last_page>949</bqs:last_page>
        </bqs:JournalArticle>
      </bqs:reference>
    </rdf:Description>
  </rdf:RDF> 
</model>