nr. 477 - 2011 - Institut for Natur, Systemer og Modeller (NSM)
nr. 477 - 2011 - Institut for Natur, Systemer og Modeller (NSM)
nr. 477 - 2011 - Institut for Natur, Systemer og Modeller (NSM)
Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
6.1 Fixed Points of the IM Model 57<br />
✻<br />
Ma(t)<br />
❅❘<br />
˙<br />
Ma < 0<br />
˙<br />
Ba > 0<br />
Ma<br />
˙ < 0<br />
Ba<br />
˙ < 0<br />
❄✛<br />
✠<br />
˙<br />
Ma = 0<br />
✲<br />
✒<br />
˙<br />
Ma > 0<br />
˙<br />
Ba > 0<br />
✻<br />
❅■<br />
Ma<br />
˙ > 0<br />
˙<br />
Ba < 0<br />
˙<br />
Ba = 0<br />
✲<br />
Asymp1 Asymp2 Ba(t)<br />
Figure 6.3 Sketch of the nullclines and some representative arrows that indicate the flow <strong>for</strong><br />
the IM with dM/dt = 0. For the healthy rest state, (Ba, Ma) = (0, 0), to be stable the initial<br />
slope of the Ma-nullcline must be less than the initial slope of the Ba-nullcline. Then by<br />
assuming that the nullcline of Ma reaches its asymptotic point be<strong>for</strong>e the nullcline of Ba<br />
introduces a nontrivial fixed point at the intersection of the nullclines characteristic <strong>for</strong> the<br />
NOD dynamics. By the knowledge of the stability of the healthy rest state it is a simple<br />
matter to determine the direction of the slope field. This reveals that the nontrivial fixed<br />
point must be a saddle.<br />
are used. We have already concluded that the initial slope of the Ma nullcline must<br />
be less than the initial slope of the Ba nullcline <strong>for</strong> the healthy rest state to be stable.<br />
Thus this must still be the case when using Balb/c parameters. However unlike<br />
the NOD configuration there should be no nontrivial fixed point which implies that<br />
B a(Asymp1) > B a(Asymp2) so no intersection of the nullclines occur. This leads to the<br />
nullclines and slope field as seen on the sketch in figure 6.4. Here it is evident that any<br />
stimulated inflammation will be nonpermanent.<br />
Be<strong>for</strong>e we continue with the next subsection we would like to make a brief remark.<br />
Usually at the point where one is conducting a phase plane analysis the model under<br />
scrutiny has been recast in a dimensionless <strong>for</strong>m, since this eases the job of deciding<br />
if some terms in the model can be disregarded, if some of the parameters are (relatively)<br />
significantly smaller than others or some parameters hold a greater impact on<br />
the behavior of the dynamic. This is often helpful if some or none of the parameters are