Stability in Systems - Differential Equations

  1. Consider the linear system x′=(α−15−1)x\mathbf{x}' = \begin{pmatrix} \alpha & -1 \\ 5 & -1 \end{pmatrix} \mathbf{x}. For which range of the parameter α\alpha is the origin an asymptotically stable spiral point?
    Answer: −1−25<α<1-1 - 2\sqrt{5} < \alpha < 1 The stability is determined by the eigenvalues of the matrix A=(α−15−1)A = \begin{pmatrix} \alpha & -1 \\ 5 & -1 \end{pmatrix}. The characteristic equation is λ2−tr(A)λ+det⁡(A)=0\lambda^2 - \mathrm{tr}(A)\lambda + \det(A) = 0, which is λ2−(α−1)λ+(5−α)=0\lambda^2 - (\alpha-1)\lambda + (5-\alpha) = 0. For a spiral point, the eigenvalues must be complex, which means the discriminant of the characteristic equation must be negative: (α−1)2−4(5−α)<0(\alpha-1)^2 - 4(5-\alpha) < 0. This simplifies to α2−2α+1−20+4α<0\alpha^2 - 2\alpha + 1 - 20 + 4\alpha < 0, or α2+2α−19<0\alpha^2 + 2\alpha - 19 < 0. The roots of α2+2α−19=0\alpha^2 + 2\alpha - 19 = 0 are \alpha = -1 \pm \sqrt{1 - (-19)}} = -1 \pm \sqrt{20} = -1 \pm 2\sqrt{5}. Thus, we need −1−25<α<−1+25-1 - 2\sqrt{5} < \alpha < -1 + 2\sqrt{5}. For the spiral to be asymptotically stable, the real part of the eigenvalues must be negative. The real part is tr(A)2=α−12\frac{\mathrm{tr}(A)}{2} = \frac{\alpha-1}{2}. We need α−12<0\frac{\alpha-1}{2} < 0, which implies α<1\alpha < 1. Combining both conditions, we need −1−25<α<1-1 - 2\sqrt{5} < \alpha < 1.
  2. A biochemical reaction network has the linearized system ddt(xyz)=(−2100−1210−3)(xyz)\frac{d}{dt}\begin{pmatrix} x \\ y \\ z \end{pmatrix} = \begin{pmatrix} -2 & 1 & 0 \\ 0 & -1 & 2 \\ 1 & 0 & -3 \end{pmatrix}\begin{pmatrix} x \\ y \\ z \end{pmatrix} near an equilibrium. What can be concluded about the stability of this equilibrium?
    Answer: The equilibrium is stable because all eigenvalues have negative real parts The matrix A = [[-2,1,0],[0,-1,2],[1,0,-3]] has characteristic polynomial det(A-λI). We can use the Routh-Hurwitz criteria: the trace is -6 < 0, and we need to check that all principal minors alternate in sign appropriately. Computing: the characteristic polynomial is -λ³ - 6λ² - 9λ - 4. By Routh-Hurwitz, since all coefficients are negative (when written as λ³ + 6λ² + 9λ + 4), we need: 6 > 0 ✓, 6·9 - 1·4 = 50 > 0 ✓, and 4 > 0 ✓. All conditions are satisfied, so all eigenvalues have negative real parts. Choice A is incorrect because negative trace is necessary but not sufficient for stability in 3D. Choice B is wrong because all eigenvalues actually have negative real parts. Choice D is wrong because the Routh-Hurwitz criteria allow us to determine stability without explicit eigenvalue computation.
  3. Consider the linear system dxdt=3x+2y\frac{dx}{dt} = 3x + 2y, dydt=−x+4y\frac{dy}{dt} = -x + 4y. After determining the eigenvalues of the coefficient matrix, what can be concluded about the stability of the equilibrium point at the origin?
    Answer: The equilibrium is unstable because both eigenvalues have positive real parts The coefficient matrix is A = [[3, 2], [-1, 4]]. The characteristic equation is det(A - λI) = (3-λ)(4-λ) - (2)(-1) = λ² - 7λ + 14 = 0. Using the quadratic formula: λ = (7 ± √(49-56))/2 = (7 ± √(-7))/2 = 7/2 ± i√7/2. Both eigenvalues have positive real part 7/2 > 0, making the equilibrium unstable. Choice A is wrong because the real parts are positive, not negative. Choice C is wrong because both eigenvalues have the same sign (positive real part). Choice D is wrong because the eigenvalues have nonzero real parts.
  4. If the eigenvalues of the real matrix AA for the linear system x′=Ax\mathbf{x}' = A\mathbf{x} are λ1=3i\lambda_1 = 3i and λ2=−3i\lambda_2 = -3i, which of the following best describes the long-term behavior of non-zero solutions?
    Answer: Trajectories are ellipses or circles, representing periodic solutions. The eigenvalues are a purely imaginary conjugate pair, λ=±bi\lambda = \pm bi with b=3b=3. For a linear system, this corresponds to a stable center at the origin. The solutions are linear combinations of cos⁡(3t)\cos(3t) and sin⁡(3t)\sin(3t), which are periodic. The trajectories are closed orbits around the origin, specifically ellipses (or circles in special cases). These solutions are stable because they remain in a bounded neighborhood of the origin, but they are not asymptotically stable because they do not approach the origin as t→∞t \to \infty.
  5. For the Van der Pol oscillator d2xdt2−μ(1−x2)dxdt+x=0\frac{d^2x}{dt^2} - \mu(1-x^2)\frac{dx}{dt} + x = 0 with μ>0\mu > 0, converted to the system dxdt=y\frac{dx}{dt} = y, dydt=μ(1−x2)y−x\frac{dy}{dt} = \mu(1-x^2)y - x, what happens to the stability of the origin as μ\mu increases from 0?
    Answer: The origin becomes unstable for any μ>0\mu > 0 and a stable limit cycle emerges When analyzing the Van der Pol oscillator, you're dealing with a classic example of how nonlinear systems can exhibit bifurcations and limit cycle behavior. The key is to examine how the linear stability of the equilibrium point changes with the parameter μ\mu. To determine stability of the origin, you linearize the system around (0,0)(0,0). The Jacobian matrix is: J=(01−1μ)J = \begin{pmatrix} 0 & 1 \\ -1 & \mu \end{pmatrix} The eigenvalues satisfy λ2−μλ+1=0\lambda^2 - \mu\lambda + 1 = 0, giving λ=μ±μ2−42\lambda = \frac{\mu \pm \sqrt{\mu^2 - 4}}{2}. For any μ>0\mu > 0, the real parts are positive (since the sum of eigenvalues is μ>0\mu > 0 and the product is 1), making the origin unstable. The Van der Pol oscillator is famous because when μ>0\mu > 0, trajectories spiral outward from the unstable origin but are eventually pulled back by the nonlinear damping term μ(1−x2)y\mu(1-x^2)y, creating a stable limit cycle. Answer A is wrong because the origin is unstable, not stable, for μ>0\mu > 0. Answer B incorrectly identifies a transcritical bifurcation at μ=1\mu = 1 - the stability actually changes at μ=0\mu = 0, and this is a Hopf bifurcation. Answer C is incorrect because the origin is unstable for any μ>0\mu > 0, not just μ>1\mu > 1. Remember: For the Van der Pol oscillator, any positive μ\mu creates an unstable spiral at the origin with an attracting limit cycle - this is the hallmark of a supercritical Hopf bifurcation at μ=0\mu = 0.
  6. A system is described in polar coordinates by the equations r˙=r(4−r2)\dot{r} = r(4 - r^2) and θ˙=−1\dot{\theta} = -1. What is the long-term behavior of a trajectory that starts at the Cartesian point (x0,y0)=(3,0)(x_0, y_0) = (3, 0)?
    Answer: The trajectory approaches a limit cycle of radius 2. The initial point (x0,y0)=(3,0)(x_0, y_0) = (3, 0) corresponds to a radius of r0=32+02=3r_0 = \sqrt{3^2 + 0^2} = 3. The equation for the radius is r˙=r(4−r2)=r(2−r)(2+r)\dot{r} = r(4 - r^2) = r(2-r)(2+r). For r>0r>0, the radial velocity r˙\dot{r} is zero when r=2r=2. If 0<r<20 < r < 2, r˙>0\dot{r} > 0, so rr increases. If r>2r > 2, r˙<0\dot{r} < 0, so rr decreases. This means that r=2r=2 is a stable limit cycle. Since the trajectory starts at r0=3r_0=3, which is greater than 2, the radius r(t)r(t) will decrease and approach 2 as t→∞t \to \infty. The equation θ˙=−1\dot{\theta}=-1 indicates that the trajectory spirals in a clockwise direction. Therefore, the trajectory starting at r=3r=3 approaches the stable limit cycle at r=2r=2.
  7. For the system x′=(−311−3)x\mathbf{x}' = \begin{pmatrix} -3 & 1 \\ 1 & -3 \end{pmatrix} \mathbf{x}, trajectories approach the origin as t→∞t \to \infty. As they approach the origin, they become tangent to which line?
    Answer: y=xy = x The asymptotic behavior of trajectories near a stable node is determined by the eigenvectors. First, we find the eigenvalues. The characteristic equation is (−3−λ)2−1=0(-3-\lambda)^2 - 1 = 0, which is λ2+6λ+9−1=0\lambda^2 + 6\lambda + 9 - 1 = 0, or λ2+6λ+8=0\lambda^2 + 6\lambda + 8 = 0. Factoring gives (λ+4)(λ+2)=0(\lambda+4)(\lambda+2) = 0, so the eigenvalues are λ1=−2\lambda_1 = -2 and λ2=−4\lambda_2 = -4. The general solution is x(t)=c1e−2tv1+c2e−4tv2\mathbf{x}(t) = c_1 e^{-2t} \mathbf{v}_1 + c_2 e^{-4t} \mathbf{v}_2. As t→∞t \to \infty, the term e−4te^{-4t} decays to zero much faster than e−2te^{-2t}. Thus, for large tt, the solution is dominated by the term with the less negative eigenvalue, λ1=−2\lambda_1 = -2. The trajectory will approach the origin tangent to the eigenvector v1\mathbf{v}_1. We find v1\mathbf{v}_1 by solving (A−(−2)I)v1=0(A - (-2)I)\mathbf{v}_1 = \mathbf{0}: (−3−(−2)11−3−(−2))(xy)=(−111−1)(xy)=(00)\begin{pmatrix} -3 - (-2) & 1 \\ 1 & -3 - (-2) \end{pmatrix} \begin{pmatrix} x \\ y \end{pmatrix} = \begin{pmatrix} -1 & 1 \\ 1 & -1 \end{pmatrix} \begin{pmatrix} x \\ y \end{pmatrix} = \begin{pmatrix} 0 \\ 0 \end{pmatrix} This gives the equation −x+y=0-x+y=0, or y=xy=x. So, trajectories approach the origin tangent to the line y=xy=x.
  8. Consider the 3D linear system x′=Ax\mathbf{x}' = A\mathbf{x} where A=(−2001−110−1−1)A = \begin{pmatrix} -2 & 0 & 0 \\ 1 & -1 & 1 \\ 0 & -1 & -1 \end{pmatrix}. Which statement best describes the long-term behavior of solutions to this system?
    Answer: All non-zero solutions approach the origin as t→∞t \to \infty. The long-term behavior of a linear system is determined by the eigenvalues of the matrix AA. Since AA is block lower triangular, its eigenvalues are the eigenvalues of its diagonal blocks. The first block is the 1×11 \times 1 matrix (−2)(-2), so one eigenvalue is λ1=−2\lambda_1 = -2. The second block is the 2×22 \times 2 matrix M=(−11−1−1)M = \begin{pmatrix} -1 & 1 \\ -1 & -1 \end{pmatrix}. The characteristic equation for MM is (−1−λ)2−(1)(−1)=0(-1-\lambda)^2 - (1)(-1) = 0, which is (λ+1)2+1=0(\lambda+1)^2 + 1 = 0. This gives λ+1=±i\lambda+1 = \pm i, so the other two eigenvalues are λ2,3=−1±i\lambda_{2,3} = -1 \pm i. The eigenvalues of AA are therefore {−2,−1+i,−1−i}\{-2, -1+i, -1-i\}. For a system to be asymptotically stable, all eigenvalues must have negative real parts. The real parts are −2-2, −1-1, and −1-1. Since all are negative, the origin is an asymptotically stable equilibrium point. This means that all solutions, regardless of their initial condition, will approach the origin as t→∞t \to \infty.
  9. Which of the following statements about the stability of equilibrium points for 2D autonomous systems is true?
    Answer: For a linear system x′=Ax\mathbf{x}' = A\mathbf{x}, if det⁡(A)<0\det(A) < 0, then some solutions will move away from the origin as t→∞t \to \infty. Let's analyze each statement: A. False. If the linearization has purely imaginary eigenvalues, it is a borderline case. The stability of the nonlinear system's equilibrium point depends on the nonlinear terms and cannot be determined from the linearization alone; it could be a stable center, an unstable spiral, or a stable spiral. B. False. If all trajectories approach x0\mathbf{x}_0, it is an asymptotically stable equilibrium. However, it could be either a stable node (real negative eigenvalues) or a stable spiral (complex eigenvalues with negative real part). C. True. For a 2D linear system, det⁡(A)=λ1λ2\det(A) = \lambda_1 \lambda_2. If det⁡(A)<0\det(A) < 0, the eigenvalues must be real and have opposite signs (e.g., one positive, one negative). This corresponds to a saddle point, which is an unstable equilibrium. Trajectories along the eigenvector corresponding to the positive eigenvalue will move away from the origin. D. False. A stable equilibrium point is not necessarily asymptotically stable. A center, for example, is stable (trajectories stay nearby) but not asymptotically stable (they do not approach the equilibrium point).
  10. A predator-prey system has the form dxdt=x(a−by)\frac{dx}{dt} = x(a-by), dydt=y(−c+dx)\frac{dy}{dt} = y(-c+dx) where a,b,c,d>0a, b, c, d > 0. The equilibrium point (cd,ab)\left(\frac{c}{d}, \frac{a}{b}\right) is analyzed for stability. Which statement correctly describes the nature of this equilibrium?
    Answer: It is neutrally stable with periodic orbits because the eigenvalues are purely imaginary The Jacobian at the equilibrium is J = [[0, -bc/d], [ad/b, 0]]. The characteristic equation is λ² + ac = 0, giving eigenvalues λ = ±i√(ac). Since a,c > 0, these are purely imaginary, indicating a center with periodic orbits (neutrally stable). Choice A is wrong because the trace is 0, not negative. Choice B is wrong because det(J) = ac > 0, not negative. Choice D is wrong because the eigenvalues don't have negative real parts; they have zero real parts.
  11. A competitive species model has the form dxdt=x(1−x−ay)\frac{dx}{dt} = x(1-x-ay), dydt=y(1−bx−y)\frac{dy}{dt} = y(1-bx-y) where a,b>0a, b > 0. For the coexistence equilibrium (1−a1−ab,1−b1−ab)\left(\frac{1-a}{1-ab}, \frac{1-b}{1-ab}\right) to be stable when ab>1ab > 1, which condition must hold?
    Answer: The coexistence equilibrium cannot be stable when ab>1ab > 1 due to competitive exclusion When analyzing competitive species models, you're examining how two populations interact when competing for the same resources. The key insight is understanding when species can coexist versus when one will drive the other to extinction. To determine stability, you need to examine the coexistence equilibrium's existence and stability conditions. For the equilibrium (1−a1−ab,1−b1−ab)\left(\frac{1-a}{1-ab}, \frac{1-b}{1-ab}\right) to exist in the positive quadrant (both species present), you need both coordinates positive. This requires 1−a>01-a > 0 and 1−b>01-b > 0 (so a<1a < 1 and b<1b < 1), plus 1−ab>01-ab > 0 (so ab<1ab < 1). However, the question states ab>1ab > 1. When ab>1ab > 1, the denominator 1−ab1-ab becomes negative, making the equilibrium coordinates negative—meaning it doesn't exist in the biologically meaningful positive quadrant. Without a valid coexistence equilibrium, stability is impossible. Choice A correctly identifies that coexistence cannot be stable when ab>1ab > 1 due to competitive exclusion—one species will always eliminate the other. Choice B is wrong because ab>1ab > 1 actually prevents stability rather than guaranteeing it. Choice C incorrectly suggests the equilibrium could exist when ab>1ab > 1 if additional conditions were met. Choice D mischaracterizes the dynamics, as there's no stable interior equilibrium when ab>1ab > 1. Study tip: In competitive models, always check equilibrium existence before analyzing stability. When competition coefficients satisfy ab>1ab > 1, expect competitive exclusion where only one species survives.
  12. For the nonlinear system dxdt=x(2−x−y)\frac{dx}{dt} = x(2-x-y), dydt=y(1−x−2y)\frac{dy}{dt} = y(1-x-2y), the equilibrium point (1,12)(1, \frac{1}{2}) has been identified. To determine its stability, what is the trace of the Jacobian matrix evaluated at this equilibrium?
    Answer: −32-\frac{3}{2} The Jacobian matrix J has entries: ∂f₁/∂x = 2-2x-y, ∂f₁/∂y = -x, ∂f₂/∂x = -y, ∂f₂/∂y = 1-x-4y. At (1, 1/2): ∂f₁/∂x = 2-2(1)-1/2 = -1/2, ∂f₁/∂y = -1, ∂f₂/∂x = -1/2, ∂f₂/∂y = 1-1-4(1/2) = -2. So J = [[-1/2, -1], [-1/2, -2]]. The trace is -1/2 + (-2) = -3/2. Choice B gives only the first diagonal entry. Choice C gives the negative of the first diagonal entry. Choice D gives the negative of the correct answer.
  13. The motion of a damped pendulum can be modeled by the equation θ¨+cθ˙+sin⁡(θ)=0\ddot{\theta} + c\dot{\theta} + \sin(\theta) = 0 with a damping constant c>0c > 0. The equilibrium point (θ,θ˙)=(π,0)(\theta, \dot{\theta}) = (\pi, 0) corresponds to the pendulum balanced in the vertical, inverted position. What is the nature of this equilibrium point?
    Answer: Saddle point (unstable) First, we convert the second-order equation into a first-order system. Let x=θx = \theta and y=θ˙y = \dot{\theta}. The system is x˙=y\dot{x} = y and y˙=−sin⁡(x)−cy\dot{y} = -\sin(x) - cy. We analyze the stability at the equilibrium point (π,0)(\pi, 0) by linearizing the system. The Jacobian matrix is J(x,y)=(01−cos⁡(x)−c)J(x,y) = \begin{pmatrix} 0 & 1 \\ -\cos(x) & -c \end{pmatrix}. Evaluating at (π,0)(\pi, 0), we get J(π,0)=(01−cos⁡(π)−c)=(011−c)J(\pi, 0) = \begin{pmatrix} 0 & 1 \\ -\cos(\pi) & -c \end{pmatrix} = \begin{pmatrix} 0 & 1 \\ 1 & -c \end{pmatrix}. The determinant of this matrix is det⁡(J)=(0)(−c)−(1)(1)=−1\det(J) = (0)(-c) - (1)(1) = -1. Since the determinant is negative, the equilibrium point is a saddle point, which is always unstable. This matches the physical intuition that it is impossible to balance a damped pendulum in the inverted position; any small perturbation will cause it to fall.
  14. Consider the system dxdt=y\frac{dx}{dt} = y, dydt=−sin⁡(x)−δy\frac{dy}{dt} = -\sin(x) - \delta y where δ>0\delta > 0. This represents a damped pendulum. For small oscillations about the stable equilibrium, what determines the transition between overdamped and underdamped motion?
    Answer: The transition occurs when δ=2\delta = 2, where the discriminant of the characteristic equation changes sign When analyzing damped oscillatory systems like the pendulum, you need to linearize around equilibrium and examine the characteristic equation to determine the nature of the motion. For small oscillations about x=0x = 0, you can approximate sin⁡(x)≈x\sin(x) \approx x, giving the linearized system: dxdt=y\frac{dx}{dt} = y and dydt=−x−δy\frac{dy}{dt} = -x - \delta y. Converting to a second-order equation: d2xdt2+δdxdt+x=0\frac{d^2x}{dt^2} + \delta \frac{dx}{dt} + x = 0. The characteristic equation is r2+δr+1=0r^2 + \delta r + 1 = 0 with discriminant Δ=δ2−4\Delta = \delta^2 - 4. The nature of the roots determines the motion type: - If Δ>0\Delta > 0 (δ>2\delta > 2): two real roots → overdamped motion (exponential decay) - If Δ=0\Delta = 0 (δ=2\delta = 2): repeated real root → critically damped - If Δ<0\Delta < 0 (δ<2\delta < 2): complex roots → underdamped motion (oscillatory decay) The transition occurs precisely when δ=2\delta = 2, where the discriminant equals zero. Answer A incorrectly suggests δ=1\delta = 1 as the transition point, confusing this with other damping criteria. Answer B is wrong because damped systems can still be underdamped—dissipation doesn't eliminate oscillations entirely, just causes them to decay. Answer C incorrectly involves initial amplitude (irrelevant for the overdamped/underdamped classification) and gives δ=2\delta = \sqrt{2}. Study tip: For any damped harmonic oscillator x¨+2γx˙+ω02x=0\ddot{x} + 2\gamma\dot{x} + \omega_0^2 x = 0, the transition occurs when the damping coefficient equals the natural frequency. Always check the discriminant of the characteristic equation to classify the motion.
  15. Consider the nonlinear system x˙=−x3−y\dot{x} = -x^3 - y and y˙=x−y3\dot{y} = x - y^3. Using the Liapunov function V(x,y)=12(x2+y2)V(x,y) = \frac{1}{2}(x^2+y^2), what can be concluded about the stability of the origin (0,0)(0,0)?
    Answer: The origin is asymptotically stable. To use the Liapunov function V(x,y)=12(x2+y2)V(x,y) = \frac{1}{2}(x^2+y^2), we first note that it is positive definite (V(0,0)=0V(0,0)=0 and V(x,y)>0V(x,y) > 0 otherwise). Next, we compute its time derivative along the trajectories of the system, V˙=∂V∂xx˙+∂V∂yy˙\dot{V} = \frac{\partial V}{\partial x}\dot{x} + \frac{\partial V}{\partial y}\dot{y}. We have ∂V∂x=x\frac{\partial V}{\partial x} = x and ∂V∂y=y\frac{\partial V}{\partial y} = y. Substituting the expressions for x˙\dot{x} and y˙\dot{y}: V˙=x(−x3−y)+y(x−y3)=−x4−xy+xy−y4=−(x4+y4)\dot{V} = x(-x^3 - y) + y(x - y^3) = -x^4 - xy + xy - y^4 = -(x^4 + y^4). Since V˙(0,0)=0\dot{V}(0,0)=0 and V˙(x,y)<0\dot{V}(x,y) < 0 for all (x,y)≠(0,0)(x,y) \neq (0,0), V˙\dot{V} is negative definite. According to Liapunov's stability theorem, if there exists a positive definite function VV whose derivative V˙\dot{V} is negative definite, then the equilibrium point at the origin is asymptotically stable.
  16. The linearization of the system x˙=y+x(x2+y2)\dot{x} = y + x(x^2+y^2), y˙=−x+y(x2+y2)\dot{y} = -x + y(x^2+y^2) at the origin (0,0)(0,0) is x˙=y\dot{x} = y, y˙=−x\dot{y} = -x, which corresponds to a center. What is the actual stability of the origin for the full nonlinear system?
    Answer: The origin is unstable. When the linearization of a nonlinear system at an equilibrium point results in a center (eigenvalues are purely imaginary), the linearization is inconclusive and the stability is determined by the nonlinear terms. To analyze the full system, we can convert to polar coordinates. Let x=rcos⁡θx = r\cos\theta and y=rsin⁡θy = r\sin\theta. The time derivative of r2=x2+y2r^2 = x^2+y^2 is 2rr˙=2xx˙+2yy˙2r\dot{r} = 2x\dot{x} + 2y\dot{y}. Substituting for x˙\dot{x} and y˙\dot{y}: rr˙=x(y+x(x2+y2))+y(−x+y(x2+y2))=xy+x2(r2)−xy+y2(r2)=(x2+y2)r2=r4r\dot{r} = x(y + x(x^2+y^2)) + y(-x + y(x^2+y^2)) = xy + x^2(r^2) - xy + y^2(r^2) = (x^2+y^2)r^2 = r^4. For r≠0r \neq 0, we have r˙=r3\dot{r} = r^3. Since r˙>0\dot{r} > 0 for all r>0r>0, the radius r(t)r(t) is strictly increasing. Trajectories move away from the origin, so the origin is an unstable equilibrium point.
  17. This system has a non-zero equilibrium point representing the coexistence of both species. What is the nature of the stability at this coexistence equilibrium point?
    Answer: Asymptotically stable node First, find the coexistence equilibrium point where x>0x > 0 and y>0y > 0. This requires solving x˙=0\dot{x}=0 and y˙=0\dot{y}=0: 1−x−y=01 - x - y = 0 and 0.75−y−0.5x=00.75 - y - 0.5x = 0. From the first equation, y=1−xy = 1 - x. Substituting into the second gives 0.75−(1−x)−0.5x=00.75 - (1 - x) - 0.5x = 0, which simplifies to −0.25+0.5x=0-0.25 + 0.5x = 0, so x=0.5x=0.5. Then y=1−0.5=0.5y = 1 - 0.5 = 0.5. The equilibrium point is (0.5,0.5)(0.5, 0.5). Next, linearize the system by finding the Jacobian matrix: J(x,y)=(1−2x−y−x−0.5y0.75−2y−0.5x)J(x,y) = \begin{pmatrix} 1-2x-y & -x \\ -0.5y & 0.75-2y-0.5x \end{pmatrix}. Evaluate the Jacobian at (0.5,0.5)(0.5, 0.5): J(0.5,0.5)=(1−1−0.5−0.5−0.250.75−1−0.25)=(−0.5−0.5−0.25−0.5)J(0.5, 0.5) = \begin{pmatrix} 1-1-0.5 & -0.5 \\ -0.25 & 0.75-1-0.25 \end{pmatrix} = \begin{pmatrix} -0.5 & -0.5 \\ -0.25 & -0.5 \end{pmatrix}. The trace is T=−0.5−0.5=−1T = -0.5 - 0.5 = -1. The determinant is D=(−0.5)(−0.5)−(−0.5)(−0.25)=0.25−0.125=0.125D = (-0.5)(-0.5) - (-0.5)(-0.25) = 0.25 - 0.125 = 0.125. Since T<0T < 0 and D>0D > 0, the equilibrium is asymptotically stable. To determine if it is a node or spiral, we check the discriminant T2−4D=(−1)2−4(0.125)=1−0.5=0.5>0T^2 - 4D = (-1)^2 - 4(0.125) = 1 - 0.5 = 0.5 > 0. Since the discriminant is positive, the eigenvalues are real and distinct, meaning the point is an asymptotically stable node.
  18. For the system dxdt=μx−y−x3\frac{dx}{dt} = \mu x - y - x^3, dydt=x+μy−y3\frac{dy}{dt} = x + \mu y - y^3 where μ\mu is a parameter, at what value of μ\mu does the origin undergo a bifurcation, and what type occurs?
    Answer: μ=0\mu = 0; Hopf bifurcation as complex eigenvalues cross the imaginary axis The Jacobian at the origin is J = [[μ, -1], [1, μ]], with eigenvalues λ = μ ± i. At μ = 0, the eigenvalues are ±i (purely imaginary), and as μ increases through 0, they cross the imaginary axis with nonzero speed (Re(λ) = μ). This is the condition for a Hopf bifurcation. The cubic terms provide the nonlinearity needed to determine supercritical vs subcritical nature. Choice A is wrong because transcritical bifurcations involve real eigenvalues passing through zero. Choice C is wrong about both the parameter value and bifurcation type. Choice D is wrong because saddle-node bifurcations involve repeated real eigenvalues.
  19. For the system x′=(1kk−3)x\mathbf{x}' = \begin{pmatrix} 1 & k \\ k & -3 \end{pmatrix} \mathbf{x}, the origin is a saddle point for which values of the real parameter kk?
    Answer: For all real kk The classification of the equilibrium point at the origin for a linear system x′=Ax\mathbf{x}' = A\mathbf{x} can be determined by the trace, T=tr(A)T = \mathrm{tr}(A), and the determinant, D=det⁡(A)D = \det(A). For the given matrix, T=1+(−3)=−2T = 1 + (-3) = -2 and D=(1)(−3)−(k)(k)=−3−k2D = (1)(-3) - (k)(k) = -3 - k^2. An equilibrium point is a saddle point if and only if D<0D < 0. In this case, we need −3−k2<0-3 - k^2 < 0. This inequality can be written as k2>−3k^2 > -3. Since kk is a real parameter, k2k^2 is always non-negative (k2≥0k^2 \ge 0). Therefore, the condition k2>−3k^2 > -3 is satisfied for all real values of kk. The trace is irrelevant for identifying a saddle point, as the condition D<0D<0 is sufficient.

Start this quiz