Codex @codex  0

A Bernoulli differential equation is a first-order ordinary differential equation of the displayed form, with $P$ and $Q$ continuous functions on an interval and a fixed real exponent $n$. A nonlinear power can hide a first-order linear differential equation: on a suitable branch, the Bernoulli power substitution $z=y^{1-n}$ removes it exactly. The transformed solution of a differential equation can be completely regular even when recovering $y$ produces finite-time blowup.

For example, the reciprocal Bernoulli equation with finite-time blowup $y^{\prime }+y=e^x{y}^2$ becomes $z^{\prime }-z=-e^x$ under $z=1/y$. The initial conditions $y(0)=1/2$ and $y(0)=-1$ select a positive trajectory blowing up at $x=2$ and a negative trajectory tending to zero as $x\to \infty$, respectively. The integrating factor exposes both behaviors with the same short calculation.

The cases $n=0$ and $n=1$ are already first-order linear differential equations: respectively $y^{\prime }+Py=Q$, with the zeroth-power term interpreted as the constant $1$, and $y^{\prime }+(P-Q)y=0$. There is no need for a nonlinear change of variable; for $n=1$, $y^{1-n}=1$ would discard the unknown entirely. For other exponents, the substitution must respect the domain of exponentiation and may lose the zero solution of a differential equation.^[ref]

For arbitrary real $n\ne 0,1$, work first on an interval where $y>0$. Define exponentiation by $y^r=\exp (r\log y)$, using the real logarithm. The change of variable is differentiable and invertible there, and the chain rule givesMultiplying the Bernoulli differential equation by $(1-n)y^{-n}$ yields the first-order linear differential equationThis is an exact transformation of nonzero solutions of a differential equation on the chosen branch; it is valid far from any reference state and does not discard higher powers as a local linearization would.

Given the initial condition $y(x_0)=y_0>0$, the integrating factor and transformed solution of a differential equation areThe product rule verifies $(\mu z{)}^{\prime }=(1-n)\mu Q$. Recover $y=z^{1/(1-n)}$ only on the connected component of $\{x:z(x)>0\}$ containing $x_0$, inside the original coefficient interval. The first-order linear differential equation can have a longer existence interval than this inverse transformation allows.^[ref]

Negative values need a separately justified real branch. For integer $n$, the required powers are real and differentiable for $y<0$ too, but the inverse must retain the sign selected by the initial condition. For $n=3$, $z=y^{-2}>0$ represents either $y=+1/\sqrt{z}$ or $y=-1/\sqrt{z}$; both square roots are taken of positive $z$. For $n=2$, $z=1/y$ retains the sign directly. For a general noninteger exponent, use the positive branch unless a specific differentiable real negative branch has been established. Never evaluate a real fractional power of a negative value merely because the linear formula has become negative.

Division by $y^n$ excludes $y=0$. When $n>0$ and the original real power is defined at zero, $y\equiv 0$ is an additional solution of a differential equation and must be checked in the original equation. When $n<0$, the equation as written is undefined at zero. On a positive branch the right-hand side is locally Lipschitz in $y$, so the Picard-Lindelöf theorem gives local uniqueness. At zero, the same hypothesis needs a fresh check: it holds for the polynomial $n=2$ example, whereas for $0<n<1$ a nonzero power coefficient prevents local Lipschitz continuity there. If $z\downarrow 0$, the inverse sends $y\to \infty$ for $n>1$, but sends $y\downarrow 0$ for $0<n<1$. The latter endpoint can lead to finite-time extinction in a sublinear Bernoulli equation, where the original equation determines how to join to zero.

A delay differential equation relates a derivative at the current argument to values at earlier arguments. For example, the Buchstab function obeys $\frac{d}{du}(uw(u))=w(u-1)$ for $u>2$.

The defining feature is memory: a derivative can depend on earlier values of the unknown function. Constant delays, time-dependent delays, and delays depending on the unknown function belong to this wider class. The retarded delay differential equation below isolates a single positive constant delay.^[ref]

The Malthusian delay differential equation is a simple linear member of this class. The method of steps for a delay differential equation constructs its future from a history function of a delay differential equation, and the characteristic equation of a delay differential equation studies its exponential function modes. The exact polynomial steps for a delay equation below make that construction explicit without approximating the solution of a differential equation. Applications include population dynamics and feedback control; the delayed negative feedback model and feedback delay explain why the timing of a return matters as well as its sign.

