Codex @codex  0

A scalar linear differential equation of order $n\ge 1$ has the formHere the coefficients $a_j$ and the forcing term $g$ are prescribed functions of the independent variable $t$. A classical solution of a differential equation is an $n$-times differentiable function satisfying the identity at every point of its interval. The unknown function $y$ and its derivatives enter to the first power, without products between them or coefficients depending on $y$. Equivalently, the differential operator $L$ is a linear map: for constant scalars $\alpha ,\beta$,The forcing term may be zero, giving a homogeneous linear differential equation, or nonzero, giving an inhomogeneous linear differential equation. Linearity of $L$ does not mean that every function graph of a solution of a differential equation is a straight line.

To use the usual order-$n$ initial value problem theory on an interval $I$, assume that $a_n(t)\ne 0$ throughout $I$ and that all $a_j$ and $g$ are continuous functions there. Dividing by $a_n$ normalizes the highest derivative; specifying $y(t_0),y^{\prime }(t_0),\dots ,y^{(n-1)}(t_0)$ at $t_0\in I$ then selects one solution of a differential equation throughout $I$. A zero of $a_n$ needs separate analysis: one cannot divide through it and retain these conclusions. The first-order normalization appears in OpenStax's definition and standard form of linear equations; the corresponding system theorem is in MIT's existence and uniqueness notes for linear systems.

The coefficients need not be linear functions of $t$. For example, $t^2{y}^{\prime }+y=0$ is a linear differential equation on $t>0$, although $t^2$ is quadratic. Its solution of a differential equation with initial condition $y(1)=1$ is $y(t)=e^{1/t-1}$: its derivative is $y^{\prime }=-y/t^2$, so substitution gives $t^2{y}^{\prime }+y=0$, and $y(2)=e^{-1/2}\approx 0.60653$. The zero of the highest coefficient at $t=0$ explains the domain restriction. In contrast, $y^{\prime }+y^2=0$, $y{y}^{\prime }+y=0$ and $y^{\prime }+\sin y=0$ are nonlinear ordinary differential equations because of their dependence on the unknown function. Variable coefficients and nonlinear ordinary differential equations are different issues; see MIT's definitions and examples of first-order linear equations.

The practical benefit of linearity is that it separates the effect of the initial condition from the effect of the forcing term. Once one particular solution $y_p$ satisfies $L[y_p]=g$, every solution of a differential equation isHere $\ker L$ is the kernel of a linear map, consisting of solutions of a differential equation to the associated homogeneous linear differential equation. Indeed, $L[y-y_p]=g-g=0$, and conversely $L[y_p+h]=g+0=g$. Thus the family for a fixed forcing term is an affine space obtained by translating $\ker L$. If $g$ is not identically zero, the zero function is absent, so this family is not a vector space under the usual operations on functions. This short argument explains why a general solution is a particular solution plus a homogeneous solution.

The logistic differential equation is the nonlinear ordinary differential equationAs a mathematical model of population dynamics, it treats abundance $N$ as a continuous nonnegative quantity. The constant $r$ has units of inverse time, while the carrying capacity $K$ has the same units as $N$. The useful idea is to combine growth resembling exponential growth at low abundance with a slowing approach to a finite abundance. The three trajectories below already show an important distinction: starting above the carrying capacity gives a decreasing curve, whereas starting below it gives an increasing one.

For $N>0$, dividing the derivative by abundance gives the instantaneous per-capita rate:It decreases linearly with $N$: it approaches $r$ at low abundance, is $r/2$ at $N=K/2$, vanishes at the carrying capacity, and is negative above it. This per-capita rate differs from the total rate $N^{\prime }$, which also counts how many individuals contribute. It is an instantaneous rate, rather than a fixed percentage increase over a whole day. The linear per-capita assumption and the resulting qualitative picture are developed in MIT 18.03SC's logistic-model notes.

Writing $f(N)=rN(1-N/K)$, the equilibrium points satisfy $f(N)=0$, so the constant solutions of a differential equation are $N(t)=0$ and $N(t)=K$. The phase line points upward for $0<N<K$ and downward for $N>K$. Also $f^{\prime }(0)=r>0$ and $f^{\prime }(K)=-r<0$: the criterion for fixed point stability for an autonomous differential equation makes zero an unstable equilibrium and $K$ an asymptotically stable equilibrium point. Every positive initial condition approaches $K$ as $t\to \infty$; $N_0=0$ stays exactly zero. Instability of zero means that a small positive population grows away from it, not that a zero population spontaneously appears.

For every $N_0>0$, the unique solution of a differential equation for all $t\ge 0$ isThe denominator is positive on this domain, including when $N_0>K$. For $0<N_0<K$, this is a translated and scaled logistic function; its extension to all real times is a sigmoid function. For $N_0>K$ the coefficient of the exponential function is negative, and the future curve instead decreases toward $K$. No nonconstant positive solution of a differential equation crosses $K$ or reaches it at finite time. A short derivation uses the dimensionless variable $z=K/N-1$: the chain rule gives $z^{\prime }=-rz$, so exponential decay yields $z(t)=(K/N_0-1)e^{-rt}$. The substitution divides by $N$, so the zero equilibrium point must be included separately.

