Rendering of the source text

<?xml version='1.0' encoding='utf-8'?>
<!--  FILE :  wang_model_2007.xml

CREATED :  10th July 2007

LAST MODIFIED : 10th July 2007

AUTHOR :  Catherine Lloyd
          Bioengineering Institute
          The University of Auckland
          
MODEL STATUS :  This model conforms to the CellML 1.1 Specification.

DESCRIPTION :  This file contains a CellML description of Wang, Li, and Kuang's 2007 mathematical model and qualitative analysis of insulin therapies.

CHANGES:  
  
--><model xmlns="http://www.cellml.org/cellml/1.0#" xmlns:cmeta="http://www.cellml.org/metadata/1.0#" 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:cellml="http://www.cellml.org/cellml/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#" name="wang_model_2007" cmeta:id="wang_model_2007">

<documentation xmlns="http://cellml.org/tmp-documentation">
<article>
  <articleinfo>
  <title>Mathematical modeling and qualitative analysis of insulin therapies</title>
  <author>
    <firstname>Catherine</firstname>
          <surname>Lloyd</surname>
    <affiliation>
      <shortaffil>Auckland Bioengineering Institute, The University of Auckland</shortaffil>
    </affiliation>
  </author>
</articleinfo>
  <section id="sec_status">
    <title>Model Status</title>
    <para>
        This CellML version of the model has been checked in COR and OpenCell.  The units are consistent and the model runs to partially recreate the published results.  The infusion of glucose and insulin, and the insulin response, are the same as in the paper, but the glucose oscillations are different - this may be because the time delays in the original model can not be included in the CellML model.
          </para>
  </section>
  <sect1 id="sec_structure">
<title>Model Structure</title>

<para>
Abstract: Several insulin therapies are widely in clinical use with the basic strategy that mimics insulin secretion in a normal glucose-nsulin endocrine metabolic regulatory system. In this paper, we model the insulin therapies using a delay differential equation model. We study the dynamics of the model both qualitatively and quantitatively. The analytical results show the existence and uniqueness of a stable periodic solution that corresponds to ultradian insulin secretion oscillations. Numerically we simulate the insulin administration based on our model. The numerical simulation results are in agreement with findings of clinical studies.
</para>

<informalfigure float="0" id="fig_reaction_diagram">
<mediaobject>
  <imageobject>
    <objectinfo>
      <title>model diagram</title>
    </objectinfo>
    <imagedata fileref="wang_2007.png"/>
  </imageobject>
</mediaobject>
<caption>Glucose (blue solid line) and insulin (red dotted line) concentrations oscillate over time, with a slight time delay between the peaks as insulin spikes in response to elevated glucose levels.</caption>
</informalfigure>

<para>
The complete original paper reference is cited below:
</para>

<para>
Mathematical modeling and qualitative analysis of insulin therapies, Haiyan Wang, Jiaxu Li, and Yang Kuang, 2007, <emphasis>Mathematical Biosciences</emphasis> 210, 17-33. <ulink url="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=pubmed&amp;cmd=Retrieve&amp;dopt=AbstractPlus&amp;list_uids=17610909&amp;query_hl=1&amp;itool=pubmed_docsum">PubMed ID: 17610909</ulink>
</para>

