Rendering of the source text

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

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

CellML 1.0 was used to generate this model
http://www.cellml.org/
-->
<model xmlns:dc="http://purl.org/dc/elements/1.1/" xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#" xmlns:bqs="http://www.cellml.org/bqs/1.0#" xmlns:dcterms="http://purl.org/dc/terms/" xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" xmlns:ns7="http://www.cellml.org/metadata/simulation/1.0#" xmlns:xlink="http://www.w3.org/1999/xlink" name="ekaterinburg_rat" cmeta:id="ekaterinburg_rat" xmlns="http://www.cellml.org/cellml/1.0#" xmlns:cellml="http://www.cellml.org/cellml/1.0#" xmlns:cmeta="http://www.cellml.org/metadata/1.0#">
	<documentation xmlns="http://cellml.org/tmp-documentation">
		<article>
			<articleinfo>
				<title>...</title>
				<author>
					<firstname>Leonid</firstname>
					<surname>Katsnelson</surname>
					<affiliation>
						<shortaffil>Institute of immunology and physiology</shortaffil>
					</affiliation>
				</author>
			</articleinfo>
			<section id="sec_status">
				<title>Model Status</title>
				<para>
This version has been curated and unit-checked by Penny Noble and is known to run in COR and PCEnv. This variant is the STRETCH model, which includes acetylcholine dependent modulation and the electromechanical processes simulated by Variant 02, and also includes simulation of stretch and fibroblasts. This model is parameterised the same as the BASIC model, with all additional variables set to zero. The values for these variables can be found in <ulink url="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?cmd=Retrieve&amp;db=PubMed&amp;list_uids=9487284&amp;dopt=Abstract">Noble et al.'s 1998 paper.</ulink> This model is also associated with a PCEnv session file.
    </para>
			</section>
			<sect1 id="sec_structure">
				<title>Model Structure</title>
				<para>In their 1998 paper, Denis Noble, Anthony Varghese, Peter Kohl and Penny Noble extended the guinea-pig ventricular cell model (originally developed in 1991 by Noble <emphasis>et al</emphasis>) to include accumulation and depletion of calcium in a diadic space between the sarcolemma and the sarcoplasmic reticulum, where the majority of calcium-induced calcium release is triggered.  The model also includes rapid (i_Kr) and slow (i_Ks) components of the delayed rectifier current, and length- and tension-dependent changes in mechanical and electrophysiological processes have been incorporated (see <xref linkend="fig_cell_diagram"/> below).  Drug-receptor interactions have also started to be modeled using the sodium channel as an example.  The drugs used have the effect of blocking the inactivation h-gate of the fast sodium channel.
</para>
				<para>
The complete original paper reference is cited below:
</para>
				<para>
					<ulink url="http://www.pulsus.com/CARDIOL/14_01/nobl_ed.htm">Improved guinea-pig ventricular cell model incorporating a diadic space, I<subscript>Kr</subscript> and I<subscript>Ks</subscript>, and length- and tension-dependent processes</ulink>, Denis Noble, Anthony Varghese, Peter Kohl and Penelope Noble, 1998, <ulink url="http://www.pulsus.com/CARDIOL/home.htm">
						<emphasis>Can J Cardiol</emphasis>
					</ulink>, 14, 123-134.  <ulink url="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?cmd=Retrieve&amp;db=PubMed&amp;list_uids=9487284&amp;dopt=Abstract">PubMed ID: 9487284</ulink>
				</para>
				<para>