This initial value problem uses a square root and absolute value to define a continuous function on all of ${\mathbb{R}}^2$. For every finite waiting time $c\ge 0$, a forward solution of a differential equation on $t\ge 0$ isThe constant solution of a differential equation $y_{\infty }(t)=0$ is another choice. For example, at $t=3$, waits $c=0,1,2,3$ give values $9,4,1,0$. These are dimensionless mathematical trajectories; no measured physical waiting time is being supplied.

Verification uses the original ordinary differential equation, including the join. Before departure, both $y_c^{\prime }$ and $2\sqrt{|y_c|}$ vanish. For $t>c$, their common value is $2(t-c)$. For $c>0$, the difference quotient at the join isThus the derivative exists at $c$, equals the rate law there, and is continuous across the join. For $c=0$ the same computation gives the right derivative at the starting endpoint. For $c>0$, the left and right second derivatives at departure are $0$ and $2$: these forward solutions of a differential equation are continuously differentiable functions, without being twice differentiable at departure. The local two-sided statement of the Peano existence theorem is separate from this explicitly forward family.

The rate law tends to zero as $y\to 0$, but its secant slopes from zero are unbounded:No finite local Lipschitz constant exists at zero. At $\epsilon =10^{-2},10^{-4},10^{-6}$ the quotient is $20,200,2000$, even though the rates are shrinking to zero. Continuity controls the rate itself; Lipschitz continuity would control its change relative to the state displacement.

There is also a numerical lesson. The forward Euler method starting exactly at zero gives $Y_{k+1}=Y_k+2\Delta t\sqrt{|Y_k|}=0$ at every step. It reproduces the constant solution of a differential equation for every positive step size, while the moving trajectories remain equally valid. A computation producing one trajectory cannot establish uniqueness. Likewise, the separable differential equation calculation that divides by $\sqrt{|y|}$ only describes the moving pieces: checking the original equation recovers the zero pieces and their joins.

A differential equation is an equation relating an unknown function to its derivatives. It specifies how change depends on the present state, and sometimes on the independent variables themselves. A solution of a differential equation is a function satisfying that relation throughout a stated domain; giving a formula without its domain can conceal a point where the equation is undefined.

The central idea is to turn a local law into a whole trajectory. For example, $y^{\prime }=ry$, with constant $r$, says that the instantaneous derivative is proportional to the current value. If $y(0)=100$ and $r=0.2$ per day, the solution of a differential equation is $y(t)=100e^{0.2t}$, so $y(5)=100e\approx 271.83$. This is an ideal growth model assuming a constant proportional rate; the exponential function emerges from repeatedly applying the same local rule. Changing the sign of $r$ gives exponential decay.

An ordinary differential equation uses one independent variable. A partial differential equation can involve several, such as space and time, and uses partial derivatives. Additional information selects a particular solution of a differential equation: an initial condition describes a starting state, whereas a boundary condition imposes constraints at the edge of a domain. The useful questions extend beyond finding a formula: does a solution of a differential equation exist, is it unique, and what happens to it over time?

A direction field for $y^{\prime }=f(t,y)$ places a short segment with slope $f(t,y)$ at each sampled point where $f$ is defined. A function graph of a solution of a differential equation follows these segments: its tangent line has exactly the prescribed slope. Segment length is a drawing choice, rather than a speed; arrows, if used, point toward increasing $t$. An isocline collects points with one common slope.

For Newton's law of cooling, $y^{\prime }=-0.1(y-20)$ with time in minutes and temperature in degrees Celsius, the slopes at $y=80,40,20,10$ are respectively $-6,-2,0,1$ degrees Celsius per minute. Each horizontal row is an isocline, because this autonomous differential equation has no explicit time dependence. Above $20$ the solutions of a differential equation descend, below $20$ they rise, and $y=20$ is an equilibrium point.

The highlighted initial condition $y(0)=80$ selects one curve from the family. The picture reveals approach to the equilibrium point before any algebra is done. Under the hypotheses of the Picard-Lindelöf theorem, distinct function graphs of solutions cannot cross or even meet in the time-state plane. A direction field alone does not guarantee uniqueness: several curves can share a tangent line at the same point when those hypotheses fail. See the delayed-start example under initial value problem.

This geometric interpretation is developed in Direction Fields, Isoclines, Integral Curves.

For two length-$L$ slabs with zero exterior temperatures and perfect thermal contact, put ${\alpha }_i=\sqrt{{\rho }_i{c}_i/k_i}$, $u=\omega {\alpha }_{1}L$, and $v=\omega {\alpha }_{2}L$. A separated solution with time factor $e^{-{\omega }^{2}t}$ and $\omega \ne 0$ exists precisely when the displayed interface determinant vanishes. Dividing by $\cos u\cos v$ gives a tangent characteristic equation only away from its poles. Common zeros of both cosines are legitimate normal modes and must be retained. The zero-frequency case gives only the trivial solution of a differential equation with these positive material parameters and exterior boundary conditions.

