<?xml version="1.0" encoding="ISO-8859-1"?><!DOCTYPE model PUBLIC "-//APACHE//DTD Documentation V1.3//EN" "document-v13.dtd">
<model xmlns="http://www.cellml.org/cellml/1.0#" xmlns:cellml="http://www.cellml.org/cellml/1.0#" xmlns:cmeta="http://www.cellml.org/metadata/1.0#" name="ICC_model" cmeta:id="ICC_model_no_IP3">
  <rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#" xmlns:bqs="http://www.cellml.org/bqs/1.0#" xmlns:dc="http://purl.org/dc/elements/1.1/" xmlns:dcterms="http://purl.org/dc/terms/" xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" xmlns:cs="http://www.cellml.org/metadata/simulation/1.0#" xmlns:cg="http://www.cellml.org/metadata/graphs/1.0#">
    
    <rdf:Description rdf:about="">
      <dc:creator rdf:parseType="Resource">
        <vCard:N rdf:parseType="Resource">
          <vCard:Family>Corrias</vCard:Family>
          <vCard:Given>Alberto</vCard:Given>
        </vCard:N>
        <vCard:EMAIL rdf:parseType="Resource">
          <rdf:value>albertocorrias@nus.edu.sg</rdf:value>
          <rdf:type rdf:resource="http://imc.org/vCard/3.0#internet"/>
        </vCard:EMAIL>
        <vCard:ORG rdf:parseType="Resource">
          <vCard:Orgname>National University of Singapore</vCard:Orgname>
          <vCard:Orgunit>Division of BioEngineering</vCard:Orgunit>
        </vCard:ORG>
      </dc:creator>
      <dcterms:created rdf:parseType="Resource">
        <dcterms:W3CDTF>2008-03-01</dcterms:W3CDTF>
      </dcterms:created>
      <dc:publisher>
        National University of Singapore
      </dc:publisher>
    </rdf:Description>
    <rdf:Description rdf:about="#ICC_model_no_IP3">
      <!--  Keyword(s) -->
      <bqs:reference rdf:parseType="Resource">
        <dc:subject rdf:parseType="Resource">
          <bqs:subject_type>keyword</bqs:subject_type>
          <rdf:value>
            <rdf:Bag>
              <rdf:li>electrophysiology</rdf:li>
              <rdf:li>gastric</rdf:li>
              <rdf:li>ICC</rdf:li>
            </rdf:Bag>
          </rdf:value>
        </dc:subject>
      </bqs:reference>
      
      <bqs:reference rdf:parseType="Resource">
        <bqs:Pubmed_id>18276830</bqs:Pubmed_id>
        <bqs:JournalArticle rdf:parseType="Resource">
          <dc:creator>
            <rdf:Seq>
              <rdf:li rdf:parseType="Resource">
                <bqs:Person rdf:parseType="Resource">
                  <vCard:N rdf:parseType="Resource">
                    <vCard:Family>Corrias</vCard:Family>
                    <vCard:Given>Alberto</vCard:Given>
                  </vCard:N>
                </bqs:Person>
              </rdf:li>
              <rdf:li rdf:parseType="Resource">
                <bqs:Person rdf:parseType="Resource">
                  <vCard:N rdf:parseType="Resource">
                    <vCard:Family>Buist</vCard:Family>
                    <vCard:Given>Martin</vCard:Given>
                    <vCard:Other>L</vCard:Other>
                  </vCard:N>
                </bqs:Person>
              </rdf:li>
            </rdf:Seq>
          </dc:creator>
          <dc:title>
            Quantitative cellular description of gastric slow wave activity
          </dc:title>
          <dcterms:issued rdf:parseType="Resource">
            <dcterms:W3CDTF>2008</dcterms:W3CDTF>
          </dcterms:issued>
          <bqs:Journal rdf:parseType="Resource">
            <dc:title>
              Am J Physiol Gastrointest Liver Physiol
            </dc:title>
          </bqs:Journal>
          <bqs:volume>294</bqs:volume>
          <bqs:first_page>G989</bqs:first_page>
          <bqs:last_page>G995</bqs:last_page>
        </bqs:JournalArticle>
      </bqs:reference>
      <cs:simulation>
        <rdf:Description rdf:ID="simulation">
          <cs:simulationName>Control ICC simulation</cs:simulationName>
          <cs:multistepMethod>bdf</cs:multistepMethod>
          <cs:iterationMethod>newton</cs:iterationMethod>
          <cs:linearSolver>dense</cs:linearSolver>
          <cs:boundIntervals rdf:parseType="Collection">
            <rdf:Description>
              <cs:boundVariable>
                <rdf:Description rdf:about="#time"/>
              </cs:boundVariable>
              <cs:maximumStepSize rdf:datatype="http://www.w3.org/2001/XMLSchema#double">
                0.1
              </cs:maximumStepSize>
              <cs:tabulationStepSize rdf:datatype="http://www.w3.org/2001/XMLSchema#double">
                0.1
              </cs:tabulationStepSize>
              <cs:startingValue rdf:datatype="http://www.w3.org/2001/XMLSchema#double">
                0.0
              </cs:startingValue>
              <cs:endingValue rdf:datatype="http://www.w3.org/2001/XMLSchema#double">
                300
              </cs:endingValue>
            </rdf:Description>
          </cs:boundIntervals>
        </rdf:Description>
      </cs:simulation>
      <cg:graph>
        <rdf:Description rdf:ID="membrane_potential_graph_no_IP3">
          <cg:x-label>Time (s)</cg:x-label>
          <cg:y-label>Membrane potential (mV)</cg:y-label>
          <cg:background-colour>#ffffff</cg:background-colour>
          <cg:traces rdf:parseType="Collection">
            <rdf:Description>
              <cg:type rdf:resource="http://www.cellml.org/metadata/graphs/1.0#line"/>
              <cg:colour>#ff0000</cg:colour>
              <cg:line-type rdf:datatype="http://www.w3.org/2001/XMLSchema#integer">1</cg:line-type>
              <cg:x-variable rdf:parseType="Resource">
                <cs:simulation rdf:resource="#simulation"/>
                <cg:variable rdf:resource="#time"/>
              </cg:x-variable>
              <cg:y-variable rdf:parseType="Resource">
                <cs:simulation rdf:resource="#simulation"/>
                <cg:variable rdf:resource="#Vm"/>
              </cg:y-variable>
            </rdf:Description>
          </cg:traces>
        </rdf:Description>
      </cg:graph>
    </rdf:Description>
  </rdf:RDF>
   <units name="conductance_units">
      <unit units="siemens" prefix="nano"/>
   </units>
   <units name="millifarads">
      <unit units="farad" prefix="milli"/>
   </units>
   <units name="voltage_units">
      <unit units="volt" prefix="milli"/>
   </units>
   <units name="volume_units">
      <unit units="litre"/>
   </units>
   <units name="Inverse_Voltage_units">
      <unit units="voltage_units" exponent="-1"/>
   </units>
   <units name="Inverse_Voltage_Units_per_time_units">
      <unit units="Inverse_Voltage_units"/>
      <unit units="time_units" exponent="-1"/>
   </units>
   <units name="time_units">
      <unit units="second"/>
   </units>
   <units name="current_units">
      <unit units="voltage_units"/>
      <unit units="conductance_units"/>
   </units>
   <units name="rate_constants_units">
      <unit units="time_units" exponent="-1"/>
   </units>
   <units name="capacitance_units">
      <unit units="farad" prefix="nano"/>
   </units>
   <units name="Temperature_units">
      <unit units="kelvin"/>
   </units>
   <units name="Inverse_Temperature_units_times_conductance">
      <unit units="Temperature_units" exponent="-1"/>
      <unit units="conductance_units"/>
   </units>
   <units name="F_units">
      <unit units="mole" prefix="nano" exponent="-1"/>
      <unit units="coulomb" prefix="micro"/>
   </units>
   <units name="R_units">
      <unit units="joule" prefix="nano"/>
      <unit units="mole" prefix="nano" exponent="-1"/>
      <unit units="kelvin" exponent="-1"/>
   </units>
   <units name="millimolar">
      <unit units="mole" prefix="milli"/>
      <unit units="litre" exponent="-1"/>
   </units>
   <units name="per_millimolar">
      <unit units="litre"/>
      <unit units="mole" prefix="milli" exponent="-1"/>
   </units>
   <units name="millimolar_per_second">
      <unit units="mole" prefix="milli"/>
      <unit units="litre" exponent="-1"/>
      <unit units="second" exponent="-1"/>
   </units>
   <units name="millimolar_per_second_per_millivolt">
      <unit units="mole" prefix="milli"/>
      <unit units="litre" exponent="-1"/>
      <unit units="second" exponent="-1"/>
      <unit units="volt" prefix="milli" exponent="-1"/>
   </units>
   <units name="millimole_to_nanomole">
      <unit units="mole" prefix="nano"/>
      <unit units="mole" prefix="milli" exponent="-1"/>
   </units>
   <units name="microcoulomb_to_picocoulomb">
      <unit units="coulomb" prefix="pico"/>
      <unit units="coulomb" prefix="micro" exponent="-1"/>
   </units>
   <component name="Time">
      <variable cmeta:id="time" name="time" units="time_units" public_interface="out" private_interface="out"/>
   </component>
   <component name="Environment">
      <variable name="T" units="Temperature_units" initial_value="310" private_interface="out"/>
      <variable name="T_exp" units="Temperature_units" initial_value="297" private_interface="out"/>
      <variable name="F" units="F_units" initial_value="96.4846" public_interface="out" private_interface="out"/>
      <variable name="R" units="R_units" initial_value="8.3144" private_interface="out"/>
      <variable name="Q10Ca" units="dimensionless" initial_value="2.1" private_interface="out"/>
      <variable name="Q10K" units="dimensionless" initial_value="1.5" private_interface="out"/>
      <variable name="Q10Na" units="dimensionless" initial_value="2.45" private_interface="out"/>
      <variable name="Ca_o" units="millimolar" initial_value="2.5" public_interface="out" private_interface="out"/>
      <variable name="Na_o" units="millimolar" initial_value="137" public_interface="out" private_interface="out"/>
      <variable name="K_o" units="millimolar" initial_value="7" public_interface="out" private_interface="out"/>
      <variable name="Cl_o" units="millimolar" initial_value="134" public_interface="out" private_interface="out"/>
      <variable name="T_correction_Na" units="dimensionless" public_interface="out" private_interface="out"/>
      <variable name="T_correction_K" units="dimensionless" public_interface="out" private_interface="out"/>
      <variable name="T_correction_Ca" units="dimensionless" public_interface="out" private_interface="out"/>
      <variable name="T_correction_BK" units="conductance_units" public_interface="out" private_interface="out"/>
      <variable name="FoRT" units="Inverse_Voltage_units" public_interface="out" private_interface="out"/>
      <variable name="RToF" units="voltage_units" public_interface="out" private_interface="out"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>FoRT</ci>
            <apply>
               <divide/>
               <ci>F</ci>
               <apply>
                  <times/>
                  <ci>R</ci>
                  <ci>T</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>RToF</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <ci>R</ci>
                  <ci>T</ci>
               </apply>
               <ci>F</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>T_correction_Ca</ci>
            <apply>
               <power/>
               <ci>Q10Ca</ci>
               <apply>
                  <divide/>
                  <apply>
                     <minus/>
                     <ci>T</ci>
                     <ci>T_exp</ci>
                  </apply>
                  <cn cellml:units="Temperature_units">10</cn>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>T_correction_K</ci>
            <apply>
               <power/>
               <ci>Q10K</ci>
               <apply>
                  <divide/>
                  <apply>
                     <minus/>
                     <ci>T</ci>
                     <ci>T_exp</ci>
                  </apply>
                  <cn cellml:units="Temperature_units">10</cn>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>T_correction_Na</ci>
            <apply>
               <power/>
               <ci>Q10Na</ci>
               <apply>
                  <divide/>
                  <apply>
                     <minus/>
                     <ci>T</ci>
                     <ci>T_exp</ci>
                  </apply>
                  <cn cellml:units="Temperature_units">10</cn>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>T_correction_BK</ci>
            <apply>
               <times/>
               <cn cellml:units="Inverse_Temperature_units_times_conductance">1.1</cn>
               <apply>
                  <minus/>
                  <ci>T</ci>
                  <ci>T_exp</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="ICC_Membrane">
      <variable name="Cm" units="capacitance_units" initial_value="0.025"/>
      <variable name="Vol" units="volume_units" initial_value="1e-12" public_interface="out" private_interface="out"/>
      <variable name="P_cyto" units="dimensionless" initial_value="0.7" public_interface="out" private_interface="out"/>
      <variable name="V_cyto" units="volume_units" private_interface="out"/>
      <variable name="fc" units="dimensionless" initial_value="0.01" private_interface="out"/>
      <variable cmeta:id="Vm" name="Vm" units="voltage_units" initial_value="-67" private_interface="out"/>
      <variable name="time" units="time_units" public_interface="in" private_interface="out"/>
      <variable name="Ca_i" units="millimolar" initial_value="0.00000993087" public_interface="out" private_interface="out"/>
      <variable name="Na_i" units="millimolar" initial_value="30" public_interface="out" private_interface="out"/>
      <variable name="K_i" units="millimolar" initial_value="120" private_interface="out"/>
      <variable name="Cl_i" units="millimolar" initial_value="88" private_interface="out"/>
      <variable name="Ca_o" units="millimolar" public_interface="in" private_interface="out"/>
      <variable name="Na_o" units="millimolar" public_interface="in" private_interface="out"/>
      <variable name="K_o" units="millimolar" public_interface="in" private_interface="out"/>
      <variable name="Cl_o" units="millimolar" public_interface="in" private_interface="out"/>
      <variable name="T_correction_Ca" units="dimensionless" public_interface="in" private_interface="out"/>
      <variable name="T_correction_Na" units="dimensionless" public_interface="in" private_interface="out"/>
      <variable name="T_correction_K" units="dimensionless" public_interface="in" private_interface="out"/>
      <variable name="T_correction_BK" units="conductance_units" public_interface="in" private_interface="out"/>
      <variable name="RToF" units="voltage_units" public_interface="in" private_interface="out"/>
      <variable name="FoRT" units="Inverse_Voltage_units" public_interface="in" private_interface="out"/>
      <variable name="F" units="F_units" public_interface="in" private_interface="out"/>
      <variable name="J_leak" units="millimolar_per_second" private_interface="in"/>
      <variable name="I_Na" units="current_units" private_interface="in"/>
      <variable name="I_Ltype" units="current_units" private_interface="in"/>
      <variable name="I_VDDR" units="current_units" private_interface="in"/>
      <variable name="I_kv11" units="current_units" private_interface="in"/>
      <variable name="I_BK" units="current_units" private_interface="in"/>
      <variable name="I_ERG" units="current_units" private_interface="in"/>
      <variable name="I_CaCl" units="current_units" private_interface="in"/>
      <variable name="I_NSCC" units="current_units" private_interface="in"/>
      <variable name="I_bk" units="current_units" private_interface="in"/>
      <variable name="J_PMCA" units="millimolar_per_second" private_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>V_cyto</ci>
            <apply>
               <times/>
               <ci>Vol</ci>
               <ci>P_cyto</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>Vm</ci>
            </apply>
            <apply>
               <times/>
               <apply>
                  <divide/>
                  <apply>
                     <times/>
                     <apply>
                        <minus/>
                        <cn cellml:units="dimensionless">1</cn>
                     </apply>
                     <cn cellml:units="dimensionless">1</cn>
                  </apply>
                  <ci>Cm</ci>
               </apply>
               <apply>
                  <plus/>
                  <ci>I_Na</ci>
                  <ci>I_Ltype</ci>
                  <ci>I_VDDR</ci>
                  <ci>I_kv11</ci>
                  <ci>I_ERG</ci>
                  <ci>I_BK</ci>
                  <ci>I_CaCl</ci>
                  <ci>I_NSCC</ci>
                  <ci>I_bk</ci>
                  <apply>
                     <times/>
                     <ci>J_PMCA</ci>
                     <cn cellml:units="dimensionless">2</cn>
                     <cn cellml:units="millimole_to_nanomole">1000000</cn>
                     <cn cellml:units="microcoulomb_to_picocoulomb">1000000</cn>
                     <ci>F</ci>
                     <ci>V_cyto</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>Ca_i</ci>
            </apply>
            <apply>
               <times/>
               <ci>fc</ci>
               <apply>
                  <plus/>
                  <apply>
                     <divide/>
                     <apply>
                        <plus/>
                        <apply>
                           <times/>
                           <apply>
                              <minus/>
                              <cn cellml:units="dimensionless">1</cn>
                           </apply>
                           <ci>I_Ltype</ci>
                        </apply>
                        <apply>
                           <times/>
                           <apply>
                              <minus/>
                              <cn cellml:units="dimensionless">1</cn>
                           </apply>
                           <ci>I_VDDR</ci>
                        </apply>
                     </apply>
                     <apply>
                        <times/>
                        <cn cellml:units="dimensionless">2</cn>
                        <cn cellml:units="millimole_to_nanomole">1000000</cn>
                        <cn cellml:units="microcoulomb_to_picocoulomb">1000000</cn>
                        <ci>F</ci>
                        <ci>V_cyto</ci>
                     </apply>
                  </apply>
                  <ci>J_leak</ci>
                  <apply>
                     <times/>
                     <apply>
                        <minus/>
                        <cn cellml:units="dimensionless">1</cn>
                     </apply>
                     <ci>J_PMCA</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="d_Ltype">
      <variable name="Vm" units="voltage_units" public_interface="in"/>
      <variable name="d_inf_Ltype" units="dimensionless"/>
      <variable name="tau_d_Ltype" units="time_units"/>
      <variable name="d_Ltype" units="dimensionless" initial_value="0" public_interface="out"/>
      <variable name="time" units="time_units" public_interface="in"/>
      <variable name="T_correction_Ca" units="dimensionless" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>d_inf_Ltype</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</ci>
                           <cn cellml:units="voltage_units">17</cn>
                        </apply>
                        <apply>
                           <minus/>
                           <cn cellml:units="voltage_units">4.3</cn>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>tau_d_Ltype</ci>
            <apply>
               <times/>
               <ci>T_correction_Ca</ci>
               <cn cellml:units="time_units">0.001</cn>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>d_Ltype</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>d_inf_Ltype</ci>
                  <ci>d_Ltype</ci>
               </apply>
               <ci>tau_d_Ltype</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="f_Ltype">
      <variable name="Vm" units="voltage_units" public_interface="in"/>
      <variable name="f_inf_Ltype" units="dimensionless"/>
      <variable name="tau_f_Ltype" units="time_units"/>
      <variable name="f_Ltype" units="dimensionless" initial_value="1" public_interface="out"/>
      <variable name="time" units="time_units" public_interface="in"/>
      <variable name="T_correction_Ca" units="dimensionless" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>f_inf_Ltype</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</ci>
                           <cn cellml:units="voltage_units">43</cn>
                        </apply>
                        <cn cellml:units="voltage_units">8.9</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>tau_f_Ltype</ci>
            <apply>
               <times/>
               <ci>T_correction_Ca</ci>
               <cn cellml:units="time_units">0.086</cn>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>f_Ltype</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>f_inf_Ltype</ci>
                  <ci>f_Ltype</ci>
               </apply>
               <ci>tau_f_Ltype</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="f_ca_Ltype">
      <variable name="Ca_i" units="millimolar" public_interface="in"/>
      <variable name="f_ca_inf_Ltype" units="dimensionless"/>
      <variable name="tau_f_ca_Ltype" units="time_units"/>
      <variable name="f_ca_Ltype" units="dimensionless" initial_value="1" public_interface="out"/>
      <variable name="time" units="time_units" public_interface="in"/>
      <variable name="T_correction_Ca" units="dimensionless" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>f_ca_inf_Ltype</ci>
            <apply>
               <minus/>
               <cn cellml:units="dimensionless">1</cn>
               <apply>
                  <divide/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <plus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <minus/>
                              <apply>
                                 <minus/>
                                 <ci>Ca_i</ci>
                                 <cn cellml:units="millimolar">0.0001</cn>
                              </apply>
                              <cn cellml:units="millimolar">0.000214</cn>
                           </apply>
                           <apply>
                              <minus/>
                              <cn cellml:units="millimolar">0.0000131</cn>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>tau_f_ca_Ltype</ci>
            <apply>
               <times/>
               <ci>T_correction_Ca</ci>
               <cn cellml:units="time_units">0.002</cn>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>f_ca_Ltype</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>f_ca_inf_Ltype</ci>
                  <ci>f_ca_Ltype</ci>
               </apply>
               <ci>tau_f_ca_Ltype</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="I_Ltype">
      <variable name="I_Ltype" units="current_units" public_interface="out"/>
      <variable name="E_Ca" units="voltage_units"/>
      <variable name="G_max_Ltype" units="conductance_units" initial_value="2"/>
      <variable name="Vm" units="voltage_units" public_interface="in" private_interface="out"/>
      <variable name="f_Ltype" units="dimensionless" private_interface="in"/>
      <variable name="d_Ltype" units="dimensionless" private_interface="in"/>
      <variable name="f_ca_Ltype" units="dimensionless" private_interface="in"/>
      <variable name="RToF" units="voltage_units" public_interface="in"/>
      <variable name="Ca_o" units="millimolar" public_interface="in"/>
      <variable name="Ca_i" units="millimolar" public_interface="in" private_interface="out"/>
      <variable name="T_correction_Ca" units="dimensionless" public_interface="in" private_interface="out"/>
      <variable name="time" units="time_units" public_interface="in" private_interface="out"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>E_Ca</ci>
            <apply>
               <times/>
               <cn cellml:units="dimensionless">0.5</cn>
               <ci>RToF</ci>
               <apply>
                  <ln/>
                  <apply>
                     <divide/>
                     <ci>Ca_o</ci>
                     <ci>Ca_i</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>I_Ltype</ci>
            <apply>
               <times/>
               <ci>G_max_Ltype</ci>
               <ci>f_Ltype</ci>
               <ci>d_Ltype</ci>
               <ci>f_ca_Ltype</ci>
               <apply>
                  <minus/>
                  <ci>Vm</ci>
                  <ci>E_Ca</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="J_PMCA">
      <variable name="J_PMCA" units="millimolar_per_second" public_interface="out"/>
      <variable name="J_max_PMCA" units="millimolar_per_second" initial_value="0.088464"/>
      <variable name="Ca_i" units="millimolar" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>J_PMCA</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <ci>J_max_PMCA</ci>
                  <cn cellml:units="dimensionless">1</cn>
               </apply>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <divide/>
                     <cn cellml:units="millimolar">0.000298</cn>
                     <ci>Ca_i</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="d_VDDR">
      <variable name="Vm" units="voltage_units" public_interface="in"/>
      <variable name="d_inf_VDDR" units="dimensionless"/>
      <variable name="tau_d_VDDR" units="time_units"/>
      <variable name="d_VDDR" units="dimensionless" initial_value="0" public_interface="out"/>
      <variable name="time" units="time_units" public_interface="in"/>
      <variable name="T_correction_Ca" units="dimensionless" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>d_inf_VDDR</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</ci>
                           <cn cellml:units="voltage_units">26</cn>
                        </apply>
                        <apply>
                           <minus/>
                           <cn cellml:units="voltage_units">6</cn>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>tau_d_VDDR</ci>
            <apply>
               <times/>
               <ci>T_correction_Ca</ci>
               <cn cellml:units="time_units">0.006</cn>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>d_VDDR</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>d_inf_VDDR</ci>
                  <ci>d_VDDR</ci>
               </apply>
               <ci>tau_d_VDDR</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="f_VDDR">
      <variable name="Vm" units="voltage_units" public_interface="in"/>
      <variable name="f_inf_VDDR" units="dimensionless"/>
      <variable name="tau_f_VDDR" units="time_units"/>
      <variable name="f_VDDR" units="dimensionless" initial_value="1" public_interface="out"/>
      <variable name="time" units="time_units" public_interface="in"/>
      <variable name="T_correction_Ca" units="dimensionless" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>f_inf_VDDR</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</ci>
                           <cn cellml:units="voltage_units">66</cn>
                        </apply>
                        <cn cellml:units="voltage_units">6</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>tau_f_VDDR</ci>
            <apply>
               <times/>
               <ci>T_correction_Ca</ci>
               <cn cellml:units="time_units">0.04</cn>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>f_VDDR</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>f_inf_VDDR</ci>
                  <ci>f_VDDR</ci>
               </apply>
               <ci>tau_f_VDDR</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="I_VDDR">
      <variable name="I_VDDR" units="current_units" public_interface="out"/>
      <variable name="E_Ca" units="voltage_units"/>
      <variable name="G_max_VDDR" units="conductance_units" initial_value="3"/>
      <variable name="Vm" units="voltage_units" public_interface="in" private_interface="out"/>
      <variable name="f_VDDR" units="dimensionless" private_interface="in"/>
      <variable name="d_VDDR" units="dimensionless" private_interface="in"/>
      <variable name="RToF" units="voltage_units" public_interface="in"/>
      <variable name="Ca_o" units="millimolar" public_interface="in"/>
      <variable name="Ca_i" units="millimolar" public_interface="in"/>
      <variable name="T_correction_Ca" units="dimensionless" public_interface="in" private_interface="out"/>
      <variable name="time" units="time_units" public_interface="in" private_interface="out"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>E_Ca</ci>
            <apply>
               <times/>
               <cn cellml:units="dimensionless">0.5</cn>
               <ci>RToF</ci>
               <apply>
                  <ln/>
                  <apply>
                     <divide/>
                     <ci>Ca_o</ci>
                     <ci>Ca_i</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>I_VDDR</ci>
            <apply>
               <times/>
               <ci>G_max_VDDR</ci>
               <ci>f_VDDR</ci>
               <ci>d_VDDR</ci>
               <apply>
                  <minus/>
                  <ci>Vm</ci>
                  <ci>E_Ca</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="d_CaCl">
      <variable name="Ca_i" units="millimolar" public_interface="in"/>
      <variable name="d_inf_CaCl" units="dimensionless"/>
      <variable name="tau_d_CaCl" units="time_units" initial_value="0.03"/>
      <variable name="d_CaCl" units="dimensionless" initial_value="0" public_interface="out"/>
      <variable name="time" units="time_units" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>d_inf_CaCl</ci>
            <apply>
               <divide/>
               <cn cellml:units="dimensionless">1</cn>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <power/>
                     <apply>
                        <divide/>
                        <cn cellml:units="millimolar">0.00014</cn>
                        <ci>Ca_i</ci>
                     </apply>
                     <cn cellml:units="dimensionless">3</cn>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>d_CaCl</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>d_inf_CaCl</ci>
                  <ci>d_CaCl</ci>
               </apply>
               <ci>tau_d_CaCl</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="I_CaCl">
      <variable name="I_CaCl" units="current_units" public_interface="out"/>
      <variable name="Ca_i" units="millimolar" public_interface="in" private_interface="out"/>
      <variable name="E_Cl" units="voltage_units"/>
      <variable name="G_max_CaCl" units="conductance_units" initial_value="10.1"/>
      <variable name="Vm" units="voltage_units" public_interface="in"/>
      <variable name="d_CaCl" units="dimensionless" private_interface="in"/>
      <variable name="RToF" units="voltage_units" public_interface="in"/>
      <variable name="Cl_o" units="millimolar" public_interface="in"/>
      <variable name="Cl_i" units="millimolar" public_interface="in"/>
      <variable name="time" units="time_units" public_interface="in" private_interface="out"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>E_Cl</ci>
            <apply>
               <times/>
               <ci>RToF</ci>
               <apply>
                  <ln/>
                  <apply>
                     <divide/>
                     <ci>Cl_i</ci>
                     <ci>Cl_o</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>I_CaCl</ci>
            <apply>
               <times/>
               <ci>G_max_CaCl</ci>
               <ci>d_CaCl</ci>
               <apply>
                  <minus/>
                  <ci>Vm</ci>
                  <ci>E_Cl</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="d_BK">
      <variable name="Vm" units="voltage_units" public_interface="in"/>
      <variable name="Ca_i" units="millimolar" public_interface="in"/>
      <variable name="d_BK" units="dimensionless" public_interface="out"/>
      <variable name="time" units="time_units" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>d_BK</ci>
            <apply>
               <divide/>
               <cn cellml:units="dimensionless">1</cn>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <minus/>
                        <apply>
                           <divide/>
                           <ci>Vm</ci>
                           <apply>
                              <minus/>
                              <cn cellml:units="voltage_units">17</cn>
                           </apply>
                        </apply>
                        <apply>
                           <times/>
                           <cn cellml:units="dimensionless">2</cn>
                           <apply>
                              <ln/>
                              <apply>
                                 <divide/>
                                 <ci>Ca_i</ci>
                                 <cn cellml:units="millimolar">0.001</cn>
                              </apply>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="I_BK">
      <variable name="I_BK" units="current_units" public_interface="out"/>
      <variable name="Ca_i" units="millimolar" public_interface="in" private_interface="out"/>
      <variable name="E_K" units="voltage_units"/>
      <variable name="G_max_BK" units="conductance_units" initial_value="23"/>
      <variable name="Vm" units="voltage_units" public_interface="in" private_interface="out"/>
      <variable name="d_BK" units="dimensionless" private_interface="in"/>
      <variable name="RToF" units="voltage_units" public_interface="in"/>
      <variable name="K_o" units="millimolar" public_interface="in"/>
      <variable name="K_i" units="millimolar" public_interface="in"/>
      <variable name="T_correction_BK" units="conductance_units" public_interface="in"/>
      <variable name="time" units="time_units" public_interface="in" private_interface="out"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>E_K</ci>
            <apply>
               <times/>
               <ci>RToF</ci>
               <apply>
                  <ln/>
                  <apply>
                     <divide/>
                     <ci>K_o</ci>
                     <ci>K_i</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>I_BK</ci>
            <apply>
               <times/>
               <apply>
                  <plus/>
                  <ci>G_max_BK</ci>
                  <ci>T_correction_BK</ci>
               </apply>
               <ci>d_BK</ci>
               <apply>
                  <minus/>
                  <ci>Vm</ci>
                  <ci>E_K</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="I_bk">
      <variable name="I_bk" units="current_units" public_interface="out"/>
      <variable name="E_K" units="voltage_units"/>
      <variable name="G_max_bk" units="conductance_units" initial_value="0.15"/>
      <variable name="Vm" units="voltage_units" public_interface="in"/>
      <variable name="RToF" units="voltage_units" public_interface="in"/>
      <variable name="K_o" units="millimolar" public_interface="in"/>
      <variable name="K_i" units="millimolar" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>E_K</ci>
            <apply>
               <times/>
               <ci>RToF</ci>
               <apply>
                  <ln/>
                  <apply>
                     <divide/>
                     <ci>K_o</ci>
                     <ci>K_i</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>I_bk</ci>
            <apply>
               <times/>
               <ci>G_max_bk</ci>
               <apply>
                  <minus/>
                  <ci>Vm</ci>
                  <ci>E_K</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="d_ERG">
      <variable name="Vm" units="voltage_units" public_interface="in"/>
      <variable name="d_inf_ERG" units="dimensionless"/>
      <variable name="tau_d_ERG" units="time_units"/>
      <variable name="d_ERG" units="dimensionless" initial_value="0" public_interface="out"/>
      <variable name="time" units="time_units" public_interface="in"/>
      <variable name="T_correction_K" units="dimensionless" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>d_inf_ERG</ci>
            <apply>
               <plus/>
               <cn cellml:units="dimensionless">0.2</cn>
               <apply>
                  <divide/>
                  <cn cellml:units="dimensionless">0.8</cn>
                  <apply>
                     <plus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <plus/>
                              <ci>Vm</ci>
                              <cn cellml:units="voltage_units">20</cn>
                           </apply>
                           <apply>
                              <minus/>
                              <cn cellml:units="voltage_units">1.8</cn>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>tau_d_ERG</ci>
            <apply>
               <times/>
               <ci>T_correction_K</ci>
               <cn cellml:units="time_units">0.003</cn>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>d_ERG</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>d_inf_ERG</ci>
                  <ci>d_ERG</ci>
               </apply>
               <ci>tau_d_ERG</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="I_ERG">
      <variable name="I_ERG" units="current_units" public_interface="out"/>
      <variable name="E_K" units="voltage_units"/>
      <variable name="G_max_ERG" units="conductance_units" initial_value="2.5"/>
      <variable name="Vm" units="voltage_units" public_interface="in" private_interface="out"/>
      <variable name="d_ERG" units="dimensionless" private_interface="in"/>
      <variable name="RToF" units="voltage_units" public_interface="in"/>
      <variable name="K_o" units="millimolar" public_interface="in"/>
      <variable name="K_i" units="millimolar" public_interface="in"/>
      <variable name="T_correction_K" units="dimensionless" public_interface="in" private_interface="out"/>
      <variable name="time" units="time_units" public_interface="in" private_interface="out"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>E_K</ci>
            <apply>
               <times/>
               <ci>RToF</ci>
               <apply>
                  <ln/>
                  <apply>
                     <divide/>
                     <ci>K_o</ci>
                     <ci>K_i</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>I_ERG</ci>
            <apply>
               <times/>
               <ci>G_max_ERG</ci>
               <ci>d_ERG</ci>
               <apply>
                  <minus/>
                  <ci>Vm</ci>
                  <ci>E_K</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="d_kv11">
      <variable name="Vm" units="voltage_units" public_interface="in"/>
      <variable name="d_inf_kv11" units="dimensionless"/>
      <variable name="tau_d_kv11" units="time_units"/>
      <variable name="d_kv11" units="dimensionless" initial_value="0" public_interface="out"/>
      <variable name="time" units="time_units" public_interface="in"/>
      <variable name="T_correction_K" units="dimensionless" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>d_inf_kv11</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</ci>
                           <cn cellml:units="voltage_units">25</cn>
                        </apply>
                        <apply>
                           <minus/>
                           <cn cellml:units="voltage_units">7.7</cn>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>tau_d_kv11</ci>
            <apply>
               <times/>
               <ci>T_correction_K</ci>
               <cn cellml:units="time_units">0.005</cn>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>d_kv11</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>d_inf_kv11</ci>
                  <ci>d_kv11</ci>
               </apply>
               <ci>tau_d_kv11</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="f_kv11">
      <variable name="Vm" units="voltage_units" public_interface="in"/>
      <variable name="f_inf_kv11" units="dimensionless"/>
      <variable name="tau_f_kv11" units="time_units"/>
      <variable name="f_kv11" units="dimensionless" initial_value="1" public_interface="out"/>
      <variable name="time" units="time_units" public_interface="in"/>
      <variable name="T_correction_K" units="dimensionless" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>f_inf_kv11</ci>
            <apply>
               <plus/>
               <cn cellml:units="dimensionless">0.5</cn>
               <apply>
                  <divide/>
                  <cn cellml:units="dimensionless">0.5</cn>
                  <apply>
                     <plus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <plus/>
                              <ci>Vm</ci>
                              <cn cellml:units="voltage_units">44.8</cn>
                           </apply>
                           <cn cellml:units="voltage_units">4.4</cn>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>tau_f_kv11</ci>
            <apply>
               <times/>
               <ci>T_correction_K</ci>
               <cn cellml:units="time_units">0.005</cn>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>f_kv11</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>f_inf_kv11</ci>
                  <ci>f_kv11</ci>
               </apply>
               <ci>tau_f_kv11</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="I_kv11">
      <variable name="I_kv11" units="current_units" public_interface="out"/>
      <variable name="E_K" units="voltage_units"/>
      <variable name="G_max_kv11" units="conductance_units" initial_value="6.3"/>
      <variable name="Vm" units="voltage_units" public_interface="in" private_interface="out"/>
      <variable name="f_kv11" units="dimensionless" private_interface="in"/>
      <variable name="d_kv11" units="dimensionless" private_interface="in"/>
      <variable name="RToF" units="voltage_units" public_interface="in"/>
      <variable name="K_o" units="millimolar" public_interface="in"/>
      <variable name="K_i" units="millimolar" public_interface="in"/>
      <variable name="T_correction_K" units="dimensionless" public_interface="in" private_interface="out"/>
      <variable name="time" units="time_units" public_interface="in" private_interface="out"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>E_K</ci>
            <apply>
               <times/>
               <ci>RToF</ci>
               <apply>
                  <ln/>
                  <apply>
                     <divide/>
                     <ci>K_o</ci>
                     <ci>K_i</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>I_kv11</ci>
            <apply>
               <times/>
               <ci>G_max_kv11</ci>
               <ci>f_kv11</ci>
               <ci>d_kv11</ci>
               <apply>
                  <minus/>
                  <ci>Vm</ci>
                  <ci>E_K</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="d_Na">
      <variable name="Vm" units="voltage_units" public_interface="in"/>
      <variable name="d_inf_Na" units="dimensionless"/>
      <variable name="tau_d_Na" units="time_units"/>
      <variable name="d_Na" units="dimensionless" initial_value="0" public_interface="out"/>
      <variable name="time" units="time_units" public_interface="in"/>
      <variable name="T_correction_Na" units="dimensionless" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>d_inf_Na</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</ci>
                           <cn cellml:units="voltage_units">47</cn>
                        </apply>
                        <apply>
                           <minus/>
                           <cn cellml:units="voltage_units">4.8</cn>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>tau_d_Na</ci>
            <apply>
               <times/>
               <ci>T_correction_Na</ci>
               <cn cellml:units="time_units">0.003</cn>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>d_Na</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>d_inf_Na</ci>
                  <ci>d_Na</ci>
               </apply>
               <ci>tau_d_Na</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="f_Na">
      <variable name="Vm" units="voltage_units" public_interface="in"/>
      <variable name="f_inf_Na" units="dimensionless"/>
      <variable name="tau_f_Na" units="time_units"/>
      <variable name="f_Na" units="dimensionless" initial_value="1" public_interface="out"/>
      <variable name="time" units="time_units" public_interface="in"/>
      <variable name="T_correction_Na" units="dimensionless" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>f_inf_Na</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</ci>
                           <cn cellml:units="voltage_units">78</cn>
                        </apply>
                        <cn cellml:units="voltage_units">7</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>tau_f_Na</ci>
            <apply>
               <times/>
               <ci>T_correction_Na</ci>
               <cn cellml:units="time_units">0.0016</cn>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>f_Na</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>f_inf_Na</ci>
                  <ci>f_Na</ci>
               </apply>
               <ci>tau_f_Na</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="I_Na">
      <variable name="I_Na" units="current_units" public_interface="out"/>
      <variable name="E_Na" units="voltage_units"/>
      <variable name="G_max_Na" units="conductance_units" initial_value="20"/>
      <variable name="Vm" units="voltage_units" public_interface="in" private_interface="out"/>
      <variable name="f_Na" units="dimensionless" private_interface="in"/>
      <variable name="d_Na" units="dimensionless" private_interface="in"/>
      <variable name="RToF" units="voltage_units" public_interface="in"/>
      <variable name="Na_o" units="millimolar" public_interface="in"/>
      <variable name="Na_i" units="millimolar" public_interface="in"/>
      <variable name="T_correction_Na" units="dimensionless" public_interface="in" private_interface="out"/>
      <variable name="time" units="time_units" public_interface="in" private_interface="out"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>E_Na</ci>
            <apply>
               <times/>
               <ci>RToF</ci>
               <apply>
                  <ln/>
                  <apply>
                     <divide/>
                     <ci>Na_o</ci>
                     <ci>Na_i</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>I_Na</ci>
            <apply>
               <times/>
               <ci>G_max_Na</ci>
               <ci>f_Na</ci>
               <ci>d_Na</ci>
               <apply>
                  <minus/>
                  <ci>Vm</ci>
                  <ci>E_Na</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="d_NSCC">
      <variable name="Ca_PU" units="millimolar" public_interface="in"/>
      <variable name="d_inf_NSCC" units="dimensionless"/>
      <variable name="tau_d_NSCC" units="time_units" initial_value="0.35"/>
      <variable name="d_NSCC" units="dimensionless" initial_value="0" public_interface="out"/>
      <variable name="time" units="time_units" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>d_inf_NSCC</ci>
            <apply>
               <divide/>
               <cn cellml:units="dimensionless">1</cn>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <power/>
                     <apply>
                        <divide/>
                        <cn cellml:units="millimolar">0.0000745</cn>
                        <ci>Ca_PU</ci>
                     </apply>
                     <apply>
                        <minus/>
                        <cn cellml:units="dimensionless">85</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>d_NSCC</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>d_inf_NSCC</ci>
                  <ci>d_NSCC</ci>
               </apply>
               <ci>tau_d_NSCC</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="I_NSCC">
      <variable name="I_NSCC" units="current_units" public_interface="out"/>
      <variable name="E_NSCC" units="voltage_units"/>
      <variable name="G_max_NSCC" units="conductance_units" initial_value="12.15"/>
      <variable name="Vm" units="voltage_units" public_interface="in"/>
      <variable name="d_NSCC" units="dimensionless" private_interface="in"/>
      <variable name="RToF" units="voltage_units" public_interface="in"/>
      <variable name="Na_o" units="millimolar" public_interface="in"/>
      <variable name="Na_i" units="millimolar" public_interface="in"/>
      <variable name="K_o" units="millimolar" public_interface="in"/>
      <variable name="K_i" units="millimolar" public_interface="in"/>
      <variable name="NaPerm_o_Kperm" units="dimensionless" initial_value="1.056075"/>
      <variable name="Ca_PU" units="millimolar" public_interface="in" private_interface="out"/>
      <variable name="time" units="time_units" public_interface="in" private_interface="out"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>E_NSCC</ci>
            <apply>
               <times/>
               <ci>RToF</ci>
               <apply>
                  <ln/>
                  <apply>
                     <divide/>
                     <apply>
                        <plus/>
                        <ci>K_o</ci>
                        <apply>
                           <times/>
                           <ci>Na_o</ci>
                           <ci>NaPerm_o_Kperm</ci>
                        </apply>
                     </apply>
                     <apply>
                        <plus/>
                        <ci>K_i</ci>
                        <apply>
                           <times/>
                           <ci>Na_i</ci>
                           <ci>NaPerm_o_Kperm</ci>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>I_NSCC</ci>
            <apply>
               <times/>
               <ci>G_max_NSCC</ci>
               <ci>d_NSCC</ci>
               <apply>
                  <minus/>
                  <ci>Vm</ci>
                  <ci>E_NSCC</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="PU_unit">
      <variable name="Vol" units="volume_units" public_interface="in"/>
      <variable name="P_PU" units="dimensionless" initial_value="0.001"/>
      <variable name="P_mito" units="dimensionless" initial_value="0.12871"/>
      <variable name="P_ER" units="dimensionless" initial_value="0.1"/>
      <variable name="V_MITO" units="volume_units"/>
      <variable name="V_ER" units="volume_units"/>
      <variable name="V_PU" units="volume_units"/>
      <variable name="fe" units="dimensionless" initial_value="0.01"/>
      <variable name="fm" units="dimensionless" initial_value="0.0003"/>
      <variable name="time" units="time_units" public_interface="in"/>
      <variable name="fc" units="dimensionless" public_interface="in"/>
      <variable name="Na_i" units="millimolar" public_interface="in"/>
      <variable name="Ca_i" units="millimolar" public_interface="in"/>
      <variable name="V_cyto" units="volume_units" public_interface="in"/>
      <variable name="FoRT" units="Inverse_Voltage_units" public_interface="in"/>
      <variable name="RToF" units="voltage_units" public_interface="in"/>
      <variable name="Ca_PU" units="millimolar" initial_value="0.0000902" public_interface="out" private_interface="out"/>
      <variable name="F" units="F_units" public_interface="in"/>
      <variable name="Ca_m" units="millimolar" initial_value="0.000136"/>
      <variable name="Ca_ER" units="millimolar" initial_value="0.007299"/>
      <variable name="ADP_m" units="millimolar" initial_value="2.60093454"/>
      <variable name="ADP_i" units="millimolar" initial_value="0.0077282"/>
      <variable name="NADH_m" units="millimolar" initial_value="0.101476"/>
      <variable name="h" units="dimensionless" initial_value="0.9397"/>
      <variable name="IP3" units="millimolar" initial_value="0.0006"/>
      <variable name="deltaPsi" units="voltage_units" initial_value="164.000044"/>
      <variable name="deltapH" units="dimensionless" initial_value="-0.4"/>
      <variable name="Cmito" units="millifarads" initial_value="0.006995"/>
      <variable name="K_res" units="dimensionless" initial_value="1.35e18"/>
      <variable name="r1" units="dimensionless" initial_value="2.077e-18"/>
      <variable name="r2" units="dimensionless" initial_value="1.728e-9"/>
      <variable name="r3" units="dimensionless" initial_value="1.059e-26"/>
      <variable name="ra" units="rate_constants_units" initial_value="6.394e-10"/>
      <variable name="rb" units="rate_constants_units" initial_value="1.762e-13"/>
      <variable name="rc1" units="rate_constants_units" initial_value="2.656e-19"/>
      <variable name="rc2" units="rate_constants_units" initial_value="8.632e-27"/>
      <variable name="deltaPsi_B" units="voltage_units" initial_value="50"/>
      <variable name="g" units="dimensionless" initial_value="0.85"/>
      <variable name="K_F1" units="millimolar" initial_value="1.71e9"/>
      <variable name="Pi_m" units="millimolar" initial_value="20"/>
      <variable name="p1" units="dimensionless" initial_value="1.346e-8"/>
      <variable name="p2" units="dimensionless" initial_value="7.739e-7"/>
      <variable name="p3" units="dimensionless" initial_value="6.65e-15"/>
      <variable name="pa" units="rate_constants_units" initial_value="1.656e-5"/>
      <variable name="pb" units="rate_constants_units" initial_value="3.373e-7"/>
      <variable name="pc1" units="rate_constants_units" initial_value="9.651e-14"/>
      <variable name="pc2" units="rate_constants_units" initial_value="4.845e-19"/>
      <variable name="frac" units="dimensionless" initial_value="0.5"/>
      <variable name="K_act" units="millimolar" initial_value="0.00038"/>
      <variable name="na" units="dimensionless" initial_value="2.8"/>
      <variable name="deltaPsi_star" units="voltage_units" initial_value="91"/>
      <variable name="K_Na" units="millimolar" initial_value="9.4"/>
      <variable name="K_Ca" units="millimolar" initial_value="0.003"/>
      <variable name="K_trans" units="millimolar" initial_value="0.006"/>
      <variable name="L" units="dimensionless" initial_value="50"/>
      <variable name="b" units="dimensionless" initial_value="0.5"/>
      <variable name="beta_max" units="rate_constants_units" initial_value="2.055"/>
      <variable name="beta1" units="per_millimolar" initial_value="1.66"/>
      <variable name="beta2" units="per_millimolar" initial_value="0.0249"/>
      <variable name="beta3" units="per_millimolar" initial_value="4"/>
      <variable name="beta4" units="per_millimolar" initial_value="2.83"/>
      <variable name="beta5" units="per_millimolar" initial_value="1.3"/>
      <variable name="beta6" units="per_millimolar" initial_value="2.66"/>
      <variable name="beta7" units="per_millimolar" initial_value="0.16"/>
      <variable name="KCa_PDH" units="millimolar" initial_value="0.00005"/>
      <variable name="u1" units="dimensionless" initial_value="15"/>
      <variable name="u2" units="dimensionless" initial_value="1.1"/>
      <variable name="n" units="dimensionless" initial_value="2"/>
      <variable name="K_Glc" units="millimolar" initial_value="8.7"/>
      <variable name="nhyd" units="dimensionless" initial_value="2.7"/>
      <variable name="K_hyd" units="rate_constants_units" initial_value="0.05125"/>
      <variable name="J_ERleak" units="rate_constants_units" initial_value="1.666667"/>
      <variable name="Jmax_IP3" units="rate_constants_units" initial_value="0"/>
      <variable name="d_IP3" units="millimolar" initial_value="0.00025"/>
      <variable name="d_ACT" units="millimolar" initial_value="0.001"/>
      <variable name="d_INH" units="millimolar" initial_value="0.0014"/>
      <variable name="tauh" units="time_units" initial_value="4"/>
      <variable name="Jmax_serca" units="millimolar_per_second" initial_value="1.8333"/>
      <variable name="k_serca" units="millimolar" initial_value="0.00042"/>
      <variable name="conc" units="millimolar" initial_value="0.001"/>
      <variable name="Jmax_uni" units="rate_constants_units" initial_value="5000"/>
      <variable name="Jmax_NaCa" units="millimolar_per_second" initial_value="0.05"/>
      <variable name="J_max_leak" units="rate_constants_units" initial_value="0.01"/>
      <variable name="rho_res" units="millimolar" initial_value="0.4"/>
      <variable name="rho_F1" units="millimolar" initial_value="0.7"/>
      <variable name="g_H" units="millimolar_per_second_per_millivolt" initial_value="0.0033333"/>
      <variable name="J_red_basal" units="millimolar_per_second" initial_value="0.3333"/>
      <variable name="Jmax_ANT" units="millimolar_per_second" initial_value="15"/>
      <variable name="J_hyd_max" units="millimolar_per_second" initial_value="0.037625"/>
      <variable name="Glc" units="millimolar" initial_value="1"/>
      <variable name="total_NAD_m" units="millimolar" initial_value="8"/>
      <variable name="total_ANP_m" units="millimolar" initial_value="12"/>
      <variable name="total_ANP_i" units="millimolar" initial_value="2"/>
      <variable name="NAD_m" units="millimolar"/>
      <variable name="ATP_m" units="millimolar"/>
      <variable name="ADP_mfree" units="millimolar"/>
      <variable name="ADP3_m" units="millimolar"/>
      <variable name="ATP4_m" units="millimolar"/>
      <variable name="ATP_i" units="millimolar"/>
      <variable name="ADP_ifree" units="millimolar"/>
      <variable name="ADP3_i" units="millimolar"/>
      <variable name="MgADP_i" units="millimolar"/>
      <variable name="ATP4_i" units="millimolar"/>
      <variable name="J_ERout" units="millimolar_per_second"/>
      <variable name="J_SERCA" units="millimolar_per_second"/>
      <variable name="MWC" units="millimolar"/>
      <variable name="J_uni" units="millimolar_per_second"/>
      <variable name="J_NaCa" units="millimolar_per_second"/>
      <variable name="J_leak" units="millimolar_per_second" public_interface="out"/>
      <variable name="J_red" units="millimolar_per_second"/>
      <variable name="J_pTCA" units="millimolar_per_second"/>
      <variable name="A_F1" units="voltage_units"/>
      <variable name="J_pF1" units="millimolar_per_second"/>
      <variable name="J_HF1" units="millimolar_per_second"/>
      <variable name="A_res" units="voltage_units"/>
      <variable name="J_o" units="millimolar_per_second"/>
      <variable name="J_glyTotal" units="millimolar_per_second"/>
      <variable name="f_PDHa" units="dimensionless"/>
      <variable name="J_Hres" units="millimolar_per_second"/>
      <variable name="J_ANT" units="millimolar_per_second"/>
      <variable name="PMF" units="voltage_units"/>
      <variable name="J_Hleak" units="millimolar_per_second"/>
      <variable name="J_pGly" units="millimolar_per_second"/>
      <variable name="J_hydSS" units="millimolar_per_second"/>
      <variable name="J_hyd" units="millimolar_per_second"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>V_MITO</ci>
            <apply>
               <times/>
               <ci>Vol</ci>
               <ci>P_mito</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>V_PU</ci>
            <apply>
               <times/>
               <ci>Vol</ci>
               <ci>P_PU</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>V_ER</ci>
            <apply>
               <times/>
               <ci>Vol</ci>
               <ci>P_ER</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>J_ERout</ci>
            <apply>
               <times/>
               <apply>
                  <plus/>
                  <apply>
                     <times/>
                     <ci>Jmax_IP3</ci>
                     <apply>
                        <power/>
                        <apply>
                           <divide/>
                           <ci>IP3</ci>
                           <apply>
                              <plus/>
                              <ci>IP3</ci>
                              <ci>d_IP3</ci>
                           </apply>
                        </apply>
                        <cn cellml:units="dimensionless">3</cn>
                     </apply>
                     <apply>
                        <power/>
                        <apply>
                           <divide/>
                           <ci>Ca_PU</ci>
                           <apply>
                              <plus/>
                              <ci>Ca_PU</ci>
                              <ci>d_ACT</ci>
                           </apply>
                        </apply>
                        <cn cellml:units="dimensionless">3</cn>
                     </apply>
                     <apply>
                        <power/>
                        <ci>h</ci>
                        <cn cellml:units="dimensionless">3</cn>
                     </apply>
                  </apply>
                  <ci>J_ERleak</ci>
               </apply>
               <apply>
                  <minus/>
                  <ci>Ca_ER</ci>
                  <ci>Ca_PU</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>J_SERCA</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <ci>Jmax_serca</ci>
                  <apply>
                     <power/>
                     <ci>Ca_PU</ci>
                     <cn cellml:units="dimensionless">2</cn>
                  </apply>
               </apply>
               <apply>
                  <plus/>
                  <apply>
                     <power/>
                     <ci>k_serca</ci>
                     <cn cellml:units="dimensionless">2</cn>
                  </apply>
                  <apply>
                     <power/>
                     <ci>Ca_PU</ci>
                     <cn cellml:units="dimensionless">2</cn>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>MWC</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <ci>conc</ci>
                        <ci>Ca_PU</ci>
                     </apply>
                     <ci>K_trans</ci>
                  </apply>
                  <apply>
                     <power/>
                     <apply>
                        <plus/>
                        <cn cellml:units="dimensionless">1</cn>
                        <apply>
                           <divide/>
                           <ci>Ca_PU</ci>
                           <ci>K_trans</ci>
                        </apply>
                     </apply>
                     <cn cellml:units="dimensionless">3</cn>
                  </apply>
               </apply>
               <apply>
                  <plus/>
                  <apply>
                     <power/>
                     <apply>
                        <plus/>
                        <cn cellml:units="dimensionless">1</cn>
                        <apply>
                           <divide/>
                           <ci>Ca_PU</ci>
                           <ci>K_trans</ci>
                        </apply>
                     </apply>
                     <cn cellml:units="dimensionless">4</cn>
                  </apply>
                  <apply>
                     <divide/>
                     <ci>L</ci>
                     <apply>
                        <power/>
                        <apply>
                           <plus/>
                           <cn cellml:units="dimensionless">1</cn>
                           <apply>
                              <divide/>
                              <ci>Ca_PU</ci>
                              <ci>K_act</ci>
                           </apply>
                        </apply>
                        <ci>na</ci>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>J_uni</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <ci>Jmax_uni</ci>
                  <apply>
                     <minus/>
                     <ci>MWC</ci>
                     <apply>
                        <times/>
                        <ci>Ca_m</ci>
                        <apply>
                           <exp/>
                           <apply>
                              <times/>
                              <apply>
                                 <minus/>
                                 <cn cellml:units="dimensionless">2</cn>
                              </apply>
                              <ci>FoRT</ci>
                              <apply>
                                 <minus/>
                                 <ci>deltaPsi</ci>
                                 <ci>deltaPsi_star</ci>
                              </apply>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
                  <cn cellml:units="dimensionless">2</cn>
                  <ci>FoRT</ci>
                  <apply>
                     <minus/>
                     <ci>deltaPsi</ci>
                     <ci>deltaPsi_star</ci>
                  </apply>
               </apply>
               <apply>
                  <minus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <times/>
                        <apply>
                           <minus/>
                           <cn cellml:units="dimensionless">2</cn>
                        </apply>
                        <ci>FoRT</ci>
                        <apply>
                           <minus/>
                           <ci>deltaPsi</ci>
                           <ci>deltaPsi_star</ci>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>J_NaCa</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <ci>Jmax_NaCa</ci>
                  <apply>
                     <exp/>
                     <apply>
                        <times/>
                        <ci>b</ci>
                        <ci>FoRT</ci>
                        <apply>
                           <minus/>
                           <ci>deltaPsi</ci>
                           <ci>deltaPsi_star</ci>
                        </apply>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <apply>
                     <plus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <apply>
                        <power/>
                        <apply>
                           <divide/>
                           <ci>K_Na</ci>
                           <ci>Na_i</ci>
                        </apply>
                        <ci>n</ci>
                     </apply>
                  </apply>
                  <apply>
                     <plus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <apply>
                        <divide/>
                        <ci>K_Ca</ci>
                        <ci>Ca_m</ci>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>J_leak</ci>
            <apply>
               <times/>
               <ci>J_max_leak</ci>
               <apply>
                  <minus/>
                  <ci>Ca_PU</ci>
                  <ci>Ca_i</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>A_res</ci>
            <apply>
               <times/>
               <ci>RToF</ci>
               <apply>
                  <ln/>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <ci>K_res</ci>
                        <apply>
                           <root/>
                           <ci>NADH_m</ci>
                        </apply>
                     </apply>
                     <apply>
                        <root/>
                        <ci>NAD_m</ci>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>J_o</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <ci>rho_res</ci>
                  <cn cellml:units="dimensionless">0.5</cn>
                  <apply>
                     <plus/>
                     <apply>
                        <times/>
                        <apply>
                           <plus/>
                           <apply>
                              <times/>
                              <ci>ra</ci>
                              <apply>
                                 <power/>
                                 <cn cellml:units="dimensionless">10</cn>
                                 <apply>
                                    <times/>
                                    <cn cellml:units="dimensionless">6</cn>
                                    <ci>deltapH</ci>
                                 </apply>
                              </apply>
                           </apply>
                           <apply>
                              <times/>
                              <ci>rc1</ci>
                              <apply>
                                 <exp/>
                                 <apply>
                                    <times/>
                                    <cn cellml:units="dimensionless">6</cn>
                                    <ci>deltaPsi_B</ci>
                                    <ci>FoRT</ci>
                                 </apply>
                              </apply>
                           </apply>
                        </apply>
                        <apply>
                           <exp/>
                           <apply>
                              <times/>
                              <ci>A_res</ci>
                              <ci>FoRT</ci>
                           </apply>
                        </apply>
                     </apply>
                     <apply>
                        <times/>
                        <apply>
                           <minus/>
                           <cn cellml:units="dimensionless">1</cn>
                        </apply>
                        <ci>ra</ci>
                        <apply>
                           <exp/>
                           <apply>
                              <times/>
                              <ci>g</ci>
                              <cn cellml:units="dimensionless">6</cn>
                              <ci>FoRT</ci>
                              <ci>deltaPsi</ci>
                           </apply>
                        </apply>
                     </apply>
                     <apply>
                        <times/>
                        <ci>rc2</ci>
                        <apply>
                           <exp/>
                           <apply>
                              <times/>
                              <ci>FoRT</ci>
                              <ci>A_res</ci>
                           </apply>
                        </apply>
                        <apply>
                           <exp/>
                           <apply>
                              <times/>
                              <ci>FoRT</ci>
                              <ci>deltaPsi</ci>
                              <cn cellml:units="dimensionless">6</cn>
                              <ci>g</ci>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <plus/>
                  <apply>
                     <times/>
                     <apply>
                        <plus/>
                        <cn cellml:units="dimensionless">1</cn>
                        <apply>
                           <times/>
                           <ci>r1</ci>
                           <apply>
                              <exp/>
                              <apply>
                                 <times/>
                                 <ci>FoRT</ci>
                                 <ci>A_res</ci>
                              </apply>
                           </apply>
                        </apply>
                     </apply>
                     <apply>
                        <exp/>
                        <apply>
                           <times/>
                           <ci>FoRT</ci>
                           <ci>deltaPsi_B</ci>
                           <cn cellml:units="dimensionless">6</cn>
                        </apply>
                     </apply>
                  </apply>
                  <apply>
                     <times/>
                     <apply>
                        <plus/>
                        <ci>r2</ci>
                        <apply>
                           <times/>
                           <ci>r3</ci>
                           <apply>
                              <exp/>
                              <apply>
                                 <times/>
                                 <ci>FoRT</ci>
                                 <ci>A_res</ci>
                              </apply>
                           </apply>
                        </apply>
                     </apply>
                     <apply>
                        <exp/>
                        <apply>
                           <times/>
                           <ci>FoRT</ci>
                           <ci>deltaPsi</ci>
                           <ci>g</ci>
                           <cn cellml:units="dimensionless">6</cn>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>J_Hres</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <ci>rho_res</ci>
                  <cn cellml:units="dimensionless">3.966</cn>
                  <apply>
                     <plus/>
                     <apply>
                        <times/>
                        <ci>ra</ci>
                        <apply>
                           <power/>
                           <cn cellml:units="dimensionless">10</cn>
                           <apply>
                              <times/>
                              <cn cellml:units="dimensionless">6</cn>
                              <ci>deltapH</ci>
                           </apply>
                        </apply>
                        <apply>
                           <exp/>
                           <apply>
                              <times/>
                              <ci>FoRT</ci>
                              <ci>A_res</ci>
                           </apply>
                        </apply>
                     </apply>
                     <apply>
                        <times/>
                        <ci>rb</ci>
                        <apply>
                           <power/>
                           <cn cellml:units="dimensionless">10</cn>
                           <apply>
                              <times/>
                              <cn cellml:units="dimensionless">6</cn>
                              <ci>deltapH</ci>
                           </apply>
                        </apply>
                     </apply>
                     <apply>
                        <times/>
                        <apply>
                           <minus/>
                           <cn cellml:units="dimensionless">1</cn>
                        </apply>
                        <apply>
                           <plus/>
                           <ci>ra</ci>
                           <ci>rb</ci>
                        </apply>
                        <apply>
                           <exp/>
                           <apply>
                              <times/>
                              <ci>g</ci>
                              <ci>FoRT</ci>
                              <ci>deltaPsi</ci>
                              <cn cellml:units="dimensionless">6</cn>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <plus/>
                  <apply>
                     <times/>
                     <apply>
                        <plus/>
                        <cn cellml:units="dimensionless">1</cn>
                        <apply>
                           <times/>
                           <ci>r1</ci>
                           <apply>
                              <exp/>
                              <apply>
                                 <times/>
                                 <ci>FoRT</ci>
                                 <ci>A_res</ci>
                              </apply>
                           </apply>
                        </apply>
                     </apply>
                     <apply>
                        <exp/>
                        <apply>
                           <times/>
                           <cn cellml:units="dimensionless">6</cn>
                           <ci>FoRT</ci>
                           <ci>deltaPsi_B</ci>
                        </apply>
                     </apply>
                  </apply>
                  <apply>
                     <times/>
                     <apply>
                        <plus/>
                        <ci>r2</ci>
                        <apply>
                           <times/>
                           <ci>r3</ci>
                           <apply>
                              <exp/>
                              <apply>
                                 <times/>
                                 <ci>FoRT</ci>
                                 <ci>A_res</ci>
                              </apply>
                           </apply>
                        </apply>
                     </apply>
                     <apply>
                        <exp/>
                        <apply>
                           <times/>
                           <ci>g</ci>
                           <cn cellml:units="dimensionless">6</cn>
                           <ci>FoRT</ci>
                           <ci>deltaPsi</ci>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>J_glyTotal</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <ci>beta_max</ci>
                  <apply>
                     <plus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <apply>
                        <times/>
                        <ci>beta1</ci>
                        <ci>Glc</ci>
                     </apply>
                  </apply>
                  <ci>beta2</ci>
                  <ci>Glc</ci>
                  <ci>ATP_i</ci>
               </apply>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <times/>
                     <ci>beta3</ci>
                     <ci>ATP_i</ci>
                  </apply>
                  <apply>
                     <times/>
                     <apply>
                        <plus/>
                        <cn cellml:units="dimensionless">1</cn>
                        <apply>
                           <times/>
                           <ci>beta4</ci>
                           <ci>ATP_i</ci>
                        </apply>
                     </apply>
                     <ci>beta5</ci>
                     <ci>Glc</ci>
                  </apply>
                  <apply>
                     <times/>
                     <apply>
                        <plus/>
                        <cn cellml:units="dimensionless">1</cn>
                        <apply>
                           <times/>
                           <ci>beta6</ci>
                           <ci>ATP_i</ci>
                        </apply>
                     </apply>
                     <ci>beta7</ci>
                     <ci>Glc</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>f_PDHa</ci>
            <apply>
               <divide/>
               <cn cellml:units="dimensionless">1</cn>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <times/>
                     <ci>u2</ci>
                     <apply>
                        <plus/>
                        <cn cellml:units="dimensionless">1</cn>
                        <apply>
                           <divide/>
                           <ci>u1</ci>
                           <apply>
                              <power/>
                              <apply>
                                 <plus/>
                                 <cn cellml:units="dimensionless">1</cn>
                                 <apply>
                                    <divide/>
                                    <ci>Ca_m</ci>
                                    <ci>KCa_PDH</ci>
                                 </apply>
                              </apply>
                              <cn cellml:units="dimensionless">2</cn>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>J_red</ci>
            <apply>
               <plus/>
               <ci>J_red_basal</ci>
               <apply>
                  <times/>
                  <cn cellml:units="dimensionless">6.3944</cn>
                  <ci>f_PDHa</ci>
                  <ci>J_glyTotal</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>J_pTCA</ci>
            <apply>
               <plus/>
               <apply>
                  <divide/>
                  <ci>J_red_basal</ci>
                  <cn cellml:units="dimensionless">3</cn>
               </apply>
               <apply>
                  <times/>
                  <cn cellml:units="dimensionless">0.84</cn>
                  <ci>f_PDHa</ci>
                  <ci>J_glyTotal</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>A_F1</ci>
            <apply>
               <times/>
               <ci>RToF</ci>
               <apply>
                  <ln/>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <ci>K_F1</ci>
                        <ci>ATP_m</ci>
                     </apply>
                     <apply>
                        <times/>
                        <ci>ADP_mfree</ci>
                        <ci>Pi_m</ci>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>J_pF1</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <ci>rho_F1</ci>
                  <apply>
                     <minus/>
                     <cn cellml:units="dimensionless">1</cn>
                  </apply>
                  <apply>
                     <plus/>
                     <apply>
                        <times/>
                        <apply>
                           <plus/>
                           <apply>
                              <times/>
                              <ci>pa</ci>
                              <apply>
                                 <power/>
                                 <cn cellml:units="dimensionless">10</cn>
                                 <apply>
                                    <times/>
                                    <cn cellml:units="dimensionless">3</cn>
                                    <ci>deltapH</ci>
                                 </apply>
                              </apply>
                           </apply>
                           <apply>
                              <times/>
                              <ci>pc1</ci>
                              <apply>
                                 <exp/>
                                 <apply>
                                    <times/>
                                    <cn cellml:units="dimensionless">3</cn>
                                    <ci>FoRT</ci>
                                    <ci>deltaPsi_B</ci>
                                 </apply>
                              </apply>
                           </apply>
                        </apply>
                        <apply>
                           <exp/>
                           <apply>
                              <times/>
                              <ci>FoRT</ci>
                              <ci>A_F1</ci>
                           </apply>
                        </apply>
                     </apply>
                     <apply>
                        <times/>
                        <apply>
                           <minus/>
                           <cn cellml:units="dimensionless">1</cn>
                        </apply>
                        <ci>pa</ci>
                        <apply>
                           <exp/>
                           <apply>
                              <times/>
                              <cn cellml:units="dimensionless">3</cn>
                              <ci>FoRT</ci>
                              <ci>deltaPsi</ci>
                           </apply>
                        </apply>
                     </apply>
                     <apply>
                        <times/>
                        <ci>pc2</ci>
                        <apply>
                           <exp/>
                           <apply>
                              <times/>
                              <ci>FoRT</ci>
                              <ci>A_F1</ci>
                           </apply>
                        </apply>
                        <apply>
                           <exp/>
                           <apply>
                              <times/>
                              <cn cellml:units="dimensionless">3</cn>
                              <ci>FoRT</ci>
                              <ci>deltaPsi</ci>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <plus/>
                  <apply>
                     <times/>
                     <apply>
                        <plus/>
                        <cn cellml:units="dimensionless">1</cn>
                        <apply>
                           <times/>
                           <ci>p1</ci>
                           <apply>
                              <exp/>
                              <apply>
                                 <times/>
                                 <ci>FoRT</ci>
                                 <ci>A_F1</ci>
                              </apply>
                           </apply>
                        </apply>
                     </apply>
                     <apply>
                        <exp/>
                        <apply>
                           <times/>
                           <cn cellml:units="dimensionless">3</cn>
                           <ci>FoRT</ci>
                           <ci>deltaPsi_B</ci>
                        </apply>
                     </apply>
                  </apply>
                  <apply>
                     <times/>
                     <apply>
                        <plus/>
                        <ci>p2</ci>
                        <apply>
                           <times/>
                           <ci>p3</ci>
                           <apply>
                              <exp/>
                              <apply>
                                 <times/>
                                 <ci>FoRT</ci>
                                 <ci>A_F1</ci>
                              </apply>
                           </apply>
                        </apply>
                     </apply>
                     <apply>
                        <exp/>
                        <apply>
                           <times/>
                           <cn cellml:units="dimensionless">3</cn>
                           <ci>FoRT</ci>
                           <ci>deltaPsi</ci>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>J_HF1</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <apply>
                     <minus/>
                     <cn cellml:units="dimensionless">1</cn>
                  </apply>
                  <ci>rho_F1</ci>
                  <cn cellml:units="dimensionless">3</cn>
                  <apply>
                     <plus/>
                     <apply>
                        <times/>
                        <ci>pa</ci>
                        <apply>
                           <power/>
                           <cn cellml:units="dimensionless">10</cn>
                           <apply>
                              <times/>
                              <cn cellml:units="dimensionless">3</cn>
                              <ci>deltapH</ci>
                           </apply>
                        </apply>
                        <apply>
                           <exp/>
                           <apply>
                              <times/>
                              <ci>FoRT</ci>
                              <ci>A_F1</ci>
                           </apply>
                        </apply>
                     </apply>
                     <apply>
                        <times/>
                        <ci>pb</ci>
                        <apply>
                           <power/>
                           <cn cellml:units="dimensionless">10</cn>
                           <apply>
                              <times/>
                              <cn cellml:units="dimensionless">3</cn>
                              <ci>deltapH</ci>
                           </apply>
                        </apply>
                     </apply>
                     <apply>
                        <times/>
                        <apply>
                           <minus/>
                           <cn cellml:units="dimensionless">1</cn>
                        </apply>
                        <apply>
                           <plus/>
                           <ci>pa</ci>
                           <ci>pb</ci>
                        </apply>
                        <apply>
                           <exp/>
                           <apply>
                              <times/>
                              <cn cellml:units="dimensionless">3</cn>
                              <ci>FoRT</ci>
                              <ci>deltaPsi</ci>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <plus/>
                  <apply>
                     <times/>
                     <apply>
                        <plus/>
                        <cn cellml:units="dimensionless">1</cn>
                        <apply>
                           <times/>
                           <ci>p1</ci>
                           <apply>
                              <exp/>
                              <apply>
                                 <times/>
                                 <ci>FoRT</ci>
                                 <ci>A_F1</ci>
                              </apply>
                           </apply>
                        </apply>
                     </apply>
                     <apply>
                        <exp/>
                        <apply>
                           <times/>
                           <cn cellml:units="dimensionless">3</cn>
                           <ci>FoRT</ci>
                           <ci>deltaPsi_B</ci>
                        </apply>
                     </apply>
                  </apply>
                  <apply>
                     <times/>
                     <apply>
                        <plus/>
                        <ci>p2</ci>
                        <apply>
                           <times/>
                           <ci>p3</ci>
                           <apply>
                              <exp/>
                              <apply>
                                 <times/>
                                 <ci>FoRT</ci>
                                 <ci>A_F1</ci>
                              </apply>
                           </apply>
                        </apply>
                     </apply>
                     <apply>
                        <exp/>
                        <apply>
                           <times/>
                           <cn cellml:units="dimensionless">3</cn>
                           <ci>FoRT</ci>
                           <ci>deltaPsi</ci>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>J_ANT</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <ci>Jmax_ANT</ci>
                  <apply>
                     <minus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <apply>
                        <times/>
                        <apply>
                           <divide/>
                           <apply>
                              <times/>
                              <ci>ATP4_i</ci>
                              <ci>ADP3_m</ci>
                           </apply>
                           <apply>
                              <times/>
                              <ci>ADP3_i</ci>
                              <ci>ATP4_m</ci>
                           </apply>
                        </apply>
                        <apply>
                           <exp/>
                           <apply>
                              <times/>
                              <apply>
                                 <minus/>
                                 <cn cellml:units="dimensionless">1</cn>
                              </apply>
                              <ci>FoRT</ci>
                              <ci>deltaPsi</ci>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <apply>
                     <plus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <apply>
                        <times/>
                        <apply>
                           <divide/>
                           <ci>ATP4_i</ci>
                           <ci>ADP3_i</ci>
                        </apply>
                        <apply>
                           <exp/>
                           <apply>
                              <times/>
                              <apply>
                                 <minus/>
                                 <cn cellml:units="dimensionless">1</cn>
                              </apply>
                              <ci>frac</ci>
                              <ci>FoRT</ci>
                              <ci>deltaPsi</ci>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
                  <apply>
                     <plus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <apply>
                        <divide/>
                        <ci>ADP3_m</ci>
                        <ci>ATP4_m</ci>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>PMF</ci>
            <apply>
               <minus/>
               <ci>deltaPsi</ci>
               <apply>
                  <times/>
                  <cn cellml:units="dimensionless">2.303</cn>
                  <ci>RToF</ci>
                  <ci>deltapH</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>J_Hleak</ci>
            <apply>
               <times/>
               <ci>g_H</ci>
               <ci>PMF</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>J_pGly</ci>
            <apply>
               <times/>
               <cn cellml:units="dimensionless">0.15</cn>
               <ci>J_glyTotal</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>J_hydSS</ci>
            <apply>
               <divide/>
               <ci>J_hyd_max</ci>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <power/>
                     <apply>
                        <divide/>
                        <ci>K_Glc</ci>
                        <ci>Glc</ci>
                     </apply>
                     <ci>nhyd</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>J_hyd</ci>
            <apply>
               <plus/>
               <apply>
                  <times/>
                  <ci>K_hyd</ci>
                  <ci>ATP_i</ci>
               </apply>
               <ci>J_hydSS</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>NADH_m</ci>
            </apply>
            <apply>
               <minus/>
               <ci>J_red</ci>
               <ci>J_o</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>NAD_m</ci>
            <apply>
               <minus/>
               <ci>total_NAD_m</ci>
               <ci>NADH_m</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>ADP_m</ci>
            </apply>
            <apply>
               <plus/>
               <ci>J_ANT</ci>
               <apply>
                  <times/>
                  <apply>
                     <minus/>
                     <cn cellml:units="dimensionless">1</cn>
                  </apply>
                  <ci>J_pTCA</ci>
               </apply>
               <apply>
                  <times/>
                  <apply>
                     <minus/>
                     <cn cellml:units="dimensionless">1</cn>
                  </apply>
                  <ci>J_pF1</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>ATP_m</ci>
            <apply>
               <minus/>
               <ci>total_ANP_m</ci>
               <ci>ADP_m</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>ADP_mfree</ci>
            <apply>
               <times/>
               <cn cellml:units="dimensionless">0.8</cn>
               <ci>ADP_m</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>ADP3_m</ci>
            <apply>
               <times/>
               <cn cellml:units="dimensionless">0.45</cn>
               <ci>ADP_mfree</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>ATP4_m</ci>
            <apply>
               <times/>
               <cn cellml:units="dimensionless">0.05</cn>
               <ci>ATP_m</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>ADP_i</ci>
            </apply>
            <apply>
               <plus/>
               <apply>
                  <divide/>
                  <apply>
                     <times/>
                     <apply>
                        <minus/>
                        <cn cellml:units="dimensionless">1</cn>
                     </apply>
                     <ci>J_ANT</ci>
                     <ci>V_MITO</ci>
                  </apply>
                  <ci>V_cyto</ci>
               </apply>
               <ci>J_hyd</ci>
               <apply>
                  <times/>
                  <apply>
                     <minus/>
                     <cn cellml:units="dimensionless">1</cn>
                  </apply>
                  <ci>J_pGly</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>ATP_i</ci>
            <apply>
               <minus/>
               <ci>total_ANP_i</ci>
               <ci>ADP_i</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>ADP_ifree</ci>
            <apply>
               <times/>
               <cn cellml:units="dimensionless">0.3</cn>
               <ci>ADP_i</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>ADP3_i</ci>
            <apply>
               <times/>
               <cn cellml:units="dimensionless">0.45</cn>
               <ci>ADP_ifree</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>MgADP_i</ci>
            <apply>
               <times/>
               <cn cellml:units="dimensionless">0.55</cn>
               <ci>ADP_ifree</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>ATP4_i</ci>
            <apply>
               <times/>
               <cn cellml:units="dimensionless">0.05</cn>
               <ci>ATP_i</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>Ca_PU</ci>
            </apply>
            <apply>
               <times/>
               <ci>fc</ci>
               <apply>
                  <plus/>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <apply>
                           <minus/>
                           <ci>J_NaCa</ci>
                           <ci>J_uni</ci>
                        </apply>
                        <ci>V_MITO</ci>
                     </apply>
                     <ci>V_PU</ci>
                  </apply>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <apply>
                           <minus/>
                           <ci>J_ERout</ci>
                           <ci>J_SERCA</ci>
                        </apply>
                        <ci>V_ER</ci>
                     </apply>
                     <ci>V_PU</ci>
                  </apply>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <apply>
                           <minus/>
                           <cn cellml:units="dimensionless">1</cn>
                        </apply>
                        <ci>J_leak</ci>
                        <ci>V_cyto</ci>
                     </apply>
                     <ci>V_PU</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>Ca_m</ci>
            </apply>
            <apply>
               <times/>
               <ci>fm</ci>
               <apply>
                  <minus/>
                  <ci>J_uni</ci>
                  <ci>J_NaCa</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>Ca_ER</ci>
            </apply>
            <apply>
               <times/>
               <ci>fe</ci>
               <apply>
                  <minus/>
                  <ci>J_SERCA</ci>
                  <ci>J_ERout</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>deltaPsi</ci>
            </apply>
            <apply>
               <times/>
               <apply>
                  <divide/>
                  <apply>
                     <times/>
                     <apply>
                        <minus/>
                        <cn cellml:units="dimensionless">1</cn>
                     </apply>
                     <ci>F</ci>
                     <ci>V_MITO</ci>
                     <cn cellml:units="millimole_to_nanomole">1000000</cn>
                     <cn cellml:units="dimensionless">1</cn>
                  </apply>
                  <ci>Cmito</ci>
               </apply>
               <apply>
                  <plus/>
                  <ci>J_Hleak</ci>
                  <apply>
                     <times/>
                     <apply>
                        <minus/>
                        <cn cellml:units="dimensionless">1</cn>
                     </apply>
                     <ci>J_Hres</ci>
                  </apply>
                  <ci>J_ANT</ci>
                  <ci>J_HF1</ci>
                  <apply>
                     <times/>
                     <cn cellml:units="dimensionless">2</cn>
                     <ci>J_uni</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>h</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <cn cellml:units="per_millimolar">1</cn>
                  <apply>
                     <minus/>
                     <ci>d_INH</ci>
                     <apply>
                        <times/>
                        <ci>h</ci>
                        <apply>
                           <plus/>
                           <ci>Ca_PU</ci>
                           <ci>d_INH</ci>
                        </apply>
                     </apply>
                  </apply>
               </apply>
               <ci>tauh</ci>
            </apply>
         </apply>
      </math>
   </component>
   <group>
      <relationship_ref relationship="encapsulation"/>
      <component_ref component="ICC_Membrane">
         <component_ref component="I_Na">
            <component_ref component="d_Na"/>
            <component_ref component="f_Na"/>
         </component_ref>
         <component_ref component="I_Ltype">
            <component_ref component="d_Ltype"/>
            <component_ref component="f_Ltype"/>
            <component_ref component="f_ca_Ltype"/>
         </component_ref>
         <component_ref component="I_VDDR">
            <component_ref component="d_VDDR"/>
            <component_ref component="f_VDDR"/>
         </component_ref>
         <component_ref component="I_BK">
            <component_ref component="d_BK"/>
         </component_ref>
         <component_ref component="I_CaCl">
            <component_ref component="d_CaCl"/>
         </component_ref>
         <component_ref component="I_kv11">
            <component_ref component="d_kv11"/>
            <component_ref component="f_kv11"/>
         </component_ref>
         <component_ref component="I_ERG">
            <component_ref component="d_ERG"/>
         </component_ref>
         <component_ref component="I_bk"/>
         <component_ref component="I_NSCC">
            <component_ref component="d_NSCC"/>
         </component_ref>
         <component_ref component="J_PMCA"/>
         <component_ref component="PU_unit"/>
      </component_ref>
   </group>
   <connection>
      <map_components component_1="d_Na" component_2="I_Na"/>
      <map_variables variable_1="Vm" variable_2="Vm"/>
      <map_variables variable_1="d_Na" variable_2="d_Na"/>
      <map_variables variable_1="time" variable_2="time"/>
      <map_variables variable_1="T_correction_Na" variable_2="T_correction_Na"/>
   </connection>
   <connection>
      <map_components component_1="f_Na" component_2="I_Na"/>
      <map_variables variable_1="Vm" variable_2="Vm"/>
      <map_variables variable_1="f_Na" variable_2="f_Na"/>
      <map_variables variable_1="time" variable_2="time"/>
      <map_variables variable_1="T_correction_Na" variable_2="T_correction_Na"/>
   </connection>
   <connection>
      <map_components component_1="I_Na" component_2="ICC_Membrane"/>
      <map_variables variable_1="Na_i" variable_2="Na_i"/>
      <map_variables variable_1="Vm" variable_2="Vm"/>
      <map_variables variable_1="I_Na" variable_2="I_Na"/>
      <map_variables variable_1="time" variable_2="time"/>
      <map_variables variable_1="RToF" variable_2="RToF"/>
      <map_variables variable_1="Na_o" variable_2="Na_o"/>
      <map_variables variable_1="T_correction_Na" variable_2="T_correction_Na"/>
   </connection>
   <connection>
      <map_components component_1="d_Ltype" component_2="I_Ltype"/>
      <map_variables variable_1="Vm" variable_2="Vm"/>
      <map_variables variable_1="d_Ltype" variable_2="d_Ltype"/>
      <map_variables variable_1="time" variable_2="time"/>
      <map_variables variable_1="T_correction_Ca" variable_2="T_correction_Ca"/>
   </connection>
   <connection>
      <map_components component_1="f_Ltype" component_2="I_Ltype"/>
      <map_variables variable_1="Vm" variable_2="Vm"/>
      <map_variables variable_1="f_Ltype" variable_2="f_Ltype"/>
      <map_variables variable_1="time" variable_2="time"/>
      <map_variables variable_1="T_correction_Ca" variable_2="T_correction_Ca"/>
   </connection>
   <connection>
      <map_components component_1="f_ca_Ltype" component_2="I_Ltype"/>
      <map_variables variable_1="f_ca_Ltype" variable_2="f_ca_Ltype"/>
      <map_variables variable_1="Ca_i" variable_2="Ca_i"/>
      <map_variables variable_1="time" variable_2="time"/>
      <map_variables variable_1="T_correction_Ca" variable_2="T_correction_Ca"/>
   </connection>
   <connection>
      <map_components component_1="I_Ltype" component_2="ICC_Membrane"/>
      <map_variables variable_1="Ca_i" variable_2="Ca_i"/>
      <map_variables variable_1="Vm" variable_2="Vm"/>
      <map_variables variable_1="I_Ltype" variable_2="I_Ltype"/>
      <map_variables variable_1="time" variable_2="time"/>
      <map_variables variable_1="RToF" variable_2="RToF"/>
      <map_variables variable_1="Ca_o" variable_2="Ca_o"/>
      <map_variables variable_1="T_correction_Ca" variable_2="T_correction_Ca"/>
   </connection>
   <connection>
      <map_components component_1="d_VDDR" component_2="I_VDDR"/>
      <map_variables variable_1="Vm" variable_2="Vm"/>
      <map_variables variable_1="d_VDDR" variable_2="d_VDDR"/>
      <map_variables variable_1="time" variable_2="time"/>
      <map_variables variable_1="T_correction_Ca" variable_2="T_correction_Ca"/>
   </connection>
   <connection>
      <map_components component_1="f_VDDR" component_2="I_VDDR"/>
      <map_variables variable_1="Vm" variable_2="Vm"/>
      <map_variables variable_1="f_VDDR" variable_2="f_VDDR"/>
      <map_variables variable_1="time" variable_2="time"/>
      <map_variables variable_1="T_correction_Ca" variable_2="T_correction_Ca"/>
   </connection>
   <connection>
      <map_components component_1="I_VDDR" component_2="ICC_Membrane"/>
      <map_variables variable_1="Ca_i" variable_2="Ca_i"/>
      <map_variables variable_1="Vm" variable_2="Vm"/>
      <map_variables variable_1="I_VDDR" variable_2="I_VDDR"/>
      <map_variables variable_1="time" variable_2="time"/>
      <map_variables variable_1="RToF" variable_2="RToF"/>
      <map_variables variable_1="Ca_o" variable_2="Ca_o"/>
      <map_variables variable_1="T_correction_Ca" variable_2="T_correction_Ca"/>
   </connection>
   <connection>
      <map_components component_1="d_kv11" component_2="I_kv11"/>
      <map_variables variable_1="Vm" variable_2="Vm"/>
      <map_variables variable_1="d_kv11" variable_2="d_kv11"/>
      <map_variables variable_1="time" variable_2="time"/>
      <map_variables variable_1="T_correction_K" variable_2="T_correction_K"/>
   </connection>
   <connection>
      <map_components component_1="f_kv11" component_2="I_kv11"/>
      <map_variables variable_1="Vm" variable_2="Vm"/>
      <map_variables variable_1="f_kv11" variable_2="f_kv11"/>
      <map_variables variable_1="time" variable_2="time"/>
      <map_variables variable_1="T_correction_K" variable_2="T_correction_K"/>
   </connection>
   <connection>
      <map_components component_1="I_kv11" component_2="ICC_Membrane"/>
      <map_variables variable_1="K_i" variable_2="K_i"/>
      <map_variables variable_1="Vm" variable_2="Vm"/>
      <map_variables variable_1="I_kv11" variable_2="I_kv11"/>
      <map_variables variable_1="time" variable_2="time"/>
      <map_variables variable_1="RToF" variable_2="RToF"/>
      <map_variables variable_1="K_o" variable_2="K_o"/>
      <map_variables variable_1="T_correction_K" variable_2="T_correction_K"/>
   </connection>
   <connection>
      <map_components component_1="d_ERG" component_2="I_ERG"/>
      <map_variables variable_1="Vm" variable_2="Vm"/>
      <map_variables variable_1="d_ERG" variable_2="d_ERG"/>
      <map_variables variable_1="time" variable_2="time"/>
      <map_variables variable_1="T_correction_K" variable_2="T_correction_K"/>
   </connection>
   <connection>
      <map_components component_1="I_ERG" component_2="ICC_Membrane"/>
      <map_variables variable_1="K_i" variable_2="K_i"/>
      <map_variables variable_1="Vm" variable_2="Vm"/>
      <map_variables variable_1="I_ERG" variable_2="I_ERG"/>
      <map_variables variable_1="time" variable_2="time"/>
      <map_variables variable_1="RToF" variable_2="RToF"/>
      <map_variables variable_1="K_o" variable_2="K_o"/>
      <map_variables variable_1="T_correction_K" variable_2="T_correction_K"/>
   </connection>
   <connection>
      <map_components component_1="d_CaCl" component_2="I_CaCl"/>
      <map_variables variable_1="Ca_i" variable_2="Ca_i"/>
      <map_variables variable_1="d_CaCl" variable_2="d_CaCl"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="I_CaCl" component_2="ICC_Membrane"/>
      <map_variables variable_1="Cl_i" variable_2="Cl_i"/>
      <map_variables variable_1="Ca_i" variable_2="Ca_i"/>
      <map_variables variable_1="Vm" variable_2="Vm"/>
      <map_variables variable_1="I_CaCl" variable_2="I_CaCl"/>
      <map_variables variable_1="time" variable_2="time"/>
      <map_variables variable_1="RToF" variable_2="RToF"/>
      <map_variables variable_1="Cl_o" variable_2="Cl_o"/>
   </connection>
   <connection>
      <map_components component_1="d_BK" component_2="I_BK"/>
      <map_variables variable_1="Vm" variable_2="Vm"/>
      <map_variables variable_1="Ca_i" variable_2="Ca_i"/>
      <map_variables variable_1="d_BK" variable_2="d_BK"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="I_BK" component_2="ICC_Membrane"/>
      <map_variables variable_1="K_i" variable_2="K_i"/>
      <map_variables variable_1="Ca_i" variable_2="Ca_i"/>
      <map_variables variable_1="Vm" variable_2="Vm"/>
      <map_variables variable_1="I_BK" variable_2="I_BK"/>
      <map_variables variable_1="time" variable_2="time"/>
      <map_variables variable_1="RToF" variable_2="RToF"/>
      <map_variables variable_1="K_o" variable_2="K_o"/>
      <map_variables variable_1="T_correction_BK" variable_2="T_correction_BK"/>
   </connection>
   <connection>
      <map_components component_1="I_bk" component_2="ICC_Membrane"/>
      <map_variables variable_1="K_i" variable_2="K_i"/>
      <map_variables variable_1="Vm" variable_2="Vm"/>
      <map_variables variable_1="I_bk" variable_2="I_bk"/>
      <map_variables variable_1="RToF" variable_2="RToF"/>
      <map_variables variable_1="K_o" variable_2="K_o"/>
   </connection>
   <connection>
      <map_components component_1="d_NSCC" component_2="I_NSCC"/>
      <map_variables variable_1="Ca_PU" variable_2="Ca_PU"/>
      <map_variables variable_1="d_NSCC" variable_2="d_NSCC"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="I_NSCC" component_2="ICC_Membrane"/>
      <map_variables variable_1="K_i" variable_2="K_i"/>
      <map_variables variable_1="Na_i" variable_2="Na_i"/>
      <map_variables variable_1="Vm" variable_2="Vm"/>
      <map_variables variable_1="I_NSCC" variable_2="I_NSCC"/>
      <map_variables variable_1="time" variable_2="time"/>
      <map_variables variable_1="RToF" variable_2="RToF"/>
      <map_variables variable_1="K_o" variable_2="K_o"/>
      <map_variables variable_1="Na_o" variable_2="Na_o"/>
   </connection>
   <connection>
      <map_components component_1="PU_unit" component_2="I_NSCC"/>
      <map_variables variable_1="Ca_PU" variable_2="Ca_PU"/>
   </connection>
   <connection>
      <map_components component_1="PU_unit" component_2="ICC_Membrane"/>
      <map_variables variable_1="Vol" variable_2="Vol"/>
      <map_variables variable_1="fc" variable_2="fc"/>
      <map_variables variable_1="V_cyto" variable_2="V_cyto"/>
      <map_variables variable_1="FoRT" variable_2="FoRT"/>
      <map_variables variable_1="RToF" variable_2="RToF"/>
      <map_variables variable_1="Na_i" variable_2="Na_i"/>
      <map_variables variable_1="Ca_i" variable_2="Ca_i"/>
      <map_variables variable_1="time" variable_2="time"/>
      <map_variables variable_1="J_leak" variable_2="J_leak"/>
      <map_variables variable_1="F" variable_2="F"/>
   </connection>
   <connection>
      <map_components component_1="J_PMCA" component_2="ICC_Membrane"/>
      <map_variables variable_1="J_PMCA" variable_2="J_PMCA"/>
      <map_variables variable_1="Ca_i" variable_2="Ca_i"/>
   </connection>
   <connection>
      <map_components component_1="ICC_Membrane" component_2="Environment"/>
      <map_variables variable_1="Ca_o" variable_2="Ca_o"/>
      <map_variables variable_1="Na_o" variable_2="Na_o"/>
      <map_variables variable_1="K_o" variable_2="K_o"/>
      <map_variables variable_1="Cl_o" variable_2="Cl_o"/>
      <map_variables variable_1="RToF" variable_2="RToF"/>
      <map_variables variable_1="FoRT" variable_2="FoRT"/>
      <map_variables variable_1="F" variable_2="F"/>
      <map_variables variable_1="T_correction_Ca" variable_2="T_correction_Ca"/>
      <map_variables variable_1="T_correction_Na" variable_2="T_correction_Na"/>
      <map_variables variable_1="T_correction_K" variable_2="T_correction_K"/>
      <map_variables variable_1="T_correction_BK" variable_2="T_correction_BK"/>
   </connection>
   <connection>
      <map_components component_1="Time" component_2="ICC_Membrane"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
</model>