The raw CellML description of the Noble'98 model can be downloaded in various formats as described in <xref linkend="sec_download_this_model"/>.  For an example of a more complete documentation for an electrophysiological model, see <ulink url="${HTML_EXMPL_HHSA_INTRO}">The Hodgkin-Huxley Squid Axon Model, 1952</ulink>. 
</para>
				<informalfigure float="0" id="fig_cell_diagram">
					<mediaobject>
						<imageobject>
							<objectinfo>
								<title>cell diagram of Noble'98 model showing ionic currents, pumps and exchangers within the sarcolemma and the sarcoplasmic reticulum</title>
							</objectinfo>
							<imagedata fileref="noble_1998.png"/>
						</imageobject>
					</mediaobject>
					<caption>A schematic diagram describing the current flows across the cell membrane that are captured in the Noble'98 model.</caption>
				</informalfigure>
			</sect1>
		</article>
	</documentation>
	<rdf:RDF>
		<rdf:Seq rdf:about="rdf:#0ce2c5c0-6262-4a24-917e-1fc336b0cabb">
			<rdf:li rdf:resource="rdf:#492a8d76-a4ab-4919-8506-fd55bb48c9e6"/>
			<rdf:li rdf:resource="rdf:#8bb3cc2e-3939-4d46-b384-50f8f57d72d5"/>
			<rdf:li rdf:resource="rdf:#8bda7cea-82ea-48a2-a242-f068440dd53b"/>
			<rdf:li rdf:resource="rdf:#1c28a87b-ecfc-4a17-9e63-1069ce345cd1"/>
		</rdf:Seq>
		<rdf:Description rdf:about="rdf:#72025205-c2ee-46ba-945e-6ad6e24afc85">
			<vCard:Given>D</vCard:Given>
			<vCard:Family>Noble</vCard:Family>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#$hm4db">
			<ns7:boundIntervals rdf:resource="rdf:#$im4db"/>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#db5ec880-ef58-4c15-8e0d-36bf8e7e497f">
			<vCard:Given>Penny</vCard:Given>
			<vCard:Family>Noble</vCard:Family>
			<vCard:Other>J</vCard:Other>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#9c09dc78-7f49-4457-93ed-c7b9955a38c8">
			<vCard:Given>A</vCard:Given>
			<vCard:Family>Varghese</vCard:Family>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#a5f2a5b3-2edd-47c1-8582-1eab7e8f752b">
			<vCard:Given>James</vCard:Given>
			<vCard:Family>Lawson</vCard:Family>
			<vCard:Other>Richard</vCard:Other>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#d7beb4e1-8910-4c28-b560-61c8d1d4a814">
			<dcterms:W3CDTF>2007-12-04T15:11:55+13:00</dcterms:W3CDTF>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#2ecb01cc-a373-40fa-9501-4aa93ccd5d8a">
			<dcterms:modified rdf:resource="rdf:#a96b2bdc-cbd5-4502-ac7c-203b452d742b"/>
			<rdf:value>Penny Noble was informed by a community member that INaCads was not included in the dCads/dt ODE and hence the system was out of balance. This has been fixed in version 07. An alteration was also made to the model to remove two differential expressions in the intracellular_calcium concentration component to allow PCEnv 0.2 to run the model.</rdf:value>
			<cmeta:modifier xmlns:cmeta="http://www.cellml.org/metadata/1.0#" rdf:resource="rdf:#dec66365-7c06-4d35-8dfb-28989da96412"/>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#6b5feb3f-e3a7-44d9-ab1c-7874f76bb698">
			<vCard:FN>James Lawson</vCard:FN>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#f7a88ef0-b25b-4f1b-8c84-925bdd5b9e1c">
			<vCard:N rdf:resource="rdf:#16333cab-2375-4624-a3e1-16842ffa5ddc"/>
		</rdf:Description>
		<rdf:Description rdf:about="#noble_model_1998_stretch">
			<ns7:simulation rdf:resource="rdf:#$hm4db"/>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#8bda7cea-82ea-48a2-a242-f068440dd53b">
			<rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
			<vCard:N rdf:resource="rdf:#7342d718-2ba4-41c7-8a1d-ff647241975f"/>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#9f557c82-27e9-4db6-93a8-a9bd42afe2a6">
			<vCard:ORG rdf:resource="rdf:#f663fe1d-2bd6-497b-94d0-eef64ccedde3"/>
			<vCard:EMAIL rdf:resource="rdf:#8c5aca23-4409-46cb-bbca-4ca7826d19e6"/>
			<vCard:N rdf:resource="rdf:#db5ec880-ef58-4c15-8e0d-36bf8e7e497f"/>
		</rdf:Description>
		<rdf:Description rdf:about="">
			<dc:publisher>Department of Physiology, Anatomy &amp; Genetics, University of Oxford</dc:publisher>
			<cmeta:comment xmlns:cmeta="http://www.cellml.org/metadata/1.0#" rdf:resource="rdf:#1b1008fc-cda4-422b-a815-3a82b5e9c32a"/>
			<dcterms:created rdf:resource="rdf:#1f496d47-a63b-480e-941e-53eb648706ec"/>
			<dc:creator rdf:resource="rdf:#9f557c82-27e9-4db6-93a8-a9bd42afe2a6"/>
			<cmeta:modification xmlns:cmeta="http://www.cellml.org/metadata/1.0#" rdf:resource="rdf:#6761e204-596b-44f6-b8cd-51afc006050e"/>
			<cmeta:modification xmlns:cmeta="http://www.cellml.org/metadata/1.0#" rdf:resource="rdf:#2daea026-24aa-4913-8020-2620ff205c30"/>
			<cmeta:modification xmlns:cmeta="http://www.cellml.org/metadata/1.0#" rdf:resource="rdf:#2ecb01cc-a373-40fa-9501-4aa93ccd5d8a"/>
			<cmeta:modification xmlns:cmeta="http://www.cellml.org/metadata/1.0#" rdf:resource="rdf:#43217e3a-8f02-4e5e-bf24-4d9a005b8777"/>
			<cmeta:modification xmlns:cmeta="http://www.cellml.org/metadata/1.0#" rdf:resource="rdf:#4fa87bf8-9c73-415d-ac3a-b97e8e3887d6"/>
			<cmeta:modification xmlns:cmeta="http://www.cellml.org/metadata/1.0#" rdf:resource="rdf:#bb802211-58b6-4242-94b0-8eec067e1049"/>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#1f496d47-a63b-480e-941e-53eb648706ec">
			<dcterms:W3CDTF>2006-05-03T00:00:00+00:00</dcterms:W3CDTF>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#f663fe1d-2bd6-497b-94d0-eef64ccedde3">
			<vCard:Orgname>Oxford University</vCard:Orgname>
			<vCard:Orgunit>Department of Physiology, Anatomy &amp; Genetics</vCard:Orgunit>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#16333cab-2375-4624-a3e1-16842ffa5ddc">
			<vCard:Given>James</vCard:Given>
			<vCard:Family>Lawson</vCard:Family>
			<vCard:Other>Richard</vCard:Other>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#1bbdcfe9-2109-4062-86de-7fc13439a418">
			<dcterms:W3CDTF>2007-07-18T12:59:52+12:00</dcterms:W3CDTF>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#6761e204-596b-44f6-b8cd-51afc006050e">
			<dcterms:modified rdf:resource="rdf:#1bca4043-64c7-44e2-97ec-092c2d8c44cb"/>
			<rdf:value>Note that the missing i_NCX_ds component was added to the dCai/dt rather than the dCads/dt equation in order to achieve balance in the model.</rdf:value>
			<cmeta:modifier xmlns:cmeta="http://www.cellml.org/metadata/1.0#" rdf:resource="rdf:#e3a54851-287a-46f2-9a61-0f77ae0c1458"/>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#$im4db">
			<rdf:rest rdf:resource="http://www.w3.org/1999/02/22-rdf-syntax-ns#nil"/>
			<rdf:first rdf:resource="rdf:#$jm4db"/>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#4fa87bf8-9c73-415d-ac3a-b97e8e3887d6">
			<dcterms:modified rdf:resource="rdf:#86d71735-52ed-4a55-a0a2-1cc257a8217f"/>
			<rdf:value>This model has been curated by Penny Noble</rdf:value>
			<cmeta:modifier xmlns:cmeta="http://www.cellml.org/metadata/1.0#" rdf:resource="rdf:#e6b753c8-3d08-4984-9f32-e577c2597f13"/>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#a96b2bdc-cbd5-4502-ac7c-203b452d742b">
			<dcterms:W3CDTF>2007-09-06T13:01:21+12:00</dcterms:W3CDTF>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#2daea026-24aa-4913-8020-2620ff205c30">
			<dcterms:modified rdf:resource="rdf:#d7beb4e1-8910-4c28-b560-61c8d1d4a814"/>
			<rdf:value>Units checked, curated.</rdf:value>
			<cmeta:modifier xmlns:cmeta="http://www.cellml.org/metadata/1.0#" rdf:resource="rdf:#19aeec40-a962-4cbe-9ab2-65f9b1438afd"/>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#1b1008fc-cda4-422b-a815-3a82b5e9c32a">
			<dc:creator rdf:resource="rdf:#6b5feb3f-e3a7-44d9-ab1c-7874f76bb698"/>
			<rdf:value>This version was created by Penny Noble and is known to read into COR. James Lawson fixed an error (19/04/07) that prevented it reading in PCEnv (unsupported predefined operator diff in component intracellular_calcium_concentration. This variant of version 06 is the STRETCH model, which includes acetylcholine dependent modulation and the electromechanical processes simulated by Variant 02, and also includes simulation of stretch and fibroblasts. This model is parameterised the same as the BASIC model, with all additional variables set to zero.
To simulate multiple action potentials with this model the following integrator parameters must be used: max step size must be 0.001 or less, max point density must be 10,000 or more.</rdf:value>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#43217e3a-8f02-4e5e-bf24-4d9a005b8777">
			<dcterms:modified rdf:resource="rdf:#1bbdcfe9-2109-4062-86de-7fc13439a418"/>
			<rdf:value>Changed stim_end in stimulus protocol from 10.1 to 1,000,000 seconds</rdf:value>
			<cmeta:modifier xmlns:cmeta="http://www.cellml.org/metadata/1.0#" rdf:resource="rdf:#e2a1ec0e-f041-4b1c-af55-110096409042"/>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#8bb3cc2e-3939-4d46-b384-50f8f57d72d5">
			<rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
			<vCard:N rdf:resource="rdf:#9c09dc78-7f49-4457-93ed-c7b9955a38c8"/>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#7f2d6d6b-344c-4a7d-a15b-a4747f3898ca">
			<dc:title>Canadian Journal of Cardiology</dc:title>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#cfbad6be-cc46-4275-a900-a78cb6fd8bd4">
			<dcterms:W3CDTF>2007-06-22T13:59:26+12:00</dcterms:W3CDTF>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#8c52a9a5-590f-4c2e-b319-a344696e33f1">
			<vCard:Given>Penny</vCard:Given>
			<vCard:Family>Noble</vCard:Family>
			<vCard:Other>J</vCard:Other>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#$6wmB13">
			<ns7:nonstandard-pointDensity>100000</ns7:nonstandard-pointDensity>
			<ns7:maximumStepSize>0.001</ns7:maximumStepSize>
			<ns7:nonstandard-algorithm>bdf15</ns7:nonstandard-algorithm>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#e3a54851-287a-46f2-9a61-0f77ae0c1458">
			<vCard:N rdf:resource="rdf:#feb02716-3dd9-4ca6-8307-90e863c534d1"/>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#86d71735-52ed-4a55-a0a2-1cc257a8217f">
			<dcterms:W3CDTF>2007-02-06T13:00:00+13:00</dcterms:W3CDTF>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#dec66365-7c06-4d35-8dfb-28989da96412">
			<vCard:N rdf:resource="rdf:#b3346d0b-66a3-4213-a4e0-321b7948aa9d"/>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#bb802211-58b6-4242-94b0-8eec067e1049">
			<dcterms:modified rdf:resource="rdf:#cfbad6be-cc46-4275-a900-a78cb6fd8bd4"/>
			<rdf:value>Several variables have been given cmeta:id's to allow creation of a PCEnv session file.</rdf:value>
			<cmeta:modifier xmlns:cmeta="http://www.cellml.org/metadata/1.0#" rdf:resource="rdf:#f7a88ef0-b25b-4f1b-8c84-925bdd5b9e1c"/>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#492a8d76-a4ab-4919-8506-fd55bb48c9e6">
			<rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
			<vCard:N rdf:resource="rdf:#72025205-c2ee-46ba-945e-6ad6e24afc85"/>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#feb02716-3dd9-4ca6-8307-90e863c534d1">
			<vCard:Given>Penny</vCard:Given>
			<vCard:Family>Noble</vCard:Family>
			<vCard:Other/>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#34f1c61b-8e98-4ed5-a473-3ba951668467">
			<vCard:Given>P</vCard:Given>
			<vCard:Family>Noble</vCard:Family>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#1c28a87b-ecfc-4a17-9e63-1069ce345cd1">
			<rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
			<vCard:N rdf:resource="rdf:#34f1c61b-8e98-4ed5-a473-3ba951668467"/>
		</rdf:Description>
		<rdf:Description rdf:about="#noble_varghese_kohl_noble_1998_stretch">
			<ns7:simulation rdf:resource="rdf:#$4wmB13"/>
			<bqs:reference rdf:resource="rdf:#89bb8cae-ce4c-427e-b18f-fe308d30dd9e"/>
			<bqs:reference rdf:parseType="Resource">
				<dc:subject rdf:parseType="Resource">
					<bqs:subject_type>keyword</bqs:subject_type>
					<rdf:value>
						<rdf:Bag>
							<rdf:li>diadic space</rdf:li>
							<rdf:li>guinea-pig</rdf:li>
							<rdf:li>excitation-contraction coupling</rdf:li>
							<rdf:li>ventricular myocyte</rdf:li>
							<rdf:li>calcium dynamics</rdf:li>
							<rdf:li>electrophysiology</rdf:li>
							<rdf:li>cardiac</rdf:li>
						</rdf:Bag>
					</rdf:value>
				</dc:subject>
			</bqs:reference>
			<cmeta:comment xmlns:cmeta="http://www.cellml.org/metadata/1.0#" rdf:resource="rdf:#56512c4d-6eaf-4d72-9205-e876c9ea5e70"/>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#56512c4d-6eaf-4d72-9205-e876c9ea5e70">
			<dc:creator rdf:resource="rdf:#5ba8104a-bf61-4827-8ecf-ec222459ce5b"/>
			<rdf:value>This file is a CellML description of Noble, Carghese, Kohl and Noble's 1998 extension of their 1991 guinea pig ventricular cell model. This model incorporates a diadic space, IKr and IKs currents, and length and tension dependent processes. Variant 3 is parameterised for the 'stretch' model.</rdf:value>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#e6b753c8-3d08-4984-9f32-e577c2597f13">
			<vCard:N rdf:resource="rdf:#8c52a9a5-590f-4c2e-b319-a344696e33f1"/>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#7342d718-2ba4-41c7-8a1d-ff647241975f">
			<vCard:Given>P</vCard:Given>
			<vCard:Family>Kohl</vCard:Family>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#89bb8cae-ce4c-427e-b18f-fe308d30dd9e">
			<bqs:Pubmed_id>9487284</bqs:Pubmed_id>
			<bqs:JournalArticle rdf:resource="rdf:#7412123a-3acd-4f20-aea4-11f746133c40"/>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#07a73155-cbd2-4f61-848e-2c11a2608cf0">
			<dcterms:W3CDTF>1998-01-00 00:00</dcterms:W3CDTF>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#8c5aca23-4409-46cb-bbca-4ca7826d19e6">
			<rdf:type rdf:resource="http://imc.org/vCard/3.0#internet"/>
			<rdf:value>penny.noble@dpag.ox.ac.uk</rdf:value>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#$jm4db">
			<ns7:nonstandard-pointDensity>50000</ns7:nonstandard-pointDensity>
			<ns7:maximumStepSize>0.001</ns7:maximumStepSize>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#1bca4043-64c7-44e2-97ec-092c2d8c44cb">
			<dcterms:W3CDTF>2008-05-14T03:16:52+12:00</dcterms:W3CDTF>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#c261f321-6315-485c-b701-bff0ae199be0">
			<vCard:Given>Penny</vCard:Given>
			<vCard:Family>Noble</vCard:Family>
			<vCard:Other>J</vCard:Other>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#e2a1ec0e-f041-4b1c-af55-110096409042">
			<vCard:N rdf:resource="rdf:#a5f2a5b3-2edd-47c1-8582-1eab7e8f752b"/>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#$4wmB13">
			<ns7:boundIntervals rdf:resource="rdf:#$5wmB13"/>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#5ba8104a-bf61-4827-8ecf-ec222459ce5b">
			<vCard:FN>James Lawson</vCard:FN>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#19aeec40-a962-4cbe-9ab2-65f9b1438afd">
			<vCard:N rdf:resource="rdf:#c261f321-6315-485c-b701-bff0ae199be0"/>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#7412123a-3acd-4f20-aea4-11f746133c40">
			<dc:creator rdf:resource="rdf:#0ce2c5c0-6262-4a24-917e-1fc336b0cabb"/>
			<dc:title>Improved guinea-pig ventricular cell model incorporating a diadic space, IKr and IKs, and length- and tension-dependent processes</dc:title>
			<bqs:volume>14(1)</bqs:volume>
			<bqs:first_page>123</bqs:first_page>
			<bqs:Journal rdf:resource="rdf:#7f2d6d6b-344c-4a7d-a15b-a4747f3898ca"/>
			<dcterms:issued rdf:resource="rdf:#07a73155-cbd2-4f61-848e-2c11a2608cf0"/>
			<bqs:last_page>134</bqs:last_page>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#$5wmB13">
			<rdf:rest rdf:resource="http://www.w3.org/1999/02/22-rdf-syntax-ns#nil"/>
			<rdf:first rdf:resource="rdf:#$6wmB13"/>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#b3346d0b-66a3-4213-a4e0-321b7948aa9d">
			<vCard:Given>James</vCard:Given>
			<vCard:Family>Lawson</vCard:Family>
			<vCard:Other>Richard</vCard:Other>
		</rdf:Description>
	</rdf:RDF>   <units name="per_millimolar">
      <unit units="millimolar" exponent="-1"/>
   </units>
   <units name="per_second2">
      <unit units="second" exponent="-2"/>
   </units>
   <units name="micrometre_per_second">
      <unit units="micrometre"/>
      <unit units="per_second"/>
   </units>
   <units name="micrometre_per_second2">
      <unit units="micrometre"/>
      <unit units="per_second2"/>
   </units>
   <units name="second_per_micrometre">
      <unit units="second"/>
      <unit units="per_micrometre"/>
   </units>
   <units name="millinewton">
      <unit units="newton" prefix="milli"/>
   </units>
   <units name="millinewton_per_second">
      <unit units="millinewton"/>
      <unit units="per_second"/>
   </units>
   <units name="millinewton_second_per_micrometre">
      <unit units="millinewton"/>
      <unit units="second_per_micrometre"/>
   </units>
   <units name="micrometre">
      <unit units="metre" prefix="micro"/>
   </units>
   <units name="per_micrometre">
      <unit units="micrometre" exponent="-1"/>
   </units>
   <units name="per_second">
      <unit units="second" exponent="-1"/>
   </units>
   <units name="millivolt">
      <unit units="volt" prefix="milli"/>
   </units>
   <units name="per_millivolt">
      <unit units="volt" prefix="milli" exponent="-1"/>
   </units>
   <units name="per_millivolt_second">
      <unit units="millivolt" exponent="-1"/>
      <unit units="second" exponent="-1"/>
   </units>
   <units name="microS">
      <unit units="siemens" prefix="micro"/>
   </units>
   <units name="microF">
      <unit units="farad" prefix="micro"/>
   </units>
   <units name="nanoA">
      <unit units="ampere" prefix="nano"/>
   </units>
   <units name="mA_nA">
      <unit units="ampere" prefix="milli"/>
      <unit units="nanoA" exponent="-1"/>
   </units>
   <units name="nanoA_per_millimolar">
      <unit units="ampere" prefix="nano"/>
      <unit units="millimolar" exponent="-1"/>
   </units>
   <units name="millimolar">
      <unit units="mole" prefix="milli"/>
      <unit units="litre" exponent="-1"/>
   </units>
   <units name="millimolar4">
      <unit units="millimolar" exponent="4"/>
   </units>
   <units name="millimolar_per_second">
      <unit units="millimolar"/>
      <unit units="second" exponent="-1"/>
   </units>
   <units name="per_millimolar_second">
      <unit units="millimolar" exponent="-1"/>
      <unit units="second" exponent="-1"/>
   </units>
   <units name="joule_per_kilomole_kelvin">
      <unit units="joule"/>
      <unit units="mole" prefix="kilo" exponent="-1"/>
      <unit units="kelvin" exponent="-1"/>
   </units>
   <units name="coulomb_per_mole">
      <unit units="coulomb"/>
      <unit units="mole" exponent="-1"/>
   </units>
   <units name="micrometre3">
      <unit units="metre" prefix="micro" exponent="3"/>
   </units>
   <units name="litre_micrometre3">
      <unit units="litre"/>
      <unit units="micrometre3" exponent="-1"/>
   </units>
   <component name="environment">
      <variable name="time" units="second" public_interface="out"/>
   </component>
   <component name="parameters">
      <variable name="isotonic" units="dimensionless" initial_value="0" public_interface="out"/>
      <variable name="alpha_1" units="per_micrometre" initial_value="14.6" public_interface="out"/>
      <variable name="beta_1" units="millinewton" initial_value="0.84" public_interface="out"/>
      <variable name="alpha_2" units="per_micrometre" initial_value="14.6" public_interface="out"/>
      <variable name="beta_2" units="millinewton" initial_value="0.0018" public_interface="out"/>
      <variable name="alpha_3" units="per_micrometre" initial_value="48" public_interface="out"/>
      <variable name="beta_3" units="millinewton" initial_value="0.015" public_interface="out"/>
      <variable name="lambda" units="millinewton" initial_value="30" public_interface="out"/>
      <variable name="mu" units="dimensionless" initial_value="3" public_interface="out"/>
      <variable name="k_mu" units="dimensionless" initial_value="0.6" public_interface="out"/>
      <variable name="kappa" units="dimensionless" initial_value="0.705" public_interface="out"/>
      <variable name="kappa_0" units="dimensionless" initial_value="3" public_interface="out"/>
      <variable name="m_0" units="dimensionless" initial_value="0.9" public_interface="out"/>
      <variable name="v_max" units="micrometre_per_second" initial_value="5.5" public_interface="out"/>
      <variable name="a" units="dimensionless" initial_value="0.25" public_interface="out"/>
      <variable name="d_h" units="dimensionless" initial_value="0.5" public_interface="out"/>
      <variable name="alpha_P" units="dimensionless" initial_value="4" public_interface="out"/>
   </component>
   <component name="isotonic">
      <variable name="isotonic" units="dimensionless" public_interface="in"/>
      <variable name="l" units="micrometre" public_interface="in"/>
      <variable name="F_muscle" units="millinewton" public_interface="in"/>
      <variable name="F_afterload" units="millinewton" initial_value="2"/>
      <variable name="isotonic_mode" units="dimensionless" public_interface="out"/>
      <variable name="l_0" units="micrometre" public_interface="in"/>
      <variable name="stim_period" units="second" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>isotonic_mode</ci>
            <piecewise>
               <piece>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <and/>
                     <apply>
                        <eq/>
                        <ci>isotonic</ci>
                        <cn cellml:units="dimensionless">1</cn>
                     </apply>
                     <apply>
                        <gt/>
                        <ci>F_muscle</ci>
                        <ci>F_afterload</ci>
                     </apply>
                     <apply>
                        <leq/>
                        <ci>l</ci>
                        <apply>
                           <times/>
                           <ci>l_0</ci>
                           <apply>
                              <plus/>
                              <cn cellml:units="dimensionless">1</cn>
                              <cn type="e-notation" cellml:units="dimensionless">1   <sep/>
                              -4</cn>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </piece>
               <otherwise>
                  <cn cellml:units="dimensionless">0</cn>
               </otherwise>
            </piecewise>
         </apply>
      </math>
   </component>
   <component name="parameters_izakov_et_al_1991">
      <variable name="S_0" units="micrometre" initial_value="1.14" public_interface="out"/>
      <variable name="q_v" units="per_second" public_interface="out"/>
      <variable name="q_1" units="per_second" initial_value="17.3"/>
      <variable name="q_2" units="per_second" initial_value="259"/>
      <variable name="q_3" units="per_second" initial_value="17.3"/>
      <variable name="q_4" units="per_second" initial_value="15" public_interface="out"/>
      <variable name="x_st" units="dimensionless" initial_value="0.964285"/>
      <variable name="v_st" units="micrometre_per_second" public_interface="out"/>
      <variable name="v_1" units="micrometre_per_second" public_interface="out"/>
      <variable name="alpha_G" units="dimensionless" initial_value="1" public_interface="out"/>
      <variable name="k_A" units="per_millimolar" initial_value="40" public_interface="out"/>
      <variable name="v" units="micrometre_per_second" public_interface="in"/>
      <variable name="v_max" units="micrometre_per_second" public_interface="in"/>
      <variable name="alpha_Q" units="dimensionless" initial_value="10"/>
      <variable name="beta_Q" units="dimensionless" initial_value="5"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>v_st</ci>
            <apply>
               <times/>
               <ci>x_st</ci>
               <ci>v_max</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>q_v</ci>
            <piecewise>
               <piece>
                  <apply>
                     <minus/>
                     <ci>q_1</ci>
                     <apply>
                        <divide/>
                        <apply>
                           <times/>
                           <ci>q_2</ci>
                           <ci>v</ci>
                        </apply>
                        <ci>v_max</ci>
                     </apply>
                  </apply>
                  <apply>
                     <leq/>
                     <ci>v</ci>
                     <cn cellml:units="micrometre_per_second">0</cn>
                  </apply>
               </piece>
               <piece>
                  <apply>
                     <plus/>
                     <apply>
                        <divide/>
                        <apply>
                           <times/>
                           <apply>
                              <minus/>
                              <ci>q_4</ci>
                              <ci>q_3</ci>
                           </apply>
                           <ci>v</ci>
                        </apply>
                        <ci>v_st</ci>
                     </apply>
                     <ci>q_3</ci>
                  </apply>
                  <apply>
                     <and/>
                     <apply>
                        <leq/>
                        <ci>v</ci>
                        <ci>v_st</ci>
                     </apply>
                     <apply>
                        <lt/>
                        <cn cellml:units="micrometre_per_second">0</cn>
                        <ci>v</ci>
                     </apply>
                  </apply>
               </piece>
               <otherwise>
                  <apply>
                     <divide/>
                     <ci>q_4</ci>
                     <apply>
                        <power/>
                        <apply>
                           <plus/>
                           <cn cellml:units="dimensionless">1</cn>
                           <apply>
                              <divide/>
                              <apply>
                                 <times/>
                                 <ci>beta_Q</ci>
                                 <apply>
                                    <minus/>
                                    <ci>v</ci>
                                    <ci>v_st</ci>
                                 </apply>
                              </apply>
                              <ci>v_max</ci>
                           </apply>
                        </apply>
                        <ci>alpha_Q</ci>
                     </apply>
                  </apply>
               </otherwise>
            </piecewise>
         </apply>
         <apply>
            <eq/>
            <ci>v_1</ci>
            <apply>
               <divide/>
               <ci>v_max</ci>
               <cn cellml:units="dimensionless">10</cn>
            </apply>
         </apply>
      </math>
   </component>
   <component name="force">
      <variable name="F_CE" units="millinewton"/>
      <variable name="F_muscle" units="millinewton" public_interface="out"/>
      <variable name="F_XSE" units="millinewton"/>
      <variable name="F_SE" units="millinewton"/>
      <variable name="F_PE" units="millinewton"/>
      <variable name="N" units="dimensionless" public_interface="in"/>
      <variable name="k_P_vis" units="millinewton_second_per_micrometre" public_interface="in"/>
      <variable name="k_S_vis" units="millinewton_second_per_micrometre" public_interface="in"/>
      <variable name="w" units="micrometre_per_second" public_interface="in"/>
      <variable name="v" units="micrometre_per_second" public_interface="in"/>
      <variable name="l_1" units="micrometre" public_interface="in"/>
      <variable name="l_2" units="micrometre" public_interface="in"/>
      <variable name="l_3" units="micrometre" public_interface="in"/>
      <variable name="p_v" units="dimensionless" public_interface="in"/>
      <variable name="alpha_1" units="per_micrometre" public_interface="in"/>
      <variable name="beta_1" units="millinewton" public_interface="in"/>
      <variable name="alpha_2" units="per_micrometre" public_interface="in"/>
      <variable name="beta_2" units="millinewton" public_interface="in"/>
      <variable name="alpha_3" units="per_micrometre" public_interface="in"/>
      <variable name="beta_3" units="millinewton" public_interface="in"/>
      <variable name="lambda" units="millinewton" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>F_CE</ci>
            <apply>
               <times/>
               <ci>lambda</ci>
               <ci>p_v</ci>
               <ci>N</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>F_muscle</ci>
            <ci>F_XSE</ci>
         </apply>
         <apply>
            <eq/>
            <ci>F_SE</ci>
            <apply>
               <times/>
               <ci>beta_1</ci>
               <apply>
                  <minus/>
                  <apply>
                     <exp/>
                     <apply>
                        <times/>
                        <ci>alpha_1</ci>
                        <apply>
                           <minus/>
                           <ci>l_2</ci>
                           <ci>l_1</ci>
                        </apply>
                     </apply>
                  </apply>
                  <cn cellml:units="dimensionless">1</cn>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>F_PE</ci>
            <apply>
               <times/>
               <ci>beta_2</ci>
               <apply>
                  <minus/>
                  <apply>
                     <exp/>
                     <apply>
                        <times/>
                        <ci>alpha_2</ci>
                        <ci>l_2</ci>
                     </apply>
                  </apply>
                  <cn cellml:units="dimensionless">1</cn>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>F_XSE</ci>
            <apply>
               <times/>
               <ci>beta_3</ci>
               <apply>
                  <minus/>
                  <apply>
                     <exp/>
                     <apply>
                        <times/>
                        <ci>alpha_3</ci>
                        <ci>l_3</ci>
                     </apply>
                  </apply>
                  <cn cellml:units="dimensionless">1</cn>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="crossbridge_kinetics">
      <variable name="N" units="dimensionless" initial_value="7.455e-8" public_interface="out"/>
      <variable name="K_kappa" units="per_second" public_interface="out"/>
      <variable name="M_A" units="dimensionless"/>
      <variable name="n_1" units="dimensionless"/>
      <variable name="L_oz" units="dimensionless" public_interface="out"/>
      <variable name="k_p_v" units="per_second"/>
      <variable name="k_m_v" units="per_second"/>
      <variable name="time" units="second" public_interface="in"/>
      <variable name="A" units="millimolar" public_interface="in"/>
      <variable name="l_1" units="micrometre" public_interface="in"/>
      <variable name="A_tot" units="millimolar" public_interface="in"/>
      <variable name="mu" units="dimensionless" public_interface="in"/>
      <variable name="k_mu" units="dimensionless" public_interface="in"/>
      <variable name="kappa" units="dimensionless" public_interface="in"/>
      <variable name="kappa_0" units="dimensionless" public_interface="in"/>
      <variable name="m_0" units="dimensionless" public_interface="in"/>
      <variable name="S_0" units="micrometre" public_interface="in"/>
      <variable name="q_v" units="per_second" public_interface="in"/>
      <variable name="v" units="micrometre_per_second" public_interface="in"/>
      <variable name="v_st" units="micrometre_per_second" public_interface="in"/>
      <variable name="G_star" units="dimensionless" public_interface="in"/>
      <variable name="P_star" units="dimensionless" public_interface="in"/>
      <variable name="q_4" units="per_second" public_interface="in"/>
      <variable name="g_1" units="per_micrometre" initial_value="0.6"/>
      <variable name="g_2" units="dimensionless" initial_value="0.52"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>M_A</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <apply>
                     <power/>
                     <apply>
                        <divide/>
                        <ci>A</ci>
                        <ci>A_tot</ci>
                     </apply>
                     <ci>mu</ci>
                  </apply>
                  <apply>
                     <plus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <apply>
                        <power/>
                        <ci>k_mu</ci>
                        <ci>mu</ci>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <plus/>
                  <apply>
                     <power/>
                     <apply>
                        <divide/>
                        <ci>A</ci>
                        <ci>A_tot</ci>
                     </apply>
                     <ci>mu</ci>
                  </apply>
                  <apply>
                     <power/>
                     <ci>k_mu</ci>
                     <ci>mu</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>n_1</ci>
            <piecewise>
               <piece>
                  <cn cellml:units="dimensionless">0</cn>
                  <apply>
                     <lt/>
                     <apply>
                        <plus/>
                        <apply>
                           <times/>
                           <ci>g_1</ci>
                           <ci>l_1</ci>
                        </apply>
                        <ci>g_2</ci>
                     </apply>
                     <cn cellml:units="dimensionless">0</cn>
                  </apply>
               </piece>
               <piece>
                  <apply>
                     <plus/>
                     <apply>
                        <times/>
                        <ci>g_1</ci>
                        <ci>l_1</ci>
                     </apply>
                     <ci>g_2</ci>
                  </apply>
                  <apply>
                     <lt/>
                     <apply>
                        <plus/>
                        <apply>
                           <times/>
                           <ci>g_1</ci>
                           <ci>l_1</ci>
                        </apply>
                        <ci>g_2</ci>
                     </apply>
                     <cn cellml:units="dimensionless">1</cn>
                  </apply>
               </piece>
               <otherwise>
                  <cn cellml:units="dimensionless">1</cn>
               </otherwise>
            </piecewise>
         </apply>
         <apply>
            <eq/>
            <ci>L_oz</ci>
            <piecewise>
               <piece>
                  <apply>
                     <divide/>
                     <apply>
                        <plus/>
                        <ci>l_1</ci>
                        <ci>S_0</ci>
                     </apply>
                     <apply>
                        <plus/>
                        <cn cellml:units="micrometre">0.46</cn>
                        <ci>S_0</ci>
                     </apply>
                  </apply>
                  <apply>
                     <gt/>
                     <ci>l_1</ci>
                     <cn cellml:units="micrometre">0.55</cn>
                  </apply>
               </piece>
               <otherwise>
                  <apply>
                     <times/>
                     <apply>
                        <plus/>
                        <ci>S_0</ci>
                        <cn cellml:units="micrometre">0.55</cn>
                     </apply>
                     <cn cellml:units="per_micrometre">1</cn>
                  </apply>
               </otherwise>
            </piecewise>
         </apply>
         <apply>
            <eq/>
            <ci>k_p_v</ci>
            <apply>
               <times/>
               <ci>kappa</ci>
               <ci>kappa_0</ci>
               <ci>q_v</ci>
               <ci>m_0</ci>
               <ci>G_star</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>k_m_v</ci>
            <apply>
               <times/>
               <ci>kappa_0</ci>
               <ci>q_v</ci>
               <apply>
                  <minus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <times/>
                     <ci>kappa</ci>
                     <ci>m_0</ci>
                     <ci>G_star</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>K_kappa</ci>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>k_p_v</ci>
                  <ci>M_A</ci>
                  <ci>n_1</ci>
                  <ci>L_oz</ci>
                  <apply>
                     <minus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <ci>N</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>k_m_v</ci>
                  <ci>N</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>N</ci>
            </apply>
            <ci>K_kappa</ci>
         </apply>
      </math>
   </component>
   <component name="length">
      <variable name="dl_1_dt" units="micrometre_per_second" public_interface="out"/>
      <variable name="l_1" units="micrometre" initial_value="0.436333342969918" public_interface="out"/>
      <variable name="dl_2_dt" units="micrometre_per_second" public_interface="out"/>
      <variable name="l_2" units="micrometre" initial_value="0.436333525334166" public_interface="out"/>
      <variable name="dl_3_dt" units="micrometre_per_second"/>
      <variable name="l_3" units="micrometre" initial_value="0.088805830771694" public_interface="out"/>
      <variable name="time" units="second" public_interface="in"/>
      <variable name="isotonic_mode" units="dimensionless" public_interface="in"/>
      <variable name="v" units="micrometre_per_second" public_interface="in"/>
      <variable name="w" units="micrometre_per_second" public_interface="in"/>
      <variable name="l" units="micrometre" public_interface="out"/>
      <variable name="k_S_vis" units="millinewton_second_per_micrometre" public_interface="in"/>
      <variable name="phi_chi_2" units="micrometre_per_second" public_interface="in"/>
      <variable name="stim_start" units="second" public_interface="in"/>
      <variable name="stim_end" units="second" public_interface="in"/>
      <variable name="stim_period" units="second" public_interface="in"/>
      <variable name="stim_duration" units="second" public_interface="in"/>
      <variable name="stim_amplitude" units="nanoA" public_interface="in"/>
      <variable name="l_0" units="micrometre" initial_value="0.525139356105856" public_interface="out"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>l</ci>
            <apply>
               <plus/>
               <ci>l_2</ci>
               <ci>l_3</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>dl_1_dt</ci>
            <ci>v</ci>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>l_1</ci>
            </apply>
            <ci>dl_1_dt</ci>
         </apply>
         <apply>
            <eq/>
            <ci>dl_2_dt</ci>
            <piecewise>
               <piece>
                  <ci>phi_chi_2</ci>
                  <apply>
                     <eq/>
                     <ci>k_S_vis</ci>
                     <cn cellml:units="millinewton_second_per_micrometre">0</cn>
                  </apply>
               </piece>
               <otherwise>
                  <ci>w</ci>
               </otherwise>
            </piecewise>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>l_2</ci>
            </apply>
            <ci>dl_2_dt</ci>
         </apply>
         <apply>
            <eq/>
            <ci>dl_3_dt</ci>
            <piecewise>
               <piece>
                  <cn cellml:units="micrometre_per_second">0</cn>
                  <apply>
                     <eq/>
                     <ci>isotonic_mode</ci>
                     <cn cellml:units="dimensionless">1</cn>
                  </apply>
               </piece>
               <piece>
                  <apply>
                     <minus/>
                     <ci>phi_chi_2</ci>
                  </apply>
                  <apply>
                     <and/>
                     <apply>
                        <eq/>
                        <ci>isotonic_mode</ci>
                        <cn cellml:units="dimensionless">0</cn>
                     </apply>
                     <apply>
                        <eq/>
                        <ci>k_S_vis</ci>
                        <cn cellml:units="millinewton_second_per_micrometre">0</cn>
                     </apply>
                  </apply>
               </piece>
               <piece>
                  <apply>
                     <minus/>
                     <ci>w</ci>
                  </apply>
                  <apply>
                     <and/>
                     <apply>
                        <eq/>
                        <ci>isotonic_mode</ci>
                        <cn cellml:units="dimensionless">0</cn>
                     </apply>
                     <apply>
                        <neq/>
                        <ci>k_S_vis</ci>
                        <cn cellml:units="millinewton_second_per_micrometre">0</cn>
                     </apply>
                  </apply>
               </piece>
            </piecewise>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>l_3</ci>
            </apply>
            <piecewise>
               <piece>
                  <apply>
                     <divide/>
                     <apply>
                        <minus/>
                        <apply>
                           <minus/>
                           <apply>
                              <plus/>
                              <ci>l_2</ci>
                              <ci>l_3</ci>
                           </apply>
                           <ci>l_0</ci>
                        </apply>
                     </apply>
                     <ci>stim_duration</ci>
                  </apply>
                  <apply>
                     <and/>
                     <apply>
                        <geq/>
                        <ci>time</ci>
                        <ci>stim_start</ci>
                     </apply>
                     <apply>
                        <leq/>
                        <ci>time</ci>
                        <ci>stim_end</ci>
                     </apply>
                     <apply>
                        <leq/>
                        <apply>
                           <minus/>
                           <apply>
                              <minus/>
                              <ci>time</ci>
                              <ci>stim_start</ci>
                           </apply>
                           <apply>
                              <times/>
                              <apply>
                                 <floor/>
                                 <apply>
                                    <divide/>
                                    <apply>
                                       <minus/>
                                       <ci>time</ci>
                                       <ci>stim_start</ci>
                                    </apply>
                                    <ci>stim_period</ci>
                                 </apply>
                              </apply>
                              <ci>stim_period</ci>
                           </apply>
                        </apply>
                        <ci>stim_duration</ci>
                     </apply>
                  </apply>
               </piece>
               <otherwise>
                  <ci>dl_3_dt</ci>
               </otherwise>
            </piecewise>
         </apply>
      </math>
   </component>
   <component name="CE_velocity">
      <variable name="v" units="micrometre_per_second" initial_value="0" public_interface="out"/>
      <variable name="phi_chi" units="micrometre_per_second2" public_interface="out"/>
      <variable name="phi_chi_2" units="micrometre_per_second" public_interface="out"/>
      <variable name="k_P_vis" units="millinewton_second_per_micrometre" public_interface="out"/>
      <variable name="k_S_vis" units="millinewton_second_per_micrometre" public_interface="in"/>
      <variable name="time" units="second" public_interface="in"/>
      <variable name="w" units="micrometre_per_second" public_interface="in"/>
      <variable name="K_kappa" units="per_second" public_interface="in"/>
      <variable name="N" units="dimensionless" public_interface="in"/>
      <variable name="l_1" units="micrometre" public_interface="in"/>
      <variable name="dl_1_dt" units="micrometre_per_second" public_interface="in"/>
      <variable name="l_2" units="micrometre" public_interface="in"/>
      <variable name="l_3" units="micrometre" public_interface="in"/>
      <variable name="alpha_1" units="per_micrometre" public_interface="in"/>
      <variable name="beta_1" units="millinewton" public_interface="in"/>
      <variable name="alpha_2" units="per_micrometre" public_interface="in"/>
      <variable name="beta_2" units="millinewton" public_interface="in"/>
      <variable name="alpha_3" units="per_micrometre" public_interface="in"/>
      <variable name="beta_3" units="millinewton" public_interface="in"/>
      <variable name="lambda" units="millinewton" public_interface="in"/>
      <variable name="isotonic_mode" units="dimensionless" public_interface="in"/>
      <variable name="p_v" units="dimensionless" public_interface="in"/>
      <variable name="p_prime_v" units="second_per_micrometre" public_interface="in"/>
      <variable name="alpha_P_lengthening" units="per_micrometre" initial_value="16"/>
      <variable name="beta_P_lengthening" units="millinewton_second_per_micrometre" initial_value="0.0015"/>
      <variable name="alpha_P_shortening" units="per_micrometre" initial_value="16"/>
      <variable name="beta_P_shortening" units="millinewton_second_per_micrometre" initial_value="0.0015"/>
      <variable name="alp_p" units="per_micrometre"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>alp_p</ci>
            <piecewise>
               <piece>
                  <ci>alpha_P_lengthening</ci>
                  <apply>
                     <leq/>
                     <ci>v</ci>
                     <cn cellml:units="micrometre_per_second">0</cn>
                  </apply>
               </piece>
               <otherwise>
                  <ci>alpha_P_shortening</ci>
               </otherwise>
            </piecewise>
         </apply>
         <apply>
            <eq/>
            <ci>k_P_vis</ci>
            <piecewise>
               <piece>
                  <apply>
                     <times/>
                     <ci>beta_P_lengthening</ci>
                     <apply>
                        <exp/>
                        <apply>
                           <times/>
                           <ci>alpha_P_lengthening</ci>
                           <ci>l_1</ci>
                        </apply>
                     </apply>
                  </apply>
                  <apply>
                     <leq/>
                     <ci>v</ci>
                     <cn cellml:units="micrometre_per_second">0</cn>
                  </apply>
               </piece>
               <otherwise>
                  <apply>
                     <times/>
                     <ci>beta_P_shortening</ci>
                     <apply>
                        <exp/>
                        <apply>
                           <times/>
                           <ci>alpha_P_shortening</ci>
                           <ci>l_1</ci>
                        </apply>
                     </apply>
                  </apply>
               </otherwise>
            </piecewise>
         </apply>
         <apply>
            <eq/>
            <ci>phi_chi</ci>
            <piecewise>
               <piece>
                  <apply>
                     <divide/>
                     <apply>
                        <minus/>
                        <apply>
                           <plus/>
                           <apply>
                              <times/>
                              <ci>lambda</ci>
                              <ci>K_kappa</ci>
                              <ci>p_v</ci>
                           </apply>
                           <apply>
                              <times/>
                              <ci>alp_p</ci>
                              <ci>k_P_vis</ci>
                              <apply>
                                 <power/>
                                 <ci>v</ci>
                                 <cn cellml:units="dimensionless">2</cn>
                              </apply>
                           </apply>
                           <apply>
                              <times/>
                              <ci>alpha_2</ci>
                              <ci>beta_2</ci>
                              <apply>
                                 <exp/>
                                 <apply>
                                    <times/>
                                    <ci>alpha_2</ci>
                                    <ci>l_2</ci>
                                 </apply>
                              </apply>
                              <ci>w</ci>
                           </apply>
                        </apply>
                     </apply>
                     <apply>
                        <plus/>
                        <apply>
                           <times/>
                           <ci>lambda</ci>
                           <ci>N</ci>
                           <ci>p_prime_v</ci>
                        </apply>
                        <ci>k_P_vis</ci>
                     </apply>
                  </apply>
                  <apply>
                     <eq/>
                     <ci>isotonic_mode</ci>
                     <cn cellml:units="dimensionless">1</cn>
                  </apply>
               </piece>
               <otherwise>
                  <apply>
                     <divide/>
                     <apply>
                        <minus/>
                        <apply>
                           <plus/>
                           <apply>
                              <times/>
                              <ci>lambda</ci>
                              <ci>K_kappa</ci>
                              <ci>p_v</ci>
                           </apply>
                           <apply>
                              <times/>
                              <ci>alp_p</ci>
                              <ci>k_P_vis</ci>
                              <apply>
                                 <power/>
                                 <ci>v</ci>
                                 <cn cellml:units="dimensionless">2</cn>
                              </apply>
                           </apply>
                           <apply>
                              <times/>
                              <apply>
                                 <plus/>
                                 <apply>
                                    <times/>
                                    <ci>alpha_2</ci>
                                    <ci>beta_2</ci>
                                    <apply>
                                       <exp/>
                                       <apply>
                                          <times/>
                                          <ci>alpha_2</ci>
                                          <ci>l_2</ci>
                                       </apply>
                                    </apply>
                                 </apply>
                                 <apply>
                                    <times/>
                                    <ci>alpha_3</ci>
                                    <ci>beta_3</ci>
                                    <apply>
                                       <exp/>
                                       <apply>
                                          <times/>
                                          <ci>alpha_3</ci>
                                          <ci>l_3</ci>
                                       </apply>
                                    </apply>
                                 </apply>
                              </apply>
                              <ci>w</ci>
                           </apply>
                        </apply>
                     </apply>
                     <apply>
                        <plus/>
                        <apply>
                           <times/>
                           <ci>lambda</ci>
                           <ci>N</ci>
                           <ci>p_prime_v</ci>
                        </apply>
                        <ci>k_P_vis</ci>
                     </apply>
                  </apply>
               </otherwise>
            </piecewise>
         </apply>
         <apply>
            <eq/>
            <ci>phi_chi_2</ci>
            <piecewise>
               <piece>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <ci>alpha_1</ci>
                        <ci>beta_1</ci>
                        <apply>
                           <exp/>
                           <apply>
                              <times/>
                              <ci>alpha_1</ci>
                              <apply>
                                 <minus/>
                                 <ci>l_2</ci>
                                 <ci>l_1</ci>
                              </apply>
                           </apply>
                        </apply>
                        <ci>v</ci>
                     </apply>
                     <apply>
                        <plus/>
                        <apply>
                           <times/>
                           <ci>alpha_1</ci>
                           <ci>beta_1</ci>
                           <apply>
                              <exp/>
                              <apply>
                                 <times/>
                                 <ci>alpha_1</ci>
                                 <apply>
                                    <minus/>
                                    <ci>l_2</ci>
                                    <ci>l_1</ci>
                                 </apply>
                              </apply>
                           </apply>
                        </apply>
                        <apply>
                           <times/>
                           <ci>alpha_2</ci>
                           <ci>beta_2</ci>
                           <apply>
                              <exp/>
                              <apply>
                                 <times/>
                                 <ci>alpha_2</ci>
                                 <ci>l_2</ci>
                              </apply>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
                  <apply>
                     <eq/>
                     <ci>isotonic_mode</ci>
                     <cn cellml:units="dimensionless">1</cn>
                  </apply>
               </piece>
               <otherwise>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <ci>alpha_1</ci>
                        <ci>beta_1</ci>
                        <apply>
                           <exp/>
                           <apply>
                              <times/>
                              <ci>alpha_1</ci>
                              <apply>
                                 <minus/>
                                 <ci>l_2</ci>
                                 <ci>l_1</ci>
                              </apply>
                           </apply>
                        </apply>
                        <ci>v</ci>
                     </apply>
                     <apply>
                        <plus/>
                        <apply>
                           <times/>
                           <ci>alpha_1</ci>
                           <ci>beta_1</ci>
                           <apply>
                              <exp/>
                              <apply>
                                 <times/>
                                 <ci>alpha_1</ci>
                                 <apply>
                                    <minus/>
                                    <ci>l_2</ci>
                                    <ci>l_1</ci>
                                 </apply>
                              </apply>
                           </apply>
                        </apply>
                        <apply>
                           <times/>
                           <ci>alpha_2</ci>
                           <ci>beta_2</ci>
                           <apply>
                              <exp/>
                              <apply>
                                 <times/>
                                 <ci>alpha_2</ci>
                                 <ci>l_2</ci>
                              </apply>
                           </apply>
                        </apply>
                        <apply>
                           <times/>
                           <ci>alpha_3</ci>
                           <ci>beta_3</ci>
                           <apply>
                              <exp/>
                              <apply>
                                 <times/>
                                 <ci>alpha_3</ci>
                                 <ci>l_3</ci>
                              </apply>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </otherwise>
            </piecewise>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>v</ci>
            </apply>
            <piecewise>
               <piece>
                  <apply>
                     <divide/>
                     <apply>
                        <minus/>
                        <apply>
                           <times/>
                           <ci>alpha_1</ci>
                           <ci>beta_1</ci>
                           <apply>
                              <exp/>
                              <apply>
                                 <times/>
                                 <ci>alpha_1</ci>
                                 <apply>
                                    <minus/>
                                    <ci>l_2</ci>
                                    <ci>l_1</ci>
                                 </apply>
                              </apply>
                           </apply>
                           <apply>
                              <minus/>
                              <ci>phi_chi_2</ci>
                              <ci>v</ci>
                           </apply>
                        </apply>
                        <apply>
                           <plus/>
                           <apply>
                              <times/>
                              <ci>lambda</ci>
                              <ci>K_kappa</ci>
                              <ci>p_v</ci>
                           </apply>
                           <apply>
                              <times/>
                              <ci>alp_p</ci>
                              <ci>k_P_vis</ci>
                              <apply>
                                 <power/>
                                 <ci>v</ci>
                                 <cn cellml:units="dimensionless">2</cn>
                              </apply>
                           </apply>
                        </apply>
                     </apply>
                     <apply>
                        <plus/>
                        <apply>
                           <times/>
                           <ci>lambda</ci>
                           <ci>N</ci>
                           <ci>p_prime_v</ci>
                        </apply>
                        <ci>k_P_vis</ci>
                     </apply>
                  </apply>
                  <apply>
                     <eq/>
                     <ci>k_S_vis</ci>
                     <cn cellml:units="millinewton_second_per_micrometre">0</cn>
                  </apply>
               </piece>
               <otherwise>
                  <ci>phi_chi</ci>
               </otherwise>
            </piecewise>
         </apply>
      </math>
   </component>
   <component name="PE_velocity">
      <variable name="w" units="micrometre_per_second" initial_value="0" public_interface="out"/>
      <variable name="k_S_vis" units="millinewton_second_per_micrometre" public_interface="out"/>
      <variable name="time" units="second" public_interface="in"/>
      <variable name="v" units="micrometre_per_second" public_interface="in"/>
      <variable name="phi_chi" units="micrometre_per_second2" public_interface="in"/>
      <variable name="l_1" units="micrometre" public_interface="in"/>
      <variable name="dl_1_dt" units="micrometre_per_second" public_interface="in"/>
      <variable name="l_2" units="micrometre" public_interface="in"/>
      <variable name="dl_2_dt" units="micrometre_per_second" public_interface="in"/>
      <variable name="l_3" units="micrometre" public_interface="in"/>
      <variable name="alpha_1" units="per_micrometre" public_interface="in"/>
      <variable name="beta_1" units="millinewton" public_interface="in"/>
      <variable name="alpha_2" units="per_micrometre" public_interface="in"/>
      <variable name="beta_2" units="millinewton" public_interface="in"/>
      <variable name="alpha_3" units="per_micrometre" public_interface="in"/>
      <variable name="beta_3" units="millinewton" public_interface="in"/>
      <variable name="lambda" units="millinewton" public_interface="in"/>
      <variable name="isotonic_mode" units="dimensionless" public_interface="in"/>
      <variable name="alpha_S_lengthening" units="per_micrometre" initial_value="46"/>
      <variable name="beta_S_lengthening" units="millinewton_second_per_micrometre" initial_value="0"/>
      <variable name="alpha_S_shortening" units="per_micrometre" initial_value="39"/>
      <variable name="beta_S_shortening" units="millinewton_second_per_micrometre" initial_value="0"/>
      <variable name="alp_s" units="per_micrometre"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>alp_s</ci>
            <piecewise>
               <piece>
                  <ci>alpha_S_lengthening</ci>
                  <apply>
                     <leq/>
                     <ci>w</ci>
                     <ci>v</ci>
                  </apply>
               </piece>
               <otherwise>
                  <ci>alpha_S_shortening</ci>
               </otherwise>
            </piecewise>
         </apply>
         <apply>
            <eq/>
            <ci>k_S_vis</ci>
            <piecewise>
               <piece>
                  <apply>
                     <times/>
                     <ci>beta_S_lengthening</ci>
                     <apply>
                        <exp/>
                        <apply>
                           <times/>
                           <ci>alpha_S_lengthening</ci>
                           <apply>
                              <minus/>
                              <ci>l_2</ci>
                              <ci>l_1</ci>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
                  <apply>
                     <leq/>
                     <ci>w</ci>
                     <ci>v</ci>
                  </apply>
               </piece>
               <otherwise>
                  <apply>
                     <times/>
                     <ci>beta_S_shortening</ci>
                     <apply>
                        <exp/>
                        <apply>
                           <times/>
                           <ci>alpha_S_shortening</ci>
                           <apply>
                              <minus/>
                              <ci>l_2</ci>
                              <ci>l_1</ci>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </otherwise>
            </piecewise>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>w</ci>
            </apply>
            <piecewise>
               <piece>
                  <apply>
                     <divide/>
                     <apply>
                        <minus/>
                        <apply>
                           <minus/>
                           <apply>
                              <times/>
                              <ci>k_S_vis</ci>
                              <apply>
                                 <minus/>
                                 <ci>phi_chi</ci>
                                 <apply>
                                    <times/>
                                    <ci>alp_s</ci>
                                    <apply>
                                       <power/>
                                       <apply>
                                          <minus/>
                                          <ci>w</ci>
                                          <ci>v</ci>
                                       </apply>
                                       <cn cellml:units="dimensionless">2</cn>
                                    </apply>
                                 </apply>
                              </apply>
                           </apply>
                           <apply>
                              <times/>
                              <ci>alpha_1</ci>
                              <ci>beta_1</ci>
                              <apply>
                                 <exp/>
                                 <apply>
                                    <times/>
                                    <ci>alpha_1</ci>
                                    <apply>
                                       <minus/>
                                       <ci>l_2</ci>
                                       <ci>l_1</ci>
                                    </apply>
                                 </apply>
                              </apply>
                              <apply>
                                 <minus/>
                                 <ci>w</ci>
                                 <ci>v</ci>
                              </apply>
                           </apply>
                        </apply>
                        <apply>
                           <times/>
                           <ci>alpha_2</ci>
                           <ci>beta_2</ci>
                           <apply>
                              <exp/>
                              <apply>
                                 <times/>
                                 <ci>alpha_2</ci>
                                 <ci>l_2</ci>
                              </apply>
                           </apply>
                           <ci>w</ci>
                        </apply>
                     </apply>
                     <ci>k_S_vis</ci>
                  </apply>
                  <apply>
                     <and/>
                     <apply>
                        <eq/>
                        <ci>isotonic_mode</ci>
                        <cn cellml:units="dimensionless">1</cn>
                     </apply>
                     <apply>
                        <neq/>
                        <ci>k_S_vis</ci>
                        <cn cellml:units="millinewton_second_per_micrometre">0</cn>
                     </apply>
                  </apply>
               </piece>
               <piece>
                  <apply>
                     <minus/>
                     <apply>
                        <minus/>
                        <ci>phi_chi</ci>
                        <apply>
                           <times/>
                           <ci>alp_s</ci>
                           <apply>
                              <power/>
                              <apply>
                                 <minus/>
                                 <ci>w</ci>
                                 <ci>v</ci>
                              </apply>
                              <cn cellml:units="dimensionless">2</cn>
                           </apply>
                        </apply>
                     </apply>
                     <apply>
                        <divide/>
                        <apply>
                           <plus/>
                           <apply>
                              <times/>
                              <ci>alpha_1</ci>
                              <ci>beta_1</ci>
                              <apply>
                                 <exp/>
                                 <apply>
                                    <times/>
                                    <ci>alpha_1</ci>
                                    <apply>
                                       <minus/>
                                       <ci>l_2</ci>
                                       <ci>l_1</ci>
                                    </apply>
                                 </apply>
                              </apply>
                              <apply>
                                 <minus/>
                                 <ci>w</ci>
                                 <ci>v</ci>
                              </apply>
                           </apply>
                           <apply>
                              <times/>
                              <apply>
                                 <plus/>
                                 <apply>
                                    <times/>
                                    <ci>alpha_2</ci>
                                    <ci>beta_2</ci>
                                    <apply>
                                       <exp/>
                                       <apply>
                                          <times/>
                                          <ci>alpha_2</ci>
                                          <ci>l_2</ci>
                                       </apply>
                                    </apply>
                                 </apply>
                                 <apply>
                                    <times/>
                                    <ci>alpha_3</ci>
                                    <ci>beta_3</ci>
                                    <apply>
                                       <exp/>
                                       <apply>
                                          <times/>
                                          <ci>alpha_3</ci>
                                          <ci>l_3</ci>
                                       </apply>
                                    </apply>
                                 </apply>
                              </apply>
                              <ci>w</ci>
                           </apply>
                        </apply>
                        <ci>k_S_vis</ci>
                     </apply>
                  </apply>
                  <apply>
                     <and/>
                     <apply>
                        <eq/>
                        <ci>isotonic_mode</ci>
                        <cn cellml:units="dimensionless">0</cn>
                     </apply>
                     <apply>
                        <neq/>
                        <ci>k_S_vis</ci>
                        <cn cellml:units="millinewton_second_per_micrometre">0</cn>
                     </apply>
                  </apply>
               </piece>
               <piece>
                  <cn cellml:units="micrometre_per_second2">0</cn>
                  <apply>
                     <eq/>
                     <ci>k_S_vis</ci>
                     <cn cellml:units="millinewton_second_per_micrometre">0</cn>
                  </apply>
               </piece>
            </piecewise>
         </apply>
      </math>
   </component>
   <component name="average_crossbridge_force">
      <variable name="p_v" units="dimensionless" public_interface="out"/>
      <variable name="p_prime_v" units="second_per_micrometre" public_interface="out"/>
      <variable name="P_star" units="dimensionless" public_interface="out"/>
      <variable name="G_star" units="dimensionless" public_interface="out"/>
      <variable name="gamma" units="dimensionless"/>
      <variable name="case_1" units="second_per_micrometre"/>
      <variable name="case_2" units="second_per_micrometre"/>
      <variable name="case_3" units="second_per_micrometre"/>
      <variable name="case_4" units="second_per_micrometre"/>
      <variable name="v" units="micrometre_per_second" public_interface="in"/>
      <variable name="phi_chi" units="micrometre_per_second2" public_interface="in"/>
      <variable name="v_1" units="micrometre_per_second" public_interface="in"/>
      <variable name="v_max" units="micrometre_per_second" public_interface="in"/>
      <variable name="a" units="dimensionless" public_interface="in"/>
      <variable name="d_h" units="dimensionless" public_interface="in"/>
      <variable name="alpha_G" units="dimensionless" public_interface="in"/>
      <variable name="alpha_P" units="dimensionless" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>gamma</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <ci>a</ci>
                  <ci>d_h</ci>
                  <apply>
                     <power/>
                     <apply>
                        <divide/>
                        <ci>v_1</ci>
                        <ci>v_max</ci>
                     </apply>
                     <cn cellml:units="dimensionless">2</cn>
                  </apply>
               </apply>
               <apply>
                  <minus/>
                  <apply>
                     <times/>
                     <cn cellml:units="dimensionless">3</cn>
                     <ci>a</ci>
                     <ci>d_h</ci>
                  </apply>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <apply>
                           <plus/>
                           <ci>a</ci>
                           <cn cellml:units="dimensionless">1</cn>
                        </apply>
                        <ci>v_1</ci>
                     </apply>
                     <ci>v_max</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>P_star</ci>
            <piecewise>
               <piece>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <ci>a</ci>
                        <apply>
                           <plus/>
                           <cn cellml:units="dimensionless">1</cn>
                           <apply>
                              <divide/>
                              <ci>v</ci>
                              <ci>v_max</ci>
                           </apply>
                        </apply>
                     </apply>
                     <apply>
                        <minus/>
                        <ci>a</ci>
                        <apply>
                           <divide/>
                           <ci>v</ci>
                           <ci>v_max</ci>
                        </apply>
                     </apply>
                  </apply>
                  <apply>
                     <leq/>
                     <ci>v</ci>
                     <cn cellml:units="micrometre_per_second">0</cn>
                  </apply>
               </piece>
               <otherwise>
                  <apply>
                     <minus/>
                     <apply>
                        <plus/>
                        <cn cellml:units="dimensionless">1</cn>
                        <ci>d_h</ci>
                     </apply>
                     <apply>
                        <divide/>
                        <apply>
                           <times/>
                           <apply>
                              <power/>
                              <ci>d_h</ci>
                              <cn cellml:units="dimensionless">2</cn>
                           </apply>
                           <ci>a</ci>
                        </apply>
                        <apply>
                           <plus/>
                           <apply>
                              <times/>
                              <apply>
                                 <divide/>
                                 <apply>
                                    <times/>
                                    <ci>a</ci>
                                    <ci>d_h</ci>
                                 </apply>
                                 <ci>gamma</ci>
                              </apply>
                              <apply>
                                 <power/>
                                 <apply>
                                    <divide/>
                                    <ci>v</ci>
                                    <ci>v_max</ci>
                                 </apply>
                                 <cn cellml:units="dimensionless">2</cn>
                              </apply>
                           </apply>
                           <apply>
                              <divide/>
                              <apply>
                                 <times/>
                                 <apply>
                                    <plus/>
                                    <ci>a</ci>
                                    <cn cellml:units="dimensionless">1</cn>
                                 </apply>
                                 <ci>v</ci>
                              </apply>
                              <ci>v_max</ci>
                           </apply>
                           <apply>
                              <times/>
                              <ci>a</ci>
                              <ci>d_h</ci>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </otherwise>
            </piecewise>
         </apply>
         <apply>
            <eq/>
            <ci>G_star</ci>
            <piecewise>
               <piece>
                  <apply>
                     <plus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <apply>
                        <divide/>
                        <apply>
                           <times/>
                           <cn cellml:units="dimensionless">0.6</cn>
                           <ci>v</ci>
                        </apply>
                        <ci>v_max</ci>
                     </apply>
                  </apply>
                  <apply>
                     <and/>
                     <apply>
                        <leq/>
                        <apply>
                           <minus/>
                           <ci>v_max</ci>
                        </apply>
                        <ci>v</ci>
                     </apply>
                     <apply>
                        <leq/>
                        <ci>v</ci>
                        <cn cellml:units="micrometre_per_second">0</cn>
                     </apply>
                  </apply>
               </piece>
               <piece>
                  <apply>
                     <divide/>
                     <ci>P_star</ci>
                     <apply>
                        <plus/>
                        <apply>
                           <divide/>
                           <apply>
                              <times/>
                              <apply>
                                 <plus/>
                                 <apply>
                                    <times/>
                                    <cn cellml:units="dimensionless">0.4</cn>
                                    <ci>a</ci>
                                 </apply>
                                 <cn cellml:units="dimensionless">1</cn>
                              </apply>
                              <ci>v</ci>
                           </apply>
                           <apply>
                              <times/>
                              <ci>a</ci>
                              <ci>v_max</ci>
                           </apply>
                        </apply>
                        <cn cellml:units="dimensionless">1</cn>
                     </apply>
                  </apply>
                  <apply>
                     <and/>
                     <apply>
                        <lt/>
                        <cn cellml:units="micrometre_per_second">0</cn>
                        <ci>v</ci>
                     </apply>
                     <apply>
                        <leq/>
                        <ci>v</ci>
                        <ci>v_1</ci>
                     </apply>
                  </apply>
               </piece>
               <otherwise>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <ci>P_star</ci>
                        <apply>
                           <exp/>
                           <apply>
                              <times/>
                              <apply>
                                 <minus/>
                                 <ci>alpha_G</ci>
                              </apply>
                              <apply>
                                 <power/>
                                 <apply>
                                    <divide/>
                                    <apply>
                                       <minus/>
                                       <ci>v</ci>
                                       <ci>v_1</ci>
                                    </apply>
                                    <ci>v_max</ci>
                                 </apply>
                                 <ci>alpha_P</ci>
                              </apply>
                           </apply>
                        </apply>
                     </apply>
                     <apply>
                        <plus/>
                        <apply>
                           <divide/>
                           <apply>
                              <times/>
                              <apply>
                                 <plus/>
                                 <apply>
                                    <times/>
                                    <cn cellml:units="dimensionless">0.4</cn>
                                    <ci>a</ci>
                                 </apply>
                                 <cn cellml:units="dimensionless">1</cn>
                              </apply>
                              <ci>v</ci>
                           </apply>
                           <apply>
                              <times/>
                              <ci>a</ci>
                              <ci>v_max</ci>
                           </apply>
                        </apply>
                        <cn cellml:units="dimensionless">1</cn>
                     </apply>
                  </apply>
               </otherwise>
            </piecewise>
         </apply>
         <apply>
            <eq/>
            <ci>case_1</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <ci>a</ci>
                  <apply>
                     <plus/>
                     <cn cellml:units="dimensionless">0.4</cn>
                     <apply>
                        <times/>
                        <cn cellml:units="dimensionless">0.4</cn>
                        <ci>a</ci>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>v_max</ci>
                  <apply>
                     <power/>
                     <apply>
                        <times/>
                        <apply>
                           <plus/>
                           <ci>a</ci>
                           <cn cellml:units="dimensionless">1</cn>
                        </apply>
                        <cn cellml:units="dimensionless">0.4</cn>
                     </apply>
                     <cn cellml:units="dimensionless">2</cn>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>case_2</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <ci>a</ci>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <plus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <apply>
                        <times/>
                        <cn cellml:units="dimensionless">0.4</cn>
                        <ci>a</ci>
                     </apply>
                     <apply>
                        <divide/>
                        <apply>
                           <times/>
                           <cn cellml:units="dimensionless">1.2</cn>
                           <ci>v</ci>
                        </apply>
                        <ci>v_max</ci>
                     </apply>
                     <apply>
                        <times/>
                        <cn cellml:units="dimensionless">0.6</cn>
                        <apply>
                           <power/>
                           <apply>
                              <divide/>
                              <ci>v</ci>
                              <ci>v_max</ci>
                           </apply>
                           <cn cellml:units="dimensionless">2</cn>
                        </apply>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>v_max</ci>
                  <apply>
                     <power/>
                     <apply>
                        <times/>
                        <apply>
                           <minus/>
                           <ci>a</ci>
                           <apply>
                              <divide/>
                              <ci>v</ci>
                              <ci>v_max</ci>
                           </apply>
                        </apply>
                        <apply>
                           <plus/>
                           <cn cellml:units="dimensionless">1</cn>
                           <apply>
                              <divide/>
                              <apply>
                                 <times/>
                                 <cn cellml:units="dimensionless">0.6</cn>
                                 <ci>v</ci>
                              </apply>
                              <ci>v_max</ci>
                           </apply>
                        </apply>
                     </apply>
                     <cn cellml:units="dimensionless">2</cn>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>case_3</ci>
            <apply>
               <divide/>
               <apply>
                  <plus/>
                  <apply>
                     <times/>
                     <cn cellml:units="dimensionless">0.4</cn>
                     <ci>a</ci>
                  </apply>
                  <cn cellml:units="dimensionless">1</cn>
               </apply>
               <apply>
                  <times/>
                  <ci>a</ci>
                  <ci>v_max</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>case_4</ci>
            <apply>
               <times/>
               <apply>
                  <divide/>
                  <cn cellml:units="dimensionless">1</cn>
                  <ci>v_max</ci>
               </apply>
               <apply>
                  <exp/>
                  <apply>
                     <times/>
                     <apply>
                        <minus/>
                        <ci>alpha_G</ci>
                     </apply>
                     <apply>
                        <power/>
                        <apply>
                           <minus/>
                           <apply>
                              <divide/>
                              <ci>v</ci>
                              <ci>v_max</ci>
                           </apply>
                           <apply>
                              <divide/>
                              <ci>v_1</ci>
                              <ci>v_max</ci>
                           </apply>
                        </apply>
                        <ci>alpha_P</ci>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <plus/>
                  <apply>
                     <divide/>
                     <apply>
                        <plus/>
                        <apply>
                           <times/>
                           <cn cellml:units="dimensionless">0.4</cn>
                           <ci>a</ci>
                        </apply>
                        <cn cellml:units="dimensionless">1</cn>
                     </apply>
                     <ci>a</ci>
                  </apply>
                  <apply>
                     <times/>
                     <ci>alpha_G</ci>
                     <ci>alpha_P</ci>
                     <apply>
                        <plus/>
                        <cn cellml:units="dimensionless">1</cn>
                        <apply>
                           <divide/>
                           <apply>
                              <times/>
                              <apply>
                                 <plus/>
                                 <apply>
                                    <times/>
                                    <cn cellml:units="dimensionless">0.4</cn>
                                    <ci>a</ci>
                                 </apply>
                                 <cn cellml:units="dimensionless">1</cn>
                              </apply>
                              <ci>v</ci>
                           </apply>
                           <apply>
                              <times/>
                              <ci>a</ci>
                              <ci>v_max</ci>
                           </apply>
                        </apply>
                     </apply>
                     <apply>
                        <power/>
                        <apply>
                           <minus/>
                           <apply>
                              <divide/>
                              <ci>v</ci>
                              <ci>v_max</ci>
                           </apply>
                           <apply>
                              <divide/>
                              <ci>v_1</ci>
                              <ci>v_max</ci>
                           </apply>
                        </apply>
                        <apply>
                           <minus/>
                           <ci>alpha_P</ci>
                           <cn cellml:units="dimensionless">1</cn>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>p_v</ci>
            <apply>
               <divide/>
               <ci>P_star</ci>
               <ci>G_star</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>p_prime_v</ci>
            <piecewise>
               <piece>
                  <ci>case_1</ci>
                  <apply>
                     <leq/>
                     <ci>v</ci>
                     <apply>
                        <minus/>
                        <ci>v_max</ci>
                     </apply>
                  </apply>
               </piece>
               <piece>
                  <ci>case_2</ci>
                  <apply>
                     <and/>
                     <apply>
                        <lt/>
                        <apply>
                           <minus/>
                           <ci>v_max</ci>
                        </apply>
                        <ci>v</ci>
                     </apply>
                     <apply>
                        <leq/>
                        <ci>v</ci>
                        <cn cellml:units="micrometre_per_second">0</cn>
                     </apply>
                  </apply>
               </piece>
               <piece>
                  <ci>case_3</ci>
                  <apply>
                     <and/>
                     <apply>
                        <lt/>
                        <cn cellml:units="micrometre_per_second">0</cn>
                        <ci>v</ci>
                     </apply>
                     <apply>
                        <leq/>
                        <ci>v</ci>
                        <ci>v_1</ci>
                     </apply>
                  </apply>
               </piece>
               <otherwise>
                  <ci>case_4</ci>
               </otherwise>
            </piecewise>
         </apply>
      </math>
   </component>
   <component name="membrane">
      <variable name="V" units="millivolt" initial_value="-93.658148" public_interface="out"/>
      <variable name="R" units="joule_per_kilomole_kelvin" initial_value="8314.472" public_interface="out"/>
      <variable name="T" units="kelvin" initial_value="310" public_interface="out"/>
      <variable name="F" units="coulomb_per_mole" initial_value="96485.3415" public_interface="out"/>
      <variable name="Cm" units="microF" initial_value="9.5e-5"/>
      <variable name="time" units="second" public_interface="in"/>
      <variable name="i_K1" units="nanoA" public_interface="in"/>
      <variable name="i_to" units="nanoA" public_interface="in"/>
      <variable name="i_Kr" units="nanoA" public_interface="in"/>
      <variable name="i_Ks" units="nanoA" public_interface="in"/>
      <variable name="i_Ca_L_K_cyt" units="nanoA" public_interface="in"/>
      <variable name="i_Ca_L_K_ds" units="nanoA" public_interface="in"/>
      <variable name="i_NaK" units="nanoA" public_interface="in"/>
      <variable name="i_Na" units="nanoA" public_interface="in"/>
      <variable name="i_b_Na" units="nanoA" public_interface="in"/>
      <variable name="i_p_Na" units="nanoA" public_interface="in"/>
      <variable name="i_Ca_L_Na_cyt" units="nanoA" public_interface="in"/>
      <variable name="i_Ca_L_Na_ds" units="nanoA" public_interface="in"/>
      <variable name="i_NaCa_cyt" units="nanoA" public_interface="in"/>
      <variable name="i_NaCa_ds" units="nanoA" public_interface="in"/>
      <variable name="i_Ca_L_Ca_cyt" units="nanoA" public_interface="in"/>
      <variable name="i_Ca_L_Ca_ds" units="nanoA" public_interface="in"/>
      <variable name="i_b_Ca" units="nanoA" public_interface="in"/>
      <variable name="i_Stim" units="nanoA"/>
      <variable name="stim_start" units="second" initial_value="0.06" public_interface="out"/>
      <variable name="stim_end" units="second" initial_value="10000" public_interface="out"/>
      <variable name="stim_period" units="second" initial_value="1" public_interface="out"/>
      <variable name="stim_duration" units="second" initial_value="0.0025" public_interface="out"/>
      <variable name="stim_amplitude" units="nanoA" initial_value="-3" public_interface="out"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_Stim</ci>
            <piecewise>
               <piece>
                  <ci>stim_amplitude</ci>
                  <apply>
                     <and/>
                     <apply>
                        <geq/>
                        <ci>time</ci>
                        <ci>stim_start</ci>
                     </apply>
                     <apply>
                        <leq/>
                        <ci>time</ci>
                        <ci>stim_end</ci>
                     </apply>
                     <apply>
                        <leq/>
                        <apply>
                           <minus/>
                           <apply>
                              <minus/>
                              <ci>time</ci>
                              <ci>stim_start</ci>
                           </apply>
                           <apply>
                              <times/>
                              <apply>
                                 <floor/>
                                 <apply>
                                    <divide/>
                                    <apply>
                                       <minus/>
                                       <ci>time</ci>
                                       <ci>stim_start</ci>
                                    </apply>
                                    <ci>stim_period</ci>
                                 </apply>
                              </apply>
                              <ci>stim_period</ci>
                           </apply>
                        </apply>
                        <ci>stim_duration</ci>
                     </apply>
                  </apply>
               </piece>
               <otherwise>
                  <cn cellml:units="nanoA">0</cn>
               </otherwise>
            </piecewise>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>V</ci>
            </apply>
            <apply>
               <times/>
               <apply>
                  <divide/>
                  <apply>
                     <minus/>
                     <cn cellml:units="dimensionless">1</cn>
                  </apply>
                  <ci>Cm</ci>
               </apply>
               <apply>
                  <plus/>
                  <ci>i_Stim</ci>
                  <ci>i_K1</ci>
                  <ci>i_to</ci>
                  <ci>i_Kr</ci>
                  <ci>i_Ks</ci>
                  <ci>i_NaK</ci>
                  <ci>i_Na</ci>
                  <ci>i_b_Na</ci>
                  <ci>i_p_Na</ci>
                  <ci>i_Ca_L_Na_cyt</ci>
                  <ci>i_Ca_L_Na_ds</ci>
                  <ci>i_NaCa_cyt</ci>
                  <ci>i_NaCa_ds</ci>
                  <ci>i_Ca_L_Ca_cyt</ci>
                  <ci>i_Ca_L_Ca_ds</ci>
                  <ci>i_Ca_L_K_cyt</ci>
                  <ci>i_Ca_L_K_ds</ci>
                  <ci>i_b_Ca</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="reversal_potentials">
      <variable name="E_Na" units="millivolt" public_interface="out"/>
      <variable name="E_K" units="millivolt" public_interface="out"/>
      <variable name="E_Ks" units="millivolt" public_interface="out"/>
      <variable name="E_Ca" units="millivolt" public_interface="out"/>
      <variable name="E_mh" units="millivolt" public_interface="out"/>
      <variable name="P_kna" units="dimensionless" initial_value="0.03"/>
      <variable name="K_o" units="millimolar" public_interface="in"/>
      <variable name="Na_o" units="millimolar" public_interface="in"/>
      <variable name="K_i" units="millimolar" public_interface="in"/>
      <variable name="Na_i" units="millimolar" public_interface="in"/>
      <variable name="Ca_o" units="millimolar" public_interface="in"/>
      <variable name="Ca_i" units="millimolar" public_interface="in"/>
      <variable name="R" units="joule_per_kilomole_kelvin" public_interface="in"/>
      <variable name="F" units="coulomb_per_mole" public_interface="in"/>
      <variable name="T" units="kelvin" public_interface="in"/>
      <variable name="time" units="second" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>E_Na</ci>
            <apply>
               <times/>
               <apply>
                  <divide/>
                  <apply>
                     <times/>
                     <ci>R</ci>
                     <ci>T</ci>
                  </apply>
                  <ci>F</ci>
               </apply>
               <apply>
                  <ln/>
                  <apply>
                     <divide/>
                     <ci>Na_o</ci>
                     <ci>Na_i</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>E_K</ci>
            <apply>
               <times/>
               <apply>
                  <divide/>
                  <apply>
                     <times/>
                     <ci>R</ci>
                     <ci>T</ci>
                  </apply>
                  <ci>F</ci>
               </apply>
               <apply>
                  <ln/>
                  <apply>
                     <divide/>
                     <ci>K_o</ci>
                     <ci>K_i</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>E_Ks</ci>
            <apply>
               <times/>
               <apply>
                  <divide/>
                  <apply>
                     <times/>
                     <ci>R</ci>
                     <ci>T</ci>
                  </apply>
                  <ci>F</ci>
               </apply>
               <apply>
                  <ln/>
                  <apply>
                     <divide/>
                     <apply>
                        <plus/>
                        <ci>K_o</ci>
                        <apply>
                           <times/>
                           <ci>P_kna</ci>
                           <ci>Na_o</ci>
                        </apply>
                     </apply>
                     <apply>
                        <plus/>
                        <ci>K_i</ci>
                        <apply>
                           <times/>
                           <ci>P_kna</ci>
                           <ci>Na_i</ci>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>E_Ca</ci>
            <apply>
               <times/>
               <apply>
                  <divide/>
                  <apply>
                     <times/>
                     <cn cellml:units="dimensionless">0.5</cn>
                     <ci>R</ci>
                     <ci>T</ci>
                  </apply>
                  <ci>F</ci>
               </apply>
               <apply>
                  <ln/>
                  <apply>
                     <divide/>
                     <ci>Ca_o</ci>
                     <ci>Ca_i</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>E_mh</ci>
            <apply>
               <times/>
               <apply>
                  <divide/>
                  <apply>
                     <times/>
                     <ci>R</ci>
                     <ci>T</ci>
                  </apply>
                  <ci>F</ci>
               </apply>
               <apply>
                  <ln/>
                  <apply>
                     <divide/>
                     <apply>
                        <plus/>
                        <ci>Na_o</ci>
                        <apply>
                           <times/>
                           <cn cellml:units="dimensionless">0.12</cn>
                           <ci>K_o</ci>
                        </apply>
                     </apply>
                     <apply>
                        <plus/>
                        <ci>Na_i</ci>
                        <apply>
                           <times/>
                           <cn cellml:units="dimensionless">0.12</cn>
                           <ci>K_i</ci>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="time_independent_potassium_current">
      <variable name="i_K1" units="nanoA" public_interface="out"/>
      <variable name="K_mk1" units="millimolar" initial_value="10" public_interface="out"/>
      <variable name="g_K1" units="microS" initial_value="0.5"/>
      <variable name="K_o" units="millimolar" public_interface="in"/>
      <variable name="V" units="millivolt" public_interface="in"/>
      <variable name="E_K" units="millivolt" public_interface="in"/>
      <variable name="R" units="joule_per_kilomole_kelvin" public_interface="in"/>
      <variable name="F" units="coulomb_per_mole" public_interface="in"/>
      <variable name="T" units="kelvin" public_interface="in"/>
      <variable name="time" units="second" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_K1</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <ci>g_K1</ci>
                        <ci>K_o</ci>
                     </apply>
                     <apply>
                        <plus/>
                        <ci>K_o</ci>
                        <ci>K_mk1</ci>
                     </apply>
                  </apply>
                  <apply>
                     <minus/>
                     <ci>V</ci>
                     <ci>E_K</ci>
                  </apply>
               </apply>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <times/>
                           <apply>
                              <minus/>
                              <apply>
                                 <minus/>
                                 <ci>V</ci>
                                 <ci>E_K</ci>
                              </apply>
                              <cn cellml:units="millivolt">10</cn>
                           </apply>
                           <ci>F</ci>
                           <cn cellml:units="dimensionless">1.25</cn>
                        </apply>
                        <apply>
                           <times/>
                           <ci>R</ci>
                           <ci>T</ci>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="rapid_delayed_rectifier_potassium_current">
      <variable name="i_Kr" units="nanoA" public_interface="out"/>
      <variable name="g_Kr1" units="microS" initial_value="0.0021"/>
      <variable name="g_Kr2" units="microS" initial_value="0.0013"/>
      <variable name="E_K" units="millivolt" public_interface="in"/>
      <variable name="time" units="second" public_interface="in" private_interface="out"/>
      <variable name="V" units="millivolt" public_interface="in" private_interface="out"/>
      <variable name="xr1" units="dimensionless" private_interface="in"/>
      <variable name="xr2" units="dimensionless" private_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_Kr</ci>
            <apply>
               <times/>
               <apply>
                  <divide/>
                  <apply>
                     <times/>
                     <apply>
                        <plus/>
                        <apply>
                           <times/>
                           <ci>g_Kr1</ci>
                           <ci>xr1</ci>
                        </apply>
                        <apply>
                           <times/>
                           <ci>g_Kr2</ci>
                           <ci>xr2</ci>
                        </apply>
                     </apply>
                     <cn cellml:units="dimensionless">1</cn>
                  </apply>
                  <apply>
                     <plus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <plus/>
                              <ci>V</ci>
                              <cn cellml:units="millivolt">9</cn>
                           </apply>
                           <cn cellml:units="millivolt">22.4</cn>
                        </apply>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>E_K</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="rapid_delayed_rectifier_potassium_current_xr1_gate">
      <variable name="xr1" units="dimensionless" initial_value="8.88859784542779e-6" public_interface="out"/>
      <variable name="V" units="millivolt" public_interface="in"/>
      <variable name="alpha_xr1" units="per_second"/>
      <variable name="beta_xr1" units="per_second"/>
      <variable name="time" units="second" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>alpha_xr1</ci>
            <apply>
               <divide/>
               <cn cellml:units="per_second">50</cn>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <minus/>
                           <apply>
                              <minus/>
                              <ci>V</ci>
                              <cn cellml:units="millivolt">5</cn>
                           </apply>
                        </apply>
                        <cn cellml:units="millivolt">9</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>beta_xr1</ci>
            <apply>
               <times/>
               <cn cellml:units="per_second">0.05</cn>
               <apply>
                  <exp/>
                  <apply>
                     <divide/>
                     <apply>
                        <minus/>
                        <apply>
                           <minus/>
                           <ci>V</ci>
                           <cn cellml:units="millivolt">20</cn>
                        </apply>
                     </apply>
                     <cn cellml:units="millivolt">15</cn>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>xr1</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>alpha_xr1</ci>
                  <apply>
                     <minus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <ci>xr1</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>beta_xr1</ci>
                  <ci>xr1</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="rapid_delayed_rectifier_potassium_current_xr2_gate">
      <variable name="xr2" units="dimensionless" initial_value="1.53745791069154e-7" public_interface="out"/>
      <variable name="V" units="millivolt" public_interface="in"/>
      <variable name="alpha_xr2" units="per_second"/>
      <variable name="beta_xr2" units="per_second"/>
      <variable name="time" units="second" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>alpha_xr2</ci>
            <apply>
               <divide/>
               <cn cellml:units="per_second">50</cn>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <minus/>
                           <apply>
                              <minus/>
                              <ci>V</ci>
                              <cn cellml:units="millivolt">5</cn>
                           </apply>
                        </apply>
                        <cn cellml:units="millivolt">9</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>beta_xr2</ci>
            <apply>
               <times/>
               <cn cellml:units="per_second">0.4</cn>
               <apply>
                  <exp/>
                  <apply>
                     <minus/>
                     <apply>
                        <power/>
                        <apply>
                           <divide/>
                           <apply>
                              <plus/>
                              <ci>V</ci>
                              <cn cellml:units="millivolt">30</cn>
                           </apply>
                           <cn cellml:units="millivolt">30</cn>
                        </apply>
                        <cn cellml:units="dimensionless">3</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>xr2</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>alpha_xr2</ci>
                  <apply>
                     <minus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <ci>xr2</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>beta_xr2</ci>
                  <ci>xr2</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="slow_delayed_rectifier_potassium_current">
      <variable name="i_Ks" units="nanoA" public_interface="out"/>
      <variable name="g_Ks" units="microS" initial_value="0.0026"/>
      <variable name="E_Ks" units="millivolt" public_interface="in"/>
      <variable name="time" units="second" public_interface="in" private_interface="out"/>
      <variable name="V" units="millivolt" public_interface="in" private_interface="out"/>
      <variable name="xs" units="dimensionless" private_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_Ks</ci>
            <apply>
               <times/>
               <ci>g_Ks</ci>
               <apply>
                  <power/>
                  <ci>xs</ci>
                  <cn cellml:units="dimensionless">2</cn>
               </apply>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>E_Ks</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="slow_delayed_rectifier_potassium_current_xs_gate">
      <variable name="xs" units="dimensionless" initial_value="0.001" public_interface="out"/>
      <variable name="V" units="millivolt" public_interface="in"/>
      <variable name="alpha_xs" units="per_second"/>
      <variable name="beta_xs" units="per_second"/>
      <variable name="time" units="second" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>alpha_xs</ci>
            <apply>
               <divide/>
               <cn cellml:units="per_second">14</cn>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <minus/>
                           <apply>
                              <minus/>
                              <ci>V</ci>
                              <cn cellml:units="millivolt">40</cn>
                           </apply>
                        </apply>
                        <cn cellml:units="millivolt">9</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>beta_xs</ci>
            <apply>
               <times/>
               <cn cellml:units="per_second">1</cn>
               <apply>
                  <exp/>
                  <apply>
                     <divide/>
                     <apply>
                        <minus/>
                        <ci>V</ci>
                     </apply>
                     <cn cellml:units="millivolt">45</cn>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>xs</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>alpha_xs</ci>
                  <apply>
                     <minus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <ci>xs</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>beta_xs</ci>
                  <ci>xs</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="sodium_activated_potassium_current">
      <variable name="i_KNa" units="nanoA" public_interface="out"/>
      <variable name="g_K_Na" units="microS" initial_value="0"/>
      <variable name="K_kna" units="millimolar" initial_value="20"/>
      <variable name="Na_i" units="millimolar" public_interface="in"/>
      <variable name="E_K" units="millivolt" public_interface="in"/>
      <variable name="V" units="millivolt" public_interface="in"/>
      <variable name="time" units="second" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_KNa</ci>
            <apply>
               <times/>
               <apply>
                  <divide/>
                  <apply>
                     <times/>
                     <ci>g_K_Na</ci>
                     <ci>Na_i</ci>
                  </apply>
                  <apply>
                     <plus/>
                     <ci>Na_i</ci>
                     <ci>K_kna</ci>
                  </apply>
               </apply>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>E_K</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="fast_sodium_current">
      <variable name="i_Na" units="nanoA" public_interface="out"/>
      <variable name="g_Na" units="microS" initial_value="2.5"/>
      <variable name="E_mh" units="millivolt" public_interface="in"/>
      <variable name="time" units="second" public_interface="in" private_interface="out"/>
      <variable name="V" units="millivolt" public_interface="in" private_interface="out"/>
      <variable name="m" units="dimensionless" private_interface="in"/>
      <variable name="h" units="dimensionless" private_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_Na</ci>
            <apply>
               <times/>
               <ci>g_Na</ci>
               <apply>
                  <power/>
                  <ci>m</ci>
                  <cn cellml:units="dimensionless">3</cn>
               </apply>
               <ci>h</ci>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>E_mh</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="fast_sodium_current_m_gate">
      <variable name="m" units="dimensionless" initial_value="0.0015" public_interface="out"/>
      <variable name="alpha_m" units="per_second"/>
      <variable name="beta_m" units="per_second"/>
      <variable name="V" units="millivolt" public_interface="in"/>
      <variable name="time" units="second" public_interface="in"/>
      <variable name="delta_m" units="millivolt" initial_value="1e-5"/>
      <variable name="E0_m" units="millivolt"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>E0_m</ci>
            <apply>
               <plus/>
               <ci>V</ci>
               <cn cellml:units="millivolt">41</cn>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>alpha_m</ci>
            <piecewise>
               <piece>
                  <cn cellml:units="per_second">2000</cn>
                  <apply>
                     <lt/>
                     <apply>
                        <abs/>
                        <ci>E0_m</ci>
                     </apply>
                     <ci>delta_m</ci>
                  </apply>
               </piece>
               <otherwise>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <cn cellml:units="per_millivolt_second">200</cn>
                        <ci>E0_m</ci>
                     </apply>
                     <apply>
                        <minus/>
                        <cn cellml:units="dimensionless">1</cn>
                        <apply>
                           <exp/>
                           <apply>
                              <times/>
                              <apply>
                                 <minus/>
                                 <cn cellml:units="per_millivolt">0.1</cn>
                              </apply>
                              <ci>E0_m</ci>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </otherwise>
            </piecewise>
         </apply>
         <apply>
            <eq/>
            <ci>beta_m</ci>
            <apply>
               <times/>
               <cn cellml:units="per_second">8000</cn>
               <apply>
                  <exp/>
                  <apply>
                     <times/>
                     <apply>
                        <minus/>
                        <cn cellml:units="per_millivolt">0.056</cn>
                     </apply>
                     <apply>
                        <plus/>
                        <ci>V</ci>
                        <cn cellml:units="millivolt">66</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>m</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>alpha_m</ci>
                  <apply>
                     <minus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <ci>m</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>beta_m</ci>
                  <ci>m</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="fast_sodium_current_h_gate">
      <variable name="h" units="dimensionless" initial_value="0.995" public_interface="out"/>
      <variable name="alpha_h" units="per_second"/>
      <variable name="beta_h" units="per_second"/>
      <variable name="V" units="millivolt" public_interface="in"/>
      <variable name="time" units="second" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>alpha_h</ci>
            <apply>
               <times/>
               <cn cellml:units="per_second">20</cn>
               <apply>
                  <exp/>
                  <apply>
                     <times/>
                     <apply>
                        <minus/>
                        <cn cellml:units="per_millivolt">0.125</cn>
                     </apply>
                     <apply>
                        <plus/>
                        <ci>V</ci>
                        <cn cellml:units="millivolt">75</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>beta_h</ci>
            <apply>
               <divide/>
               <cn cellml:units="per_second">2000</cn>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <times/>
                     <cn cellml:units="dimensionless">320</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <times/>
                           <apply>
                              <minus/>
                              <cn cellml:units="per_millivolt">0.1</cn>
                           </apply>
                           <apply>
                              <plus/>
                              <ci>V</ci>
                              <cn cellml:units="millivolt">75</cn>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>h</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>alpha_h</ci>
                  <apply>
                     <minus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <ci>h</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>beta_h</ci>
                  <ci>h</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="persistent_sodium_current">
      <variable name="i_p_Na" units="nanoA" public_interface="out"/>
      <variable name="g_pna" units="microS" initial_value="0.004"/>
      <variable name="E_Na" units="millivolt" public_interface="in"/>
      <variable name="time" units="second" public_interface="in"/>
      <variable name="V" units="millivolt" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_p_Na</ci>
            <apply>
               <times/>
               <apply>
                  <divide/>
                  <apply>
                     <times/>
                     <ci>g_pna</ci>
                     <cn cellml:units="dimensionless">1</cn>
                  </apply>
                  <apply>
                     <plus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <minus/>
                              <apply>
                                 <plus/>
                                 <ci>V</ci>
                                 <cn cellml:units="millivolt">52</cn>
                              </apply>
                           </apply>
                           <cn cellml:units="millivolt">8</cn>
                        </apply>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>E_Na</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="sodium_background_current">
      <variable name="i_b_Na" units="nanoA" public_interface="out"/>
      <variable name="g_bna" units="microS" initial_value="0.0006"/>
      <variable name="Na_i" units="millimolar" public_interface="in"/>
      <variable name="E_Na" units="millivolt" public_interface="in"/>
      <variable name="time" units="second" public_interface="in"/>
      <variable name="V" units="millivolt" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_b_Na</ci>
            <apply>
               <times/>
               <ci>g_bna</ci>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>E_Na</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="L_type_Ca_channel">
      <variable name="i_Ca_L_Ca_cyt" units="nanoA" public_interface="out"/>
      <variable name="i_Ca_L_K_cyt" units="nanoA" public_interface="out"/>
      <variable name="i_Ca_L_Na_cyt" units="nanoA" public_interface="out"/>
      <variable name="i_Ca_L_Ca_ds" units="nanoA" public_interface="out"/>
      <variable name="i_Ca_L_K_ds" units="nanoA" public_interface="out"/>
      <variable name="i_Ca_L_Na_ds" units="nanoA" public_interface="out"/>
      <variable name="i_Ca_L" units="nanoA" public_interface="out"/>
      <variable name="P_Ca_L" units="nanoA_per_millimolar" initial_value="0.1"/>
      <variable name="P_CaK" units="dimensionless" initial_value="0.002"/>
      <variable name="P_CaNa" units="dimensionless" initial_value="0.01"/>
      <variable name="Ca_o" units="millimolar" public_interface="in"/>
      <variable name="Ca_i" units="millimolar" public_interface="in" private_interface="out"/>
      <variable name="Ca_ds" units="millimolar" public_interface="in" private_interface="out"/>
      <variable name="K_o" units="millimolar" public_interface="in"/>
      <variable name="K_i" units="millimolar" public_interface="in"/>
      <variable name="Na_o" units="millimolar" public_interface="in"/>
      <variable name="Na_i" units="millimolar" public_interface="in"/>
      <variable name="R" units="joule_per_kilomole_kelvin" public_interface="in"/>
      <variable name="F" units="coulomb_per_mole" public_interface="in"/>
      <variable name="T" units="kelvin" public_interface="in"/>
      <variable name="time" units="second" public_interface="in" private_interface="out"/>
      <variable name="V" units="millivolt" public_interface="in" private_interface="out"/>
      <variable name="d" units="dimensionless" private_interface="in"/>
      <variable name="f" units="dimensionless" private_interface="in"/>
      <variable name="f2" units="dimensionless" private_interface="in"/>
      <variable name="f2ds" units="dimensionless" private_interface="in"/>
      <variable name="Km_f2" units="millimolar" initial_value="100000" private_interface="out"/>
      <variable name="Km_f2ds" units="millimolar" initial_value="0.001" private_interface="out"/>
      <variable name="R_decay" units="per_second" initial_value="20" public_interface="out" private_interface="out"/>
      <variable name="FrICa" units="dimensionless" initial_value="1"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_Ca_L_Ca_cyt</ci>
            <apply>
               <times/>
               <apply>
                  <divide/>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <apply>
                           <minus/>
                           <cn cellml:units="dimensionless">1</cn>
                           <ci>FrICa</ci>
                        </apply>
                        <cn cellml:units="dimensionless">4</cn>
                        <ci>P_Ca_L</ci>
                        <ci>d</ci>
                        <ci>f</ci>
                        <ci>f2</ci>
                        <apply>
                           <minus/>
                           <ci>V</ci>
                           <cn cellml:units="millivolt">50</cn>
                        </apply>
                        <ci>F</ci>
                     </apply>
                     <apply>
                        <times/>
                        <ci>R</ci>
                        <ci>T</ci>
                     </apply>
                  </apply>
                  <apply>
                     <minus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <times/>
                              <apply>
                                 <minus/>
                                 <apply>
                                    <minus/>
                                    <ci>V</ci>
                                    <cn cellml:units="millivolt">50</cn>
                                 </apply>
                              </apply>
                              <ci>F</ci>
                              <cn cellml:units="dimensionless">2</cn>
                           </apply>
                           <apply>
                              <times/>
                              <ci>R</ci>
                              <ci>T</ci>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <minus/>
                  <apply>
                     <times/>
                     <ci>Ca_i</ci>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <times/>
                              <cn cellml:units="millivolt">100</cn>
                              <ci>F</ci>
                           </apply>
                           <apply>
                              <times/>
                              <ci>R</ci>
                              <ci>T</ci>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
                  <apply>
                     <times/>
                     <ci>Ca_o</ci>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <times/>
                              <apply>
                                 <minus/>
                                 <apply>
                                    <minus/>
                                    <ci>V</ci>
                                    <cn cellml:units="millivolt">50</cn>
                                 </apply>
                              </apply>
                              <ci>F</ci>
                              <cn cellml:units="dimensionless">2</cn>
                           </apply>
                           <apply>
                              <times/>
                              <ci>R</ci>
                              <ci>T</ci>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>i_Ca_L_K_cyt</ci>
            <apply>
               <times/>
               <apply>
                  <divide/>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <apply>
                           <minus/>
                           <cn cellml:units="dimensionless">1</cn>
                           <ci>FrICa</ci>
                        </apply>
                        <ci>P_CaK</ci>
                        <ci>P_Ca_L</ci>
                        <ci>d</ci>
                        <ci>f</ci>
                        <ci>f2</ci>
                        <apply>
                           <minus/>
                           <ci>V</ci>
                           <cn cellml:units="millivolt">50</cn>
                        </apply>
                        <ci>F</ci>
                     </apply>
                     <apply>
                        <times/>
                        <ci>R</ci>
                        <ci>T</ci>
                     </apply>
                  </apply>
                  <apply>
                     <minus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <times/>
                              <apply>
                                 <minus/>
                                 <apply>
                                    <minus/>
                                    <ci>V</ci>
                                    <cn cellml:units="millivolt">50</cn>
                                 </apply>
                              </apply>
                              <ci>F</ci>
                           </apply>
                           <apply>
                              <times/>
                              <ci>R</ci>
                              <ci>T</ci>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <minus/>
                  <apply>
                     <times/>
                     <ci>K_i</ci>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <times/>
                              <cn cellml:units="millivolt">50</cn>
                              <ci>F</ci>
                           </apply>
                           <apply>
                              <times/>
                              <ci>R</ci>
                              <ci>T</ci>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
                  <apply>
                     <times/>
                     <ci>K_o</ci>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <times/>
                              <apply>
                                 <minus/>
                                 <apply>
                                    <minus/>
                                    <ci>V</ci>
                                    <cn cellml:units="millivolt">50</cn>
                                 </apply>
                              </apply>
                              <ci>F</ci>
                           </apply>
                           <apply>
                              <times/>
                              <ci>R</ci>
                              <ci>T</ci>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>i_Ca_L_Na_cyt</ci>
            <apply>
               <times/>
               <apply>
                  <divide/>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <apply>
                           <minus/>
                           <cn cellml:units="dimensionless">1</cn>
                           <ci>FrICa</ci>
                        </apply>
                        <ci>P_CaNa</ci>
                        <ci>P_Ca_L</ci>
                        <ci>d</ci>
                        <ci>f</ci>
                        <ci>f2</ci>
                        <apply>
                           <minus/>
                           <ci>V</ci>
                           <cn cellml:units="millivolt">50</cn>
                        </apply>
                        <ci>F</ci>
                     </apply>
                     <apply>
                        <times/>
                        <ci>R</ci>
                        <ci>T</ci>
                     </apply>
                  </apply>
                  <apply>
                     <minus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <times/>
                              <apply>
                                 <minus/>
                                 <apply>
                                    <minus/>
                                    <ci>V</ci>
                                    <cn cellml:units="millivolt">50</cn>
                                 </apply>
                              </apply>
                              <ci>F</ci>
                           </apply>
                           <apply>
                              <times/>
                              <ci>R</ci>
                              <ci>T</ci>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <minus/>
                  <apply>
                     <times/>
                     <ci>Na_i</ci>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <times/>
                              <cn cellml:units="millivolt">50</cn>
                              <ci>F</ci>
                           </apply>
                           <apply>
                              <times/>
                              <ci>R</ci>
                              <ci>T</ci>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
                  <apply>
                     <times/>
                     <ci>Na_o</ci>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <times/>
                              <apply>
                                 <minus/>
                                 <apply>
                                    <minus/>
                                    <ci>V</ci>
                                    <cn cellml:units="millivolt">50</cn>
                                 </apply>
                              </apply>
                              <ci>F</ci>
                           </apply>
                           <apply>
                              <times/>
                              <ci>R</ci>
                              <ci>T</ci>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>i_Ca_L_Ca_ds</ci>
            <apply>
               <times/>
               <apply>
                  <divide/>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <ci>FrICa</ci>
                        <cn cellml:units="dimensionless">4</cn>
                        <ci>P_Ca_L</ci>
                        <ci>d</ci>
                        <ci>f</ci>
                        <ci>f2ds</ci>
                        <apply>
                           <minus/>
                           <ci>V</ci>
                           <cn cellml:units="millivolt">50</cn>
                        </apply>
                        <ci>F</ci>
                     </apply>
                     <apply>
                        <times/>
                        <ci>R</ci>
                        <ci>T</ci>
                     </apply>
                  </apply>
                  <apply>
                     <minus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <times/>
                              <apply>
                                 <minus/>
                                 <apply>
                                    <minus/>
                                    <ci>V</ci>
                                    <cn cellml:units="millivolt">50</cn>
                                 </apply>
                              </apply>
                              <ci>F</ci>
                              <cn cellml:units="dimensionless">2</cn>
                           </apply>
                           <apply>
                              <times/>
                              <ci>R</ci>
                              <ci>T</ci>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <minus/>
                  <apply>
                     <times/>
                     <ci>Ca_i</ci>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <times/>
                              <cn cellml:units="millivolt">100</cn>
                              <ci>F</ci>
                           </apply>
                           <apply>
                              <times/>
                              <ci>R</ci>
                              <ci>T</ci>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
                  <apply>
                     <times/>
                     <ci>Ca_o</ci>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <times/>
                              <apply>
                                 <minus/>
                                 <apply>
                                    <minus/>
                                    <ci>V</ci>
                                    <cn cellml:units="millivolt">50</cn>
                                 </apply>
                              </apply>
                              <ci>F</ci>
                              <cn cellml:units="dimensionless">2</cn>
                           </apply>
                           <apply>
                              <times/>
                              <ci>R</ci>
                              <ci>T</ci>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>i_Ca_L_K_ds</ci>
            <apply>
               <times/>
               <apply>
                  <divide/>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <ci>FrICa</ci>
                        <ci>P_CaK</ci>
                        <ci>P_Ca_L</ci>
                        <ci>d</ci>
                        <ci>f</ci>
                        <ci>f2ds</ci>
                        <apply>
                           <minus/>
                           <ci>V</ci>
                           <cn cellml:units="millivolt">50</cn>
                        </apply>
                        <ci>F</ci>
                     </apply>
                     <apply>
                        <times/>
                        <ci>R</ci>
                        <ci>T</ci>
                     </apply>
                  </apply>
                  <apply>
                     <minus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <times/>
                              <apply>
                                 <minus/>
                                 <apply>
                                    <minus/>
                                    <ci>V</ci>
                                    <cn cellml:units="millivolt">50</cn>
                                 </apply>
                              </apply>
                              <ci>F</ci>
                           </apply>
                           <apply>
                              <times/>
                              <ci>R</ci>
                              <ci>T</ci>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <minus/>
                  <apply>
                     <times/>
                     <ci>K_i</ci>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <times/>
                              <cn cellml:units="millivolt">50</cn>
                              <ci>F</ci>
                           </apply>
                           <apply>
                              <times/>
                              <ci>R</ci>
                              <ci>T</ci>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
                  <apply>
                     <times/>
                     <ci>K_o</ci>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <times/>
                              <apply>
                                 <minus/>
                                 <apply>
                                    <minus/>
                                    <ci>V</ci>
                                    <cn cellml:units="millivolt">50</cn>
                                 </apply>
                              </apply>
                              <ci>F</ci>
                           </apply>
                           <apply>
                              <times/>
                              <ci>R</ci>
                              <ci>T</ci>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>i_Ca_L_Na_ds</ci>
            <apply>
               <times/>
               <apply>
                  <divide/>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <ci>FrICa</ci>
                        <ci>P_CaNa</ci>
                        <ci>P_Ca_L</ci>
                        <ci>d</ci>
                        <ci>f</ci>
                        <ci>f2ds</ci>
                        <apply>
                           <minus/>
                           <ci>V</ci>
                           <cn cellml:units="millivolt">50</cn>
                        </apply>
                        <ci>F</ci>
                     </apply>
                     <apply>
                        <times/>
                        <ci>R</ci>
                        <ci>T</ci>
                     </apply>
                  </apply>
                  <apply>
                     <minus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <times/>
                              <apply>
                                 <minus/>
                                 <apply>
                                    <minus/>
                                    <ci>V</ci>
                                    <cn cellml:units="millivolt">50</cn>
                                 </apply>
                              </apply>
                              <ci>F</ci>
                           </apply>
                           <apply>
                              <times/>
                              <ci>R</ci>
                              <ci>T</ci>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <minus/>
                  <apply>
                     <times/>
                     <ci>Na_i</ci>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <times/>
                              <cn cellml:units="millivolt">50</cn>
                              <ci>F</ci>
                           </apply>
                           <apply>
                              <times/>
                              <ci>R</ci>
                              <ci>T</ci>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
                  <apply>
                     <times/>
                     <ci>Na_o</ci>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <times/>
                              <apply>
                                 <minus/>
                                 <apply>
                                    <minus/>
                                    <ci>V</ci>
                                    <cn cellml:units="millivolt">50</cn>
                                 </apply>
                              </apply>
                              <ci>F</ci>
                           </apply>
                           <apply>
                              <times/>
                              <ci>R</ci>
                              <ci>T</ci>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>i_Ca_L</ci>
            <apply>
               <plus/>
               <ci>i_Ca_L_Ca_cyt</ci>
               <ci>i_Ca_L_K_cyt</ci>
               <ci>i_Ca_L_Na_cyt</ci>
               <ci>i_Ca_L_Ca_ds</ci>
               <ci>i_Ca_L_K_ds</ci>
               <ci>i_Ca_L_Na_ds</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="L_type_Ca_channel_d_gate">
      <variable name="d" units="dimensionless" initial_value="0" public_interface="out"/>
      <variable name="alpha_d" units="per_second"/>
      <variable name="beta_d" units="per_second"/>
      <variable name="V" units="millivolt" public_interface="in"/>
      <variable name="time" units="second" public_interface="in"/>
      <variable name="E0_d" units="millivolt"/>
      <variable name="speed_d" units="dimensionless" initial_value="3"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>E0_d</ci>
            <apply>
               <minus/>
               <apply>
                  <plus/>
                  <ci>V</ci>
                  <cn cellml:units="millivolt">24</cn>
               </apply>
               <cn cellml:units="millivolt">5</cn>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>alpha_d</ci>
            <piecewise>
               <piece>
                  <cn cellml:units="per_second">120</cn>
                  <apply>
                     <lt/>
                     <apply>
                        <abs/>
                        <ci>E0_d</ci>
                     </apply>
                     <cn cellml:units="millivolt">0.0001</cn>
                  </apply>
               </piece>
               <otherwise>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <cn cellml:units="per_millivolt_second">30</cn>
                        <ci>E0_d</ci>
                     </apply>
                     <apply>
                        <minus/>
                        <cn cellml:units="dimensionless">1</cn>
                        <apply>
                           <exp/>
                           <apply>
                              <divide/>
                              <apply>
                                 <minus/>
                                 <ci>E0_d</ci>
                              </apply>
                              <cn cellml:units="millivolt">4</cn>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </otherwise>
            </piecewise>
         </apply>
         <apply>
            <eq/>
            <ci>beta_d</ci>
            <piecewise>
               <piece>
                  <cn cellml:units="per_second">120</cn>
                  <apply>
                     <lt/>
                     <apply>
                        <abs/>
                        <ci>E0_d</ci>
                     </apply>
                     <cn cellml:units="millivolt">0.0001</cn>
                  </apply>
               </piece>
               <otherwise>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <cn cellml:units="per_millivolt_second">12</cn>
                        <ci>E0_d</ci>
                     </apply>
                     <apply>
                        <minus/>
                        <apply>
                           <exp/>
                           <apply>
                              <divide/>
                              <ci>E0_d</ci>
                              <cn cellml:units="millivolt">10</cn>
                           </apply>
                        </apply>
                        <cn cellml:units="dimensionless">1</cn>
                     </apply>
                  </apply>
               </otherwise>
            </piecewise>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>d</ci>
            </apply>
            <apply>
               <times/>
               <ci>speed_d</ci>
               <apply>
                  <minus/>
                  <apply>
                     <times/>
                     <ci>alpha_d</ci>
                     <apply>
                        <minus/>
                        <cn cellml:units="dimensionless">1</cn>
                        <ci>d</ci>
                     </apply>
                  </apply>
                  <apply>
                     <times/>
                     <ci>beta_d</ci>
                     <ci>d</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="L_type_Ca_channel_f_gate">
      <variable name="f" units="dimensionless" initial_value="1" public_interface="out"/>
      <variable name="alpha_f" units="per_second"/>
      <variable name="beta_f" units="per_second"/>
      <variable name="V" units="millivolt" public_interface="in"/>
      <variable name="time" units="second" public_interface="in"/>
      <variable name="speed_f" units="dimensionless" initial_value="0.3"/>
      <variable name="delta_f" units="millivolt" initial_value="0.0001"/>
      <variable name="E0_f" units="millivolt"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>E0_f</ci>
            <apply>
               <plus/>
               <ci>V</ci>
               <cn cellml:units="millivolt">34</cn>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>alpha_f</ci>
            <piecewise>
               <piece>
                  <cn cellml:units="per_second">25</cn>
                  <apply>
                     <lt/>
                     <apply>
                        <abs/>
                        <ci>E0_f</ci>
                     </apply>
                     <ci>delta_f</ci>
                  </apply>
               </piece>
               <otherwise>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <cn cellml:units="per_millivolt_second">6.25</cn>
                        <ci>E0_f</ci>
                     </apply>
                     <apply>
                        <minus/>
                        <apply>
                           <exp/>
                           <apply>
                              <divide/>
                              <ci>E0_f</ci>
                              <cn cellml:units="millivolt">4</cn>
                           </apply>
                        </apply>
                        <cn cellml:units="dimensionless">1</cn>
                     </apply>
                  </apply>
               </otherwise>
            </piecewise>
         </apply>
         <apply>
            <eq/>
            <ci>beta_f</ci>
            <apply>
               <divide/>
               <cn cellml:units="per_second">12</cn>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <times/>
                           <apply>
                              <minus/>
                              <cn cellml:units="dimensionless">1</cn>
                           </apply>
                           <apply>
                              <plus/>
                              <ci>V</ci>
                              <cn cellml:units="millivolt">34</cn>
                           </apply>
                        </apply>
                        <cn cellml:units="millivolt">4</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>f</ci>
            </apply>
            <apply>
               <times/>
               <ci>speed_f</ci>
               <apply>
                  <minus/>
                  <apply>
                     <times/>
                     <ci>alpha_f</ci>
                     <apply>
                        <minus/>
                        <cn cellml:units="dimensionless">1</cn>
                        <ci>f</ci>
                     </apply>
                  </apply>
                  <apply>
                     <times/>
                     <ci>beta_f</ci>
                     <ci>f</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="L_type_Ca_channel_f2_gate">
      <variable name="f2" units="dimensionless" initial_value="1" public_interface="out"/>
      <variable name="Km_f2" units="millimolar" public_interface="in"/>
      <variable name="Ca_i" units="millimolar" public_interface="in"/>
      <variable name="V" units="millivolt" public_interface="in"/>
      <variable name="time" units="second" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>f2</ci>
            </apply>
            <apply>
               <minus/>
               <cn cellml:units="per_second">1</cn>
               <apply>
                  <times/>
                  <cn cellml:units="per_second">1</cn>
                  <apply>
                     <plus/>
                     <apply>
                        <divide/>
                        <ci>Ca_i</ci>
                        <apply>
                           <plus/>
                           <ci>Km_f2</ci>
                           <ci>Ca_i</ci>
                        </apply>
                     </apply>
                     <ci>f2</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="L_type_Ca_channel_f2ds_gate">
      <variable name="f2ds" units="dimensionless" initial_value="0.997" public_interface="out"/>
      <variable name="Km_f2ds" units="millimolar" public_interface="in"/>
      <variable name="R_decay" units="per_second" public_interface="in"/>
      <variable name="Ca_ds" units="millimolar" public_interface="in"/>
      <variable name="V" units="millivolt" public_interface="in"/>
      <variable name="time" units="second" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>f2ds</ci>
            </apply>
            <apply>
               <times/>
               <ci>R_decay</ci>
               <apply>
                  <minus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <plus/>
                     <apply>
                        <divide/>
                        <ci>Ca_ds</ci>
                        <apply>
                           <plus/>
                           <ci>Km_f2ds</ci>
                           <ci>Ca_ds</ci>
                        </apply>
                     </apply>
                     <ci>f2ds</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="calcium_background_current">
      <variable name="i_b_Ca" units="nanoA" public_interface="out"/>
      <variable name="g_bca" units="microS" initial_value="0.00025"/>
      <variable name="E_Ca" units="millivolt" public_interface="in"/>
      <variable name="time" units="second" public_interface="in"/>
      <variable name="V" units="millivolt" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_b_Ca</ci>
            <apply>
               <times/>
               <ci>g_bca</ci>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>E_Ca</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="transient_outward_current">
      <variable name="i_to" units="nanoA" public_interface="out"/>
      <variable name="g_to" units="microS" initial_value="0.006"/>
      <variable name="g_tos" units="dimensionless" initial_value="0"/>
      <variable name="E_K" units="millivolt" public_interface="in"/>
      <variable name="V" units="millivolt" public_interface="in" private_interface="out"/>
      <variable name="time" units="second" public_interface="in" private_interface="out"/>
      <variable name="s" units="dimensionless" private_interface="in"/>
      <variable name="r" units="dimensionless" private_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_to</ci>
            <apply>
               <times/>
               <ci>g_to</ci>
               <apply>
                  <plus/>
                  <ci>g_tos</ci>
                  <apply>
                     <times/>
                     <ci>s</ci>
                     <apply>
                        <minus/>
                        <cn cellml:units="dimensionless">1</cn>
                        <ci>g_tos</ci>
                     </apply>
                  </apply>
               </apply>
               <ci>r</ci>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>E_K</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="transient_outward_current_s_gate">
      <variable name="s" units="dimensionless" initial_value="0.997044616031121" public_interface="out"/>
      <variable name="alpha_s" units="per_second"/>
      <variable name="beta_s" units="per_second"/>
      <variable name="V" units="millivolt" public_interface="in"/>
      <variable name="time" units="second" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>alpha_s</ci>
            <apply>
               <times/>
               <cn cellml:units="per_second">0.033</cn>
               <apply>
                  <exp/>
                  <apply>
                     <divide/>
                     <apply>
                        <minus/>
                        <ci>V</ci>
                     </apply>
                     <cn cellml:units="millivolt">17</cn>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>beta_s</ci>
            <apply>
               <divide/>
               <cn cellml:units="per_second">33</cn>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <times/>
                        <apply>
                           <minus/>
                           <cn cellml:units="per_millivolt">0.125</cn>
                        </apply>
                        <apply>
                           <plus/>
                           <ci>V</ci>
                           <cn cellml:units="millivolt">10</cn>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>s</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>alpha_s</ci>
                  <apply>
                     <minus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <ci>s</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>beta_s</ci>
                  <ci>s</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="transient_outward_current_r_gate">
      <variable name="r" units="dimensionless" initial_value="1.63117173173398e-8" public_interface="out"/>
      <variable name="V" units="millivolt" public_interface="in"/>
      <variable name="time" units="second" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>r</ci>
            </apply>
            <apply>
               <times/>
               <cn cellml:units="per_second">333</cn>
               <apply>
                  <minus/>
                  <apply>
                     <divide/>
                     <cn cellml:units="dimensionless">1</cn>
                     <apply>
                        <plus/>
                        <cn cellml:units="dimensionless">1</cn>
                        <apply>
                           <exp/>
                           <apply>
                              <divide/>
                              <apply>
                                 <minus/>
                                 <apply>
                                    <plus/>
                                    <ci>V</ci>
                                    <cn cellml:units="millivolt">4</cn>
                                 </apply>
                              </apply>
                              <cn cellml:units="millivolt">5</cn>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
                  <ci>r</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="sodium_potassium_pump">
      <variable name="i_NaK" units="nanoA" public_interface="out"/>
      <variable name="i_NaK_max" units="nanoA" initial_value="0.7"/>
      <variable name="K_mK" units="millimolar" initial_value="1"/>
      <variable name="K_mNa" units="millimolar" initial_value="24.2"/>
      <variable name="K_o" units="millimolar" public_interface="in"/>
      <variable name="Na_i" units="millimolar" public_interface="in"/>
      <variable name="time" units="second" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_NaK</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <ci>i_NaK_max</ci>
                        <ci>K_o</ci>
                     </apply>
                     <apply>
                        <plus/>
                        <ci>K_mK</ci>
                        <ci>K_o</ci>
                     </apply>
                  </apply>
                  <ci>Na_i</ci>
               </apply>
               <apply>
                  <plus/>
                  <ci>K_mNa</ci>
                  <ci>Na_i</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="sodium_calcium_exchanger">
      <variable name="i_NaCa" units="nanoA" public_interface="out"/>
      <variable name="i_NaCa_cyt" units="nanoA" public_interface="out"/>
      <variable name="i_NaCa_ds" units="nanoA" public_interface="out"/>
      <variable name="k_NaCa" units="nanoA" initial_value="0.0005"/>
      <variable name="n_NaCa" units="dimensionless" initial_value="3" public_interface="out"/>
      <variable name="d_NaCa" units="dimensionless" initial_value="0"/>
      <variable name="gamma" units="dimensionless" initial_value="0.5"/>
      <variable name="FRiNaCa" units="dimensionless" initial_value="0.001"/>
      <variable name="Na_i" units="millimolar" public_interface="in"/>
      <variable name="Na_o" units="millimolar" public_interface="in"/>
      <variable name="Ca_i" units="millimolar" public_interface="in"/>
      <variable name="Ca_ds" units="millimolar" public_interface="in"/>
      <variable name="Ca_o" units="millimolar" public_interface="in"/>
      <variable name="R" units="joule_per_kilomole_kelvin" public_interface="in"/>
      <variable name="F" units="coulomb_per_mole" public_interface="in"/>
      <variable name="T" units="kelvin" public_interface="in"/>
      <variable name="V" units="millivolt" public_interface="in"/>
      <variable name="time" units="second" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_NaCa_cyt</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <apply>
                     <minus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <ci>FRiNaCa</ci>
                  </apply>
                  <ci>k_NaCa</ci>
                  <apply>
                     <minus/>
                     <apply>
                        <times/>
                        <apply>
                           <exp/>
                           <apply>
                              <divide/>
                              <apply>
                                 <times/>
                                 <ci>gamma</ci>
                                 <apply>
                                    <minus/>
                                    <ci>n_NaCa</ci>
                                    <cn cellml:units="dimensionless">2</cn>
                                 </apply>
                                 <ci>V</ci>
                                 <ci>F</ci>
                              </apply>
                              <apply>
                                 <times/>
                                 <ci>R</ci>
                                 <ci>T</ci>
                              </apply>
                           </apply>
                        </apply>
                        <apply>
                           <power/>
                           <ci>Na_i</ci>
                           <ci>n_NaCa</ci>
                        </apply>
                        <ci>Ca_o</ci>
                     </apply>
                     <apply>
                        <times/>
                        <apply>
                           <exp/>
                           <apply>
                              <divide/>
                              <apply>
                                 <times/>
                                 <apply>
                                    <minus/>
                                    <ci>gamma</ci>
                                    <cn cellml:units="dimensionless">1</cn>
                                 </apply>
                                 <apply>
                                    <minus/>
                                    <ci>n_NaCa</ci>
                                    <cn cellml:units="dimensionless">2</cn>
                                 </apply>
                                 <ci>V</ci>
                                 <ci>F</ci>
                              </apply>
                              <apply>
                                 <times/>
                                 <ci>R</ci>
                                 <ci>T</ci>
                              </apply>
                           </apply>
                        </apply>
                        <apply>
                           <power/>
                           <ci>Na_o</ci>
                           <ci>n_NaCa</ci>
                        </apply>
                        <ci>Ca_i</ci>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <apply>
                     <plus/>
                     <cn cellml:units="millimolar4">1</cn>
                     <apply>
                        <times/>
                        <ci>d_NaCa</ci>
                        <apply>
                           <plus/>
                           <apply>
                              <times/>
                              <ci>Ca_i</ci>
                              <apply>
                                 <power/>
                                 <ci>Na_o</ci>
                                 <ci>n_NaCa</ci>
                              </apply>
                           </apply>
                           <apply>
                              <times/>
                              <ci>Ca_o</ci>
                              <apply>
                                 <power/>
                                 <ci>Na_i</ci>
                                 <ci>n_NaCa</ci>
                              </apply>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
                  <apply>
                     <plus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <apply>
                        <divide/>
                        <ci>Ca_i</ci>
                        <cn cellml:units="millimolar">0.0069</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>i_NaCa_ds</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <ci>FRiNaCa</ci>
                  <ci>k_NaCa</ci>
                  <apply>
                     <minus/>
                     <apply>
                        <times/>
                        <apply>
                           <exp/>
                           <apply>
                              <divide/>
                              <apply>
                                 <times/>
                                 <ci>gamma</ci>
                                 <apply>
                                    <minus/>
                                    <ci>n_NaCa</ci>
                                    <cn cellml:units="dimensionless">2</cn>
                                 </apply>
                                 <ci>V</ci>
                                 <ci>F</ci>
                              </apply>
                              <apply>
                                 <times/>
                                 <ci>R</ci>
                                 <ci>T</ci>
                              </apply>
                           </apply>
                        </apply>
                        <apply>
                           <power/>
                           <ci>Na_i</ci>
                           <ci>n_NaCa</ci>
                        </apply>
                        <ci>Ca_o</ci>
                     </apply>
                     <apply>
                        <times/>
                        <apply>
                           <exp/>
                           <apply>
                              <divide/>
                              <apply>
                                 <times/>
                                 <apply>
                                    <minus/>
                                    <ci>gamma</ci>
                                    <cn cellml:units="dimensionless">1</cn>
                                 </apply>
                                 <apply>
                                    <minus/>
                                    <ci>n_NaCa</ci>
                                    <cn cellml:units="dimensionless">2</cn>
                                 </apply>
                                 <ci>V</ci>
                                 <ci>F</ci>
                              </apply>
                              <apply>
                                 <times/>
                                 <ci>R</ci>
                                 <ci>T</ci>
                              </apply>
                           </apply>
                        </apply>
                        <apply>
                           <power/>
                           <ci>Na_o</ci>
                           <ci>n_NaCa</ci>
                        </apply>
                        <ci>Ca_ds</ci>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <apply>
                     <plus/>
                     <cn cellml:units="millimolar4">1</cn>
                     <apply>
                        <times/>
                        <ci>d_NaCa</ci>
                        <apply>
                           <plus/>
                           <apply>
                              <times/>
                              <ci>Ca_ds</ci>
                              <apply>
                                 <power/>
                                 <ci>Na_o</ci>
                                 <ci>n_NaCa</ci>
                              </apply>
                           </apply>
                           <apply>
                              <times/>
                              <ci>Ca_o</ci>
                              <apply>
                                 <power/>
                                 <ci>Na_i</ci>
                                 <ci>n_NaCa</ci>
                              </apply>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
                  <apply>
                     <plus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <apply>
                        <divide/>
                        <ci>Ca_ds</ci>
                        <cn cellml:units="millimolar">0.0069</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>i_NaCa</ci>
            <apply>
               <plus/>
               <ci>i_NaCa_cyt</ci>
               <ci>i_NaCa_ds</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="sarcoplasmic_reticulum_calcium_pump">
      <variable name="i_up" units="millimolar_per_second" public_interface="out"/>
      <variable name="K_1" units="dimensionless"/>
      <variable name="K_2" units="millimolar"/>
      <variable name="K_cyca" units="millimolar" initial_value="0.00015"/>
      <variable name="K_xcs" units="dimensionless" initial_value="0.4"/>
      <variable name="K_srca" units="millimolar" initial_value="0.5"/>
      <variable name="alpha_up" units="millimolar_per_second" initial_value="1"/>
      <variable name="beta_up" units="millimolar_per_second" initial_value="0.03"/>
      <variable name="Ca_i" units="millimolar" public_interface="in"/>
      <variable name="Ca_up" units="millimolar" public_interface="in"/>
      <variable name="time" units="second" public_interface="in"/>
      <variable name="flag_ingib" units="dimensionless" initial_value="0"/>
      <variable name="K_inh" units="millimolar" initial_value="4"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>K_1</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <ci>K_cyca</ci>
                  <ci>K_xcs</ci>
               </apply>
               <ci>K_srca</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>K_2</ci>
            <apply>
               <plus/>
               <ci>Ca_i</ci>
               <apply>
                  <times/>
                  <ci>Ca_up</ci>
                  <ci>K_1</ci>
               </apply>
               <apply>
                  <times/>
                  <ci>K_cyca</ci>
                  <ci>K_xcs</ci>
               </apply>
               <ci>K_cyca</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>i_up</ci>
            <piecewise>
               <piece>
                  <apply>
                     <minus/>
                     <apply>
                        <times/>
                        <apply>
                           <divide/>
                           <ci>Ca_i</ci>
                           <ci>K_2</ci>
                        </apply>
                        <ci>alpha_up</ci>
                     </apply>
                     <apply>
                        <times/>
                        <apply>
                           <divide/>
                           <apply>
                              <times/>
                              <ci>Ca_up</ci>
                              <ci>K_1</ci>
                           </apply>
                           <ci>K_2</ci>
                        </apply>
                        <ci>beta_up</ci>
                     </apply>
                  </apply>
                  <apply>
                     <eq/>
                     <ci>flag_ingib</ci>
                     <cn cellml:units="dimensionless">0</cn>
                  </apply>
               </piece>
               <otherwise>
                  <apply>
                     <minus/>
                     <apply>
                        <divide/>
                        <apply>
                           <times/>
                           <apply>
                              <divide/>
                              <ci>Ca_i</ci>
                              <ci>K_2</ci>
                           </apply>
                           <ci>alpha_up</ci>
                        </apply>
                        <apply>
                           <plus/>
                           <cn cellml:units="dimensionless">1</cn>
                           <apply>
                              <divide/>
                              <ci>Ca_up</ci>
                              <ci>K_inh</ci>
                           </apply>
                        </apply>
                     </apply>
                     <apply>
                        <times/>
                        <apply>
                           <divide/>
                           <apply>
                              <times/>
                              <ci>Ca_up</ci>
                              <ci>K_1</ci>
                           </apply>
                           <ci>K_2</ci>
                        </apply>
                        <ci>beta_up</ci>
                     </apply>
                  </apply>
               </otherwise>
            </piecewise>
         </apply>
      </math>
   </component>
   <component name="calcium_translocation">
      <variable name="i_trans" units="millimolar_per_second" public_interface="out"/>
      <variable name="Ca_rel" units="millimolar" public_interface="in"/>
      <variable name="Ca_up" units="millimolar" public_interface="in"/>
      <variable name="time" units="second" public_interface="in"/>
      <variable name="a_tr" units="per_second" initial_value="15"/>
      <variable name="alpha_CaS" units="per_millimolar_second" initial_value="50000"/>
      <variable name="beta_CaS" units="per_second" initial_value="32500"/>
      <variable name="CaS_tot" units="millimolar" initial_value="40" public_interface="out"/>
      <variable name="beta" units="millimolar" public_interface="out"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_trans</ci>
            <apply>
               <times/>
               <ci>a_tr</ci>
               <apply>
                  <minus/>
                  <ci>Ca_up</ci>
                  <ci>Ca_rel</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>beta</ci>
            <apply>
               <divide/>
               <ci>beta_CaS</ci>
               <ci>alpha_CaS</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="calcium_release">
      <variable name="i_rel" units="millimolar_per_second" public_interface="out"/>
      <variable name="VoltDep" units="dimensionless"/>
      <variable name="RegBindSite" units="dimensionless"/>
      <variable name="CaiReg" units="dimensionless"/>
      <variable name="CadsReg" units="dimensionless"/>
      <variable name="ActRate" units="per_second"/>
      <variable name="InactRate" units="per_second"/>
      <variable name="SRLeak" units="per_second" initial_value="0.05"/>
      <variable name="K_m_rel" units="per_second" initial_value="10000"/>
      <variable name="K_m_Ca_cyt" units="millimolar" initial_value="0.0005"/>
      <variable name="K_m_Ca_ds" units="millimolar" initial_value="0.01"/>
      <variable name="PrecFrac" units="dimensionless"/>
      <variable name="ActFrac" units="dimensionless" initial_value="0.001914"/>
      <variable name="ProdFrac" units="dimensionless" initial_value="0.2854569"/>
      <variable name="SpeedRel" units="dimensionless"/>
      <variable name="time" units="second" public_interface="in"/>
      <variable name="V" units="millivolt" public_interface="in"/>
      <variable name="Ca_i" units="millimolar" public_interface="in"/>
      <variable name="Ca_ds" units="millimolar" public_interface="in"/>
      <variable name="Ca_rel" units="millimolar" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>VoltDep</ci>
            <apply>
               <exp/>
               <apply>
                  <times/>
                  <cn cellml:units="per_millivolt">0.08</cn>
                  <apply>
                     <minus/>
                     <ci>V</ci>
                     <cn cellml:units="millivolt">40</cn>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>CaiReg</ci>
            <apply>
               <divide/>
               <ci>Ca_i</ci>
               <apply>
                  <plus/>
                  <ci>Ca_i</ci>
                  <ci>K_m_Ca_cyt</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>CadsReg</ci>
            <apply>
               <divide/>
               <ci>Ca_ds</ci>
               <apply>
                  <plus/>
                  <ci>Ca_ds</ci>
                  <ci>K_m_Ca_ds</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>RegBindSite</ci>
            <apply>
               <plus/>
               <ci>CaiReg</ci>
               <apply>
                  <times/>
                  <apply>
                     <minus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <ci>CaiReg</ci>
                  </apply>
                  <ci>CadsReg</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>ActRate</ci>
            <apply>
               <times/>
               <cn cellml:units="per_second">500</cn>
               <apply>
                  <power/>
                  <ci>RegBindSite</ci>
                  <cn cellml:units="dimensionless">2</cn>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>InactRate</ci>
            <apply>
               <plus/>
               <cn cellml:units="per_second">60</cn>
               <apply>
                  <times/>
                  <cn cellml:units="per_second">500</cn>
                  <apply>
                     <power/>
                     <ci>RegBindSite</ci>
                     <cn cellml:units="dimensionless">2</cn>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>SpeedRel</ci>
            <piecewise>
               <piece>
                  <cn cellml:units="dimensionless">5</cn>
                  <apply>
                     <lt/>
                     <ci>V</ci>
                     <apply>
                        <minus/>
                        <cn cellml:units="millivolt">50</cn>
                     </apply>
                  </apply>
               </piece>
               <otherwise>
                  <cn cellml:units="dimensionless">1</cn>
               </otherwise>
            </piecewise>
         </apply>
         <apply>
            <eq/>
            <ci>PrecFrac</ci>
            <apply>
               <minus/>
               <apply>
                  <minus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <ci>ActFrac</ci>
               </apply>
               <ci>ProdFrac</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>ActFrac</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>PrecFrac</ci>
                  <ci>SpeedRel</ci>
                  <ci>ActRate</ci>
               </apply>
               <apply>
                  <times/>
                  <ci>ActFrac</ci>
                  <ci>SpeedRel</ci>
                  <ci>InactRate</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>ProdFrac</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>ActFrac</ci>
                  <ci>SpeedRel</ci>
                  <ci>InactRate</ci>
               </apply>
               <apply>
                  <times/>
                  <ci>SpeedRel</ci>
                  <cn cellml:units="per_second">0.6</cn>
                  <ci>ProdFrac</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>i_rel</ci>
            <apply>
               <times/>
               <apply>
                  <plus/>
                  <apply>
                     <times/>
                     <apply>
                        <power/>
                        <apply>
                           <divide/>
                           <ci>ActFrac</ci>
                           <apply>
                              <plus/>
                              <ci>ActFrac</ci>
                              <cn cellml:units="dimensionless">0.25</cn>
                           </apply>
                        </apply>
                        <cn cellml:units="dimensionless">2</cn>
                     </apply>
                     <ci>K_m_rel</ci>
                  </apply>
                  <ci>SRLeak</ci>
               </apply>
               <ci>Ca_rel</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="extracellular_sodium_concentration">
      <variable name="Na_o" units="millimolar" initial_value="140" public_interface="out"/>
   </component>
   <component name="intracellular_sodium_concentration">
      <variable name="Na_i" units="millimolar" initial_value="5.18787513289509" public_interface="out"/>
      <variable name="V_i" units="micrometre3" public_interface="in"/>
      <variable name="time" units="second" public_interface="in"/>
      <variable name="F" units="coulomb_per_mole" public_interface="in"/>
      <variable name="i_Na" units="nanoA" public_interface="in"/>
      <variable name="i_b_Na" units="nanoA" public_interface="in"/>
      <variable name="i_p_Na" units="nanoA" public_interface="in"/>
      <variable name="i_Ca_L_Na_cyt" units="nanoA" public_interface="in"/>
      <variable name="i_Ca_L_Na_ds" units="nanoA" public_interface="in"/>
      <variable name="i_NaK" units="nanoA" public_interface="in"/>
      <variable name="i_NaCa_cyt" units="nanoA" public_interface="in"/>
      <variable name="n_NaK" units="dimensionless" initial_value="1.5" public_interface="out"/>
      <variable name="n_NaCa" units="dimensionless" public_interface="in"/>
      <variable name="i_NaCa" units="nanoA" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>Na_i</ci>
            </apply>
            <apply>
               <times/>
               <apply>
                  <divide/>
                  <apply>
                     <minus/>
                     <cn cellml:units="mA_nA">1</cn>
                  </apply>
                  <apply>
                     <times/>
                     <cn cellml:units="litre_micrometre3">1</cn>
                     <ci>V_i</ci>
                     <ci>F</ci>
                  </apply>
               </apply>
               <apply>
                  <plus/>
                  <ci>i_Na</ci>
                  <ci>i_p_Na</ci>
                  <ci>i_b_Na</ci>
                  <ci>i_Ca_L_Na_cyt</ci>
                  <ci>i_Ca_L_Na_ds</ci>
                  <apply>
                     <times/>
                     <apply>
                        <divide/>
                        <ci>n_NaK</ci>
                        <apply>
                           <minus/>
                           <ci>n_NaK</ci>
                           <cn cellml:units="dimensionless">1</cn>
                        </apply>
                     </apply>
                     <ci>i_NaK</ci>
                  </apply>
                  <apply>
                     <times/>
                     <apply>
                        <divide/>
                        <ci>n_NaCa</ci>
                        <apply>
                           <minus/>
                           <ci>n_NaCa</ci>
                           <cn cellml:units="dimensionless">2</cn>
                        </apply>
                     </apply>
                     <ci>i_NaCa</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="extracellular_calcium_concentration">
      <variable name="Ca_o" units="millimolar" initial_value="2" public_interface="out"/>
   </component>
   <component name="extracellular_potassium_concentration">
      <variable name="K_o" units="millimolar" initial_value="3.988" public_interface="out"/>
      <variable name="K_b" units="millimolar" initial_value="4"/>
      <variable name="i_Kr" units="nanoA" public_interface="in"/>
      <variable name="i_Ks" units="nanoA" public_interface="in"/>
      <variable name="i_K1" units="nanoA" public_interface="in"/>
      <variable name="i_to" units="nanoA" public_interface="in"/>
      <variable name="i_NaK" units="nanoA" public_interface="in"/>
      <variable name="i_Ca_L_K_cyt" units="nanoA" public_interface="in"/>
      <variable name="i_Ca_L_K_ds" units="nanoA" public_interface="in"/>
      <variable name="time" units="second" public_interface="in"/>
      <variable name="n_NaK" units="dimensionless" public_interface="in"/>
      <variable name="F" units="coulomb_per_mole" public_interface="in"/>
      <variable name="V_e" units="micrometre3" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>K_o</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <divide/>
                  <apply>
                     <times/>
                     <cn cellml:units="mA_nA">1</cn>
                     <apply>
                        <plus/>
                        <apply>
                           <minus/>
                           <apply>
                              <plus/>
                              <ci>i_Kr</ci>
                              <ci>i_Ks</ci>
                              <ci>i_K1</ci>
                              <ci>i_to</ci>
                           </apply>
                           <apply>
                              <times/>
                              <apply>
                                 <divide/>
                                 <cn cellml:units="dimensionless">1</cn>
                                 <apply>
                                    <minus/>
                                    <ci>n_NaK</ci>
                                    <cn cellml:units="dimensionless">1</cn>
                                 </apply>
                              </apply>
                              <ci>i_NaK</ci>
                           </apply>
                        </apply>
                        <ci>i_Ca_L_K_cyt</ci>
                        <ci>i_Ca_L_K_ds</ci>
                     </apply>
                  </apply>
                  <apply>
                     <times/>
                     <cn cellml:units="litre_micrometre3">1</cn>
                     <ci>V_e</ci>
                     <ci>F</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <cn cellml:units="per_second">0.7</cn>
                  <apply>
                     <minus/>
                     <ci>K_o</ci>
                     <ci>K_b</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="intracellular_potassium_concentration">
      <variable name="K_i" units="millimolar" initial_value="139.054" public_interface="out"/>
      <variable name="V_i" units="micrometre3" public_interface="in"/>
      <variable name="i_K1" units="nanoA" public_interface="in"/>
      <variable name="i_Kr" units="nanoA" public_interface="in"/>
      <variable name="i_Ks" units="nanoA" public_interface="in"/>
      <variable name="i_Ca_L_K_cyt" units="nanoA" public_interface="in"/>
      <variable name="i_Ca_L_K_ds" units="nanoA" public_interface="in"/>
      <variable name="i_to" units="nanoA" public_interface="in"/>
      <variable name="i_NaK" units="nanoA" public_interface="in"/>
      <variable name="time" units="second" public_interface="in"/>
      <variable name="F" units="coulomb_per_mole" public_interface="in"/>
      <variable name="n_NaK" units="dimensionless" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>K_i</ci>
            </apply>
            <apply>
               <times/>
               <apply>
                  <divide/>
                  <apply>
                     <minus/>
                     <cn cellml:units="mA_nA">1</cn>
                  </apply>
                  <apply>
                     <times/>
                     <cn cellml:units="litre_micrometre3">1</cn>
                     <ci>V_i</ci>
                     <ci>F</ci>
                  </apply>
               </apply>
               <apply>
                  <minus/>
                  <apply>
                     <plus/>
                     <ci>i_K1</ci>
                     <ci>i_Kr</ci>
                     <ci>i_Ks</ci>
                     <ci>i_Ca_L_K_cyt</ci>
                     <ci>i_Ca_L_K_ds</ci>
                     <ci>i_to</ci>
                  </apply>
                  <apply>
                     <times/>
                     <apply>
                        <divide/>
                        <cn cellml:units="dimensionless">1</cn>
                        <apply>
                           <minus/>
                           <ci>n_NaK</ci>
                           <cn cellml:units="dimensionless">1</cn>
                        </apply>
                     </apply>
                     <ci>i_NaK</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="intracellular_calcium_concentration">
      <variable name="Ca_i" units="millimolar" initial_value="6.15e-6" public_interface="out"/>
      <variable name="Ca_ds" units="millimolar" initial_value="2.55e-6" public_interface="out"/>
      <variable name="Ca_up" units="millimolar" initial_value="0.994579" public_interface="out"/>
      <variable name="Ca_rel" units="millimolar" initial_value="0.989665" public_interface="out"/>
      <variable name="radius" units="micrometre" initial_value="12"/>
      <variable name="length" units="micrometre" initial_value="74"/>
      <variable name="V_i" units="micrometre3" public_interface="out"/>
      <variable name="V_e" units="micrometre3" public_interface="out"/>
      <variable name="V_Cell" units="micrometre3"/>
      <variable name="V_i_ratio" units="dimensionless" initial_value="0.49"/>
      <variable name="V_ds_ratio" units="dimensionless" initial_value="0.1"/>
      <variable name="V_rel_ratio" units="dimensionless" initial_value="0.003"/>
      <variable name="V_e_ratio" units="dimensionless" initial_value="0.4"/>
      <variable name="V_up_ratio" units="dimensionless" initial_value="0.03"/>
      <variable name="Kdecay" units="per_second" initial_value="10"/>
      <variable name="i_up" units="millimolar_per_second" public_interface="in"/>
      <variable name="i_trans" units="millimolar_per_second" public_interface="in"/>
      <variable name="i_rel" units="millimolar_per_second" public_interface="in"/>
      <variable name="i_NaCa_cyt" units="nanoA" public_interface="in"/>
      <variable name="i_NaCa_ds" units="nanoA" public_interface="in"/>
      <variable name="i_Ca_L_Ca_cyt" units="nanoA" public_interface="in"/>
      <variable name="i_Ca_L_Ca_ds" units="nanoA" public_interface="in"/>
      <variable name="i_b_Ca" units="nanoA" public_interface="in"/>
      <variable name="F" units="coulomb_per_mole" public_interface="in"/>
      <variable name="time" units="second" public_interface="in"/>
      <variable name="A" units="millimolar" initial_value="0.00015" public_interface="out"/>
      <variable name="N_A" units="dimensionless"/>
      <variable name="pi_N_A" units="dimensionless"/>
      <variable name="pi_min" units="dimensionless" initial_value="0.03"/>
      <variable name="N" units="dimensionless" public_interface="in"/>
      <variable name="a_on" units="per_millimolar_second" initial_value="70000"/>
      <variable name="a_off" units="per_second" initial_value="200"/>
      <variable name="k_A" units="per_millimolar" public_interface="in"/>
      <variable name="L_oz" units="dimensionless" public_interface="in"/>
      <variable name="A_tot" units="millimolar" initial_value="0.07" public_interface="out"/>
      <variable name="B_1" units="millimolar" initial_value="0" public_interface="out"/>
      <variable name="B_1_tot" units="millimolar" initial_value="0.08"/>
      <variable name="b_1_on" units="per_millimolar_second" initial_value="100000"/>
      <variable name="b_1_off" units="per_second" initial_value="182"/>
      <variable name="B_2" units="millimolar" initial_value="0" public_interface="out"/>
      <variable name="B_2_tot" units="millimolar" initial_value="0.1"/>
      <variable name="b_2_on" units="per_millimolar_second" initial_value="1000"/>
      <variable name="b_2_off" units="per_second" initial_value="3"/>
      <variable name="CaS_tot" units="millimolar" public_interface="in"/>
      <variable name="beta" units="millimolar" public_interface="in"/>
      <variable name="n_NaCa" units="dimensionless" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>N_A</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <ci>A_tot</ci>
                  <ci>N</ci>
               </apply>
               <apply>
                  <times/>
                  <ci>L_oz</ci>
                  <ci>A</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>pi_N_A</ci>
            <piecewise>
               <piece>
                  <ci>pi_min</ci>
                  <apply>
                     <leq/>
                     <ci>N_A</ci>
                     <cn cellml:units="dimensionless">0</cn>
                  </apply>
               </piece>
               <piece>
                  <apply>
                     <power/>
                     <ci>pi_min</ci>
                     <ci>N_A</ci>
                  </apply>
                  <apply>
                     <leq/>
                     <ci>N_A</ci>
                     <cn cellml:units="dimensionless">1</cn>
                  </apply>
               </piece>
               <otherwise>
                  <cn cellml:units="dimensionless">1</cn>
               </otherwise>
            </piecewise>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>A</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>a_on</ci>
                  <apply>
                     <minus/>
                     <ci>A_tot</ci>
                     <ci>A</ci>
                  </apply>
                  <ci>Ca_i</ci>
               </apply>
               <apply>
                  <times/>
                  <ci>a_off</ci>
                  <apply>
                     <exp/>
                     <apply>
                        <times/>
                        <apply>
                           <minus/>
                           <ci>k_A</ci>
                        </apply>
                        <ci>A</ci>
                     </apply>
                  </apply>
                  <ci>pi_N_A</ci>
                  <ci>A</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>B_1</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>b_1_on</ci>
                  <apply>
                     <minus/>
                     <ci>B_1_tot</ci>
                     <ci>B_1</ci>
                  </apply>
                  <ci>Ca_i</ci>
               </apply>
               <apply>
                  <times/>
                  <ci>b_1_off</ci>
                  <ci>B_1</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>B_2</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>b_2_on</ci>
                  <apply>
                     <minus/>
                     <ci>B_2_tot</ci>
                     <ci>B_2</ci>
                  </apply>
                  <ci>Ca_i</ci>
               </apply>
               <apply>
                  <times/>
                  <ci>b_2_off</ci>
                  <ci>B_2</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>V_Cell</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <cn cellml:units="dimensionless">3.141592654</cn>
                  <apply>
                     <power/>
                     <apply>
                        <divide/>
                        <ci>radius</ci>
                        <cn cellml:units="dimensionless">1000</cn>
                     </apply>
                     <cn cellml:units="dimensionless">2</cn>
                  </apply>
                  <ci>length</ci>
               </apply>
               <cn cellml:units="dimensionless">1000</cn>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>V_i</ci>
            <apply>
               <times/>
               <ci>V_Cell</ci>
               <ci>V_i_ratio</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>V_e</ci>
            <apply>
               <times/>
               <ci>V_Cell</ci>
               <ci>V_e_ratio</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>Ca_i</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <minus/>
                  <apply>
                     <minus/>
                     <apply>
                        <minus/>
                        <apply>
                           <plus/>
                           <apply>
                              <times/>
                              <apply>
                                 <divide/>
                                 <apply>
                                    <minus/>
                                    <cn cellml:units="mA_nA">1</cn>
                                 </apply>
                                 <apply>
                                    <times/>
                                    <cn cellml:units="dimensionless">2</cn>
                                    <cn cellml:units="litre_micrometre3">1</cn>
                                    <ci>V_i</ci>
                                    <ci>F</ci>
                                 </apply>
                              </apply>
                              <apply>
                                 <minus/>
                                 <apply>
                                    <plus/>
                                    <ci>i_Ca_L_Ca_cyt</ci>
                                    <ci>i_b_Ca</ci>
                                 </apply>
                                 <apply>
                                    <times/>
                                    <apply>
                                       <divide/>
                                       <cn cellml:units="dimensionless">2</cn>
                                       <apply>
                                          <minus/>
                                          <ci>n_NaCa</ci>
                                          <cn cellml:units="dimensionless">2</cn>
                                       </apply>
                                    </apply>
                                    <ci>i_NaCa_cyt</ci>
                                 </apply>
                              </apply>
                           </apply>
                           <apply>
                              <times/>
                              <ci>Ca_ds</ci>
                              <ci>V_ds_ratio</ci>
                              <ci>Kdecay</ci>
                           </apply>
                           <apply>
                              <divide/>
                              <apply>
                                 <times/>
                                 <ci>i_rel</ci>
                                 <ci>V_rel_ratio</ci>
                              </apply>
                              <ci>V_i_ratio</ci>
                           </apply>
                        </apply>
                        <apply>
                           <diff/>
                           <bvar>
                              <ci>time</ci>
                           </bvar>
                           <ci>A</ci>
                        </apply>
                     </apply>
                     <apply>
                        <diff/>
                        <bvar>
                           <ci>time</ci>
                        </bvar>
                        <ci>B_1</ci>
                     </apply>
                  </apply>
                  <apply>
                     <diff/>
                     <bvar>
                        <ci>time</ci>
                     </bvar>
                     <ci>B_2</ci>
                  </apply>
               </apply>
               <ci>i_up</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>Ca_ds</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <divide/>
                  <apply>
                     <plus/>
                     <apply>
                        <times/>
                        <apply>
                           <minus/>
                           <cn cellml:units="mA_nA">1</cn>
                        </apply>
                        <ci>i_Ca_L_Ca_ds</ci>
                     </apply>
                     <apply>
                        <divide/>
                        <apply>
                           <times/>
                           <cn cellml:units="mA_nA">2</cn>
                           <ci>i_NaCa_ds</ci>
                        </apply>
                        <apply>
                           <minus/>
                           <ci>n_NaCa</ci>
                           <cn cellml:units="dimensionless">2</cn>
                        </apply>
                     </apply>
                  </apply>
                  <apply>
                     <times/>
                     <cn cellml:units="dimensionless">2</cn>
                     <cn cellml:units="litre_micrometre3">1</cn>
                     <ci>V_ds_ratio</ci>
                     <ci>V_i</ci>
                     <ci>F</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>Ca_ds</ci>
                  <ci>Kdecay</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>Ca_up</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <apply>
                     <divide/>
                     <ci>V_i_ratio</ci>
                     <ci>V_up_ratio</ci>
                  </apply>
                  <ci>i_up</ci>
               </apply>
               <ci>i_trans</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>Ca_rel</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <apply>
                     <times/>
                     <apply>
                        <divide/>
                        <ci>V_up_ratio</ci>
                        <ci>V_rel_ratio</ci>
                     </apply>
                     <ci>i_trans</ci>
                  </apply>
                  <ci>i_rel</ci>
               </apply>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <ci>beta</ci>
                        <ci>CaS_tot</ci>
                     </apply>
                     <apply>
                        <power/>
                        <apply>
                           <plus/>
                           <ci>Ca_rel</ci>
                           <ci>beta</ci>
                        </apply>
                        <cn cellml:units="dimensionless">2</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="fibroblast">
      <variable name="g_fibro_junct" units="microS" initial_value="2.9e-4"/>
      <variable name="g_fibro" units="microS" initial_value="2e-4"/>
      <variable name="c_fibro" units="microF" initial_value="1e-5"/>
      <variable name="g_fibro_stretch" units="microS" initial_value="0"/>
      <variable name="E_fibro_stretch" units="millivolt" initial_value="0"/>
      <variable name="V_fibro" units="millivolt" initial_value="-20"/>
      <variable name="i_fibro" units="nanoA"/>
      <variable name="i_fibro_junct" units="nanoA"/>
      <variable name="V" units="millivolt" public_interface="in"/>
      <variable name="time" units="second" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_fibro</ci>
            <apply>
               <plus/>
               <apply>
                  <times/>
                  <ci>g_fibro</ci>
                  <apply>
                     <plus/>
                     <ci>V_fibro</ci>
                     <cn cellml:units="millivolt">20</cn>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>g_fibro_stretch</ci>
                  <apply>
                     <minus/>
                     <ci>V_fibro</ci>
                     <ci>E_fibro_stretch</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>i_fibro_junct</ci>
            <apply>
               <times/>
               <apply>
                  <minus/>
                  <ci>g_fibro_junct</ci>
               </apply>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>V_fibro</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>V_fibro</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <apply>
                     <plus/>
                     <ci>i_fibro</ci>
                     <ci>i_fibro_junct</ci>
                  </apply>
               </apply>
               <ci>c_fibro</ci>
            </apply>
         </apply>
      </math>
   </component>
   <group>
      <relationship_ref relationship="containment"/>
      <component_ref component="membrane">
         <component_ref component="fast_sodium_current">
            <component_ref component="fast_sodium_current_m_gate"/>
            <component_ref component="fast_sodium_current_h_gate"/>
         </component_ref>
         <component_ref component="L_type_Ca_channel">
            <component_ref component="L_type_Ca_channel_d_gate"/>
            <component_ref component="L_type_Ca_channel_f_gate"/>
            <component_ref component="L_type_Ca_channel_f2_gate"/>
            <component_ref component="L_type_Ca_channel_f2ds_gate"/>
         </component_ref>
         <component_ref component="rapid_delayed_rectifier_potassium_current">
            <component_ref component="rapid_delayed_rectifier_potassium_current_xr1_gate"/>
            <component_ref component="rapid_delayed_rectifier_potassium_current_xr2_gate"/>
         </component_ref>
         <component_ref component="slow_delayed_rectifier_potassium_current">
            <component_ref component="slow_delayed_rectifier_potassium_current_xs_gate"/>
         </component_ref>
         <component_ref component="time_independent_potassium_current"/>
         <component_ref component="sodium_activated_potassium_current"/>
         <component_ref component="persistent_sodium_current"/>
         <component_ref component="sodium_background_current"/>
         <component_ref component="calcium_background_current"/>
         <component_ref component="transient_outward_current">
            <component_ref component="transient_outward_current_s_gate"/>
            <component_ref component="transient_outward_current_r_gate"/>
         </component_ref>
         <component_ref component="sodium_potassium_pump"/>
         <component_ref component="sodium_calcium_exchanger"/>
         <component_ref component="sarcoplasmic_reticulum_calcium_pump"/>
         <component_ref component="calcium_translocation"/>
         <component_ref component="calcium_release"/>
         <component_ref component="extracellular_sodium_concentration"/>
         <component_ref component="intracellular_sodium_concentration"/>
         <component_ref component="extracellular_potassium_concentration"/>
         <component_ref component="intracellular_potassium_concentration"/>
         <component_ref component="extracellular_calcium_concentration"/>
         <component_ref component="intracellular_calcium_concentration"/>
         <component_ref component="reversal_potentials"/>
         <component_ref component="fibroblast"/>
      </component_ref>
   </group>
   <group>
      <relationship_ref relationship="encapsulation"/>
      <component_ref component="fast_sodium_current">
         <component_ref component="fast_sodium_current_m_gate"/>
         <component_ref component="fast_sodium_current_h_gate"/>
      </component_ref>
      <component_ref component="L_type_Ca_channel">
         <component_ref component="L_type_Ca_channel_d_gate"/>
         <component_ref component="L_type_Ca_channel_f_gate"/>
         <component_ref component="L_type_Ca_channel_f2_gate"/>
         <component_ref component="L_type_Ca_channel_f2ds_gate"/>
      </component_ref>
      <component_ref component="transient_outward_current">
         <component_ref component="transient_outward_current_s_gate"/>
         <component_ref component="transient_outward_current_r_gate"/>
      </component_ref>
      <component_ref component="rapid_delayed_rectifier_potassium_current">
         <component_ref component="rapid_delayed_rectifier_potassium_current_xr1_gate"/>
         <component_ref component="rapid_delayed_rectifier_potassium_current_xr2_gate"/>
      </component_ref>
      <component_ref component="slow_delayed_rectifier_potassium_current">
         <component_ref component="slow_delayed_rectifier_potassium_current_xs_gate"/>
      </component_ref>
   </group>
   <connection>
      <map_components component_1="membrane" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="reversal_potentials" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="fast_sodium_current" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="L_type_Ca_channel" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="rapid_delayed_rectifier_potassium_current" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="slow_delayed_rectifier_potassium_current" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="sodium_activated_potassium_current" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="time_independent_potassium_current" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="persistent_sodium_current" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="sodium_background_current" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="calcium_background_current" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="transient_outward_current" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="sodium_potassium_pump" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="sodium_calcium_exchanger" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="intracellular_sodium_concentration" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="intracellular_potassium_concentration" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="intracellular_calcium_concentration" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="sarcoplasmic_reticulum_calcium_pump" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="calcium_translocation" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="calcium_release" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="fibroblast" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="calcium_release" component_2="membrane"/>
      <map_variables variable_1="V" variable_2="V"/>
   </connection>
   <connection>
      <map_components component_1="membrane" component_2="reversal_potentials"/>
      <map_variables variable_1="R" variable_2="R"/>
      <map_variables variable_1="T" variable_2="T"/>
      <map_variables variable_1="F" variable_2="F"/>
   </connection>
   <connection>
      <map_components component_1="membrane" component_2="fast_sodium_current"/>
      <map_variables variable_1="V" variable_2="V"/>
      <map_variables variable_1="i_Na" variable_2="i_Na"/>
   </connection>
   <connection>
      <map_components component_1="membrane" component_2="L_type_Ca_channel"/>
      <map_variables variable_1="V" variable_2="V"/>
      <map_variables variable_1="i_Ca_L_Ca_cyt" variable_2="i_Ca_L_Ca_cyt"/>
      <map_variables variable_1="i_Ca_L_Ca_ds" variable_2="i_Ca_L_Ca_ds"/>
      <map_variables variable_1="i_Ca_L_Na_cyt" variable_2="i_Ca_L_Na_cyt"/>
      <map_variables variable_1="i_Ca_L_Na_ds" variable_2="i_Ca_L_Na_ds"/>
      <map_variables variable_1="i_Ca_L_K_cyt" variable_2="i_Ca_L_K_cyt"/>
      <map_variables variable_1="i_Ca_L_K_ds" variable_2="i_Ca_L_K_ds"/>
      <map_variables variable_1="R" variable_2="R"/>
      <map_variables variable_1="T" variable_2="T"/>
      <map_variables variable_1="F" variable_2="F"/>
   </connection>
   <connection>
      <map_components component_1="membrane" component_2="rapid_delayed_rectifier_potassium_current"/>
      <map_variables variable_1="V" variable_2="V"/>
      <map_variables variable_1="i_Kr" variable_2="i_Kr"/>
   </connection>
   <connection>
      <map_components component_1="membrane" component_2="slow_delayed_rectifier_potassium_current"/>
      <map_variables variable_1="V" variable_2="V"/>
      <map_variables variable_1="i_Ks" variable_2="i_Ks"/>
   </connection>
   <connection>
      <map_components component_1="membrane" component_2="transient_outward_current"/>
      <map_variables variable_1="V" variable_2="V"/>
      <map_variables variable_1="i_to" variable_2="i_to"/>
   </connection>
   <connection>
      <map_components component_1="membrane" component_2="time_independent_potassium_current"/>
      <map_variables variable_1="V" variable_2="V"/>
      <map_variables variable_1="i_K1" variable_2="i_K1"/>
      <map_variables variable_1="R" variable_2="R"/>
      <map_variables variable_1="T" variable_2="T"/>
      <map_variables variable_1="F" variable_2="F"/>
   </connection>
   <connection>
      <map_components component_1="membrane" component_2="sodium_activated_potassium_current"/>
      <map_variables variable_1="V" variable_2="V"/>
   </connection>
   <connection>
      <map_components component_1="membrane" component_2="persistent_sodium_current"/>
      <map_variables variable_1="V" variable_2="V"/>
      <map_variables variable_1="i_p_Na" variable_2="i_p_Na"/>
   </connection>
   <connection>
      <map_components component_1="membrane" component_2="sodium_background_current"/>
      <map_variables variable_1="V" variable_2="V"/>
      <map_variables variable_1="i_b_Na" variable_2="i_b_Na"/>
   </connection>
   <connection>
      <map_components component_1="membrane" component_2="calcium_background_current"/>
      <map_variables variable_1="V" variable_2="V"/>
      <map_variables variable_1="i_b_Ca" variable_2="i_b_Ca"/>
   </connection>
   <connection>
      <map_components component_1="membrane" component_2="sodium_potassium_pump"/>
      <map_variables variable_1="i_NaK" variable_2="i_NaK"/>
   </connection>
   <connection>
      <map_components component_1="membrane" component_2="sodium_calcium_exchanger"/>
      <map_variables variable_1="V" variable_2="V"/>
      <map_variables variable_1="i_NaCa_cyt" variable_2="i_NaCa_cyt"/>
      <map_variables variable_1="i_NaCa_ds" variable_2="i_NaCa_ds"/>
      <map_variables variable_1="R" variable_2="R"/>
      <map_variables variable_1="T" variable_2="T"/>
      <map_variables variable_1="F" variable_2="F"/>
   </connection>
   <connection>
      <map_components component_1="membrane" component_2="fibroblast"/>
      <map_variables variable_1="V" variable_2="V"/>
   </connection>
   <connection>
      <map_components component_1="membrane" component_2="intracellular_sodium_concentration"/>
      <map_variables variable_1="F" variable_2="F"/>
   </connection>
   <connection>
      <map_components component_1="membrane" component_2="intracellular_potassium_concentration"/>
      <map_variables variable_1="F" variable_2="F"/>
   </connection>
   <connection>
      <map_components component_1="membrane" component_2="intracellular_calcium_concentration"/>
      <map_variables variable_1="F" variable_2="F"/>
   </connection>
   <connection>
      <map_components component_1="reversal_potentials" component_2="sodium_background_current"/>
      <map_variables variable_1="E_Na" variable_2="E_Na"/>
   </connection>
   <connection>
      <map_components component_1="reversal_potentials" component_2="persistent_sodium_current"/>
      <map_variables variable_1="E_Na" variable_2="E_Na"/>
   </connection>
   <connection>
      <map_components component_1="reversal_potentials" component_2="time_independent_potassium_current"/>
      <map_variables variable_1="E_K" variable_2="E_K"/>
   </connection>
   <connection>
      <map_components component_1="reversal_potentials" component_2="rapid_delayed_rectifier_potassium_current"/>
      <map_variables variable_1="E_K" variable_2="E_K"/>
   </connection>
   <connection>
      <map_components component_1="reversal_potentials" component_2="slow_delayed_rectifier_potassium_current"/>
      <map_variables variable_1="E_Ks" variable_2="E_Ks"/>
   </connection>
   <connection>
      <map_components component_1="reversal_potentials" component_2="sodium_activated_potassium_current"/>
      <map_variables variable_1="E_K" variable_2="E_K"/>
   </connection>
   <connection>
      <map_components component_1="reversal_potentials" component_2="transient_outward_current"/>
      <map_variables variable_1="E_K" variable_2="E_K"/>
   </connection>
   <connection>
      <map_components component_1="reversal_potentials" component_2="calcium_background_current"/>
      <map_variables variable_1="E_Ca" variable_2="E_Ca"/>
   </connection>
   <connection>
      <map_components component_1="reversal_potentials" component_2="fast_sodium_current"/>
      <map_variables variable_1="E_mh" variable_2="E_mh"/>
   </connection>
   <connection>
      <map_components component_1="reversal_potentials" component_2="extracellular_potassium_concentration"/>
      <map_variables variable_1="K_o" variable_2="K_o"/>
   </connection>
   <connection>
      <map_components component_1="reversal_potentials" component_2="extracellular_sodium_concentration"/>
      <map_variables variable_1="Na_o" variable_2="Na_o"/>
   </connection>
   <connection>
      <map_components component_1="reversal_potentials" component_2="intracellular_potassium_concentration"/>
      <map_variables variable_1="K_i" variable_2="K_i"/>
   </connection>
   <connection>
      <map_components component_1="reversal_potentials" component_2="intracellular_sodium_concentration"/>
      <map_variables variable_1="Na_i" variable_2="Na_i"/>
   </connection>
   <connection>
      <map_components component_1="reversal_potentials" component_2="extracellular_calcium_concentration"/>
      <map_variables variable_1="Ca_o" variable_2="Ca_o"/>
   </connection>
   <connection>
      <map_components component_1="reversal_potentials" component_2="intracellular_calcium_concentration"/>
      <map_variables variable_1="Ca_i" variable_2="Ca_i"/>
   </connection>
   <connection>
      <map_components component_1="time_independent_potassium_current" component_2="intracellular_potassium_concentration"/>
      <map_variables variable_1="i_K1" variable_2="i_K1"/>
   </connection>
   <connection>
      <map_components component_1="rapid_delayed_rectifier_potassium_current" component_2="intracellular_potassium_concentration"/>
      <map_variables variable_1="i_Kr" variable_2="i_Kr"/>
   </connection>
   <connection>
      <map_components component_1="slow_delayed_rectifier_potassium_current" component_2="intracellular_potassium_concentration"/>
      <map_variables variable_1="i_Ks" variable_2="i_Ks"/>
   </connection>
   <connection>
      <map_components component_1="fast_sodium_current" component_2="intracellular_sodium_concentration"/>
      <map_variables variable_1="i_Na" variable_2="i_Na"/>
   </connection>
   <connection>
      <map_components component_1="sodium_background_current" component_2="intracellular_sodium_concentration"/>
      <map_variables variable_1="i_b_Na" variable_2="i_b_Na"/>
      <map_variables variable_1="Na_i" variable_2="Na_i"/>
   </connection>
   <connection>
      <map_components component_1="persistent_sodium_current" component_2="intracellular_sodium_concentration"/>
      <map_variables variable_1="i_p_Na" variable_2="i_p_Na"/>
   </connection>
   <connection>
      <map_components component_1="L_type_Ca_channel" component_2="intracellular_sodium_concentration"/>
      <map_variables variable_1="i_Ca_L_Na_cyt" variable_2="i_Ca_L_Na_cyt"/>
      <map_variables variable_1="i_Ca_L_Na_ds" variable_2="i_Ca_L_Na_ds"/>
      <map_variables variable_1="Na_i" variable_2="Na_i"/>
   </connection>
   <connection>
      <map_components component_1="L_type_Ca_channel" component_2="intracellular_calcium_concentration"/>
      <map_variables variable_1="i_Ca_L_Ca_cyt" variable_2="i_Ca_L_Ca_cyt"/>
      <map_variables variable_1="i_Ca_L_Ca_ds" variable_2="i_Ca_L_Ca_ds"/>
      <map_variables variable_1="Ca_i" variable_2="Ca_i"/>
      <map_variables variable_1="Ca_ds" variable_2="Ca_ds"/>
   </connection>
   <connection>
      <map_components component_1="L_type_Ca_channel" component_2="intracellular_potassium_concentration"/>
      <map_variables variable_1="i_Ca_L_K_cyt" variable_2="i_Ca_L_K_cyt"/>
      <map_variables variable_1="i_Ca_L_K_ds" variable_2="i_Ca_L_K_ds"/>
      <map_variables variable_1="K_i" variable_2="K_i"/>
   </connection>
   <connection>
      <map_components component_1="L_type_Ca_channel" component_2="extracellular_sodium_concentration"/>
      <map_variables variable_1="Na_o" variable_2="Na_o"/>
   </connection>
   <connection>
      <map_components component_1="L_type_Ca_channel" component_2="extracellular_calcium_concentration"/>
      <map_variables variable_1="Ca_o" variable_2="Ca_o"/>
   </connection>
   <connection>
      <map_components component_1="calcium_background_current" component_2="intracellular_calcium_concentration"/>
      <map_variables variable_1="i_b_Ca" variable_2="i_b_Ca"/>
   </connection>
   <connection>
      <map_components component_1="sodium_activated_potassium_current" component_2="intracellular_sodium_concentration"/>
      <map_variables variable_1="Na_i" variable_2="Na_i"/>
   </connection>
   <connection>
      <map_components component_1="sodium_potassium_pump" component_2="intracellular_sodium_concentration"/>
      <map_variables variable_1="i_NaK" variable_2="i_NaK"/>
      <map_variables variable_1="Na_i" variable_2="Na_i"/>
   </connection>
   <connection>
      <map_components component_1="sodium_potassium_pump" component_2="intracellular_potassium_concentration"/>
      <map_variables variable_1="i_NaK" variable_2="i_NaK"/>
   </connection>
   <connection>
      <map_components component_1="transient_outward_current" component_2="intracellular_potassium_concentration"/>
      <map_variables variable_1="i_to" variable_2="i_to"/>
   </connection>
   <connection>
      <map_components component_1="sodium_calcium_exchanger" component_2="intracellular_calcium_concentration"/>
      <map_variables variable_1="i_NaCa_cyt" variable_2="i_NaCa_cyt"/>
      <map_variables variable_1="i_NaCa_ds" variable_2="i_NaCa_ds"/>
      <map_variables variable_1="Ca_i" variable_2="Ca_i"/>
      <map_variables variable_1="Ca_ds" variable_2="Ca_ds"/>
      <map_variables variable_1="n_NaCa" variable_2="n_NaCa"/>
   </connection>
   <connection>
      <map_components component_1="sodium_calcium_exchanger" component_2="intracellular_sodium_concentration"/>
      <map_variables variable_1="i_NaCa_cyt" variable_2="i_NaCa_cyt"/>
      <map_variables variable_1="i_NaCa" variable_2="i_NaCa"/>
      <map_variables variable_1="n_NaCa" variable_2="n_NaCa"/>
      <map_variables variable_1="Na_i" variable_2="Na_i"/>
   </connection>
   <connection>
      <map_components component_1="sodium_calcium_exchanger" component_2="extracellular_sodium_concentration"/>
      <map_variables variable_1="Na_o" variable_2="Na_o"/>
   </connection>
   <connection>
      <map_components component_1="sodium_calcium_exchanger" component_2="extracellular_calcium_concentration"/>
      <map_variables variable_1="Ca_o" variable_2="Ca_o"/>
   </connection>
   <connection>
      <map_components component_1="sarcoplasmic_reticulum_calcium_pump" component_2="intracellular_calcium_concentration"/>
      <map_variables variable_1="i_up" variable_2="i_up"/>
      <map_variables variable_1="Ca_i" variable_2="Ca_i"/>
      <map_variables variable_1="Ca_up" variable_2="Ca_up"/>
   </connection>
   <connection>
      <map_components component_1="calcium_translocation" component_2="intracellular_calcium_concentration"/>
      <map_variables variable_1="i_trans" variable_2="i_trans"/>
      <map_variables variable_1="Ca_rel" variable_2="Ca_rel"/>
      <map_variables variable_1="Ca_up" variable_2="Ca_up"/>
      <map_variables variable_1="CaS_tot" variable_2="CaS_tot"/>
      <map_variables variable_1="beta" variable_2="beta"/>
   </connection>
   <connection>
      <map_components component_1="calcium_release" component_2="intracellular_calcium_concentration"/>
      <map_variables variable_1="i_rel" variable_2="i_rel"/>
      <map_variables variable_1="Ca_rel" variable_2="Ca_rel"/>
      <map_variables variable_1="Ca_i" variable_2="Ca_i"/>
      <map_variables variable_1="Ca_ds" variable_2="Ca_ds"/>
   </connection>
   <connection>
      <map_components component_1="rapid_delayed_rectifier_potassium_current" component_2="rapid_delayed_rectifier_potassium_current_xr1_gate"/>
      <map_variables variable_1="xr1" variable_2="xr1"/>
      <map_variables variable_1="time" variable_2="time"/>
      <map_variables variable_1="V" variable_2="V"/>
   </connection>
   <connection>
      <map_components component_1="rapid_delayed_rectifier_potassium_current" component_2="rapid_delayed_rectifier_potassium_current_xr2_gate"/>
      <map_variables variable_1="xr2" variable_2="xr2"/>
      <map_variables variable_1="time" variable_2="time"/>
      <map_variables variable_1="V" variable_2="V"/>
   </connection>
   <connection>
      <map_components component_1="slow_delayed_rectifier_potassium_current" component_2="slow_delayed_rectifier_potassium_current_xs_gate"/>
      <map_variables variable_1="xs" variable_2="xs"/>
      <map_variables variable_1="time" variable_2="time"/>
      <map_variables variable_1="V" variable_2="V"/>
   </connection>
   <connection>
      <map_components component_1="fast_sodium_current" component_2="fast_sodium_current_m_gate"/>
      <map_variables variable_1="m" variable_2="m"/>
      <map_variables variable_1="time" variable_2="time"/>
      <map_variables variable_1="V" variable_2="V"/>
   </connection>
   <connection>
      <map_components component_1="fast_sodium_current" component_2="fast_sodium_current_h_gate"/>
      <map_variables variable_1="h" variable_2="h"/>
      <map_variables variable_1="time" variable_2="time"/>
      <map_variables variable_1="V" variable_2="V"/>
   </connection>
   <connection>
      <map_components component_1="L_type_Ca_channel" component_2="L_type_Ca_channel_d_gate"/>
      <map_variables variable_1="d" variable_2="d"/>
      <map_variables variable_1="time" variable_2="time"/>
      <map_variables variable_1="V" variable_2="V"/>
   </connection>
   <connection>
      <map_components component_1="L_type_Ca_channel" component_2="L_type_Ca_channel_f_gate"/>
      <map_variables variable_1="f" variable_2="f"/>
      <map_variables variable_1="time" variable_2="time"/>
      <map_variables variable_1="V" variable_2="V"/>
   </connection>
   <connection>
      <map_components component_1="L_type_Ca_channel" component_2="L_type_Ca_channel_f2_gate"/>
      <map_variables variable_1="f2" variable_2="f2"/>
      <map_variables variable_1="time" variable_2="time"/>
      <map_variables variable_1="V" variable_2="V"/>
      <map_variables variable_1="Ca_i" variable_2="Ca_i"/>
      <map_variables variable_1="Km_f2" variable_2="Km_f2"/>
   </connection>
   <connection>
      <map_components component_1="L_type_Ca_channel" component_2="L_type_Ca_channel_f2ds_gate"/>
      <map_variables variable_1="f2ds" variable_2="f2ds"/>
      <map_variables variable_1="time" variable_2="time"/>
      <map_variables variable_1="V" variable_2="V"/>
      <map_variables variable_1="Ca_ds" variable_2="Ca_ds"/>
      <map_variables variable_1="Km_f2ds" variable_2="Km_f2ds"/>
      <map_variables variable_1="R_decay" variable_2="R_decay"/>
   </connection>
   <connection>
      <map_components component_1="transient_outward_current" component_2="transient_outward_current_s_gate"/>
      <map_variables variable_1="s" variable_2="s"/>
      <map_variables variable_1="time" variable_2="time"/>
      <map_variables variable_1="V" variable_2="V"/>
   </connection>
   <connection>
      <map_components component_1="transient_outward_current" component_2="transient_outward_current_r_gate"/>
      <map_variables variable_1="r" variable_2="r"/>
      <map_variables variable_1="time" variable_2="time"/>
      <map_variables variable_1="V" variable_2="V"/>
   </connection>
   <connection>
      <map_components component_1="intracellular_sodium_concentration" component_2="intracellular_calcium_concentration"/>
      <map_variables variable_1="V_i" variable_2="V_i"/>
   </connection>
   <connection>
      <map_components component_1="intracellular_potassium_concentration" component_2="intracellular_calcium_concentration"/>
      <map_variables variable_1="V_i" variable_2="V_i"/>
   </connection>
   <connection>
      <map_components component_1="parameters" component_2="isotonic"/>
      <map_variables variable_1="isotonic" variable_2="isotonic"/>
   </connection>
   <connection>
      <map_components component_1="force" component_2="isotonic"/>
      <map_variables variable_1="F_muscle" variable_2="F_muscle"/>
   </connection>
   <connection>
      <map_components component_1="length" component_2="isotonic"/>
      <map_variables variable_1="isotonic_mode" variable_2="isotonic_mode"/>
      <map_variables variable_1="l" variable_2="l"/>
      <map_variables variable_1="l_0" variable_2="l_0"/>
   </connection>
   <connection>
      <map_components component_1="PE_velocity" component_2="isotonic"/>
      <map_variables variable_1="isotonic_mode" variable_2="isotonic_mode"/>
   </connection>
   <connection>
      <map_components component_1="CE_velocity" component_2="isotonic"/>
      <map_variables variable_1="isotonic_mode" variable_2="isotonic_mode"/>
   </connection>
   <connection>
      <map_components component_1="membrane" component_2="isotonic"/>
      <map_variables variable_1="stim_period" variable_2="stim_period"/>
   </connection>
   <connection>
      <map_components component_1="environment" component_2="crossbridge_kinetics"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="environment" component_2="length"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="environment" component_2="CE_velocity"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="environment" component_2="PE_velocity"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="parameters" component_2="parameters_izakov_et_al_1991"/>
      <map_variables variable_1="v_max" variable_2="v_max"/>
   </connection>
   <connection>
      <map_components component_1="force" component_2="parameters"/>
      <map_variables variable_1="alpha_1" variable_2="alpha_1"/>
      <map_variables variable_1="beta_1" variable_2="beta_1"/>
      <map_variables variable_1="alpha_2" variable_2="alpha_2"/>
      <map_variables variable_1="beta_2" variable_2="beta_2"/>
      <map_variables variable_1="alpha_3" variable_2="alpha_3"/>
      <map_variables variable_1="beta_3" variable_2="beta_3"/>
      <map_variables variable_1="lambda" variable_2="lambda"/>
   </connection>
   <connection>
      <map_components component_1="force" component_2="crossbridge_kinetics"/>
      <map_variables variable_1="N" variable_2="N"/>
   </connection>
   <connection>
      <map_components component_1="force" component_2="length"/>
      <map_variables variable_1="l_1" variable_2="l_1"/>
      <map_variables variable_1="l_2" variable_2="l_2"/>
      <map_variables variable_1="l_3" variable_2="l_3"/>
   </connection>
   <connection>
      <map_components component_1="force" component_2="CE_velocity"/>
      <map_variables variable_1="v" variable_2="v"/>
      <map_variables variable_1="k_P_vis" variable_2="k_P_vis"/>
   </connection>
   <connection>
      <map_components component_1="force" component_2="PE_velocity"/>
      <map_variables variable_1="w" variable_2="w"/>
      <map_variables variable_1="k_S_vis" variable_2="k_S_vis"/>
   </connection>
   <connection>
      <map_components component_1="force" component_2="average_crossbridge_force"/>
      <map_variables variable_1="p_v" variable_2="p_v"/>
   </connection>
   <connection>
      <map_components component_1="crossbridge_kinetics" component_2="parameters"/>
      <map_variables variable_1="mu" variable_2="mu"/>
      <map_variables variable_1="k_mu" variable_2="k_mu"/>
      <map_variables variable_1="kappa" variable_2="kappa"/>
      <map_variables variable_1="kappa_0" variable_2="kappa_0"/>
      <map_variables variable_1="m_0" variable_2="m_0"/>
   </connection>
   <connection>
      <map_components component_1="crossbridge_kinetics" component_2="intracellular_calcium_concentration"/>
      <map_variables variable_1="A" variable_2="A"/>
      <map_variables variable_1="A_tot" variable_2="A_tot"/>
      <map_variables variable_1="N" variable_2="N"/>
      <map_variables variable_1="L_oz" variable_2="L_oz"/>
   </connection>
   <connection>
      <map_components component_1="crossbridge_kinetics" component_2="average_crossbridge_force"/>
      <map_variables variable_1="G_star" variable_2="G_star"/>
      <map_variables variable_1="P_star" variable_2="P_star"/>
   </connection>
   <connection>
      <map_components component_1="crossbridge_kinetics" component_2="length"/>
      <map_variables variable_1="l_1" variable_2="l_1"/>
   </connection>
   <connection>
      <map_components component_1="crossbridge_kinetics" component_2="CE_velocity"/>
      <map_variables variable_1="K_kappa" variable_2="K_kappa"/>
      <map_variables variable_1="v" variable_2="v"/>
      <map_variables variable_1="N" variable_2="N"/>
   </connection>
   <connection>
      <map_components component_1="crossbridge_kinetics" component_2="parameters_izakov_et_al_1991"/>
      <map_variables variable_1="S_0" variable_2="S_0"/>
      <map_variables variable_1="q_v" variable_2="q_v"/>
      <map_variables variable_1="v_st" variable_2="v_st"/>
      <map_variables variable_1="q_4" variable_2="q_4"/>
   </connection>
   <connection>
      <map_components component_1="length" component_2="CE_velocity"/>
      <map_variables variable_1="v" variable_2="v"/>
      <map_variables variable_1="phi_chi_2" variable_2="phi_chi_2"/>
      <map_variables variable_1="dl_1_dt" variable_2="dl_1_dt"/>
      <map_variables variable_1="l_1" variable_2="l_1"/>
      <map_variables variable_1="l_2" variable_2="l_2"/>
      <map_variables variable_1="l_3" variable_2="l_3"/>
   </connection>
   <connection>
      <map_components component_1="length" component_2="PE_velocity"/>
      <map_variables variable_1="w" variable_2="w"/>
      <map_variables variable_1="k_S_vis" variable_2="k_S_vis"/>
      <map_variables variable_1="dl_1_dt" variable_2="dl_1_dt"/>
      <map_variables variable_1="l_1" variable_2="l_1"/>
      <map_variables variable_1="dl_2_dt" variable_2="dl_2_dt"/>
      <map_variables variable_1="l_2" variable_2="l_2"/>
      <map_variables variable_1="l_3" variable_2="l_3"/>
   </connection>
   <connection>
      <map_components component_1="length" component_2="membrane"/>
      <map_variables variable_1="stim_start" variable_2="stim_start"/>
      <map_variables variable_1="stim_end" variable_2="stim_end"/>
      <map_variables variable_1="stim_period" variable_2="stim_period"/>
      <map_variables variable_1="stim_duration" variable_2="stim_duration"/>
      <map_variables variable_1="stim_amplitude" variable_2="stim_amplitude"/>
   </connection>
   <connection>
      <map_components component_1="CE_velocity" component_2="parameters"/>
      <map_variables variable_1="alpha_1" variable_2="alpha_1"/>
      <map_variables variable_1="beta_1" variable_2="beta_1"/>
      <map_variables variable_1="alpha_2" variable_2="alpha_2"/>
      <map_variables variable_1="beta_2" variable_2="beta_2"/>
      <map_variables variable_1="alpha_3" variable_2="alpha_3"/>
      <map_variables variable_1="beta_3" variable_2="beta_3"/>
      <map_variables variable_1="lambda" variable_2="lambda"/>
   </connection>
   <connection>
      <map_components component_1="CE_velocity" component_2="parameters_izakov_et_al_1991"/>
      <map_variables variable_1="v" variable_2="v"/>
   </connection>
   <connection>
      <map_components component_1="CE_velocity" component_2="PE_velocity"/>
      <map_variables variable_1="v" variable_2="v"/>
      <map_variables variable_1="w" variable_2="w"/>
      <map_variables variable_1="k_S_vis" variable_2="k_S_vis"/>
      <map_variables variable_1="phi_chi" variable_2="phi_chi"/>
   </connection>
   <connection>
      <map_components component_1="PE_velocity" component_2="parameters"/>
      <map_variables variable_1="alpha_1" variable_2="alpha_1"/>
      <map_variables variable_1="beta_1" variable_2="beta_1"/>
      <map_variables variable_1="alpha_2" variable_2="alpha_2"/>
      <map_variables variable_1="beta_2" variable_2="beta_2"/>
      <map_variables variable_1="alpha_3" variable_2="alpha_3"/>
      <map_variables variable_1="beta_3" variable_2="beta_3"/>
      <map_variables variable_1="lambda" variable_2="lambda"/>
   </connection>
   <connection>
      <map_components component_1="average_crossbridge_force" component_2="parameters"/>
      <map_variables variable_1="v_max" variable_2="v_max"/>
      <map_variables variable_1="a" variable_2="a"/>
      <map_variables variable_1="d_h" variable_2="d_h"/>
      <map_variables variable_1="alpha_P" variable_2="alpha_P"/>
   </connection>
   <connection>
      <map_components component_1="average_crossbridge_force" component_2="parameters_izakov_et_al_1991"/>
      <map_variables variable_1="v_1" variable_2="v_1"/>
      <map_variables variable_1="alpha_G" variable_2="alpha_G"/>
   </connection>
   <connection>
      <map_components component_1="average_crossbridge_force" component_2="CE_velocity"/>
      <map_variables variable_1="v" variable_2="v"/>
      <map_variables variable_1="phi_chi" variable_2="phi_chi"/>
      <map_variables variable_1="p_v" variable_2="p_v"/>
      <map_variables variable_1="p_prime_v" variable_2="p_prime_v"/>
   </connection>
   <connection>
      <map_components component_1="intracellular_calcium_concentration" component_2="parameters_izakov_et_al_1991"/>
      <map_variables variable_1="k_A" variable_2="k_A"/>
   </connection>
   <connection>
      <map_components component_1="extracellular_potassium_concentration" component_2="rapid_delayed_rectifier_potassium_current"/>
      <map_variables variable_1="i_Kr" variable_2="i_Kr"/>
   </connection>
   <connection>
      <map_components component_1="extracellular_potassium_concentration" component_2="slow_delayed_rectifier_potassium_current"/>
      <map_variables variable_1="i_Ks" variable_2="i_Ks"/>
   </connection>
   <connection>
      <map_components component_1="extracellular_potassium_concentration" component_2="time_independent_potassium_current"/>
      <map_variables variable_1="i_K1" variable_2="i_K1"/>
      <map_variables variable_1="K_o" variable_2="K_o"/>
   </connection>
   <connection>
      <map_components component_1="extracellular_potassium_concentration" component_2="transient_outward_current"/>
      <map_variables variable_1="i_to" variable_2="i_to"/>
   </connection>
   <connection>
      <map_components component_1="extracellular_potassium_concentration" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="extracellular_potassium_concentration" component_2="sodium_potassium_pump"/>
      <map_variables variable_1="i_NaK" variable_2="i_NaK"/>
      <map_variables variable_1="K_o" variable_2="K_o"/>
   </connection>
   <connection>
      <map_components component_1="extracellular_potassium_concentration" component_2="L_type_Ca_channel"/>
      <map_variables variable_1="i_Ca_L_K_ds" variable_2="i_Ca_L_K_ds"/>
      <map_variables variable_1="i_Ca_L_K_cyt" variable_2="i_Ca_L_K_cyt"/>
      <map_variables variable_1="K_o" variable_2="K_o"/>
   </connection>
   <connection>
      <map_components component_1="extracellular_potassium_concentration" component_2="membrane"/>
      <map_variables variable_1="F" variable_2="F"/>
   </connection>
   <connection>
      <map_components component_1="extracellular_potassium_concentration" component_2="intracellular_sodium_concentration"/>
      <map_variables variable_1="n_NaK" variable_2="n_NaK"/>
   </connection>
   <connection>
      <map_components component_1="intracellular_potassium_concentration" component_2="intracellular_sodium_concentration"/>
      <map_variables variable_1="n_NaK" variable_2="n_NaK"/>
   </connection>
   <connection>
      <map_components component_1="extracellular_potassium_concentration" component_2="intracellular_calcium_concentration"/>
      <map_variables variable_1="V_e" variable_2="V_e"/>
   </connection>
</model>