An energy identity is an exact balance obtained by differentiating a chosen functional along a solution of a differential equation. The functional often is a squared L2 norm, $D$ collects dissipative terms, and $F$ collects forcing or reaction terms. Dropping nonnegative terms gives an energy estimate.

For $u_t=u_{xx}+g$ on $(0,1)$ with $u_x(0,t)=\alpha u(0,t)$ and $u(1,t)=0$, integration by parts provesThe boundary term is $[u{u}_x{]}_0^1=-\alpha u(0,t{)}^2$; the minus sign at the left endpoint is essential. For $\alpha \ge 0$ it contributes dissipation. If $ug\le a(t)u^2$ and the initial data vanish, the Gronwall inequality forces $u\equiv 0$.

Use the Malthusian delay differential equation with time measured in units of its constant delay:Here $a$ is a real constant. On the first interval $y(t-1)=1$, hence $y(t)=1+at$. On the next interval the delayed input is $1+a(t-1)$, already known, so the method of steps for a delay differential equation givesEach polynomial integrates the preceding interval's delayed polynomial, and its integration constant ensures continuity. At $a=1$, the step endpoints are $y(1)=2$, $y(2)=7/2$, and $y(3)=37/6$. At $t=3/2$, the value is $21/8$, but its derivative is the earlier value $y(1/2)=3/2$. This is the geometric distinction from the ordinary differential equation $y^{\prime }=y$, whose derivative would equal its current value.

More generally, for any integer $n\ge 0$,At an integer join the newly introduced term vanishes. Inside a step, differentiation lowers each power by one and gives exactly $a$ times the preceding step evaluated at $t-1$. Together with $y(0)=1$, this verifies the formula by mathematical induction; it is an exact construction, not a truncated power series.

The history and the new equation need not prescribe matching derivatives at the starting time; their mismatch can reappear later in higher derivatives.^[ref] In this example, $y^{\prime }(0^{-})=0$ and $y^{\prime }(0^{+})=a$. For $a\ne 0$, $y^{\prime }$ is continuous at $t=1$ but $y^{\prime \prime }$ jumps there by $a^2$; at $t=2$, the first jump is in $y^{\prime \prime \prime }$, of size $a^3$. The solution of a differential equation is therefore continuously differentiable for positive time, although it need not have continuous derivatives of every order at step boundaries.

The real exponential function $\exp (x)=e^x$ converts addition into multiplication:Its values are positive, and it defines a bijection from the real numbers to the positive real numbers. The natural logarithm is its inverse function. These properties connect arithmetic with calculus: adding an amount to the input multiplies the output by a fixed factor, and the local rate of change equals the current value. The exponential series gives one construction; an equivalent characterization is the unique solution of a differential equation $y^{\prime }=y$, $y(0)=1$.

For a concrete example, $y(t)=10e^{0.2t}$ starts at $10$ and satisfies $y^{\prime }=0.2y$. The absolute derivative increases with $y$, while the relative rate $y^{\prime }/y$ stays $0.2$. At $t=5$, $y=10e\approx 27.1828$ and $y^{\prime }\approx 5.43656$. The exponential function therefore models accumulation proportional to what is already present, rather than addition of a fixed amount at every step. A positive constant rate gives exponential growth; a negative one gives exponential decay.

Every positive real base $b\ne 1$ fits the same construction through $b^x=e^{x\ln b}$. Thus its derivative is $(\ln b)b^x$: base $e$ is special because $\ln e=1$. The exponential function is also unchanged in form by differentiation and integration, helping explain its repeated appearance in linear ordinary differential equations. Its real definition, inverse relationship, and derivative rules are developed in OpenStax, integrals, exponential functions, and logarithms.

A positive quantity undergoes exponential growth at constant rate $r>0$ when $y^{\prime }=ry$. For $y_0>0$, the solution of a differential equation is $y(t)=y_0{e}^{rt}$. Every interval of duration $\Delta t$ multiplies the value by the same factor $e^{r\Delta t}$, independently of its starting value.

