The study of Ordinary Differential Equations (ODEs) is central to understanding [continuous](/page/Continuity) dynamical systems. Whether modeling the motion of planets, chemical reaction rates, or population dynamics, the fundamental object is a [function](/page/Function) describing the rate of change of a state vector. While explicit solutions are rare, we can often deduce the qualitative behavior of solutions—such as stability, periodicity, and dependence on parameters—directly from the structure of the differential equation. This article establishes the formal framework for autonomous systems and investigates the fundamental questions of existence and uniqueness.

## Formal Definition

We begin by defining the structure of the autonomous ordinary differential equation. We consider the state of a system as a vector in Euclidean space and its time evolution as a curve governed by a vector field.

[definition: Autonomous ODE] Let $n \in \mathbb{N}$ and let $U \subseteq \mathbb{R}^n$ be an [open set](/page/Open%20Set). Let $f: U \to \mathbb{R}^n$ be a continuous function. An **autonomous ODE system** is defined by the equation: \begin{align*} \frac{d}{dt} X(t) = f(X(t)) \end{align*} where the solution is a [differentiable](/page/Derivative) map: \begin{align*} X: I &\to U \\ t &\mapsto X(t) \end{align*} defined on an open interval $I \subseteq \mathbb{R}$. We also study parameter-dependent families of solutions. Let $\epsilon \in \mathbb{R}$ be a parameter. We consider the map: \begin{align*} X: I \times \mathbb{R} &\to U \\ (t, \epsilon) &\mapsto X(t, \epsilon) \end{align*} which satisfies the differential equation for each fixed $\epsilon$: \begin{align*} \frac{\partial}{\partial t} X(t, \epsilon) = f(X(t, \epsilon)). \end{align*} In this context, we analyze how the trajectory $X(\cdot, \epsilon)$ varies as the parameter $\epsilon$ changes, typically manifesting through variations in the initial condition $X(0, \epsilon)$. [/definition]

---

## Examples

[example: Harmonic Oscillator] Consider the simple harmonic oscillator with state vector $X = (x, v)^T \in \mathbb{R}^2$. The function $f: \mathbb{R}^2 \to \mathbb{R}^2$ is defined by: \begin{align*} f(x, v) = \begin{pmatrix} v \\ -x \end{pmatrix}. \end{align*} The ODE system is: \begin{align*} \frac{d}{dt} \begin{pmatrix} x \\ v \end{pmatrix} = \begin{pmatrix} v \\ -x \end{pmatrix}. \end{align*} If we introduce a parameter $\epsilon$ scaling the initial amplitude, say $X(0, \epsilon) = (\epsilon, 0)^T$, the solution is $X(t, \epsilon) = (\epsilon \cos t, -\epsilon \sin t)^T$. This illustrates a linear dependence on the parameter. [/example]

---

## Existence and Uniqueness

The fundamental theoretical basis for the study of ODEs begins with the question of existence.

[quotetheorem:68]

Paeno's theorem is a good existence tool; however, it does not guarantee uniqueness.

[example] to be added [/example]

To guarantee existence, a sufficient assumption is to have a lebntz function

[quotetheorem:69]

## Stability Results

## Problems

[problem] Consider the vector field $f(x) = |x|^{1/2}$ for $x \in \mathbb{R}$. 1. Does Peano's Theorem apply to the initial value problem $\dot{x} = f(x), x(0) = 0$? 2. Find a non-zero solution to this problem. [/problem]

[solution] 1. **Applicability:** The function $f(x) = |x|^{1/2}$ is continuous on $\mathbb{R}$. Therefore, Peano's Theorem applies and guarantees the existence of a solution. 2. **Solution:** We separate variables for $x > 0$: \begin{align*} \frac{dx}{\sqrt{x}} &= dt \\ \int x^{-1/2} \, dx &= \int dt \\ 2 x^{1/2} &= t + C. \end{align*} Using $x(0)=0$, we get $C=0$. Thus $2\sqrt{x} = t$, which implies $x(t) = \frac{t^2}{4}$ for $t \ge 0$. We can extend this to $\mathbb{R}$ by setting $x(t) = \frac{t^2}{4}$ for $t \ge 0$ and $x(t) = -\frac{t^2}{4}$ for $t < 0$ (signed solution), or simply $x(t) = \frac{1}{4} t |t|$. The function $x(t) = \frac{t^2}{4}$ is a valid non-zero solution for $t>0$. [/solution]

---