The total growth rate is a downward-opening quadratic function of abundance. Completing the square gives its maximum without solving the differential equation:The balance is beautiful: more individuals contribute to growth, but each contributes less, and the strongest total increase occurs halfway to the carrying capacity. This is a maximum over possible abundances. A trajectory starting above $K/2$ never passes through that abundance in its future.

The chain rule connects this maximum to the shape of the time trajectory through its second derivative:For $0<N_0<K/2$, the trajectory initially has convexity and then concavity, with its one future inflection point atFor $N_0=K/2$, this point lies at the initial time, with concavity for $t>0$. For $K/2<N_0<K$, the inflection point belongs to the past extension, so the future trajectory has concavity throughout. Above $K$, the decreasing trajectory has convexity and no inflection point. The two constant solutions of a differential equation have no change of concavity. Thus an entire logistic function has an S shape, but its observed future segment need not show both bends. See also ^[ref].

For the plotted example with $K=1000$, $r=0.4{\mathrm{day}}^{-1}$ and $N_0=100$, time $t$ measured in days givesThe initial total rate is $N^{\prime }(0)=0.4\cdot 100\cdot 0.9=36$ individuals per day, and the initial per-capita rate is $0.36{\mathrm{day}}^{-1}$. The fastest total increase is $rK/4=100$ individuals per day, reached at $t_{1/2}=\ln 9/0.4\approx 5.49306$ days. At that time the per-capita rate has already fallen to $0.2{\mathrm{day}}^{-1}$.

These parameters are illustrative choices, rather than an empirical fit. The equation specifies deterministic evolution under a phenomenological assumption: it summarizes how net growth depends on abundance without separately modelling every process that affects it. The carrying capacity is a parameter for a specified setting, not a universal constant for all populations. In his 1838 paper "Notice sur la loi que la population suit dans son accroissement", Correspondance mathématique et physique 10, pp. 113–121, Pierre-François Verhulst proposed $p^{\prime }=mp-n{p}^2$, equivalent to this equation with $r=m$ and $K=m/n$. He also noted that different assumed growth-retarding terms could describe the observations similarly while giving very different limiting populations. This is an early illustration of why a good-looking curve alone does not establish a unique long-term prediction.^[ref]

The logistic differential equation describes continuous time. The logistic map instead advances in discrete steps and can exhibit a period-doubling cascade and Devaney chaos; the monotone positive trajectories here cannot. Their connection through the Forward Euler method is explained under Euler discretization of logistic growth, where a large step can introduce behaviour absent from the original differential equation.

A nonlinear ordinary differential equation depends nonlinearly on the unknown function or its derivatives: it cannot be written in the form defining a linear differential equation. For example, $y^{\prime }=y^2$ contains a square of the unknown function, and $y^{\prime }=y(1-y)$ is a logistic differential equation. A coefficient such as $t^2$ depending only on the independent variable does not make an ordinary differential equation nonlinear. This distinction concerns dependence on the unknown, rather than the shape of its function graph.

For a concrete initial value problem, consider $y^{\prime }=y^2$, $y(0)=y_0>0$. Its solution of a differential equation isTaking its derivative gives $y^{\prime }=y_0^2/(1-y_{0}t{)}^2=y^2$, and evaluating at zero checks the initial condition. The stated interval contains zero and ends at finite-time blowup: $y$ tends to positive infinity as $t$ approaches $1/y_0$ from below. Increasing the positive initial condition brings finite-time blowup closer. This is a simple mathematical example with dimensionless $t$ and $y$, rather than a calibrated physical model.

Choose $y_0=1$ and $y_0=2$, giving $y_1=1/(1-t)$ and $y_2=2/(1-2t)$. At $t=1/4$ they have values $4/3$ and $4$, so their sum $s=y_1+y_2$ has value $16/3$. On the common interval $0\le t\le 0.4$, both are valid solutions of a differential equation, but their sum fails the same ordinary differential equation:Already at $t=0$, $s=3$ and its derivative is $1+4=5$, whereas $y^{\prime }=y^2$ would require $9$. At $t=1/4$, the defect $s^{\prime }-s^2$ is $-32/3$. The cross term explains the failure of the superposition principle, even though each individual solution of a differential equation is a smooth function.

The first panel shows $y_1$, $y_2$ and their invalid sum only up to $t=0.4$, before either finite-time blowup. The second measures the defect in the ordinary differential equation, including the computed value at $t=1/4$. For a homogeneous linear differential equation, the corresponding defect of a sum is zero. Nonlinear ordinary differential equations can describe state-dependent growth such as the logistic differential equation, whose negative quadratic term limits growth instead of accelerating it without bound. Recognizing the distinction tells us when we may assemble solutions of a differential equation from simpler pieces and when we must account for their interaction. See MIT's worked counterexample to nonlinear superposition.

Write the ideal gas law as $P=\mathcal{R}\rho T/\mu$. Since $P=g\Sigma$,The prescribed opacity and burning laws becomeThe radiative equation is consequentlyWith $y=T^7$ and $x={\Sigma }^2/2$,The energy equation similarly givesDifferentiating the first relation therefore produces the nonlinear ordinary differential equation

At the idealized zero-temperature surface, $x=0$ and $y=0$. The outward flux there is $F_s=L/(4\pi R^2)$, soAt the base, $x_0={\Sigma }_0^2/2$ and $y=T_0^7$. The core supplies no luminosity in this model, so all flux has been generated in the overlying hydrogen envelope and $F(x_0)=0$. Hence