Under the Picard-Lindelöf theorem's hypotheses on $U\subset \mathbb{R}\times {\mathbb{R}}^d$, extend a solution of a differential equation as far as possible, to a maximal interval $(T_{-},T_{+})$. If $T_{+}<\infty$, the time-state trajectory $(t,y(t))$ must eventually leave every compact subset of $U$. Equivalently, confinement near $T_{+}$ to one compact subset of $U$ permits continuation past $T_{+}$. See Teschl, Section 2.6, Lemma 2.14 and Corollaries 2.15–2.16, pp. 52–53.

The domain restriction matters even for a bounded trajectory. Restrict $y^{\prime }=1$, $y(0)=0$, to $U=\mathbb{R}\times (-1,1)$. Then $y(t)=t$ exists in this domain only for $-1<t<1$; at $t=1/2,3/4,9/10$ its values are $1/2,3/4,9/10$. As $t\uparrow 1$, it approaches the excluded boundary point $(1,1)$ while remaining bounded. Its trajectory cannot fit into any compact subset of $U$ near that endpoint. This is escape from a finite domain, without finite-time blowup in norm. Enlarging the domain to $\mathbb{R}\times \mathbb{R}$ gives the same formula $y(t)=t$ for every real $t$.

Conversely, $y^{\prime }=y$, $y(0)=1$, stays in the compact set $[0,T]\times [1,e^T]$ for every fixed $T>0$. Its growing exponential function is unbounded over the entire future, but has no finite-time blowup. At $T=2$, its largest value on $[0,2]$ is only $e^2\approx 7.38905610$. The relevant question is boundedness before each finite endpoint, rather than one bound valid forever. The tangent function in the diagram under the Picard-Lindelöf theorem illustrates the other possibility: its norm escapes every bound before $T=\pi /2$.

This criterion is useful because an estimate can secure continuation without an explicit formula. In a finite-dimensional vector space, a bounded closed ball is compact; a bounded state trajectory therefore supplies the needed compact set when its time-state closure is inside the domain. That final proviso prevents mistaking a domain boundary for finite-time blowup.

A common sufficient estimate comes from global Lipschitz continuity in the state. Assume $U=I\times {\mathbb{R}}^d$ for an open time interval $I$, with $f$ continuous and a Lipschitz constant $L\ge 0$ valid for all states on each compact time interval under consideration. For $[t_0,T]\subset I$, let $B={\sup }_{t_0\le s\le T}\Vert f(s,0)\Vert$, finite by continuity. Then $\Vert f(t,y)\Vert \le B+L\Vert y\Vert$, and the Gronwall inequality gives, for $t_0\le t\le T$ wherever the solution already exists,For instance, $y^{\prime }=y+1$, $y(0)=0$, has $L=B=1$ and the exact solution $y(t)=e^t-1$, attaining this estimate; at $t=1$ the value is $e-1\approx 1.71828183$. Finite bounds of this kind give continuation throughout $I$, with the analogous estimate backwards in time. When $I=\mathbb{R}$, this gives existence for all real times. If $I$ has a finite endpoint, it gives no continuation beyond that prescribed domain. The underlying linear-growth criterion is ^[ref].

A brief explanation of the continuation mechanism suffices: within the protecting compact set, the continuous function $f$ is bounded, say by $M$. The integral equation yields $\Vert y(t)-y(s)\Vert \le M|t-s|$ near the finite endpoint, so $y(t)$ has a limit $y_{\ast }$. The protecting compact set puts $(T_{+},y_{\ast })$ inside $U$. Restarting the initial value problem there and using local uniqueness joins the trajectories.

A solution of a differential equation has finite-time blowup in norm when $\Vert y(t)\Vert \to \infty$ as $t$ approaches a finite endpoint $T$ of its maximal interval of existence. This is stronger than merely approaching an excluded boundary of an open state domain while staying bounded. The finite-dimensional continuation criterion distinguishes these obstructions; for a locally Lipschitz continuous ordinary differential equation on $I\times {\mathbb{R}}^d$, a finite maximal endpoint inside $I$ necessarily entails finite-time blowup.

A simple example with a smooth right-hand side isThe derivative identity $(\tan t{)}^{\prime }=1+{\tan }^{2}t$ verifies the equation and the initial condition. Also $|u^2-v^2|=|u-v||u+v|\le 2R|u-v|$ for $|u|,|v|\le R$, giving local Lipschitz continuity and therefore uniqueness. Yet the positive blowup time is $T=\pi /2\approx 1.57079633$. Computed values show the accelerating growth:The right panel under the Picard-Lindelöf theorem deliberately clips its vertical axis at $10$: the trajectory has already left the picture at $t=\arctan 10\approx 1.47112767$, before reaching $T$. A plot ceiling is a display limit, not a finite value at blowup.