</sect1>
</article>
</documentation> 




  <units base_units="yes" name="U"/>

  <units name="minute">
    <unit units="second" multiplier="60.0"/>
  </units>

  <units name="mg">
    <unit units="gram" prefix="milli"/>
  </units>

  <units name="mU">
    <unit units="U" prefix="milli"/>
  </units>
  
  <units name="per_litre">
    <unit units="litre" exponent="-1"/>
  </units>

  <units name="mg_per_minute">
    <unit units="gram" prefix="milli"/>
    <unit units="minute" exponent="-1"/>
  </units>

  <units name="mg_per_litre">
    <unit units="gram" prefix="milli"/>
    <unit units="litre" exponent="-1"/>
  </units>
  
  <units name="litre_per_minute">
    <unit units="litre"/>
    <unit units="minute" exponent="-1"/>
  </units>
  
  <units name="mg_per_litre_minute">
    <unit units="gram" prefix="milli"/>
    <unit units="litre" exponent="-1"/>
    <unit units="minute" exponent="-1"/>
  </units>

  <units name="mU_per_litre_minute">
    <unit units="U" prefix="milli"/>
    <unit units="litre" exponent="-1"/>
    <unit units="minute" exponent="-1"/>
  </units>
  
  <units name="mU_per_litre">
    <unit units="U" prefix="milli"/>
    <unit units="litre" exponent="-1"/>
  </units>
  
  <units name="mU_per_litre2">
    <unit units="U" prefix="milli"/>
    <unit units="litre" exponent="-2"/>
  </units>
  
  <units name="litre_per_mU">
    <unit units="litre"/>
    <unit units="U" prefix="milli" exponent="-1"/>
  </units>

  <units name="first_order_rate_constant">
    <unit units="minute" exponent="-1"/>
  </units>
  
  
  <component name="environment">
    <variable units="minute" public_interface="out" name="time"/>
  </component>
  
  <component name="G" cmeta:id="G">
    <rdf:RDF>
      <rdf:Description rdf:about="#G">
        <dc:title>G</dc:title>
        <dcterms:alternative>glucose</dcterms:alternative>  
      </rdf:Description>
    </rdf:RDF>  
    <variable units="mg_per_litre" public_interface="out" name="G" initial_value="150.0"/>
    
    <variable units="mg_per_litre_minute" name="Gin"/>
    
    <variable units="mg_per_minute" public_interface="in" name="f2_G"/>
    <variable units="per_litre" public_interface="in" name="f3_G"/>
    <variable units="mg_per_minute" public_interface="in" name="f4_I"/>
    <variable units="mg_per_minute" public_interface="in" name="f5_I"/>
    <variable units="minute" public_interface="in" name="time"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci>time</ci>
          </bvar> 
          <ci> G </ci> 
        </apply>
        <apply>
          <minus/>
          <apply>
            <plus/>
            <ci> Gin </ci>
            <apply>
              <times/>
			  <cn cellml:units="per_litre"> 1 </cn>
			  <ci> f5_I </ci>
			</apply>
          </apply>
          <apply>
            <plus/>
			<apply>
              <times/>
			  <cn cellml:units="per_litre"> 1 </cn>
              <ci> f2_G </ci>
			</apply>
            <apply>
              <times/>
              <ci> f3_G </ci>
              <ci> f4_I </ci>
            </apply>
          </apply>
        </apply>
      </apply>

      <apply>
        <eq/> 
        <ci> Gin </ci> 
        <piecewise>
          <piece>
            <apply>
              <plus/>
              <cn cellml:units="mg_per_litre_minute"> 0.05 </cn>
              <apply>
                <times/>
                <apply>
                  <divide/>
                  <cn cellml:units="mg_per_litre_minute"> 5.0 </cn>
                  <cn cellml:units="minute"> 15.0 </cn>
                </apply>
                <ci> time </ci>
              </apply>
            </apply>
            <apply>
              <and/>
              <apply>
                <lt/>
                <ci> time </ci>
                <cn cellml:units="minute"> 15.0 </cn>
              </apply>
              <apply>
                <geq/>
                <ci> time </ci>
                <cn cellml:units="minute"> 0.0 </cn>
              </apply>
            </apply>
          </piece>
          <piece>
            <apply>
              <plus/>
              <cn cellml:units="mg_per_litre_minute"> 0.05 </cn>
              <apply>
                <times/>
                <cn cellml:units="mg_per_litre_minute"> 5.0 </cn>
                <apply>
                  <divide/>
                  <apply>
                    <minus/>
                    <cn cellml:units="minute"> 45.0 </cn>
                    <ci> time </ci>
                  </apply>
                  <apply>
                    <minus/>
                    <cn cellml:units="minute"> 45.0 </cn>
                    <cn cellml:units="minute"> 15.0 </cn>
                  </apply>
                </apply>
              </apply>
            </apply>
            <apply>
              <and/>
              <apply>
                <lt/>
                <ci> time </ci>
                <cn cellml:units="minute"> 45.0 </cn>
              </apply>
              <apply>
                <geq/>
                <ci> time </ci>
                <cn cellml:units="minute"> 15.0 </cn>
              </apply>
            </apply>
          </piece>
          <piece>
            <cn cellml:units="mg_per_litre_minute"> 0.05 </cn>
            <apply>
              <and/>
              <apply>
                <leq/>
                <ci> time </ci>
                <cn cellml:units="minute"> 240.0 </cn>
              </apply>
              <apply>
                <geq/>
                <ci> time </ci>
                <cn cellml:units="minute"> 45.0 </cn>
              </apply>
            </apply>
          </piece>
          <piece>
            <apply>
              <plus/>
              <cn cellml:units="mg_per_litre_minute"> 0.05 </cn>
              <apply>
                <times/>
                <apply>
                  <divide/>
                  <cn cellml:units="mg_per_litre_minute"> 5.0 </cn>
                  <cn cellml:units="minute"> 255.0 </cn>
                </apply>
                <ci> time </ci>
              </apply>
            </apply>
            <apply>
              <and/>
              <apply>
                <gt/>
                <ci> time </ci>
                <cn cellml:units="minute"> 240.0 </cn>
              </apply>
              <apply>
                <lt/>
                <ci> time </ci>
                <cn cellml:units="minute"> 255.0 </cn>
              </apply>
            </apply>
          </piece>
          <piece>
            <apply>
              <plus/>
              <cn cellml:units="mg_per_litre_minute"> 0.05 </cn>
              <apply>
                <times/>
                <cn cellml:units="mg_per_litre_minute"> 5.0 </cn>
                <apply>
                  <divide/>
                  <apply>
                    <minus/>
                    <cn cellml:units="minute"> 285.0 </cn>
                    <ci> time </ci>
                  </apply>
                  <apply>
                    <minus/>
                    <cn cellml:units="minute"> 285.0 </cn>
                    <cn cellml:units="minute"> 255.0 </cn>
                  </apply>
                </apply>
              </apply>
            </apply>
            <apply>
              <and/>
              <apply>
                <lt/>
                <ci> time </ci>
                <cn cellml:units="minute"> 285.0 </cn>
              </apply>
              <apply>
                <geq/>
                <ci> time </ci>
                <cn cellml:units="minute"> 255.0 </cn>
              </apply>
            </apply>
          </piece>
          <piece>
            <cn cellml:units="mg_per_litre_minute"> 0.05 </cn>
            <apply>
              <and/>
              <apply>
                <leq/>
                <ci> time </ci>
                <cn cellml:units="minute"> 480.0 </cn>
              </apply>
              <apply>
                <geq/>
                <ci> time </ci>
                <cn cellml:units="minute"> 285.0 </cn>
              </apply>
            </apply>
          </piece>
          <piece>
            <apply>
              <plus/>
              <cn cellml:units="mg_per_litre_minute"> 0.05 </cn>
              <apply>
                <times/>
                <apply>
                  <divide/>
                  <cn cellml:units="mg_per_litre_minute"> 5.0 </cn>
                  <cn cellml:units="minute"> 495.0 </cn>
                </apply>
                <ci> time </ci>
              </apply>
            </apply>
            <apply>
              <and/>
              <apply>
                <lt/>
                <ci> time </ci>
                <cn cellml:units="minute"> 495.0 </cn>
              </apply>
              <apply>
                <gt/>
                <ci> time </ci>
                <cn cellml:units="minute"> 480.0 </cn>
              </apply>
            </apply>
          </piece>
          <piece>
            <apply>
              <plus/>
              <cn cellml:units="mg_per_litre_minute"> 0.05 </cn>
              <apply>
                <times/>
                <cn cellml:units="mg_per_litre_minute"> 5.0 </cn>
                <apply>
                  <divide/>
                  <apply>
                    <minus/>
                    <cn cellml:units="minute"> 525.0 </cn>
                    <ci> time </ci>
                  </apply>
                  <apply>
                    <minus/>
                    <cn cellml:units="minute"> 525.0 </cn>
                    <cn cellml:units="minute"> 495.0 </cn>
                  </apply>
                </apply>
              </apply>
            </apply>
            <apply>
              <and/>
              <apply>
                <lt/>
                <ci> time </ci>
                <cn cellml:units="minute"> 525.0 </cn>
              </apply>
              <apply>
                <geq/>
                <ci> time </ci>
                <cn cellml:units="minute"> 495.0 </cn>
              </apply>
            </apply>
          </piece>
          <piece>
            <cn cellml:units="mg_per_litre_minute"> 0.05 </cn>
            <apply>
              <and/>
              <apply>
                <leq/>
                <ci> time </ci>
                <cn cellml:units="minute"> 720.0 </cn>
              </apply>
              <apply>
                <geq/>
                <ci> time </ci>
                <cn cellml:units="minute"> 525.0 </cn>
              </apply>
            </apply>
          </piece>
          <piece>
            <apply>
              <plus/>
              <cn cellml:units="mg_per_litre_minute"> 0.05 </cn>
              <apply>
                <times/>
                <apply>
                  <divide/>
                  <cn cellml:units="mg_per_litre_minute"> 5.0 </cn>
                  <cn cellml:units="minute"> 735.0 </cn>
                </apply>
                <ci> time </ci>
              </apply>
            </apply>
            <apply>
              <and/>
              <apply>
                <gt/>
                <ci> time </ci>
                <cn cellml:units="minute"> 720.0 </cn>
              </apply>
              <apply>
                <lt/>
                <ci> time </ci>
                <cn cellml:units="minute"> 735.0 </cn>
              </apply>
            </apply>
          </piece>
          <piece>
            <apply>
              <plus/>
              <cn cellml:units="mg_per_litre_minute"> 0.05 </cn>
              <apply>
                <times/>
                <cn cellml:units="mg_per_litre_minute"> 5.0 </cn>
                <apply>
                  <divide/>
                  <apply>
                    <minus/>
                    <cn cellml:units="minute"> 765.0 </cn>
                    <ci> time </ci>
                  </apply>
                  <apply>
                    <minus/>
                    <cn cellml:units="minute"> 765.0 </cn>
                    <cn cellml:units="minute"> 735.0 </cn>
                  </apply>
                </apply>
              </apply>
            </apply>
            <apply>
              <and/>
              <apply>
                <lt/>
                <ci> time </ci>
                <cn cellml:units="minute"> 765.0 </cn>
              </apply>
              <apply>
                <geq/>
                <ci> time </ci>
                <cn cellml:units="minute"> 735.0 </cn>
              </apply>
            </apply>
          </piece>
          <piece>
            <cn cellml:units="mg_per_litre_minute"> 0.05 </cn>
            <apply>
              <and/>
              <apply>
                <leq/>
                <ci> time </ci>
                <cn cellml:units="minute"> 960.0 </cn>
              </apply>
              <apply>
                <geq/>
                <ci> time </ci>
                <cn cellml:units="minute"> 765.0 </cn>
              </apply>
            </apply>
          </piece>
        </piecewise>
      </apply>    
    </math>
  </component>
  
  <component name="I" cmeta:id="I">
    <rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
      <rdf:Description rdf:about="#I">
        <dc:title>I</dc:title>
        <dcterms:alternative>insulin</dcterms:alternative>  
      </rdf:Description>
    </rdf:RDF>  
    <variable units="mU_per_litre" public_interface="out" name="I" initial_value="5.0"/>
    
    <variable units="mU_per_litre_minute" name="Iin"/>
    
    <variable units="first_order_rate_constant" public_interface="in" name="di"/>
    <variable units="minute" public_interface="in" name="time"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci>time</ci>
          </bvar> 
          <ci> I </ci> 
        </apply>
        <apply>
          <minus/>
          <ci> Iin </ci>
          <apply>
            <times/>
            <ci> di </ci>
            <ci> I </ci>
          </apply>
        </apply>
      </apply>

      <apply>
        <eq/> 
        <ci> Iin </ci> 
        <piecewise>
          <piece>
            <cn cellml:units="mU_per_litre_minute"> 0.25 </cn>
            <apply>
              <and/>
              <apply>
                <leq/>
                <ci> time </ci>
                <cn cellml:units="minute"> 30.0 </cn>
              </apply>
              <apply>
                <geq/>
                <ci> time </ci>
                <cn cellml:units="minute"> 0.0 </cn>
              </apply>
            </apply>
          </piece>
          <piece>
            <apply>
              <plus/>
              <cn cellml:units="mU_per_litre_minute"> 0.25 </cn>
              <apply>
                <times/>
                <cn cellml:units="mU_per_litre_minute"> 1.0 </cn>
                <apply>
                  <plus/>
                  <cn cellml:units="dimensionless"> 1.0 </cn>
                  <apply>
                    <divide/>
                    <apply>
                      <minus/>
                      <ci> time </ci>
                      <cn cellml:units="minute"> 120.0 </cn>
                    </apply>
                    <cn cellml:units="minute"> 90.0 </cn>
                  </apply>
                </apply>
              </apply>
            </apply>
            <apply>
              <and/>
              <apply>
                <lt/>
                <ci> time </ci>
                <cn cellml:units="minute"> 120.0 </cn>
              </apply>
              <apply>
                <geq/>
                <ci> time </ci>
                <cn cellml:units="minute"> 30.0 </cn>
              </apply>
            </apply>
          </piece>
          <piece>
            <apply>
              <plus/>
              <cn cellml:units="mU_per_litre_minute"> 0.25 </cn>
              <apply>
                <times/>
                <cn cellml:units="mU_per_litre_minute"> 1.0 </cn>
                <apply>
                  <minus/>
                  <cn cellml:units="dimensionless"> 1.0 </cn>
                  <apply>
                    <times/>
                    <cn cellml:units="dimensionless"> 0.5 </cn>
                    <apply>
                      <divide/>
                      <apply>
                        <minus/>
                        <ci> time </ci>
                        <cn cellml:units="minute"> 120.0 </cn>
                      </apply>
                      <cn cellml:units="minute"> 120.0 </cn>
                    </apply>
                  </apply>
                </apply>
              </apply>
            </apply>
            <apply>
              <and/>
              <apply>
                <lt/>
                <ci> time </ci>
                <cn cellml:units="minute"> 240.0 </cn>
              </apply>
              <apply>
                <geq/>
                <ci> time </ci>
                <cn cellml:units="minute"> 120.0 </cn>
              </apply>
            </apply>
          </piece>
          <piece>
            <apply>
              <plus/>
              <cn cellml:units="mU_per_litre_minute"> 0.25 </cn>
              <apply>
                <times/>
                <cn cellml:units="mU_per_litre_minute"> 0.5 </cn>
                <apply>
                  <minus/>
                  <cn cellml:units="dimensionless"> 1.0 </cn>
                  <apply>
                    <divide/>
                    <apply>
                      <minus/>
                      <ci> time </ci>
                      <cn cellml:units="minute"> 240.0 </cn>
                    </apply>
                    <cn cellml:units="minute"> 240.0 </cn>
                  </apply>
                </apply>
              </apply>
            </apply>
            <apply>
              <and/>
              <apply>
                <leq/>
                <ci> time </ci>
                <cn cellml:units="minute"> 480.0 </cn>
              </apply>
              <apply>
                <geq/>
                <ci> time </ci>
                <cn cellml:units="minute"> 240.0 </cn>
              </apply>
            </apply>
          </piece>
          <piece>
            <cn cellml:units="mU_per_litre_minute"> 0.25 </cn>
            <apply>
              <and/>
              <apply>
                <lt/>
                <ci> time </ci>
                <cn cellml:units="minute"> 530.0 </cn>
              </apply>
              <apply>
                <gt/>
                <ci> time </ci>
                <cn cellml:units="minute"> 480.0 </cn>
              </apply>
            </apply>
          </piece>
          <piece>
            <apply>
              <plus/>
              <cn cellml:units="mU_per_litre_minute"> 0.25 </cn>
              <apply>
                <times/>
                <cn cellml:units="mU_per_litre_minute"> 1.0 </cn>
                <apply>
                  <plus/>
                  <cn cellml:units="dimensionless"> 1.0 </cn>
                  <apply>
                    <divide/>
                    <apply>
                      <minus/>
                      <ci> time </ci>
                      <cn cellml:units="minute"> 620.0 </cn>
                    </apply>
                    <cn cellml:units="minute"> 90.0 </cn>
                  </apply>
                </apply>
              </apply>
            </apply>
            <apply>
              <and/>
              <apply>
                <lt/>
                <ci> time </ci>
                <cn cellml:units="minute"> 620.0 </cn>
              </apply>
              <apply>
                <geq/>
                <ci> time </ci>
                <cn cellml:units="minute"> 530.0 </cn>
              </apply>
            </apply>
          </piece>
          <piece>
            <apply>
              <plus/>
              <cn cellml:units="mU_per_litre_minute"> 0.25 </cn>
              <apply>
                <times/>
                <cn cellml:units="mU_per_litre_minute"> 1.0 </cn>
                <apply>
                  <minus/>
                  <cn cellml:units="dimensionless"> 1.0 </cn>
                  <apply>
                    <times/>
                    <cn cellml:units="dimensionless"> 0.5 </cn>
                    <apply>
                      <divide/>
                      <apply>
                        <minus/>
                        <ci> time </ci>
                        <cn cellml:units="minute"> 620.0 </cn>
                      </apply>
                      <cn cellml:units="minute"> 120.0 </cn>
                    </apply>
                  </apply>
                </apply>
              </apply>
            </apply>
            <apply>
              <and/>
              <apply>
                <lt/>
                <ci> time </ci>
                <cn cellml:units="minute"> 740.0 </cn>
              </apply>
              <apply>
                <geq/>
                <ci> time </ci>
                <cn cellml:units="minute"> 620.0 </cn>
              </apply>
            </apply>
          </piece>
          <piece>
            <apply>
              <plus/>
              <cn cellml:units="mU_per_litre_minute"> 0.25 </cn>
              <apply>
                <times/>
                <cn cellml:units="mU_per_litre_minute"> 0.5 </cn>
                <apply>
                  <minus/>
                  <cn cellml:units="dimensionless"> 1.0 </cn>
                  <apply>
                    <divide/>
                    <apply>
                      <minus/>
                      <ci> time </ci>
                      <cn cellml:units="minute"> 740.0 </cn>
                    </apply>
                    <cn cellml:units="minute"> 240.0 </cn>
                  </apply>
                </apply>
              </apply>
            </apply>
            <apply>
              <and/>
              <apply>
                <leq/>
                <ci> time </ci>
                <cn cellml:units="minute"> 980.0 </cn>
              </apply>
              <apply>
                <geq/>
                <ci> time </ci>
                <cn cellml:units="minute"> 740.0 </cn>
              </apply>
            </apply>
          </piece>
        </piecewise>
      </apply>        
    </math>
  </component>

  <component name="f2_G">  
    <variable units="mg_per_minute" public_interface="out" name="f2_G"/>
    
    <variable units="mg_per_litre" public_interface="in" name="G"/>
    <variable units="mg_per_litre" public_interface="in" name="C2"/>
    <variable units="mg_per_minute" public_interface="in" name="Ub"/>
	<variable units="litre" public_interface="in" name="Vg"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <ci> f2_G </ci> 
        <apply>
          <times/>
          <ci> Ub </ci>
          <apply>
            <minus/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
              <exp/>
              <apply>
                <divide/>
                <apply>
                  <minus/>
                  <ci> G </ci>
                </apply>
                <apply>
                  <times/>
                  <ci> C2 </ci>
                  <ci> Vg </ci>
				  <cn cellml:units="per_litre"> 1.0 </cn>
                </apply>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>    
    </math>
  </component>

  <component name="f3_G">  
    <variable units="per_litre" public_interface="out" name="f3_G"/>
    
    <variable units="mg_per_litre" public_interface="in" name="G"/>
    <variable units="mg_per_litre" public_interface="in" name="C3"/>
    <variable units="litre" public_interface="in" name="Vg"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/>
        <ci> f3_G </ci> 
        <apply>
          <divide/>
          <ci> G </ci>
          <apply>
            <times/>
            <ci> C3 </ci>
            <ci> Vg </ci>
          </apply>
        </apply>
      </apply>    
    </math>
  </component>

  <component name="f4_I">  
    <variable units="mg_per_minute" public_interface="out" name="f4_I"/>
    
    <variable units="mU_per_litre" public_interface="in" name="I"/>
    <variable units="mU_per_litre" public_interface="in" name="C4"/>
    <variable units="litre" public_interface="in" name="Vi"/>
    <variable units="mg_per_minute" public_interface="in" name="U0"/>
    <variable units="mg_per_minute" public_interface="in" name="Um"/>
    <variable units="dimensionless" public_interface="in" name="beta"/>
    <variable units="minute" public_interface="in" name="ti"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <ci> f4_I </ci> 
        <apply>
          <plus/>
          <ci> U0 </ci>
          <apply>
            <divide/>
            <apply>
              <minus/>
              <ci> Um </ci>
              <ci> U0 </ci>
            </apply>
            <apply>
              <plus/>
              <cn cellml:units="dimensionless"> 1.0 </cn>
              <apply>
                <exp/>
                <apply>
                  <times/>
                  <apply>
                    <minus/>
                    <ci> beta </ci>
                  </apply>
                  <apply>
                    <ln/>
                    <apply>
                      <divide/>
                      <ci> I </ci>
                      <apply>
                        <times/>
                        <ci> C4 </ci>
                        <apply>
                          <plus/>
						  <apply>
                            <divide/>
                            <cn cellml:units="litre"> 1.0 </cn>
                            <ci> Vi </ci>
					      </apply>
                          <apply>
                            <divide/>
                            <cn cellml:units="litre"> 1.0 </cn>
                            <apply>
                              <times/>
                              <cn cellml:units="litre_per_minute"> 0.2 </cn>
                              <ci> ti </ci>
                            </apply>
                          </apply>
                        </apply>
                      </apply>
                    </apply>
                  </apply>
                </apply>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>    
    </math>
  </component>

  <component name="f5_I">  
    <variable units="mg_per_minute" public_interface="out" name="f5_I"/>
    
	<variable units="litre" public_interface="in" name="Vp"/>
    <variable units="mU_per_litre" public_interface="in" name="I"/>
    <variable units="mU_per_litre" public_interface="in" name="C5"/>
    <variable units="mg_per_minute" public_interface="in" name="Rg"/>
    <variable units="litre_per_mU" public_interface="in" name="alpha"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <ci> f5_I </ci> 
        <apply>
          <divide/>
          <ci> Rg </ci>
          <apply>
            <plus/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
              <exp/>
              <apply>
                <times/>
                <ci> alpha </ci>
		        <apply>
                  <divide/>
                  <ci> I </ci>
				  <apply>
                    <times/>
					<cn cellml:units="litre_per_mU"> 1.0 </cn>
					<apply>
                      <minus/>
                      <apply>
                        <times/>
					    <cn cellml:units="mU_per_litre2"> 1.0 </cn>
					    <ci> Vp </ci>
					  </apply>
	                  <ci> C5 </ci>
					</apply>
				  </apply>
                </apply>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>    
    </math>
  </component>
 
  <component name="model_parameters"> 
    <variable units="litre" public_interface="out" name="Vg" initial_value="10.0"/>
	<variable units="litre" public_interface="out" name="Vp" initial_value="3.0"/>
    <variable units="litre" public_interface="out" name="Vi" initial_value="11.0"/>
    <variable units="minute" public_interface="out" name="ti" initial_value="100.0"/>
    <variable units="mg_per_minute" public_interface="out" name="U0" initial_value="40.0"/>
    <variable units="mg_per_minute" public_interface="out" name="Um" initial_value="940.0"/>
	<variable units="litre_per_mU" public_interface="out" name="alpha" initial_value="0.29"/>
    <variable units="dimensionless" public_interface="out" name="beta" initial_value="1.77"/>
    <variable units="mg_per_litre" public_interface="out" name="C2" initial_value="144"/>
    <variable units="mg_per_litre" public_interface="out" name="C3" initial_value="1000.0"/>
    <variable units="mU_per_litre" public_interface="out" name="C4" initial_value="80.0"/>
    <variable units="mU_per_litre" public_interface="out" name="C5" initial_value="26.0"/>
    <variable units="mg_per_minute" public_interface="out" name="Ub" initial_value="72.0"/>
    <variable units="mg_per_minute" public_interface="out" name="Rg" initial_value="180.0"/>
    <variable units="first_order_rate_constant" public_interface="out" name="di" initial_value="0.0107"/>
  </component>
  
  <connection>
    <map_components component_2="environment" component_1="G"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="environment" component_1="I"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>

  <connection>
    <map_components component_2="f2_G" component_1="G"/>
    <map_variables variable_2="G" variable_1="G"/>
    <map_variables variable_2="f2_G" variable_1="f2_G"/>
  </connection> 

  <connection>
    <map_components component_2="f3_G" component_1="G"/>
    <map_variables variable_2="G" variable_1="G"/>
    <map_variables variable_2="f3_G" variable_1="f3_G"/>
  </connection>

  <connection>
    <map_components component_2="f4_I" component_1="G"/>
    <map_variables variable_2="f4_I" variable_1="f4_I"/>
  </connection> 

  <connection>
    <map_components component_2="f5_I" component_1="G"/>
    <map_variables variable_2="f5_I" variable_1="f5_I"/>
  </connection>

  <connection>
    <map_components component_2="model_parameters" component_1="I"/>
    <map_variables variable_2="di" variable_1="di"/>
  </connection>

  <connection>
    <map_components component_2="model_parameters" component_1="f2_G"/>
    <map_variables variable_2="C2" variable_1="C2"/>
    <map_variables variable_2="Ub" variable_1="Ub"/>
	<map_variables variable_2="Vg" variable_1="Vg"/>
  </connection>

  <connection>
    <map_components component_2="model_parameters" component_1="f3_G"/>
    <map_variables variable_2="C3" variable_1="C3"/>
    <map_variables variable_2="Vg" variable_1="Vg"/>
  </connection>

  <connection>
    <map_components component_2="f4_I" component_1="I"/>
    <map_variables variable_2="I" variable_1="I"/>
  </connection>

  <connection>
    <map_components component_2="model_parameters" component_1="f4_I"/>
    <map_variables variable_2="C4" variable_1="C4"/>
    <map_variables variable_2="Vi" variable_1="Vi"/>
    <map_variables variable_2="Um" variable_1="Um"/>
    <map_variables variable_2="U0" variable_1="U0"/>
    <map_variables variable_2="beta" variable_1="beta"/>
    <map_variables variable_2="ti" variable_1="ti"/>
  </connection> 
 
  <connection>
    <map_components component_2="f5_I" component_1="I"/>
    <map_variables variable_2="I" variable_1="I"/>
  </connection>

  <connection>
    <map_components component_2="model_parameters" component_1="f5_I"/>
    <map_variables variable_2="C5" variable_1="C5"/>
    <map_variables variable_2="Rg" variable_1="Rg"/>
    <map_variables variable_2="alpha" variable_1="alpha"/>
	<map_variables variable_2="Vp" variable_1="Vp"/>
  </connection> 

<rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
  <rdf:Bag rdf:about="rdf:#$KN9jO1">
    <rdf:li>endocrine</rdf:li>
    <rdf:li>glucose</rdf:li>
    <rdf:li>insulin</rdf:li>
    <rdf:li>diabetes</rdf:li>
  </rdf:Bag>
  <rdf:Seq rdf:about="rdf:#8c81b56f-e8a4-4897-9a43-8f4659b706f0">
    <rdf:li rdf:resource="rdf:#63551c37-5b7e-45da-9486-d2926af2a9aa"/>
    <rdf:li rdf:resource="rdf:#c5067023-1971-43e0-a76b-1ce7339b0873"/>
    <rdf:li rdf:resource="rdf:#7d81d12f-74fa-4a87-9586-b8cc777615e4"/>
  </rdf:Seq>
  <rdf:Description rdf:about="rdf:#7d81d12f-74fa-4a87-9586-b8cc777615e4">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#aeb3bc7a-5650-44f8-8ef7-7faddea5ea71"/>
  </rdf:Description>
  <rdf:Description rdf:about="#wang_model_2007">
    <ns7:simulation rdf:resource="rdf:#$ZN9jO1"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#140579b6-5f11-47f2-ab5a-e0319e805f81">
    <dc:creator rdf:resource="rdf:#8c81b56f-e8a4-4897-9a43-8f4659b706f0"/>
    <dc:title>Mathematical modeling and qualitative analysis of insulin therapies</dc:title>
    <bqs:volume/>
    <bqs:first_page/>
    <bqs:Journal rdf:resource="rdf:#3b38a479-df71-4fad-bba2-051b275d9016"/>
    <dcterms:issued rdf:resource="rdf:#58baaf52-e9b0-4c90-8f19-8804ed7c7cf8"/>
    <bqs:last_page/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#ad1cc13d-bab8-4a73-b8e0-1163b54ecf10">
    <dcterms:modified rdf:resource="rdf:#67f7dafd-4db9-4141-9acc-fbac4c3a82c4"/>
    <rdf:value>I've made the units consistent and have fixed the piecewise equations.</rdf:value>
    <cmeta:modifier rdf:resource="rdf:#e22387f3-0f80-4a1a-bd60-76397a8d52fe"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#3aecf87b-50fd-4090-8bb9-55ec60b59c67">
    <vCard:Given>Jiaxu</vCard:Given>
    <vCard:Family>Li</vCard:Family>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#$+N9jO1">
    <ns7:endingValue>1000</ns7:endingValue>
    <ns7:maximumStepSize>0.001</ns7:maximumStepSize>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#3ea0e661-7509-487a-848d-7c6e6a3b89b3">
    <vCard:Given>Catherine</vCard:Given>
    <vCard:Family>Lloyd</vCard:Family>
    <vCard:Other>May</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#cae8c6a5-a7a9-4a56-b04d-c2a5c8d3494a">
    <dcterms:modified rdf:resource="rdf:#955ff128-e109-42c7-96bc-f553021d4b0b"/>
    <rdf:value>I've made the units consistent and have fixed the piecewise equations.</rdf:value>
    <cmeta:modifier rdf:resource="rdf:#4ebba1d0-927f-4c2a-bd27-aa32568cf8d8"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#63551c37-5b7e-45da-9486-d2926af2a9aa">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#f498a2ad-3b2c-42bd-a55d-b9569dde42de"/>
  </rdf:Description>
  <rdf:Description rdf:about="#wang_model_2007">
    <dc:title>Wang, Li, and Kuang's 2007 mathematical model and qualitative analysis of insulin therapies.</dc:title>
    <cmeta:comment rdf:resource="rdf:#40096c26-98e9-4120-ab7b-8e3f8a393f30"/>
    <bqs:reference rdf:resource="rdf:#$IN9jO1"/>
    <bqs:reference rdf:resource="rdf:#d86ffac5-606e-4ec4-929b-983864dfc1af"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#3b38a479-df71-4fad-bba2-051b275d9016">
    <dc:title>Mathematical Biosciences</dc:title>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#b58424e2-cb8e-4f93-b4ab-55c7db6a0e2f">
    <vCard:Given>Catherine</vCard:Given>
    <vCard:Family>Lloyd</vCard:Family>
    <vCard:Other>May</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#cd11322d-0ee4-4088-abc5-4336e491e267">
    <rdf:type rdf:resource="http://imc.org/vCard/3.0#internet"/>
    <rdf:value>c.lloyd@auckland.ac.nz</rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#$ZN9jO1">
    <ns7:boundIntervals rdf:resource="rdf:#$.N9jO1"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#58baaf52-e9b0-4c90-8f19-8804ed7c7cf8">
    <dcterms:W3CDTF>2007-05-29</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#3b75e81a-cf16-4586-ac81-37419d07105b">
    <vCard:Given>Catherine</vCard:Given>
    <vCard:Family>Lloyd</vCard:Family>
    <vCard:Other>May</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#aeb3bc7a-5650-44f8-8ef7-7faddea5ea71">
    <vCard:Given>Yang</vCard:Given>
    <vCard:Family>Kuang</vCard:Family>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#2be7ea16-285e-443b-beb9-3d7fc73457fa">
    <vCard:Orgname>The University of Auckland</vCard:Orgname>
    <vCard:Orgunit>The Bioengineering Institute</vCard:Orgunit>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#$IN9jO1">
    <dc:subject rdf:resource="rdf:#$JN9jO1"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#7b94ea6e-f0cd-46c3-ab19-4eb845f865b5">
    <dc:creator rdf:resource="rdf:#97438361-466c-47e6-8fd4-87d577553a72"/>
    <rdf:value>The model runs in PCEnv to partially recreate the published results.  The inability of the CellML model to recreate the glucose oscillations by the original model may be due the lack of time delays in the CellML model.</rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#e22387f3-0f80-4a1a-bd60-76397a8d52fe">
    <vCard:N rdf:resource="rdf:#b58424e2-cb8e-4f93-b4ab-55c7db6a0e2f"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#955ff128-e109-42c7-96bc-f553021d4b0b">
    <dcterms:W3CDTF>2007-09-06T12:36:12+12:00</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#$JN9jO1">
    <bqs:subject_type>keyword</bqs:subject_type>
    <rdf:value rdf:resource="rdf:#$KN9jO1"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#97438361-466c-47e6-8fd4-87d577553a72">
    <vCard:FN>Catherine Lloyd</vCard:FN>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#7e801b9e-1b67-4167-ac23-0328a14a2df2">
    <dcterms:W3CDTF>2007-07-10T00:00:00+00:00</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#$.N9jO1">
    <rdf:rest rdf:resource="http://www.w3.org/1999/02/22-rdf-syntax-ns#nil"/>
    <rdf:first rdf:resource="rdf:#$+N9jO1"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#67f7dafd-4db9-4141-9acc-fbac4c3a82c4">
    <dcterms:W3CDTF>2007-09-06T12:36:12+12:00</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#c5067023-1971-43e0-a76b-1ce7339b0873">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#3aecf87b-50fd-4090-8bb9-55ec60b59c67"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#40096c26-98e9-4120-ab7b-8e3f8a393f30">
    <dc:creator rdf:resource="rdf:#af3c82b9-0361-48c0-a017-024bcaace00b"/>
    <rdf:value>This is a CellML description of Wang, Li, and Kuang's 2007 mathematical model and qualitative analysis of insulin therapies.</rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#af3c82b9-0361-48c0-a017-024bcaace00b">
    <vCard:FN>Catherine Lloyd</vCard:FN>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#f498a2ad-3b2c-42bd-a55d-b9569dde42de">
    <vCard:Given>Haiyan</vCard:Given>
    <vCard:Family>Wang</vCard:Family>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#032b8a27-172e-409a-bccb-39998585e301">
    <vCard:ORG rdf:resource="rdf:#2be7ea16-285e-443b-beb9-3d7fc73457fa"/>
    <vCard:EMAIL rdf:resource="rdf:#cd11322d-0ee4-4088-abc5-4336e491e267"/>
    <vCard:N rdf:resource="rdf:#3b75e81a-cf16-4586-ac81-37419d07105b"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#4ebba1d0-927f-4c2a-bd27-aa32568cf8d8">
    <vCard:N rdf:resource="rdf:#3ea0e661-7509-487a-848d-7c6e6a3b89b3"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#d86ffac5-606e-4ec4-929b-983864dfc1af">
    <bqs:Pubmed_id>17610909</bqs:Pubmed_id>
    <bqs:JournalArticle rdf:resource="rdf:#140579b6-5f11-47f2-ab5a-e0319e805f81"/>
  </rdf:Description>
  <rdf:Description rdf:about="">
    <dc:publisher>The University of Auckland, Bioengineering Institute</dc:publisher>
    <cmeta:comment rdf:resource="rdf:#7b94ea6e-f0cd-46c3-ab19-4eb845f865b5"/>
    <dcterms:created rdf:resource="rdf:#7e801b9e-1b67-4167-ac23-0328a14a2df2"/>
    <dc:creator rdf:resource="rdf:#032b8a27-172e-409a-bccb-39998585e301"/>
    <cmeta:modification rdf:resource="rdf:#ad1cc13d-bab8-4a73-b8e0-1163b54ecf10"/>
    <cmeta:modification rdf:resource="rdf:#cae8c6a5-a7a9-4a56-b04d-c2a5c8d3494a"/>
  </rdf:Description>
</rdf:RDF>
</model>