Worked Problems for Critical Point Analysis in Nonlinear Systems of ODEs

📚 Understanding Critical Point Analysis in Nonlinear Systems

Critical point analysis is a fundamental technique for understanding the qualitative behavior of nonlinear systems of ordinary differential equations (ODEs). It involves finding the equilibrium points of the system and analyzing their stability. These equilibrium points, also known as critical points or stationary points, represent states where the system remains constant over time. Let's explore this concept in detail.

📜 History and Background

The study of dynamical systems and their critical points dates back to the late 19th century with the work of Henri Poincaré. He emphasized the importance of understanding the qualitative behavior of solutions, rather than finding explicit solutions, which are often impossible to obtain for nonlinear systems. This approach led to the development of stability theory and the analysis of critical points.

📌 Key Principles

• 📍 Definition of a Critical Point: A critical point $(\bar{x}, \bar{y})$ of a system of ODEs $\frac{dx}{dt} = f(x, y)$ and $\frac{dy}{dt} = g(x, y)$ is a point where $f(\bar{x}, \bar{y}) = 0$ and $g(\bar{x}, \bar{y}) = 0$.
• 🧮 Finding Critical Points: To find the critical points, solve the system of equations $f(x, y) = 0$ and $g(x, y) = 0$ simultaneously.
• 📈 Linearization: Linearize the system around each critical point by computing the Jacobian matrix: $$J(x, y) = \begin{bmatrix} \frac{\partial f}{\partial x} & \frac{\partial f}{\partial y} \\ \frac{\partial g}{\partial x} & \frac{\partial g}{\partial y} \end{bmatrix}$$ Evaluate the Jacobian at each critical point $(\bar{x}, \bar{y})$ to obtain $J(\bar{x}, \bar{y})$.
• 🔑 Stability Analysis: Determine the eigenvalues $\lambda$ of the Jacobian matrix $J(\bar{x}, \bar{y})$. The stability of the critical point is determined by the nature of these eigenvalues:
  • If both eigenvalues have negative real parts, the critical point is a stable node or stable spiral.
  • If both eigenvalues have positive real parts, the critical point is an unstable node or unstable spiral.
  • If the eigenvalues are real and have opposite signs, the critical point is a saddle point.
  • If the eigenvalues are purely imaginary, the critical point is a center (stable, but not asymptotically stable).

⚙️ Worked Examples

Example 1:

Consider the system:

$$\frac{dx}{dt} = x - x^2 - xy$$ $$\frac{dy}{dt} = -y + xy$$

1. Find Critical Points:

Set $\frac{dx}{dt} = 0$ and $\frac{dy}{dt} = 0$:

$$x(1 - x - y) = 0$$ $$y(-1 + x) = 0$$

From the second equation, either $y = 0$ or $x = 1$.

• If $y = 0$, then from the first equation, $x(1 - x) = 0$, so $x = 0$ or $x = 1$. This gives critical points $(0, 0)$ and $(1, 0)$.
• If $x = 1$, then from the first equation, $1 - 1 - y = 0$, so $y = 0$. This confirms the critical point $(1, 0)$.
• Also, if $x=1$, then $1-1-y=0 \implies y=0$ and $-y+xy=0 \implies -y+y=0$. So, $(1,0)$ is a critical point.

2. Compute the Jacobian Matrix:

$$J(x, y) = \begin{bmatrix} 1 - 2x - y & -x \\ y & -1 + x \end{bmatrix}$$

3. Evaluate the Jacobian at each critical point:

• At $(0, 0)$: $J(0, 0) = \begin{bmatrix} 1 & 0 \\ 0 & -1 \end{bmatrix}$. Eigenvalues are $\lambda_1 = 1$ and $\lambda_2 = -1$. Saddle point (unstable).
• At $(1, 0)$: $J(1, 0) = \begin{bmatrix} -1 & -1 \\ 0 & 0 \end{bmatrix}$. Eigenvalues are $\lambda_1 = -1$ and $\lambda_2 = 0$. This requires further analysis (e.g., center manifold theory) to determine stability.

Example 2:

Consider the system:

$$\frac{dx}{dt} = x(2 - x - y)$$ $$\frac{dy}{dt} = y(y - x)$$

1. Find Critical Points:

Set $\frac{dx}{dt} = 0$ and $\frac{dy}{dt} = 0$:

$$x(2 - x - y) = 0$$ $$y(y - x) = 0$$

From the first equation, either $x = 0$ or $2 - x - y = 0$. From the second equation, either $y = 0$ or $y = x$.

• If $x = 0$, then from the second equation, $y(y - 0) = 0$, so $y = 0$. This gives the critical point $(0, 0)$.
• If $y = 0$, then from the first equation, $x(2 - x) = 0$, so $x = 0$ or $x = 2$. This gives critical points $(0, 0)$ and $(2, 0)$.
• If $y = x$, then from the first equation, $x(2 - x - x) = 0$, so $x(2 - 2x) = 0$, which gives $x = 0$ or $x = 1$. This gives critical points $(0, 0)$ and $(1, 1)$.

Thus, the critical points are $(0, 0)$, $(2, 0)$, and $(1, 1)$.

2. Compute the Jacobian Matrix:

$$J(x, y) = \begin{bmatrix} 2 - 2x - y & -x \\ -y & 2y - x \end{bmatrix}$$

3. Evaluate the Jacobian at each critical point:

• At $(0, 0)$: $J(0, 0) = \begin{bmatrix} 2 & 0 \\ 0 & 0 \end{bmatrix}$. Eigenvalues are $\lambda_1 = 2$ and $\lambda_2 = 0$. Requires further analysis.
• At $(2, 0)$: $J(2, 0) = \begin{bmatrix} -2 & -2 \\ 0 & -2 \end{bmatrix}$. Eigenvalues are $\lambda_1 = -2$ and $\lambda_2 = -2$. Stable node.
• At $(1, 1)$: $J(1, 1) = \begin{bmatrix} -1 & -1 \\ -1 & 1 \end{bmatrix}$. The characteristic equation is $(-1 - \lambda)(1 - \lambda) - 1 = 0$, which simplifies to $\lambda^2 = 2$, so $\lambda = \pm \sqrt{2}$. Saddle point (unstable).

💡 Conclusion

Critical point analysis provides valuable insights into the behavior of nonlinear systems. By finding and classifying critical points, we can understand the stability and long-term dynamics of these systems. This method is widely used in various fields, including physics, engineering, and economics, to model and analyze complex phenomena.