The tangent function also quantifies the last part of the trajectory. With $\delta =T-t>0$,Thus the height is approximately the reciprocal of the remaining time. At remaining times $\delta =0.1,0.01,0.001$, the heights are approximately $9.96664442,99.99666664,999.99966667$. These follow from the sine and cosine limits at zero and explain the steep rise without relying on the appearance of a clipped function graph.

The distinction is useful when interpreting a solution of a differential equation: uniqueness certifies a trajectory while it exists, while growth estimates decide how far it can be followed. For $y^{\prime }=1+y^2$, the Lipschitz constant $2R$ on a bounded state interval grows with $R$, and the right-hand side grows quadratically. There is no global Lipschitz constant or linear-growth bound of the kind used by the finite-dimensional continuation criterion. Smoothness at every finite state is entirely compatible with finite-time blowup.

A Floquet multiplier is an eigenvalue of a one-period linear evolution map. For a homogeneous linear differential equation with continuous coefficients of period $T>0$, this map is the monodromy matrix of a periodic linear system $N$. If $Nv=mv$ for a nonzero eigenvector $v$, the associated solution of a differential equation with $Y(0)=v$ satisfies $Y(t+T)=mY(t)$. Its exponential rate is $T^{-1}\log |m|$, the corresponding Floquet growth rate. A Floquet multiplier equal to one gives a $T$-periodic solution in that eigenvector direction; a nontrivial Jordan block at a unit-modulus Floquet multiplier can still produce polynomial growth.

For a nonlinear periodic orbit, the Floquet multipliers are eigenvalues of the one-period map of its variational equation. A nonconstant periodic orbit of a smooth autonomous differential equation has a Floquet multiplier equal to one in its tangent direction: differentiating the solution of a differential equation with respect to its starting phase gives a periodic solution of the variational equation. The derivative of a transverse Poincaré map removes this phase direction and retains the transverse Floquet multipliers. Transverse Floquet multipliers strictly inside the unit circle give local attraction, while one outside gives instability; a transverse Floquet multiplier on the unit circle requires further analysis. For a planar periodic orbit, the transverse Floquet multiplier is given by the divergence test for a planar periodic orbit.

A history function of a delay differential equation supplies the past interval needed to start a retarded delay differential equation. A single initial condition $y(0)=\varphi (0)$ generally leaves earlier delayed values unspecified.^[ref]

For a continuous function history, the natural state at time $t$ is the entire segmentIt belongs to the function space $C([-\tau ,0],\mathbb{R})$, an infinite-dimensional vector space: the functions $1,\theta ,{\theta }^2,\dots$ are linearly independent. Shifting $y_t$ forward retains the part of the past that the equation can still read. Thus even one scalar unknown can require a state in an infinite-dimensional vector space.

The insufficiency of one value is visible without any approximation. Consider $y^{\prime }(t)=y(t-1)$ with the two history functions of a delay differential equationBoth give $y(0)=1$, yet on the first step the respective derivatives are $1$ and $t$. Direct integration givesIn particular $y_A(1)=2$ while $y_B(1)=3/2$. On the next step,At $t=2$ they are $7/2$ and $8/3$, a gap of $5/6$. These different futures reflect different prescribed histories; each fully specified problem has a unique solution of a differential equation. More strongly, for $y^{\prime }(t)=ay(t-1)$ with $a\ne 0$, observing the entire first-step derivative recovers every history value through $\varphi (s)=y^{\prime }(s+1)/a$ for $-1\le s\le 0$, using the right derivative at $0$. The future can depend on the whole function, rather than a finite list of its samples.

A homogeneous linear differential equation is a linear differential equation with identically zero forcing term, $L[y]=0$. On a common interval, its solutions of a differential equation form the vector space $\ker L$, the kernel of a linear map. They include the zero function and are closed under addition and scalar multiplication. For any two solutions of a differential equation $u,v$ and constant scalars $\alpha ,\beta$, linearity gives $L[\alpha u+\beta v]=\alpha 0+\beta 0=0$: this is the superposition principle. The scalars must be constant; variable multipliers introduce extra derivatives.

For example, $y^{\prime }+2y=0$ has general solution $y=C{e}^{-2t}$, where $C$ is any constant real number. A short verification uses the product rule: $(e^{2t}y{)}^{\prime }=e^{2t}(y^{\prime }+2y)=0$, so $e^{2t}y$ is constant. For $y_1=e^{-2t}$ and $y_2=2e^{-2t}$, the sum is $3e^{-2t}$ and also solves the same homogeneous linear differential equation. At $t=1$, their values are approximately $0.13534$, $0.27067$ and $0.40601$. Every curve is a multiple of one exponential function; the solution space of a homogeneous linear differential equation has dimension one. Both $t$ and $y$ in this example and the following diagram are dimensionless mathematical variables.

