Author Topic: Night Sections  (Read 14460 times)

Victor Ivrii

  • Administrator
  • Elder Member
  • *****
  • Posts: 2607
  • Karma: 0
    • View Profile
    • Personal website of Victor Ivrii
Night Sections
« on: April 03, 2013, 08:07:49 PM »
\begin{equation*}
\left\{\begin{aligned}
&\frac{dx}{dt} = x(1.5-0.5x-y),\\
&\frac{dy}{dt} = y(2-y-1.125x).
\end{aligned}\right.
\end{equation*}


This problem can be interpreted as describing the interaction of two species with population densities $x$ and $y$.


(b)  Find the critical points.

(c)  For each critical point find the corresponding linear system.  Find the eigenvalues and eigenvectors of the linear system;   classify each critical point as to type, and determine whether  it is asymptotically stable, stable, or unstable.

(d) Sketch the trajectories in the neighborhood of each critical point.

(e) [Bonus]  Find, if possible, solution in the form $H(x,y)=C$ and sketch the phase portrait.


Frank (Yi) Gao

  • Newbie
  • *
  • Posts: 1
  • Karma: 1
    • View Profile
Re: Night Sections
« Reply #1 on: April 03, 2013, 08:49:26 PM »
(b) Finding the critical points means solving:
 \begin{equation*}
  0 = x(1.5 - 0.5x - y),
  0 = y(2 - x - 1.125y)
 \end{equation*}
There are four possibilities then,
 \begin{equation*}
 (x,y) = (0,0),(0,2),(3,0),(\frac{4}{5},\frac{11}{10})
 \end{equation*}
(c) The Jacobian for this question is then,
 \begin{equation*}
 J = \left( \begin{array}{cc} 1.5-1.5x-y & -x \\ -1.125y & 2-2y-1.125x \end{array} \right).
 \end{equation*}
For the point (0,0):
 \begin{equation*}
 J= \left( \begin{array}{cc} 1.5 & 0 \\ 0 & 2 \end{array} \right).
 \end{equation*}
So,
 \begin{equation*}
 \lambda_{1} = 1.5, \lambda_{2} = 2 \\
 \xi_{1} = \left(\begin{array}{cc} 1 \\ 0 \end{array} \right), \xi_{2} = \left(\begin{array}{cc} 0 \\ 1 \end{array} \right) 
 \end{equation*}
and thus, (0,0) represents a unstable node. For (0,2),
 \begin{equation*}
  J= \left( \begin{array}{cc} -0.5 & 0 \\ -2.25 & -2 \end{array} \right).
 \end{equation*}
So,
  \begin{equation*}
  \lambda_{1} = -0.5, \lambda_{2} = -2 \\
  \xi_{1} = \left(\begin{array}{cc} 3 \\ 1 \end{array} \right), \xi_{2} = \left(\begin{array}{cc} 0 \\ 1 \end{array} \right)
 \end{equation*}
and thus (0,2) its a node which is asymptotically stable. For (3,0),
 \begin{equation*}
  J= \left( \begin{array}{cc} -1.5 & -3 \\  0 & -1.375 \end{array} \right).
 \end{equation*}
So,
  \begin{equation*}
  \lambda_{1} = -1.5, \lambda_{2} = -1.375 \\
  \xi_{1} = \left(\begin{array}{cc} 1 \\ 0 \end{array} \right), \xi_{2} = \left(\begin{array}{cc} 24 \\ -1 \end{array} \right)
  \end{equation*}
and thus (3,0) is also an asymptotically stable node. Finally, for (4/5,11/10),
\begin{equation*}
  J= \left( \begin{array}{cc} -\frac{2}{5} & -\frac{4}{5} \\ -\frac{99}{80} & -\frac{11}{10} \end{array} \right).
 \end{equation*}
So,
  \begin{equation*}
  \lambda^{2} + \frac{3}{2}\lambda - \frac{11}{20} = 0 \\   
  \lambda_{1} = -1.8, \lambda_{2} = 0.3 \\
  \xi_{1} = \left(\begin{array}{cc} 1 \\ -1.4 \end{array} \right), \xi_{2} = \left(\begin{array}{cc} 1.57 \\ -1 \end{array} \right)
  \end{equation*}
and thus (4/5,11/10) is a saddle point (unstable).
« Last Edit: April 03, 2013, 10:26:42 PM by Frank (Yi) Gao »

Benny Ho

  • Newbie
  • *
  • Posts: 4
  • Karma: 2
    • View Profile
Re: Night Sections
« Reply #2 on: April 03, 2013, 08:56:47 PM »
solution

Devangi Vaghela

  • Jr. Member
  • **
  • Posts: 6
  • Karma: 3
    • View Profile
Re: Night Sections
« Reply #3 on: April 03, 2013, 09:19:55 PM »
Part One of Answer

Devangi Vaghela

  • Jr. Member
  • **
  • Posts: 6
  • Karma: 3
    • View Profile
Re: Night Sections
« Reply #4 on: April 03, 2013, 09:45:23 PM »
Part Two of Answer

Victor Ivrii

  • Administrator
  • Elder Member
  • *****
  • Posts: 2607
  • Karma: 0
    • View Profile
    • Personal website of Victor Ivrii
Re: Night Sections
« Reply #5 on: April 04, 2013, 12:04:28 AM »
Benny, nice try, but ...

Devangi – good, Frank—even better. But phase portrait is still missing (and bonus question too)

Sabrina (Man) Luo

  • Jr. Member
  • **
  • Posts: 9
  • Karma: 4
    • View Profile
Re: Night Sections
« Reply #6 on: April 04, 2013, 08:20:40 AM »
Benny, nice try, but ...

Devangi – good, Frank—even better. But phase portrait is still missing (and bonus question too)

Sabrina (Man) Luo

  • Jr. Member
  • **
  • Posts: 9
  • Karma: 4
    • View Profile
Re: Night Sections
« Reply #7 on: April 04, 2013, 08:46:34 AM »
Benny, nice try, but ...

Devangi – good, Frank—even better. But phase portrait is still missing (and bonus question too)
« Last Edit: April 04, 2013, 08:52:27 AM by Sabrina (Man) Luo »

Alexander Jankowski

  • Full Member
  • ***
  • Posts: 23
  • Karma: 19
    • View Profile
Re: Night Sections
« Reply #8 on: April 04, 2013, 11:17:02 AM »
To complete this thread, here is a stream plot. Overall, the analysis tells us that in this case, one of the competing species must eventually die out--unless, by chance, the initial populations lie on the separatrix that passes through $(4/5,11/10)$.

Victor Ivrii

  • Administrator
  • Elder Member
  • *****
  • Posts: 2607
  • Karma: 0
    • View Profile
    • Personal website of Victor Ivrii
Re: Night Sections
« Reply #9 on: April 04, 2013, 12:16:56 PM »
To complete this thread, here is a stream plot. Overall, the analysis tells us that in this case, one of the competing species must eventually die out--unless, by chance, the initial populations lie on the separatrix that passes through $(4/5,11/10)$.

But it exactly "one of" depending on initial conditions.