Location: Pasek, Simurda, Christe, 2006 @ 49d6ba3b0e16 / pasek_model_2006_surface_membrane.cellml

Author:
Catherine Lloyd <c.lloyd@auckland.ac.nz>
Date:
2009-11-10 13:47:18+13:00
Desc:
Deleted variant 01 as it is a duplicate of b). Also corrected the link to the diagram in the main pasek 1.1 file. Note that this folder contains 3 1.1 models and 2 variants of the model in 1.0.
Permanent Source URI:
https://models.cellml.org/workspace/pasek_simurda_christe_2006/rawfile/49d6ba3b0e16b16b8a9f890041a007b1d540b05f/pasek_model_2006_surface_membrane.cellml

<?xml version="1.0" encoding="utf-8"?>
<model name="pasek_model_2006_surface_membrane" cmeta:id="pasek_model_2006_surface_membrane" xmlns="http://www.cellml.org/cellml/1.1#" xmlns:cellml="http://www.cellml.org/cellml/1.1#" xmlns:cmeta="http://www.cellml.org/metadata/1.0#" xmlns:xlink="http://www.w3.org/1999/xlink" 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#">
	
   <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="per_millivolt2">
      <unit units="volt" prefix="milli" exponent="-2"/>
   </units>
   <units name="milliS">
      <unit units="siemens" prefix="milli"/>
   </units>
   <units name="milliS_per_cm2">
      <unit units="siemens" prefix="milli"/>
	  <unit units="metre" prefix="centi" exponent="-2"/>
   </units>
   <units name="microA_per_cm2">
      <unit units="ampere" prefix="micro"/>
	  <unit units="metre" prefix="centi" exponent="-2"/>
   </units>
   <units name="cm2">
      <unit units="metre" prefix="centi" exponent="2"/>
   </units>
   <units name="cm_per_second">
      <unit units="metre" prefix="centi"/>
      <unit units="second" exponent="-1"/>
   </units>
   <units name="cm3_per_second">
      <unit units="metre" prefix="centi" exponent="3"/>
      <unit units="second" exponent="-1"/>
   </units>
   <units name="microF">
      <unit units="farad" prefix="micro"/>
   </units>
   <units name="microA">
      <unit units="ampere" prefix="micro"/>
   </units>
   <units name="millimolar">
      <unit units="mole" prefix="milli"/>
      <unit units="litre" exponent="-1"/>
   </units>
   <units name="millimolar_per_second">
      <unit units="millimolar"/>
      <unit units="second" exponent="-1"/>
   </units>
   <units name="per_millimolar_per_second">
      <unit units="millimolar" exponent="-1"/>
      <unit units="second" exponent="-1"/>
   </units>
   <units name="per_millimolar3_per_second">
      <unit units="millimolar" exponent="-3"/>
      <unit units="second" exponent="-1"/>
   </units>
   <units name="millimolar4">
      <unit units="millimolar" exponent="4"/>
   </units>
   <units name="microA_per_millimolar">
      <unit units="microA"/>
      <unit units="millimolar" exponent="-1"/>
   </units>
   <units name="joule_per_kilomole_kelvin">
      <unit units="joule"/>
      <unit units="kelvin" exponent="-1"/>
      <unit units="mole" prefix="kilo" exponent="-1"/>
   </units>
   <units name="coulomb_per_mole">
      <unit units="coulomb"/>
      <unit units="mole" exponent="-1"/>
   </units>
   <units name="mole_per_coulomb">
      <unit units="coulomb_per_mole" exponent="-1"/>
   </units>
   <units name="micro_litre">
      <unit units="litre" prefix="micro"/>
   </units>
   
   <cmeta:comment rdf:parseType="Resource">
        <rdf:value>
          In order to facilitate the import of this surface membrane model by the main Pasek 2008 model I will encapsulate the current model in a top component called &quot;surface_membrane&quot;.
        </rdf:value>
        <dc:creator rdf:parseType="Resource">
          <vCard:FN>Catherine Lloyd</vCard:FN>
        </dc:creator>
      </cmeta:comment>

  <component name="surface_membrane" cmeta:id="surface_membrane">
    <rdf:RDF xmlns:rdf="http://www.w3.org/1999/0P/PP-rdf-syntax-ns#" xmlns:dc="http://purl.org/dc/elements/1.1/" xmlns:dcterms="http://purl.org/dc/terms/">
      <rdf:Description rdf:about="surface_membrane">
        <dc:title>surface_membrane</dc:title>
        <dcterms:alternative>top component which encapsulates the entire surface membrane model</dcterms:alternative>  
      </rdf:Description>
    </rdf:RDF>
    <variable name="Vm_s" units="millivolt" public_interface="out" private_interface="in"/>
	<variable name="i_Na_s" units="microA" public_interface="out" private_interface="in"/>
    <variable name="i_Ca_s" units="microA" public_interface="out" private_interface="in"/>
    <variable name="i_Kto_s" units="microA" public_interface="out" private_interface="in"/>
    <variable name="i_Kss_s" units="microA" public_interface="out" private_interface="in"/>
    <variable name="i_fNa_s" units="microA" public_interface="out" private_interface="in"/>
	<variable name="i_fK_s" units="microA" public_interface="out" private_interface="in"/>
    <variable name="i_K1_s" units="microA" public_interface="out" private_interface="in"/>
    <variable name="i_BNa_s" units="microA" public_interface="out" private_interface="in"/>
	<variable name="i_BK_s" units="microA" public_interface="out" private_interface="in"/>
	<variable name="i_BCa_s" units="microA" public_interface="out" private_interface="in"/>
    <variable name="i_NaK_s" units="microA" public_interface="out" private_interface="in"/>
    <variable name="i_NaCa_s" units="microA" public_interface="out" private_interface="in"/>
    <variable name="i_CaP_s" units="microA" public_interface="out" private_interface="in"/>
	
    <variable name="time" public_interface="in" private_interface="out" units="second"/>
    <variable name="R" units="joule_per_kilomole_kelvin" public_interface="in" private_interface="out"/>
    <variable name="T" units="kelvin" public_interface="in" private_interface="out"/>
    <variable name="F" units="coulomb_per_mole" public_interface="in" private_interface="out"/>
	<variable name="Na_e" units="millimolar" public_interface="in" private_interface="out"/>
    <variable name="Ca_e" units="millimolar" public_interface="in" private_interface="out"/>
    <variable name="K_e" units="millimolar" public_interface="in" private_interface="out"/>
	<variable name="Na_i" units="millimolar" public_interface="in" private_interface="out"/>
    <variable name="Ca_i" units="millimolar" public_interface="in" private_interface="out"/>
	<variable name="Ca_ss" units="millimolar" public_interface="in" private_interface="out"/>
    <variable name="K_i" units="millimolar" public_interface="in" private_interface="out"/>
    <variable name="i_circ" units="microA" public_interface="in" private_interface="out"/>	
	<variable name="i_Kext" units="microA" public_interface="in" private_interface="out"/>
	<variable name="Cms" units="microF" public_interface="in" private_interface="out"/>
	<variable name="Sms" units="cm2" public_interface="in" private_interface="out"/>
	<variable name="Smt" units="cm2" public_interface="in" private_interface="out"/>
  </component>
  
  <!-- 1.1 connections -->
  
  <connection>
      <map_components component_1="membrane" component_2="surface_membrane"/>
      <map_variables variable_1="time" variable_2="time"/>
	  <map_variables variable_2="Vm_s" variable_1="Vm_s"/>
	  <map_variables variable_1="i_circ" variable_2="i_circ"/>
	  <map_variables variable_1="i_Kext" variable_2="i_Kext"/>
	  <map_variables variable_2="Cms" variable_1="Cms"/>
   </connection>
   <connection>
      <map_components component_1="sodium_current" component_2="surface_membrane"/>
      <map_variables variable_1="time" variable_2="time"/>
	  <map_variables variable_2="R" variable_1="R"/>
      <map_variables variable_2="T" variable_1="T"/>
      <map_variables variable_2="F" variable_1="F"/>
	  <map_variables variable_2="Na_e" variable_1="Na_e"/>
	  <map_variables variable_1="Na_i" variable_2="Na_i"/>
	  <map_variables variable_2="K_e" variable_1="K_e"/>
	  <map_variables variable_1="K_i" variable_2="K_i"/>
	  <map_variables variable_1="i_Na_s" variable_2="i_Na_s"/>
	  <map_variables variable_2="Sms" variable_1="Sms"/>
      <map_variables variable_2="Smt" variable_1="Smt"/>
   </connection>
   <connection>
      <map_components component_1="calcium_current" component_2="surface_membrane"/>
      <map_variables variable_1="time" variable_2="time"/>
	  <map_variables variable_2="R" variable_1="R"/>
      <map_variables variable_2="T" variable_1="T"/>
      <map_variables variable_2="F" variable_1="F"/>
	  <map_variables variable_2="Ca_e" variable_1="Ca_e"/>
	  <map_variables variable_2="Ca_ss" variable_1="Ca_ss"/>
	  <map_variables variable_1="i_Ca_s" variable_2="i_Ca_s"/>
	  <map_variables variable_2="Sms" variable_1="Sms"/>
      <map_variables variable_2="Smt" variable_1="Smt"/>
   </connection>
   <connection>
      <map_components component_1="Ca_independent_transient_outward_K_current" component_2="surface_membrane"/>
      <map_variables variable_1="time" variable_2="time"/>
	  <map_variables variable_2="R" variable_1="R"/>
      <map_variables variable_2="T" variable_1="T"/>
      <map_variables variable_2="F" variable_1="F"/>
	  <map_variables variable_2="K_e" variable_1="K_e"/>
	  <map_variables variable_2="K_i" variable_1="K_i"/>
	  <map_variables variable_1="i_Kto_s" variable_2="i_Kto_s"/>
	  <map_variables variable_2="Sms" variable_1="Sms"/>
      <map_variables variable_2="Smt" variable_1="Smt"/>
   </connection>
   <connection>
      <map_components component_1="steady_state_outward_K_current" component_2="surface_membrane"/>
      <map_variables variable_1="time" variable_2="time"/>
	  <map_variables variable_1="i_Kss_s" variable_2="i_Kss_s"/>
	  <map_variables variable_2="Sms" variable_1="Sms"/>
      <map_variables variable_2="Smt" variable_1="Smt"/>
   </connection>
   <connection>
      <map_components component_1="inward_rectifier" component_2="surface_membrane"/>
      <map_variables variable_1="time" variable_2="time"/>
	  <map_variables variable_2="R" variable_1="R"/>
      <map_variables variable_2="T" variable_1="T"/>
      <map_variables variable_2="F" variable_1="F"/>
	  <map_variables variable_2="K_e" variable_1="K_e"/>
	  <map_variables variable_1="i_K1_s" variable_2="i_K1_s"/>
	  <map_variables variable_2="Sms" variable_1="Sms"/>
      <map_variables variable_2="Smt" variable_1="Smt"/>
   </connection>
   <connection>
      <map_components component_1="hyperpolarisation_activated_current" component_2="surface_membrane"/>
      <map_variables variable_1="time" variable_2="time"/>
	  <map_variables variable_1="i_fNa_s" variable_2="i_fNa_s"/>
	  <map_variables variable_1="i_fK_s" variable_2="i_fK_s"/>
	  <map_variables variable_2="Sms" variable_1="Sms"/>
      <map_variables variable_2="Smt" variable_1="Smt"/>
   </connection>
   <connection>
      <map_components component_1="sarcolemmal_calcium_pump_current" component_2="surface_membrane"/>
      <map_variables variable_1="time" variable_2="time"/>
	  <map_variables variable_1="i_CaP_s" variable_2="i_CaP_s"/>
	  <map_variables variable_2="Ca_i" variable_1="Ca_i"/>
	  <map_variables variable_2="Sms" variable_1="Sms"/>
      <map_variables variable_2="Smt" variable_1="Smt"/>
   </connection>
   <connection>
      <map_components component_1="Na_Ca_ion_exchanger_current" component_2="surface_membrane"/>
      <map_variables variable_1="time" variable_2="time"/>
	  <map_variables variable_2="Na_e" variable_1="Na_e"/>
	  <map_variables variable_2="Ca_e" variable_1="Ca_e"/>
	  <map_variables variable_2="Na_i" variable_1="Na_i"/>
	  <map_variables variable_2="Ca_i" variable_1="Ca_i"/>
	  <map_variables variable_1="i_NaCa_s" variable_2="i_NaCa_s"/>
	  <map_variables variable_2="Sms" variable_1="Sms"/>
      <map_variables variable_2="Smt" variable_1="Smt"/>
   </connection>
   <connection>
      <map_components component_1="sodium_potassium_pump" component_2="surface_membrane"/>
      <map_variables variable_2="R" variable_1="R"/>
      <map_variables variable_2="T" variable_1="T"/>
      <map_variables variable_2="F" variable_1="F"/>
	  <map_variables variable_2="Na_e" variable_1="Na_e"/>
	  <map_variables variable_2="Na_i" variable_1="Na_i"/>
	  <map_variables variable_2="K_e" variable_1="K_e"/>
	  <map_variables variable_1="i_NaK_s" variable_2="i_NaK_s"/>
	  <map_variables variable_2="Sms" variable_1="Sms"/>
      <map_variables variable_2="Smt" variable_1="Smt"/>
   </connection>
   <connection>
      <map_components component_1="background_currents" component_2="surface_membrane"/>
      <map_variables variable_2="R" variable_1="R"/>
      <map_variables variable_2="T" variable_1="T"/>
      <map_variables variable_2="F" variable_1="F"/>
	  <map_variables variable_2="Ca_e" variable_1="Ca_e"/>
	  <map_variables variable_1="i_BNa_s" variable_2="i_BNa_s"/>
      <map_variables variable_1="i_BK_s" variable_2="i_BK_s"/>
      <map_variables variable_1="i_BCa_s" variable_2="i_BCa_s"/>
      <map_variables variable_1="Ca_i" variable_2="Ca_i"/>
	  <map_variables variable_2="Sms" variable_1="Sms"/>
      <map_variables variable_2="Smt" variable_1="Smt"/>
   </connection>

   <!-- surface membrane model -->   
   
   <component name="membrane">
      <variable name="Vm_s" units="millivolt" initial_value="-7.9709547e1" public_interface="out"/>
      
	  <variable name="Cms" units="microF" public_interface="in"/>
	  <variable name="i_circ" units="microA" public_interface="in"/>
      <variable name="i_Kext" units="microA" public_interface="in"/>
	  <variable name="time" units="second" public_interface="in"/>
      <variable name="i_Na_s" units="microA" public_interface="in"/>
      <variable name="i_Ca_s" units="microA" public_interface="in"/>
      <variable name="i_Kto_s" units="microA" public_interface="in"/>
      <variable name="i_Kss_s" units="microA" public_interface="in"/>
      <variable name="i_f_s" units="microA" public_interface="in"/>
      <variable name="i_K1_s" units="microA" public_interface="in"/>
      <variable name="i_B_s" units="microA" public_interface="in"/>
      <variable name="i_NaK_s" units="microA" public_interface="in"/>
      <variable name="i_NaCa_s" units="microA" public_interface="in"/>
      <variable name="i_CaP_s" units="microA" public_interface="in"/>
      
	  <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>Vm_s</ci>
            </apply>
            <apply><times/>
			   <cn cellml:units="dimensionless">1000</cn>
			   <apply><divide/>
                  <apply><minus/>
                     <apply><plus/>
					    <ci>i_Kext</ci>
                        <ci>i_circ</ci>
					 </apply>
					 <apply><plus/>
                        <ci>i_Na_s</ci>
                        <ci>i_Ca_s</ci>
                        <ci>i_Kto_s</ci>
                        <ci>i_Kss_s</ci>
                        <ci>i_f_s</ci>
                        <ci>i_K1_s</ci>
                        <ci>i_B_s</ci>
                        <ci>i_NaK_s</ci>
                        <ci>i_NaCa_s</ci>
                        <ci>i_CaP_s</ci>
				     </apply>
			      </apply>
				  <ci>Cms</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   
   <component name="sodium_current">
      <variable name="i_Na_s" units="microA" public_interface="out"/>
      <variable name="E_Na" units="millivolt" public_interface="out"/>
      
	  <variable name="g_Na" units="milliS_per_cm2" initial_value="10.0"/>
	  <variable name="g_Nas" units="milliS"/>
	  <variable name="fNat" units="dimensionless" initial_value="0.56"/>
	  
	  <variable name="Sms" units="cm2" public_interface="in"/>
	  <variable name="Smt" units="cm2" public_interface="in"/>
      <variable name="Na_i" units="millimolar" public_interface="in"/>
      <variable name="Na_e" units="millimolar" public_interface="in"/>
	  <variable name="K_i" units="millimolar" public_interface="in"/>
      <variable name="K_e" 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="Vm_s" 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>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_e</ci>
                     <ci>Na_i</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         
		 <apply>
            <eq/>
            <ci>i_Na_s</ci>
            <apply>
               <times/>
               <ci>g_Nas</ci>
               <apply>
                  <power/>
                  <ci>m</ci>
                  <cn cellml:units="dimensionless">3</cn>
               </apply>
               <ci>h</ci>
               <apply>
                  <minus/>
                  <ci>Vm_s</ci>
				  <apply>
                     <times/>
					 <apply>
                        <log/>
						<apply>
                           <divide/>
                           <apply>
                              <plus/>
                              <ci>Na_e</ci>
							  <apply>
                                 <times/>
								 <cn cellml:units="dimensionless">0.12</cn>
								 <ci>K_e</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>
                        <divide/>
                        <apply>
                           <times/>
                           <ci>R</ci>
                           <ci>T</ci>
                        </apply>
                        <ci>F</ci>
				     </apply>
				  </apply>
               </apply>
            </apply>
         </apply>
		 
		 <apply><eq/>
            <ci>g_Nas</ci>
            <apply><times/>
               <ci>g_Na</ci>
               <apply><plus/>
                  <ci>Sms</ci>
				  <ci>Smt</ci>
               </apply>
               <apply><minus/>
                  <cn cellml:units="dimensionless">1</cn>
				  <ci>fNat</ci>
			   </apply>
			</apply>
         </apply>
      </math>
   </component>
   
   <component name="sodium_current_m_gate">
      <variable name="m" units="dimensionless" initial_value="2.1832558e-6" public_interface="out"/>
      <variable name="alpha_m" units="per_second"/>
      <variable name="beta_m" units="per_second"/>
      <variable name="Vm_s" 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_m</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <cn cellml:units="per_millivolt_second">117.26</cn>
                  <apply>
                     <plus/>
                     <ci>Vm_s</ci>
                     <cn cellml:units="millivolt">59.3</cn>
                  </apply>
               </apply>
               <apply>
                  <minus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <times/>
                        <apply>
                           <minus/>
                           <cn cellml:units="per_millivolt">0.55</cn>
                        </apply>
                        <apply>
                           <plus/>
                           <ci>Vm_s</ci>
                           <cn cellml:units="millivolt">59.3</cn>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>beta_m</ci>
            <apply>
               <times/>
               <cn cellml:units="per_second">3800</cn>
               <apply>
                  <exp/>
                  <apply>
                     <times/>
                     <apply>
                        <minus/>
                        <cn cellml:units="per_millivolt">0.072</cn>
                     </apply>
                     <apply>
                        <plus/>
                        <ci>Vm_s</ci>
                        <cn cellml:units="millivolt">61</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="sodium_current_h_gate">
      <variable name="h" units="dimensionless" initial_value="5.8232781e-1" public_interface="out"/>
      <variable name="alpha_h" units="per_second"/>
      <variable name="beta_h" units="per_second"/>
      <variable name="Vm_s" 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>
               <divide/>
               <cn cellml:units="per_second">284.4</cn>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <times/>
                        <cn cellml:units="per_millivolt">0.0812</cn>
                        <apply>
                           <plus/>
                           <ci>Vm_s</ci>
                           <cn cellml:units="millivolt">115.9</cn>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>beta_h</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <cn cellml:units="per_millivolt_second">18.77</cn>
                  <apply>
                     <plus/>
                     <ci>Vm_s</ci>
                     <cn cellml:units="millivolt">64.4</cn>
                  </apply>
               </apply>
               <apply>
                  <minus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <times/>
                        <apply>
                           <minus/>
                           <cn cellml:units="per_millivolt">0.22</cn>
                        </apply>
                        <apply>
                           <plus/>
                           <ci>Vm_s</ci>
                           <cn cellml:units="millivolt">64.4</cn>
                        </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="calcium_current">
      <variable name="i_Ca_s" units="microA" public_interface="out"/>
      
	  <variable name="Co" units="dimensionless" initial_value="2.5436465e-6"/>
	  <variable name="C_cast" units="dimensionless" initial_value="9.9570063e-1"/>
	  <variable name="C1" units="dimensionless" initial_value="9.9994546e-1"/>
	  <variable name="C2" units="dimensionless" initial_value="9.9994546e-1"/>
	  <variable name="tauCo" units="second"/>
	  <variable name="tauC_cast" units="second" initial_value="0.009"/>
	  <variable name="tauC1" units="second"/>
	  <variable name="tauC2" units="second"/>
	  <variable name="Co_inf" units="dimensionless"/>
	  <variable name="C_cast_inf" units="dimensionless"/>
	  <variable name="C1_inf" units="dimensionless"/>
	  <variable name="C2_inf" units="dimensionless"/>
      <variable name="g_Ca" units="cm_per_second" initial_value="0.00021"/>
	  <variable name="g_Cas" units="cm3_per_second"/>
	  <variable name="fCat" units="dimensionless" initial_value="0.87"/>
	  
	  <variable name="Sms" units="cm2" public_interface="in"/>
	  <variable name="Smt" units="cm2" public_interface="in"/>
	  <variable name="time" units="second" public_interface="in" private_interface="out"/>
      <variable name="R" units="joule_per_kilomole_kelvin" public_interface="in"/>
      <variable name="T" units="kelvin" public_interface="in"/>
      <variable name="F" units="coulomb_per_mole" public_interface="in"/>
      <variable name="Vm_s" units="millivolt" public_interface="in"/>
      <variable name="Ca_ss" units="millimolar" public_interface="in"/>
      <variable name="Ca_e" units="millimolar" public_interface="in"/>
      
	  <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>Co_inf</ci>
            <apply>
               <divide/>
               <cn cellml:units="dimensionless">1</cn>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <plus/>
                           <ci>Vm_s</ci>
                           <cn cellml:units="millivolt">15.3</cn>
                        </apply>
                        <apply>
                           <minus/>
                           <cn cellml:units="millivolt">5</cn>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
		 
		 <apply>
            <eq/>
            <ci>tauCo</ci>
            <apply>
               <plus/>
               <apply>
                  <times/>
                  <cn cellml:units="second">0.00305</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <times/>
                        <apply>
                           <minus/>
                           <cn cellml:units="per_millivolt2">0.0045</cn>
                        </apply>
                        <apply>
                           <power/>
                           <apply>
                              <plus/>
                              <ci>Vm_s</ci>
                              <cn cellml:units="millivolt">7</cn>
                           </apply>
                           <cn cellml:units="dimensionless">2</cn>
                        </apply>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <cn cellml:units="second">0.00105</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <times/>
                        <apply>
                           <minus/>
                           <cn cellml:units="per_millivolt2">0.002</cn>
                        </apply>
                        <apply>
                           <power/>
                           <apply>
                              <minus/>
                              <ci>Vm_s</ci>
                              <cn cellml:units="millivolt">18</cn>
                           </apply>
                           <cn cellml:units="dimensionless">2</cn>
                        </apply>
                     </apply>
                  </apply>
               </apply>
               <cn cellml:units="second">0.00025</cn>
            </apply>
         </apply>
		 
		 <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>Co</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>Co_inf</ci>
                  <ci>Co</ci>
               </apply>
               <ci>tauCo</ci>
            </apply>
         </apply>
		 
		 <apply>
            <eq/>
            <ci>C1_inf</ci>
            <apply>
               <divide/>
               <cn cellml:units="dimensionless">1</cn>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <plus/>
                           <ci>Vm_s</ci>
                           <cn cellml:units="millivolt">26.7</cn>
                        </apply>
                        <cn cellml:units="millivolt">5.4</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         
		 <apply>
            <eq/>
            <ci>tauC1</ci>
            <apply>
               <plus/>
               <apply>
                  <times/>
                  <cn cellml:units="second">0.105</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <minus/>
                        <apply>
                           <power/>
                           <apply>
                              <divide/>
                              <apply>
                                 <plus/>
                                 <ci>Vm_s</ci>
                                 <cn cellml:units="millivolt">45</cn>
                              </apply>
                              <cn cellml:units="millivolt">12</cn>
                           </apply>
                           <cn cellml:units="dimensionless">2</cn>
                        </apply>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <divide/>
                  <cn cellml:units="second">0.04</cn>
                  <apply>
                     <plus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <plus/>
                              <apply>
                                 <minus/>
                                 <ci>Vm_s</ci>
                              </apply>
                              <cn cellml:units="millivolt">25</cn>
                           </apply>
                           <cn cellml:units="millivolt">25</cn>
                        </apply>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <divide/>
                  <cn cellml:units="second">0.015</cn>
                  <apply>
                     <plus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <plus/>
                              <ci>Vm_s</ci>
                              <cn cellml:units="millivolt">75</cn>
                           </apply>
                           <cn cellml:units="millivolt">25</cn>
                        </apply>
                     </apply>
                  </apply>
               </apply>
               <cn cellml:units="second">0.0017</cn>
            </apply>
         </apply>
         
		 <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>C1</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>C1_inf</ci>
                  <ci>C1</ci>
               </apply>
               <ci>tauC1</ci>
            </apply>
         </apply>
		 
		 <apply>
            <eq/>
            <ci>C2_inf</ci>
            <apply>
               <divide/>
               <cn cellml:units="dimensionless">1</cn>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <plus/>
                           <ci>Vm_s</ci>
                           <cn cellml:units="millivolt">26.7</cn>
                        </apply>
                        <cn cellml:units="millivolt">5.4</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>tauC2</ci>
            <apply>
               <plus/>
               <apply>
                  <times/>
                  <cn cellml:units="second">0.041</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <minus/>
                        <apply>
                           <power/>
                           <apply>
                              <divide/>
                              <apply>
                                 <plus/>
                                 <ci>Vm_s</ci>
                                 <cn cellml:units="millivolt">47</cn>
                              </apply>
                              <cn cellml:units="millivolt">12</cn>
                           </apply>
                           <cn cellml:units="dimensionless">2</cn>
                        </apply>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <divide/>
                  <cn cellml:units="second">0.08</cn>
                  <apply>
                     <plus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <plus/>
                              <ci>Vm_s</ci>
                              <cn cellml:units="millivolt">55</cn>
                           </apply>
                           <apply>
                              <minus/>
                              <cn cellml:units="millivolt">5</cn>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <divide/>
                  <cn cellml:units="second">0.015</cn>
                  <apply>
                     <plus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <plus/>
                              <ci>Vm_s</ci>
                              <cn cellml:units="millivolt">75</cn>
                           </apply>
                           <cn cellml:units="millivolt">25</cn>
                        </apply>
                     </apply>
                  </apply>
               </apply>
               <cn cellml:units="second">0.0017</cn>
            </apply>
         </apply>
         
		 <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>C2</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>C2_inf</ci>
                  <ci>C2</ci>
               </apply>
               <ci>tauC2</ci>
            </apply>
         </apply>
		 
		 <apply>
            <eq/>
            <ci>C_cast_inf</ci>
            <apply>
               <divide/>
               <cn cellml:units="dimensionless">1</cn>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <divide/>
                     <ci>Ca_ss</ci>
                     <cn cellml:units="millimolar">0.01</cn>
                  </apply>
               </apply>
            </apply>
         </apply>
		 
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>C_cast</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>C_cast_inf</ci>
                  <ci>C_cast</ci>
               </apply>
               <ci>tauC_cast</ci>
            </apply>
         </apply>
		 
		 <apply><eq/>
            <ci>g_Cas</ci>
            <apply><times/>
               <ci>g_Ca</ci>
               <apply><plus/>
                  <ci>Sms</ci>
				  <ci>Smt</ci>
               </apply>
               <apply><minus/>
                  <cn cellml:units="dimensionless">1</cn>
				  <ci>fCat</ci>
			   </apply>
			</apply>
         </apply>
		 
		 <apply><eq/>
            <ci>i_Ca_s</ci>
            <apply><times/>
			   <ci>g_Cas</ci>
               <cn cellml:units="dimensionless">4</cn>
			   <apply><divide/>
				  <apply><times/>
				     <ci>Vm_s</ci>
                     <apply><power/>
                        <ci>F</ci>
                        <cn cellml:units="dimensionless">2</cn>
                     </apply>
                  </apply>
				  <apply><times/>
                     <ci>R</ci>
                     <ci>T</ci>
				  </apply>
			   </apply>
               <apply><divide/>
			      <apply><minus/>
                     <apply><times/>
                        <ci>Ca_ss</ci>
                        <apply><exp/>
                           <apply><divide/>
                              <apply><times/>
                                 <cn cellml:units="mole_per_coulomb">2</cn>
                                 <ci>Vm_s</ci>
                                 <apply><power/>
                                    <ci>F</ci>
                                    <cn cellml:units="dimensionless">2</cn>
                                 </apply>
                              </apply>
							  <apply><times/>
                                 <ci>R</ci>
                                 <ci>T</ci>
						      </apply>
                           </apply>
                        </apply>
                     </apply>
                     <apply><times/>
                        <cn cellml:units="dimensionless">0.341</cn>
                        <ci>Ca_e</ci>
                     </apply>
                  </apply>
				  <apply><minus/>
                     <apply><exp/>
                        <apply><divide/>
                           <apply><times/>
                              <cn cellml:units="mole_per_coulomb">2</cn>
                              <ci>Vm_s</ci>
                              <apply><power/>
                                 <ci>F</ci>
                                 <cn cellml:units="dimensionless">2</cn>
                              </apply>
                           </apply>
                           <apply><times/>
                              <ci>R</ci>
                              <ci>T</ci>
						   </apply>
                        </apply>
                     </apply>
					 <cn cellml:units="dimensionless">1</cn>
                  </apply>
               </apply>
               <ci>Co</ci>
               <apply><plus/>
                  <apply><times/>
				     <apply><plus/>
					    <cn cellml:units="dimensionless">0.9</cn>
						<apply><divide/>
						   <ci>C_cast</ci>
						   <cn cellml:units="dimensionless">10.0</cn>
						</apply>
				     </apply>
					 <ci>C1</ci>
				  </apply>
				  <apply><times/>
				     <apply><plus/>
					    <cn cellml:units="dimensionless">0.1</cn>
						<apply><divide/>
						   <ci>C_cast</ci>
						   <cn cellml:units="dimensionless">10.0</cn>
						</apply>
				     </apply>
					 <ci>C2</ci>
				  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   
   <component name="Ca_independent_transient_outward_K_current">
      <variable name="i_Kto_s" units="microA" public_interface="out"/>
      <variable name="E_K" units="millivolt" public_interface="out"/>
      
	  <variable name="a" units="dimensionless" initial_value="0.886"/>
      <variable name="b" units="dimensionless" initial_value="0.114"/>
      <variable name="g_Kto" units="milliS_per_cm2" initial_value="0.35"/>
	  <variable name="g_Ktos" units="milliS"/>
	  <variable name="fKtot" units="dimensionless" initial_value="0.56"/>
	  
	  <variable name="Sms" units="cm2" public_interface="in"/>
	  <variable name="Smt" units="cm2" public_interface="in"/>
	  <variable name="time" units="second" public_interface="in" private_interface="out"/>
      <variable name="Vm_s" units="millivolt" public_interface="in" private_interface="out"/>
      <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="K_e" units="millimolar" public_interface="in"/>
      <variable name="K_i" units="millimolar" public_interface="in"/>
      
	  <variable name="r" units="dimensionless" private_interface="in"/>
      <variable name="s" units="dimensionless" private_interface="in"/>
      <variable name="s_slow" units="dimensionless" private_interface="in"/>
      
	  <math xmlns="http://www.w3.org/1998/Math/MathML">
         <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_e</ci>
                     <ci>K_i</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         
		 <apply>
            <eq/>
            <ci>i_Kto_s</ci>
            <apply>
               <times/>
               <ci>g_Ktos</ci>
               <ci>r</ci>
               <apply>
                  <plus/>
                  <apply>
                     <times/>
                     <ci>a</ci>
                     <ci>s</ci>
                  </apply>
                  <apply>
                     <times/>
                     <ci>b</ci>
                     <ci>s_slow</ci>
                  </apply>
               </apply>
               <apply>
                  <minus/>
                  <ci>Vm_s</ci>
                  <ci>E_K</ci>
               </apply>
            </apply>
         </apply>
		 
		 <apply><eq/>
            <ci>g_Ktos</ci>
            <apply><times/>
               <ci>g_Kto</ci>
               <apply><plus/>
                  <ci>Sms</ci>
				  <ci>Smt</ci>
               </apply>
               <apply><minus/>
                  <cn cellml:units="dimensionless">1</cn>
				  <ci>fKtot</ci>
			   </apply>
			</apply>
         </apply>
      </math>
   </component>
   
   <component name="Ca_independent_transient_outward_K_current_r_gate">
      <variable name="r" units="dimensionless" initial_value="2.3485072e-3" public_interface="out"/>
      <variable name="tau_r" units="second"/>
      <variable name="r_max" units="dimensionless"/>
      <variable name="Vm_s" 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>r_max</ci>
            <apply>
               <divide/>
               <cn cellml:units="dimensionless">1</cn>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <plus/>
                           <ci>Vm_s</ci>
                           <cn cellml:units="millivolt">10.6</cn>
                        </apply>
                        <apply>
                           <minus/>
                           <cn cellml:units="millivolt">11.42</cn>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>tau_r</ci>
            <apply>
               <divide/>
               <cn cellml:units="second">1</cn>
               <apply>
                  <plus/>
                  <apply>
                     <times/>
                     <cn cellml:units="dimensionless">45.16</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <times/>
                           <cn cellml:units="per_millivolt">0.03577</cn>
                           <apply>
                              <plus/>
                              <ci>Vm_s</ci>
                              <cn cellml:units="millivolt">50</cn>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
                  <apply>
                     <times/>
                     <cn cellml:units="dimensionless">98.9</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <times/>
                           <apply>
                              <minus/>
                              <cn cellml:units="per_millivolt">0.1</cn>
                           </apply>
                           <apply>
                              <plus/>
                              <ci>Vm_s</ci>
                              <cn cellml:units="millivolt">38</cn>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>r</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>r_max</ci>
                  <ci>r</ci>
               </apply>
               <ci>tau_r</ci>
            </apply>
         </apply>
      </math>
   </component>
   
   <component name="Ca_independent_transient_outward_K_current_s_gate">
      <variable name="s" units="dimensionless" initial_value="9.8315843e-1" public_interface="out"/>
      <variable name="tau_s" units="second"/>
      <variable name="s_max" units="dimensionless"/>
      <variable name="Vm_s" 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>s_max</ci>
            <apply>
               <divide/>
               <cn cellml:units="dimensionless">1</cn>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <plus/>
                           <ci>Vm_s</ci>
                           <cn cellml:units="millivolt">45.3</cn>
                        </apply>
                        <cn cellml:units="millivolt">6.8841</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>tau_s</ci>
            <apply>
               <plus/>
               <apply>
                  <times/>
                  <cn cellml:units="second">0.35</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <minus/>
                        <apply>
                           <power/>
                           <apply>
                              <times/>
							  <cn cellml:units="per_millivolt">0.0666</cn>
                              <apply>
                                 <plus/>
                                 <ci>Vm_s</ci>
                                 <cn cellml:units="millivolt">70</cn>
                              </apply>
                           </apply>
                           <cn cellml:units="dimensionless">2</cn>
                        </apply>
                     </apply>
                  </apply>
               </apply>
               <cn cellml:units="second">0.035</cn>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>s</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>s_max</ci>
                  <ci>s</ci>
               </apply>
               <ci>tau_s</ci>
            </apply>
         </apply>
      </math>
   </component>
   
   <component name="Ca_independent_transient_outward_K_current_s_slow_gate">
      <variable name="s_slow" units="dimensionless" initial_value="6.7483370e-1" public_interface="out"/>
      <variable name="tau_s_slow" units="second"/>
      <variable name="s_slow_max" units="dimensionless"/>
      <variable name="Vm_s" 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>s_slow_max</ci>
            <apply>
               <divide/>
               <cn cellml:units="dimensionless">1</cn>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <plus/>
                           <ci>Vm_s</ci>
                           <cn cellml:units="millivolt">45.3</cn>
                        </apply>
                        <cn cellml:units="millivolt">6.8841</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>tau_s_slow</ci>
            <apply>
               <plus/>
               <apply>
                  <times/>
                  <cn cellml:units="second">3.7</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <minus/>
                        <apply>
                           <power/>
                           <apply>
                              <times/>
							  <cn cellml:units="per_millivolt">0.0333</cn>
                              <apply>
                                 <plus/>
                                 <ci>Vm_s</ci>
                                 <cn cellml:units="millivolt">70</cn>
                              </apply>
                           </apply>
                           <cn cellml:units="dimensionless">2</cn>
                        </apply>
                     </apply>
                  </apply>
               </apply>
               <cn cellml:units="second">0.035</cn>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>s_slow</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>s_slow_max</ci>
                  <ci>s_slow</ci>
               </apply>
               <ci>tau_s_slow</ci>
            </apply>
         </apply>
      </math>
   </component>
   
   <component name="steady_state_outward_K_current">
      <variable name="i_Kss_s" units="microA" public_interface="out"/>
      
	  <variable name="g_Kss" units="milliS_per_cm2" initial_value="0.07"/>
	  <variable name="g_Ksss" units="milliS"/>
	  <variable name="fKsst" units="dimensionless" initial_value="0.76"/>
	  
	  <variable name="Sms" units="cm2" public_interface="in"/>
	  <variable name="Smt" units="cm2" public_interface="in"/>
      <variable name="time" units="second" public_interface="in" private_interface="out"/>
      <variable name="Vm_s" units="millivolt" public_interface="in" private_interface="out"/>
      <variable name="E_K" units="millivolt" public_interface="in"/>
      <variable name="r_ss" units="dimensionless" private_interface="in"/>
      <variable name="s_ss" units="dimensionless" private_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_Kss_s</ci>
            <apply>
               <times/>
               <ci>g_Ksss</ci>
               <ci>r_ss</ci>
               <ci>s_ss</ci>
               <apply>
                  <minus/>
                  <ci>Vm_s</ci>
                  <ci>E_K</ci>
               </apply>
            </apply>
         </apply>
		 
		 <apply><eq/>
            <ci>g_Ksss</ci>
            <apply><times/>
               <ci>g_Kss</ci>
               <apply><plus/>
                  <ci>Sms</ci>
				  <ci>Smt</ci>
               </apply>
               <apply><minus/>
                  <cn cellml:units="dimensionless">1</cn>
				  <ci>fKsst</ci>
			   </apply>
			</apply>
         </apply>
      </math>
   </component>
  
   <component name="steady_state_outward_K_current_r_ss_gate">
      <variable name="r_ss" units="dimensionless" initial_value="3.1079213e-3" public_interface="out"/>
      <variable name="tau_r_ss" units="second"/>
      <variable name="r_ss_max" units="dimensionless"/>
      <variable name="Vm_s" 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>r_ss_max</ci>
            <apply>
               <divide/>
               <cn cellml:units="dimensionless">1</cn>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <minus/>
						   <apply>
                              <plus/>
                              <ci>Vm_s</ci>
                              <cn cellml:units="millivolt">11.5</cn>
                           </apply>
						</apply>
                        <cn cellml:units="millivolt">11.82</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>tau_r_ss</ci>
            <apply>
               <divide/>
               <cn cellml:units="second">10</cn>
               <apply>
                  <plus/>
                  <apply>
                     <times/>
                     <cn cellml:units="dimensionless">45.16</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <times/>
                           <cn cellml:units="per_millivolt">0.03577</cn>
                           <apply>
                              <plus/>
                              <ci>Vm_s</ci>
                              <cn cellml:units="millivolt">50</cn>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
                  <apply>
                     <times/>
                     <cn cellml:units="dimensionless">98.9</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <times/>
                           <apply>
                              <minus/>
                              <cn cellml:units="per_millivolt">0.1</cn>
                           </apply>
                           <apply>
                              <plus/>
                              <ci>Vm_s</ci>
                              <cn cellml:units="millivolt">38</cn>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>r_ss</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>r_ss_max</ci>
                  <ci>r_ss</ci>
               </apply>
               <ci>tau_r_ss</ci>
            </apply>
         </apply>
      </math>
   </component>
   
   <component name="steady_state_outward_K_current_s_ss_gate">
      <variable name="s_ss" units="dimensionless" initial_value="3.0773393e-1" public_interface="out"/>
      <variable name="tau_s_ss" units="second"/>
      <variable name="s_ss_max" units="dimensionless"/>
      <variable name="Vm_s" 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>s_ss_max</ci>
            <apply>
               <divide/>
               <cn cellml:units="dimensionless">1</cn>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <plus/>
                           <ci>Vm_s</ci>
                           <cn cellml:units="millivolt">87.5</cn>
                        </apply>
                        <cn cellml:units="millivolt">10.3</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>tau_s_ss</ci>
            <cn cellml:units="second">2.1</cn>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>s_ss</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>s_ss_max</ci>
                  <ci>s_ss</ci>
               </apply>
               <ci>tau_s_ss</ci>
            </apply>
         </apply>
      </math>
   </component>
   
   <component name="inward_rectifier">
      <variable name="i_K1_s" units="microA" public_interface="out"/>
      
	  <variable name="g_K1" units="milliS_per_cm2" initial_value="0.24"/>
	  <variable name="g_K1s" units="milliS"/>
	  <variable name="fK1t" units="dimensionless" initial_value="0.56"/>
	  
	  <variable name="Sms" units="cm2" public_interface="in"/>
	  <variable name="Smt" units="cm2" public_interface="in"/>
      <variable name="time" units="second" public_interface="in"/>
      <variable name="Vm_s" 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="K_e" units="millimolar" public_interface="in"/>
      <variable name="E_K" units="millivolt" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_K1_s</ci>
            <apply>
               <times/>
               <apply>
                  <divide/>
                  <apply>
                     <times/>
                     <ci>g_K1s</ci>
                     <apply>
                        <minus/>
                        <ci>Vm_s</ci>
                        <apply>
                           <plus/>
                           <ci>E_K</ci>
                           <cn cellml:units="millivolt">1.73</cn>
                        </apply>
                     </apply>
                  </apply>
                  <apply>
                     <plus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <times/>
                           <apply>
                              <divide/>
                              <apply>
                                 <times/>
                                 <cn cellml:units="dimensionless">1.613</cn>
                                 <ci>F</ci>
                              </apply>
                              <apply>
                                 <times/>
                                 <ci>R</ci>
                                 <ci>T</ci>
                              </apply>
                           </apply>
                           <apply>
                              <minus/>
                              <ci>Vm_s</ci>
                              <apply>
                                 <plus/>
                                 <ci>E_K</ci>
                                 <cn cellml:units="millivolt">1.73</cn>
                              </apply>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <minus/>
                           <apply>
                              <minus/>
                              <ci>K_e</ci>
                              <cn cellml:units="millimolar">0.9988</cn>
                           </apply>
                        </apply>
                        <cn cellml:units="millimolar">0.124</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
		 
		 <apply><eq/>
            <ci>g_K1s</ci>
            <apply><times/>
               <ci>g_K1</ci>
               <apply><plus/>
                  <ci>Sms</ci>
				  <ci>Smt</ci>
               </apply>
               <apply><minus/>
                  <cn cellml:units="dimensionless">1</cn>
				  <ci>fK1t</ci>
			   </apply>
			</apply>
         </apply>
      </math>
   </component>
   
   <component name="hyperpolarisation_activated_current">
      <variable name="i_f_s" units="microA" public_interface="out"/>
      <variable name="i_fNa_s" units="microA" public_interface="out"/>
      <variable name="i_fK_s" units="microA" public_interface="out"/>
      
	  <variable name="g_f" units="milliS_per_cm2" initial_value="0.0145"/>
	  <variable name="g_fs" units="milliS"/>
	  <variable name="fft" units="dimensionless" initial_value="0.56"/>
	  
	  <variable name="Sms" units="cm2" public_interface="in"/>
	  <variable name="Smt" units="cm2" public_interface="in"/>
      <variable name="f_Na" units="dimensionless" initial_value="0.2"/>
      <variable name="f_K" units="dimensionless"/>
      <variable name="time" units="second" public_interface="in" private_interface="out"/>
      <variable name="Vm_s" units="millivolt" public_interface="in" private_interface="out"/>
      <variable name="E_K" units="millivolt" public_interface="in"/>
      <variable name="E_Na" units="millivolt" public_interface="in"/>
      <variable name="y" units="dimensionless" private_interface="in"/>
      
	  <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>f_K</ci>
            <apply>
               <minus/>
               <cn cellml:units="dimensionless">1</cn>
               <ci>f_Na</ci>
            </apply>
         </apply>
         
		 <apply>
            <eq/>
            <ci>i_fNa_s</ci>
            <apply>
               <times/>
               <ci>g_fs</ci>
               <ci>y</ci>
               <ci>f_Na</ci>
               <apply>
                  <minus/>
                  <ci>Vm_s</ci>
                  <ci>E_Na</ci>
               </apply>
            </apply>
         </apply>
         
		 <apply>
            <eq/>
            <ci>i_fK_s</ci>
            <apply>
               <times/>
               <ci>g_fs</ci>
               <ci>y</ci>
               <ci>f_K</ci>
               <apply>
                  <minus/>
                  <ci>Vm_s</ci>
                  <ci>E_K</ci>
               </apply>
            </apply>
         </apply>
         
		 <apply>
            <eq/>
            <ci>i_f_s</ci>
            <apply>
               <times/>
               <ci>g_fs</ci>
               <ci>y</ci>
               <apply>
                  <plus/>
				  <apply>
                     <times/>
					 <cn cellml:units="dimensionless">0.2</cn>
				     <apply>
                        <minus/>
                        <ci>Vm_s</ci>
                        <ci>E_Na</ci>
					 </apply>
                  </apply>
				  <apply>
                     <times/>
					 <cn cellml:units="dimensionless">0.8</cn>
				     <apply>
                        <minus/>
                        <ci>Vm_s</ci>
                        <ci>E_K</ci>
					 </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
		 
		 <apply><eq/>
            <ci>g_fs</ci>
            <apply><times/>
               <ci>g_f</ci>
               <apply><plus/>
                  <ci>Sms</ci>
				  <ci>Smt</ci>
               </apply>
               <apply><minus/>
                  <cn cellml:units="dimensionless">1</cn>
				  <ci>fft</ci>
			   </apply>
			</apply>
         </apply>
      </math>
   </component>
   
   <component name="hyperpolarisation_activated_current_y_gate">
      <variable name="y" units="dimensionless" initial_value="3.4357356e-3" public_interface="out"/>
      <variable name="tau_y" units="second"/>
      <variable name="y_infinity" units="dimensionless"/>
      <variable name="Vm_s" 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>y_infinity</ci>
            <apply>
               <divide/>
               <cn cellml:units="dimensionless">1</cn>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <plus/>
                           <ci>Vm_s</ci>
                           <cn cellml:units="millivolt">138.6</cn>
                        </apply>
                        <cn cellml:units="millivolt">10.48</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>tau_y</ci>
            <apply>
               <divide/>
               <cn cellml:units="second">1</cn>
               <apply>
                  <plus/>
                  <apply>
                     <times/>
                     <cn cellml:units="dimensionless">0.11885</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <plus/>
                              <ci>Vm_s</ci>
                              <cn cellml:units="millivolt">80</cn>
                           </apply>
                           <cn cellml:units="millivolt">28.37</cn>
                        </apply>
                     </apply>
                  </apply>
                  <apply>
                     <times/>
                     <cn cellml:units="dimensionless">0.5623</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <plus/>
                              <ci>Vm_s</ci>
                              <cn cellml:units="millivolt">80</cn>
                           </apply>
                           <apply>
                              <minus/>
                              <cn cellml:units="millivolt">14.19</cn>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>y</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>y_infinity</ci>
                  <ci>y</ci>
               </apply>
               <ci>tau_y</ci>
            </apply>
         </apply>
      </math>
   </component>
   
   <component name="background_currents">
      <variable name="i_B_s" units="microA" public_interface="out"/>
      <variable name="i_BNa_s" units="microA" public_interface="out"/>
      <variable name="i_BCa_s" units="microA" public_interface="out"/>
      <variable name="i_BK_s" units="microA" public_interface="out"/>
      
	  <variable name="g_B_Ca" units="milliS_per_cm2" initial_value="0.000648"/>
	  <variable name="g_B_Na" units="milliS_per_cm2" initial_value="0.0008015"/>
	  <variable name="g_B_K" units="milliS_per_cm2" initial_value="0.00138"/>
	  <variable name="g_B_Cas" units="milliS"/>
	  <variable name="g_B_Nas" units="milliS"/>
	  <variable name="g_B_Ks" units="milliS"/>
	  <variable name="fCabt" units="dimensionless" initial_value="0.56"/>
	  <variable name="fNabt" units="dimensionless" initial_value="0.56"/>
	  <variable name="fKbt" units="dimensionless" initial_value="0.56"/>
	  
	  <variable name="Sms" units="cm2" public_interface="in"/>
	  <variable name="Smt" units="cm2" public_interface="in"/>
	  <variable name="E_Ca" units="millivolt" initial_value="65"/>
      <variable name="E_Na" units="millivolt" public_interface="in"/>
      <variable name="E_K" units="millivolt" public_interface="in"/>
      <variable name="Ca_e" 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="T" units="kelvin" public_interface="in"/>
      <variable name="F" units="coulomb_per_mole" public_interface="in"/>
      <variable name="Vm_s" units="millivolt" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_BNa_s</ci>
            <apply>
               <times/>
               <ci>g_B_Nas</ci>
               <apply>
                  <minus/>
                  <ci>Vm_s</ci>
                  <ci>E_Na</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>i_BCa_s</ci>
            <apply>
               <times/>
               <ci>g_B_Cas</ci>
               <apply>
                  <minus/>
                  <ci>Vm_s</ci>
                  <ci>E_Ca</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>i_BK_s</ci>
            <apply>
               <times/>
               <ci>g_B_Ks</ci>
               <apply>
                  <minus/>
                  <ci>Vm_s</ci>
                  <ci>E_K</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>i_B_s</ci>
            <apply>
               <plus/>
               <ci>i_BNa_s</ci>
               <ci>i_BCa_s</ci>
               <ci>i_BK_s</ci>
            </apply>
         </apply>
		 
		 <apply><eq/>
            <ci>g_B_Cas</ci>
            <apply><times/>
               <ci>g_B_Ca</ci>
               <apply><plus/>
                  <ci>Sms</ci>
				  <ci>Smt</ci>
               </apply>
               <apply><minus/>
                  <cn cellml:units="dimensionless">1</cn>
				  <ci>fCabt</ci>
			   </apply>
			</apply>
         </apply>
		 
		 <apply><eq/>
            <ci>g_B_Nas</ci>
            <apply><times/>
               <ci>g_B_Na</ci>
               <apply><plus/>
                  <ci>Sms</ci>
				  <ci>Smt</ci>
               </apply>
               <apply><minus/>
                  <cn cellml:units="dimensionless">1</cn>
				  <ci>fNabt</ci>
			   </apply>
			</apply>
         </apply>
		 
		 <apply><eq/>
            <ci>g_B_Ks</ci>
            <apply><times/>
               <ci>g_B_K</ci>
               <apply><plus/>
                  <ci>Sms</ci>
				  <ci>Smt</ci>
               </apply>
               <apply><minus/>
                  <cn cellml:units="dimensionless">1</cn>
				  <ci>fKbt</ci>
			   </apply>
			</apply>
         </apply>
      </math>
   </component>
  
   <component name="sodium_potassium_pump">
      <variable name="i_NaK_s" units="microA" public_interface="out"/>
      
	  <variable name="i_NaK_max" units="microA_per_cm2" initial_value="1"/>
	  <variable name="i_NaK_max_s" units="microA"/>
	  <variable name="fNaKt" units="dimensionless" initial_value="0.59"/>
	  
	  <variable name="Sms" units="cm2" public_interface="in"/>
	  <variable name="Smt" units="cm2" public_interface="in"/>
	  <variable name="K_m_K" units="millimolar" initial_value="1.5"/>
      <variable name="K_m_Na" units="millimolar" initial_value="10"/>
      <variable name="sigma" units="dimensionless"/>
      <variable name="K_e" units="millimolar" public_interface="in"/>
      <variable name="Na_e" 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="Vm_s" units="millivolt" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>sigma</ci>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <ci>Na_e</ci>
                        <cn cellml:units="millimolar">67.3</cn>
                     </apply>
                  </apply>
                  <cn cellml:units="dimensionless">1</cn>
               </apply>
               <cn cellml:units="dimensionless">7</cn>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>i_NaK_s</ci>
            <apply>
               <divide/>
               <apply>
                  <divide/>
                  <apply>
                     <times/>
                     <apply>
                        <divide/>
                        <ci>i_NaK_max_s</ci>
                        <apply>
                           <plus/>
                           <cn cellml:units="dimensionless">1</cn>
                           <apply>
                              <times/>
                              <cn cellml:units="dimensionless">0.1245</cn>
                              <apply>
                                 <exp/>
                                 <apply>
                                    <divide/>
                                    <apply>
                                       <times/>
                                       <apply>
                                          <minus/>
                                          <cn cellml:units="dimensionless">0.1</cn>
                                       </apply>
                                       <ci>Vm_s</ci>
                                       <ci>F</ci>
                                    </apply>
                                    <apply>
                                       <times/>
                                       <ci>R</ci>
                                       <ci>T</ci>
                                    </apply>
                                 </apply>
                              </apply>
                           </apply>
                           <apply>
                              <times/>
                              <cn cellml:units="dimensionless">0.0365</cn>
                              <ci>sigma</ci>
                              <apply>
                                 <exp/>
                                 <apply>
                                    <divide/>
                                    <apply>
                                       <times/>
                                       <apply>
                                          <minus/>
                                          <ci>Vm_s</ci>
                                       </apply>
                                       <ci>F</ci>
                                    </apply>
                                    <apply>
                                       <times/>
                                       <ci>R</ci>
                                       <ci>T</ci>
                                    </apply>
                                 </apply>
                              </apply>
                           </apply>
                        </apply>
                     </apply>
                     <ci>K_e</ci>
                  </apply>
                  <apply>
                     <plus/>
                     <ci>K_e</ci>
                     <ci>K_m_K</ci>
                  </apply>
               </apply>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <power/>
                     <apply>
                        <divide/>
                        <ci>K_m_Na</ci>
                        <ci>Na_i</ci>
                     </apply>
                     <cn cellml:units="dimensionless">1.5</cn>
                  </apply>
               </apply>
            </apply>
         </apply>
		 
		 <apply><eq/>
            <ci>i_NaK_max_s</ci>
            <apply><times/>
               <ci>i_NaK_max</ci>
               <apply><plus/>
                  <ci>Sms</ci>
				  <ci>Smt</ci>
               </apply>
               <apply><minus/>
                  <cn cellml:units="dimensionless">1</cn>
				  <ci>fNaKt</ci>
			   </apply>
			</apply>
         </apply>
      </math>
   </component>
   
   <component name="sarcolemmal_calcium_pump_current">
      <variable name="i_CaP_s" units="microA" public_interface="out"/>
      
	  <variable name="i_CaP_max" units="microA_per_cm2" initial_value="0.85"/>
	  <variable name="i_CaP_max_s" units="microA"/>
	  <variable name="fCaPt" units="dimensionless" initial_value="0.56"/>
	  
	  <variable name="Sms" units="cm2" public_interface="in"/>
	  <variable name="Smt" units="cm2" public_interface="in"/>
      <variable name="Ca_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_CaP_s</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <ci>i_CaP_max_s</ci>
                  <ci>Ca_i</ci>
               </apply>
               <apply>
                  <plus/>
                  <ci>Ca_i</ci>
                  <cn cellml:units="millimolar">0.0004</cn>
               </apply>
            </apply>
         </apply>
		 
		 <apply><eq/>
            <ci>i_CaP_max_s</ci>
            <apply><times/>
               <ci>i_CaP_max</ci>
               <apply><plus/>
                  <ci>Sms</ci>
				  <ci>Smt</ci>
               </apply>
               <apply><minus/>
                  <cn cellml:units="dimensionless">1</cn>
				  <ci>fCaPt</ci>
			   </apply>
			</apply>
         </apply>
      </math>
   </component>
   
   <component name="Na_Ca_ion_exchanger_current">
      <variable name="i_NaCa_s" units="microA" public_interface="out"/>
      
	  <variable name="i_NaCa_max" units="microA_per_cm2" initial_value="1.8e-4"/>
	  <variable name="i_NaCa_max_s" units="microA"/>
	  <variable name="fNaCat" units="dimensionless" initial_value="0.81"/>
	  
	  <variable name="Sms" units="cm2" public_interface="in"/>
	  <variable name="Smt" units="cm2" public_interface="in"/>
	  <variable name="Na_i" units="millimolar" public_interface="in"/>
      <variable name="Na_e" units="millimolar" public_interface="in"/>
      <variable name="Ca_i" units="millimolar" public_interface="in"/>
      <variable name="Ca_e" units="millimolar" public_interface="in"/>
      <variable name="Vm_s" 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_s</ci>
            <apply><divide/>
               <apply><times/>
                  <ci>i_NaCa_max_s</ci>
				  <apply><minus/>
				     <apply><times/>
				        <apply><exp/>
                           <apply><times/>
                              <cn cellml:units="per_millivolt">0.0187</cn>
                              <ci>Vm_s</ci>
                           </apply>
                        </apply>
				        <apply><power/>
                           <ci>Na_i</ci>
                           <cn cellml:units="dimensionless">3</cn>
                        </apply>
                        <ci>Ca_e</ci>
					 </apply>
                     <apply><times/>
				        <apply><exp/>
                           <apply><times/>
                              <cn cellml:units="per_millivolt">-0.0187</cn>
                              <ci>Vm_s</ci>
                           </apply>
                        </apply>
				        <apply><power/>
                           <ci>Na_e</ci>
                           <cn cellml:units="dimensionless">3</cn>
                        </apply>
                        <ci>Ca_i</ci>
					 </apply>
				  </apply>
			   </apply>
			   <apply><plus/>
                  <cn cellml:units="millimolar4">1</cn>
                  <apply><times/>
                     <cn cellml:units="dimensionless">0.0001</cn>
                     <apply><plus/>
                        <apply><times/>
                           <ci>Ca_i</ci>
                           <apply><power/>
                              <ci>Na_e</ci>
                              <cn cellml:units="dimensionless">3</cn>
                           </apply>
                        </apply>
                        <apply><times/>
                           <ci>Ca_e</ci>
                           <apply><power/>
                              <ci>Na_i</ci>
                              <cn cellml:units="dimensionless">3</cn>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
		 
		 <apply><eq/>
            <ci>i_NaCa_max_s</ci>
            <apply><times/>
               <ci>i_NaCa_max</ci>
               <apply><plus/>
                  <ci>Sms</ci>
				  <ci>Smt</ci>
               </apply>
               <apply><minus/>
                  <cn cellml:units="dimensionless">1</cn>
				  <ci>fNaCat</ci>
			   </apply>
			</apply>
         </apply>
      </math>
   </component>
   
   <group>
      <relationship_ref relationship="containment"/>
      <component_ref component="membrane">
         <component_ref component="sodium_current">
            <component_ref component="sodium_current_m_gate"/>
            <component_ref component="sodium_current_h_gate"/>
         </component_ref>
         <component_ref component="calcium_current"/>
         <component_ref component="Ca_independent_transient_outward_K_current">
            <component_ref component="Ca_independent_transient_outward_K_current_r_gate"/>
            <component_ref component="Ca_independent_transient_outward_K_current_s_gate"/>
            <component_ref component="Ca_independent_transient_outward_K_current_s_slow_gate"/>
         </component_ref>
         <component_ref component="steady_state_outward_K_current">
            <component_ref component="steady_state_outward_K_current_r_ss_gate"/>
            <component_ref component="steady_state_outward_K_current_s_ss_gate"/>
         </component_ref>
         <component_ref component="inward_rectifier"/>
         <component_ref component="hyperpolarisation_activated_current">
            <component_ref component="hyperpolarisation_activated_current_y_gate"/>
         </component_ref>
         <component_ref component="background_currents"/>
         <component_ref component="sodium_potassium_pump"/>
         <component_ref component="sarcolemmal_calcium_pump_current"/>
         <component_ref component="Na_Ca_ion_exchanger_current"/>
      </component_ref>
   </group>
   
   <group>
      <relationship_ref relationship="encapsulation"/>
      <component_ref component="surface_membrane">
	     <component_ref component="membrane"/>
		 <component_ref component="sodium_current">
            <component_ref component="sodium_current_m_gate"/>
            <component_ref component="sodium_current_h_gate"/>
         </component_ref>
		 <component_ref component="calcium_current"/>
         <component_ref component="Ca_independent_transient_outward_K_current">
            <component_ref component="Ca_independent_transient_outward_K_current_r_gate"/>
            <component_ref component="Ca_independent_transient_outward_K_current_s_gate"/>
            <component_ref component="Ca_independent_transient_outward_K_current_s_slow_gate"/>
         </component_ref>
         <component_ref component="steady_state_outward_K_current">
            <component_ref component="steady_state_outward_K_current_r_ss_gate"/>
            <component_ref component="steady_state_outward_K_current_s_ss_gate"/>
         </component_ref>
		 <component_ref component="inward_rectifier"/>
         <component_ref component="hyperpolarisation_activated_current">
            <component_ref component="hyperpolarisation_activated_current_y_gate"/>
         </component_ref>
		 <component_ref component="background_currents"/>
         <component_ref component="sodium_potassium_pump"/>
         <component_ref component="sarcolemmal_calcium_pump_current"/>
         <component_ref component="Na_Ca_ion_exchanger_current"/>
      </component_ref>
   </group>
   
   
   <connection>
      <map_components component_1="membrane" component_2="sodium_current"/>
      <map_variables variable_1="Vm_s" variable_2="Vm_s"/>
      <map_variables variable_1="i_Na_s" variable_2="i_Na_s"/>
   </connection>
   <connection>
      <map_components component_1="membrane" component_2="calcium_current"/>
      <map_variables variable_1="Vm_s" variable_2="Vm_s"/>
      <map_variables variable_1="i_Ca_s" variable_2="i_Ca_s"/>
   </connection>
   <connection>
      <map_components component_1="membrane" component_2="Ca_independent_transient_outward_K_current"/>
      <map_variables variable_1="Vm_s" variable_2="Vm_s"/>
      <map_variables variable_1="i_Kto_s" variable_2="i_Kto_s"/>
   </connection>
   <connection>
      <map_components component_1="membrane" component_2="steady_state_outward_K_current"/>
      <map_variables variable_1="Vm_s" variable_2="Vm_s"/>
      <map_variables variable_1="i_Kss_s" variable_2="i_Kss_s"/>
   </connection>
   <connection>
      <map_components component_1="membrane" component_2="inward_rectifier"/>
      <map_variables variable_1="Vm_s" variable_2="Vm_s"/>
      <map_variables variable_1="i_K1_s" variable_2="i_K1_s"/>
   </connection>
   <connection>
      <map_components component_1="membrane" component_2="hyperpolarisation_activated_current"/>
      <map_variables variable_1="Vm_s" variable_2="Vm_s"/>
      <map_variables variable_1="i_f_s" variable_2="i_f_s"/>
   </connection>
   <connection>
      <map_components component_1="membrane" component_2="background_currents"/>
      <map_variables variable_1="Vm_s" variable_2="Vm_s"/>
      <map_variables variable_1="i_B_s" variable_2="i_B_s"/>
   </connection>
   <connection>
      <map_components component_1="membrane" component_2="sodium_potassium_pump"/>
      <map_variables variable_1="Vm_s" variable_2="Vm_s"/>
      <map_variables variable_1="i_NaK_s" variable_2="i_NaK_s"/>
   </connection>
   <connection>
      <map_components component_1="membrane" component_2="sarcolemmal_calcium_pump_current"/>
      <map_variables variable_1="i_CaP_s" variable_2="i_CaP_s"/>
   </connection>
   <connection>
      <map_components component_1="membrane" component_2="Na_Ca_ion_exchanger_current"/>
      <map_variables variable_1="Vm_s" variable_2="Vm_s"/>
      <map_variables variable_1="i_NaCa_s" variable_2="i_NaCa_s"/>
   </connection>
   <connection>
      <map_components component_1="sodium_current" component_2="background_currents"/>
      <map_variables variable_1="E_Na" variable_2="E_Na"/>
   </connection>
   <connection>
      <map_components component_1="sodium_current" component_2="hyperpolarisation_activated_current"/>
      <map_variables variable_1="E_Na" variable_2="E_Na"/>
   </connection>
   <connection>
      <map_components component_1="Ca_independent_transient_outward_K_current" component_2="steady_state_outward_K_current"/>
      <map_variables variable_1="E_K" variable_2="E_K"/>
   </connection>
   <connection>
      <map_components component_1="Ca_independent_transient_outward_K_current" component_2="inward_rectifier"/>
      <map_variables variable_1="E_K" variable_2="E_K"/>
   </connection>
   <connection>
      <map_components component_1="Ca_independent_transient_outward_K_current" component_2="hyperpolarisation_activated_current"/>
      <map_variables variable_1="E_K" variable_2="E_K"/>
   </connection>
   <connection>
      <map_components component_1="Ca_independent_transient_outward_K_current" component_2="background_currents"/>
      <map_variables variable_1="E_K" variable_2="E_K"/>
   </connection>
   <connection>
      <map_components component_1="sodium_current" component_2="sodium_current_m_gate"/>
      <map_variables variable_1="m" variable_2="m"/>
      <map_variables variable_1="time" variable_2="time"/>
      <map_variables variable_1="Vm_s" variable_2="Vm_s"/>
   </connection>
   <connection>
      <map_components component_1="sodium_current" component_2="sodium_current_h_gate"/>
      <map_variables variable_1="h" variable_2="h"/>
      <map_variables variable_1="time" variable_2="time"/>
      <map_variables variable_1="Vm_s" variable_2="Vm_s"/>
   </connection>
   <connection>
      <map_components component_1="Ca_independent_transient_outward_K_current" component_2="Ca_independent_transient_outward_K_current_r_gate"/>
      <map_variables variable_1="r" variable_2="r"/>
      <map_variables variable_1="time" variable_2="time"/>
      <map_variables variable_1="Vm_s" variable_2="Vm_s"/>
   </connection>
   <connection>
      <map_components component_1="Ca_independent_transient_outward_K_current" component_2="Ca_independent_transient_outward_K_current_s_gate"/>
      <map_variables variable_1="s" variable_2="s"/>
      <map_variables variable_1="time" variable_2="time"/>
      <map_variables variable_1="Vm_s" variable_2="Vm_s"/>
   </connection>
   <connection>
      <map_components component_1="Ca_independent_transient_outward_K_current" component_2="Ca_independent_transient_outward_K_current_s_slow_gate"/>
      <map_variables variable_1="s_slow" variable_2="s_slow"/>
      <map_variables variable_1="time" variable_2="time"/>
      <map_variables variable_1="Vm_s" variable_2="Vm_s"/>
   </connection>
   <connection>
      <map_components component_1="steady_state_outward_K_current" component_2="steady_state_outward_K_current_r_ss_gate"/>
      <map_variables variable_1="r_ss" variable_2="r_ss"/>
      <map_variables variable_1="time" variable_2="time"/>
      <map_variables variable_1="Vm_s" variable_2="Vm_s"/>
   </connection>
   <connection>
      <map_components component_1="steady_state_outward_K_current" component_2="steady_state_outward_K_current_s_ss_gate"/>
      <map_variables variable_1="s_ss" variable_2="s_ss"/>
      <map_variables variable_1="time" variable_2="time"/>
      <map_variables variable_1="Vm_s" variable_2="Vm_s"/>
   </connection>
   <connection>
      <map_components component_1="hyperpolarisation_activated_current" component_2="hyperpolarisation_activated_current_y_gate"/>
      <map_variables variable_1="y" variable_2="y"/>
      <map_variables variable_1="time" variable_2="time"/>
      <map_variables variable_1="Vm_s" variable_2="Vm_s"/>
   </connection>
</model>