Now keep the same differential operator $L[y]=y^{\prime }+2y$ but prescribe the forcing term $1$. The inhomogeneous linear differential equation $y^{\prime }+2y=1$ has particular solution $y_p=1/2$ and general solution $y=1/2+C{e}^{-2t}$. The initial conditions $y(0)=0,1,2$ select, respectively,At $t=1$, their values are approximately $0.43233$, $0.56767$ and $0.70300$. The midpoint $\frac{1}{2}(u_0+u_2)=u_1$ solves the same inhomogeneous linear differential equation, because $L[\frac{1}{2}(u_0+u_2)]=\frac{1}{2}(1+1)=1$. The sum $s=u_0+u_2=1+e^{-2t}$ instead satisfies $s^{\prime }+2s=2$, so it fails the equation with forcing term $1$. More generally, a linear combination of solutions of a differential equation for the same fixed nonzero forcing term retains that forcing term precisely when its scalar weights add to one. This is the affine space structure, rather than unrestricted superposition within that fixed family.

On the left, all three curves solve the same homogeneous linear differential equation. On the right, $u_0,u_1,u_2$ share the particular solution $1/2$, and the shaded difference $u_2-1/2$ is the homogeneous solution $\frac{3}{2}{e}^{-2t}$. The dashed sum has twice the forcing term, even though its starting value coincides with that of $u_2$. Adding the two initial conditions also adds their inputs; taking their midpoint averages both. This distinction makes the superposition principle useful for assembling responses to different inputs without mistaking that assembly for a solution of a differential equation with the original input. See MIT's superposition notes, including the particular solution one half.

An implicit ordinary differential equation relates the derivative to the independent and dependent variables without requiring it to be explicitly isolated. For example $y{y}^{\prime }=1$ is meaningful at every pair $(x,y)$ as an equation, but no differentiable solution of a differential equation can pass through $y=0$: its left side would vanish there. Dividing by $y$ gives $y^{\prime }=1/y$ only on the domain $y\ne 0$. Conversely multiplying an explicit equation by a vanishing denominator can introduce points where the original equation was undefined. A solution of the multiplied equation at such a point need not solve the original initial value problem.

An implicit solution describes a solution of a differential equation through a relation $F(x,y)=C$, instead of isolating $y$ in an explicit formula. The relation describes a level set; an actual scalar solution of a differential equation must be a differentiable function $y(x)$ on an interval whose function graph lies in the equation's domain. A whole level set need not be one such function graph.

For example, considerImplicit differentiation of $x^2+y^2=C$ gives $2x+2y{y}^{\prime }=0$. For $C>0$, the circle contains two solutions of a differential equation,Both satisfy the equation on the stated open interval, and their function graphs stay in $D$. An initial condition $(x_0,y_0)\in D$ fixes $C=x_0^2+y_0^2$ and selects the sign of $y_0$. In particular, $y(0)=2$ gives $C=4$ and the upper branch on $(-2,2)$. The constant alone does not select the branch.^[ref]

The precise local criterion uses the implicit function theorem. Let $F$ be continuously differentiable on an open subset of the plane, with $F(x_0,y_0)=C$. If its partial derivative $F_y(x_0,y_0)\ne 0$, the nearby level set is a unique continuously differentiable function graph $y=g(x)$ through that point. The chain rule givesThis proves that the local function graph solves $y^{\prime }=-F_x/F_y$ there. For a different prescribed equation $y^{\prime }=f(x,y)$, one must also check $F_x+F_{y}f=0$ along that branch and that it stays in the domain of $f$. A level set by itself does not verify an arbitrary ordinary differential equation.^[ref]

For $F=x^2+y^2$, the upper branch through $(1,\sqrt{3})$ has $F_y=2\sqrt{3}$ and tangent line slope $-1/\sqrt{3}$. At $(2,0)$, however, $F_y=0$ and $F_x=4$. The circle remains a smooth regular level set with a vertical tangent line: it can be written locally as $x=\sqrt{4-y^2}$, but not as a single function graph $y=g(x)$ over an open interval containing $2$. Moreover, $-x/y$ is undefined there, so the scalar solution of a differential equation cannot include this point in its domain. Changing the parameter draws the curve through the point but does not extend that scalar initial value problem.

When both partial derivatives vanish, even the geometry needs separate inspection. The level $x^2+y^2=0$ is an isolated point, not a function graph on any nonempty open interval. The zero level set of $G(x,y)=y^2-x^2$ instead consists of the two lines $y=x$ and $y=-x$ crossing at the origin, where $\nabla G=0$. Each line satisfies the implicit ordinary differential equation $y{y}^{\prime }=x$, including at the origin; the explicit form $y^{\prime }=x/y$ excludes the origin. Thus algebraically dividing by a vanishing coefficient changes the domain and can change which solutions of a differential equation are admissible. Vanishing $F_y$ alone is only failure of a sufficient test: $F(x,y)=y^3$ has the smooth zero function graph $y=0$ despite $F_y=0$ everywhere on it.

Implicit solutions are especially useful for exact first-order ordinary differential equations and for finding a first integral of an ordinary differential equation: a conserved relation can reveal allowed states without an explicit formula. For the normalized harmonic oscillator written as $\dot{x}=y$, $\dot{y}=-x$, the energy $H=(x^2+y^2)/2$ has $\dot{H}=x\dot{x}+y\dot{y}=0$. Its positive level sets are the same circles. Where $y\ne 0$, eliminating time gives $dy/dx=-x/y$ and the local implicit solution above. The harmonic oscillator's parametric curve continues through $y=0$ because its vector field governing evolution in time is defined there; treating $x$ as the independent variable fails at that turn. The origin is an equilibrium point of the harmonic oscillator, although it is excluded from the scalar equation $y^{\prime }=-x/y$.

An initial value problem combines an ordinary differential equation with an initial condition, for example $y^{\prime }=f(t,y)$ and $y(t_0)=y_0$. The task is to find a solution of a differential equation on an interval containing $t_0$ which has that starting value. For a continuous linear differential equation $u^{\prime }=A(t)u$, every prescribed finite vector $u_0$ determines a unique solution on the coefficient interval. A boundary value problem can instead prescribe values at distinct endpoints.

For the cooling initial value problem $y^{\prime }=-k(y-T_a)$, $y(t_0)=y_0$, where $k>0$ and $T_a$ are constant, the unique solution isTake $T_a=20$ degrees Celsius, $y_0=80$ degrees Celsius, $t_0=0$, and $k=0.1$ per minute. Then $y(10)=20+60e^{-1}\approx 42.07$ degrees Celsius. Reaching $30$ degrees Celsius requires $60e^{-0.1t}=10$, so $t=10\log 6\approx 17.92$ minutes. The time constant $1/k=10$ minutes reduces the excess temperature by a factor $e^{-1}$; it does not reduce the absolute temperature by that factor.

Uniqueness needs assumptions. If $f$ is a continuous function on an open neighborhood of $(t_0,y_0)$ and satisfies a local Lipschitz condition in $y$, uniformly in $t$, the Picard-Lindelöf theorem gives a unique local solution. Concretely, a finite $L$ bounds $|f(t,u)-f(t,v)|\le L|u-v|$ in a sufficiently small rectangle. Cooling satisfies this with $L=k$. Local existence does not promise existence for all future time: $y^{\prime }=y^2$, $y(0)=1$, has $y(t)=1/(1-t)$ on $(-\infty ,1)$ and finite-time blowup at $t=1$.

Continuity alone is insufficient for uniqueness. Consider the dimensionless initial value problem involving a square root and an absolute value:Besides $y(t)\equiv 0$, each delay $a\ge 0$ gives a solution for $t\ge 0$:At $t=4$, delays $a=0,1,2$ give $y=4,9/4,1$, despite identical starting data. These are valid solutions of a differential equation: for $a>0$, both one-sided derivatives at the join are zero; for $a=0$, the right derivative at the starting time is zero. After the join, $y_a^{\prime }=(t-a)/2=\sqrt{|y_a|}$ because $t>a$. Squaring alone would also suggest a parabola before the join, but its negative derivative there would violate the equation. Since $y^{\prime }\ge 0$, a forward solution starting from zero stays nonnegative.

The right-hand side is a continuous function but fails the Lipschitz condition at zero: for $h>0$, $|\sqrt{|h|}-0|/|h|=1/\sqrt{h}$ is unbounded. Thus the Peano existence theorem still supplies local existence, while the Picard-Lindelöf theorem does not supply uniqueness. Related waiting-time examples and a sufficient uniqueness condition appear in MIT 18.034, Lecture 2: Fundamental Principles.

A brief proof outline explains the mechanism: rewrite the initial value problem as the integral equation $y(t)=y_0+{\int }_{t_0}^{t}f(s,y(s))ds$. On a short enough interval, the Lipschitz condition makes the right-hand side a contraction mapping on a closed set of continuous functions near $y_0$. The contraction mapping theorem gives a unique fixed point; the fundamental theorem of calculus turns it back into a solution.

An integrating factor is a nonzero multiplier that makes a differential equation directly integrable. For a first-order linear differential equation, it gathers the unknown function and its derivative into one product rule expression. Suppose $a$ and $b$ are continuous functions on an interval containing $t_0$, and consider the initial value problemDefine $A(t)={\int }_{t_0}^{t}a(s)ds$ and $\mu (t)=e^{A(t)}$. The exact solution of a differential equation throughout that interval isThe definite integral fixes $\mu (t_0)=1$, so the initial condition appears without an extra constant to determine. If the equation originally has a coefficient multiplying $y^{\prime }$, first divide by that coefficient on an interval where it is nonzero; a zero of that coefficient needs separate attention. See OpenStax's first-order linear equations chapter for this normalization and the integrating factor procedure.

For a concrete balance, solve $y^{\prime }+2y=3$ with $y(0)=0$. Here $\mu =e^{2t}$, so multiplying the entire ordinary differential equation gives $(e^{2t}y{)}^{\prime }=3e^{2t}$. Accumulating from zero to $t$ givesThe computed values are $y(0.5)\approx 0.94818$, $y(1)\approx 1.29700$ and $y(2)\approx 1.47253$. Its derivative is $3e^{-2t}$, which verifies $y^{\prime }+2y=3$ exactly. Interpreting $3$ as a steady input and $2y$ as proportional removal explains the shape: initially all input increases $y$, then removal catches up. The constant particular solution $y_p=3/2$ balances input and removal, and the homogeneous solution $y_h=-\frac{3}{2}{e}^{-2t}$ is the transient response that corrects it to the initial condition. Because this autonomous differential equation has a time-independent right-hand side, $3/2$ is also an equilibrium point.

The first panel adds the particular solution and homogeneous solution; the second shows the same balance in rates, $3=y^{\prime }+2y$. This links the algebra to a physical accounting rule. Newton's law of cooling has this form with $a=k$ and $b=k{T}_a$, while an electrical circuit can have a varying input. In both cases the integrating factor systematically handles the input rather than requiring a guess for a particular solution.

Variable coefficients work just as directly. For $y^{\prime }+2ty=t$, $y(0)=0$, the integrating factor is $e^{t^2}$, andThus $y(0.5)\approx 0.11060$, $y(1)\approx 0.31606$ and $y(2)\approx 0.49084$. The derivative $y^{\prime }=t{e}^{-t^2}$ verifies $y^{\prime }+2ty=t$, including at zero, where both coefficient and input vanish. No division by $t$ is necessary. The general solution is $y=1/2+C{e}^{-t^2}$: $1/2$ is a constant particular solution, and $C{e}^{-t^2}$ solves the associated homogeneous linear differential equation. Although this special forcing term admits a constant trajectory, this is a nonautonomous differential equation: it depends explicitly on time. Merely setting its right-hand side to zero at one instant would not establish a constant solution of a differential equation; at $t=0$ every state has zero slope.

The direction field shows the selected trajectory leaving its initially horizontal tangent line. For $t\ge 0$, the difference between any two solutions of a differential equation is multiplied by $e^{-t^2}$, rather than $e^{-2t}$ in the constant-rate example. This dependence on elapsed history is clearer if the exact formula is rewritten asFor $t\ge t_0$, each contribution from the input at $s$ is weighted by the accumulated coefficient between $s$ and $t$. The homogeneous solution remembers the initial condition; the particular solution given by the second term accumulates the forcing term with zero initial value. This decomposition differs from choosing a constant particular solution, but gives the same total. Initial-data memory tends to zero when ${\int }_{t_0}^{t}a\to +\infty$; a negative accumulated coefficient instead amplifies it. The reciprocal relation between integrating factor and homogeneous solution is developed in MIT's notes on solutions to linear first-order equations.

Why does the trick work? The fundamental theorem of calculus and chain rule give ${\mu }^{\prime }=a\mu$. Consequently the product rule givesIntegrating this identity and dividing by the positive exponential function $\mu$ produces the stated formula. Any nonzero constant multiple of $\mu$ works equally well; normalization is a convenience. This is the beauty of the method: a sum involving an unknown function becomes one exact derivative, leaving only an integral of known functions. See Oregon State's product-rule derivation. The related integrating factor for a differential one-form extends the idea of making an expression exact beyond this linear setting.