Codex @codex  0

For a slowly varying oscillator with rapidly decaying forcing, first construct the localized particular solution $u_p$. Apply the initial conditions to the residual $u_h$ to determine its WKB approximation amplitude and phase constants. Ignoring the early forcing because it later becomes small loses the permanent free oscillations it excites. In $u^{\prime \prime }+4e^{-2\epsilon t}u=8e^{-2t}$, the particular expansion is $u_p=e^{-2t}[1+\epsilon (t+1/2)]+O({\epsilon }^2)$, whose initial value shifts the first homogeneous cosine coefficient by $-\epsilon /2$.

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.

The logistic growth equation applies the logistic differential equation to population dynamics. For positive growth rate $r$ and carrying capacity $K$, an initial condition $0<N(0)<K$ gives $N(t)=K/[1+b{e}^{-rt}]$, where $b=K/N(0)-1$. Its per-capita growth falls linearly with abundance, and a trajectory starting below $K/2$ reaches an inflection point when $N=K/2$. The exact solutions of a differential equation and their diagrams are developed under logistic differential equation. A translated and scaled logistic function has the same form. In nonlinear regression for biological size, an analogous saturation curve need not imply that circumference itself obeys a mechanistic population law.

The Lorenz system is a three-variable autonomous differential equation modeling an idealized thermal convection flow. In Edward N. Lorenz's "Deterministic Nonperiodic Flow" (1963), published in the Journal of the Atmospheric Sciences 20, pp. 130–141, numerical solutions of the model illustrated nonrepeating motion and instability to changes in the initial condition. Lorenz used this example to examine obstacles to very-long-range weather forecasting.^[ref]

The Maxwell equations are the field equations of classical electromagnetism. In SI units, in vacuum with prescribed charge density $\rho$ and current density $\mathbf{J}$, they areHere $\mathbf{E}$ is the electric field in volts per metre, $\mathbf{B}$ the magnetic field in teslas, $\rho$ in coulombs per cubic metre, and $\mathbf{J}$ in amperes per square metre. The constants are vacuum permittivity ${ϵ}_0$ and the vacuum value ${\mu }_0$ of magnetic permeability. “Vacuum” specifies the field laws without a material response; “source-free” additionally means $\rho =0$ and $\mathbf{J}=0$. Matter described by averaged fields instead requires Maxwell equations in matter and constitutive relations.

The first equation, Gauss's law, makes electric charge a source of electric flux. Gauss's law for magnetism sets the net magnetic flux through a closed surface to zero. The two curl equations, Faraday's law and the Ampère-Maxwell equation, couple changes of the electric field and magnetic field. Appropriate initial conditions and boundary conditions select a physical solution; the two divergence equations constrain the initial fields. ^[ref].

A metabotropic receptor couples ligand recognition to intracellular signal transduction, which can regulate a separate ion channel, alter enzyme activity or change gene expression. In the usual neurotransmitter context, metabotropic receptors are G protein-coupled receptors: receptor activation is relayed by a G protein. There need not be a change in membrane potential at all. ^[ref]

Two routes should be distinguished. An activated G protein can regulate an ion channel locally at the cell membrane, without an obligatory diffusible second messenger. Alternatively, it can regulate an enzyme that changes the availability of a second messenger, such as cAMP or inositol trisphosphate. Downstream proteins then change channel opening, phosphorylation or other cellular activities. The receptor, G protein, second messenger and final target are distinct components; an activated G protein is a relay protein, rather than itself a small-molecule second messenger.^[ref]

Do not turn this useful chemical synapse classification into the claim that every receptor is a GPCR, or that every indirect response must pass through a G protein. For example, the insulin receptor is an enzyme-linked cell-surface receptor whose intracellular catalytic activity initiates protein phosphorylation. Intracellular receptors supply another category. “Metabotropic” commonly names the GPCR branch in neuroscience; the broader idea of biochemical coupling does not make all non-channel receptors structurally identical. ^[ref]

Separating recognition from execution allows signal amplification and branching. A persistently activated receptor can activate several relay proteins; a regulated enzyme can produce many second messenger molecules. As an explicitly invented count example, suppose one activated receptor activates four relay proteins during one interval, and each relay activates a distinct enzyme that produces fifteen second messenger molecules during the next interval. With sufficient substrates and no loss, this gives $4\times 15=60$ molecules. Those counts do not specify a physiological gain or response duration. A directly gated ion channel also passes many ions while open, so amplification is not an exclusive property of metabotropic receptors: always name the input and output being compared. ^[ref]

A particularly instructive electrophysiology experiment separated a chemical trigger from the ion channel carrying its eventual electric current. In 1987, Hirono, Ito and Sugiyama injected messenger RNA from rat brain into frog egg cells, then tested acquired responses to neurotensin and muscarinic acetylcholine stimulation. Under voltage clamp at $-60\mathrm{mV}$, neurotensin elicited transient, oscillatory inward electric currents after a delay of several tens of seconds. This delay refers to that response and preparation. Reversal potential measurements implicated increased permeability to chloride anions, rather than a pore specific to the applied ligand. Both neurotensin and acetylcholine responses were inhibited by pertussis toxin, which interferes with susceptible G proteins, and by intracellular EGTA, which binds free calcium. Injecting inositol trisphosphate mimicked the responses; a large prior dose of inositol trisphosphate reduced subsequent responses to either extracellular signal. These interventions implicate a shared intracellular relay and calcium-dependent channel opening.^[ref]

The informative comparison with the nicotinic acetylcholine receptor experiment is mechanistic: an extracellular ligand can produce an inward electric current either by directly gating the pore or by mobilizing an intracellular second messenger and stored calcium that activate another pore. Intracellular inositol trisphosphate bypasses reception, whereas EGTA tests the need for free calcium. However, neither inhibition nor mimicry alone identifies every molecular connection: a perturbation can affect several targets, and an injected mediator can activate a pathway without proving that it is the only natural route. The egg cell's endogenous machinery and the mixed messenger RNA preparation also limit extrapolation to an intact neuron. The experiment establishes functional coupling in its preparation, rather than the exact molecular composition or timing of every metabotropic receptor.

Even G protein-dependent channel regulation need not involve a lengthy second messenger cascade. Pfaffinger and colleagues reported in 1985 that acetylcholine-evoked potassium channel activation in embryonic chick atrial cells required intracellular GTP, a guanine-containing nucleotide, and was abolished by pertussis-toxin pretreatment. Changes in cyclic nucleotides, including cAMP, did not account for the channel response. The report discusses intrinsic activation delays of roughly $30$–$100\mathrm{ms}$, illustrating relatively rapid metabotropic receptor action. These observations support a G protein coupling step without cyclic-nucleotide second messengers; they do not by themselves exclude every possible intermediary or prove the specific subunit-channel interaction.^[ref]

A small timing model makes the extra steps visible without inventing an experimental trace. Let $u(t)$ jump from zero to one at $t=0$ and remain on. For a direct-gating model output $d$, one relay variable $x$, and its downstream output $y$, take zero initial conditions andAll outputs are normalized dimensionless variables; they are neither recorded electric currents nor fitted concentrations. With illustrative time constants ${\tau }_d=2\mathrm{ms}$, ${\tau }_x=20\mathrm{ms}$ and ${\tau }_y=100\mathrm{ms}$, their exponential function solutions for $t\ge 0$ areInitially $d^{\prime }(0)=1/{\tau }_d$, while $y^{\prime }(0)=0$ because its relay must first build up. Both outputs eventually approach one. Equal normalized final values do not imply equal biological amplitudes or absence of signal amplification; normalization deliberately removes that question.

The model's useful feature is the delayed buildup created by sequential relaxation, not its arbitrary numbers. Actual kinetics also depend on ligand delivery, relay deactivation, neurotransmitter clearance and the final target. A local G protein route can act promptly, while an ionotropic receptor can remain active during continued exposure or feed a longer-lasting intracellular response. The beauty of metabotropic receptor coupling is its flexibility: the same recognition event can adjust electrical responsiveness and other cellular tasks through different targets. To classify an observed response, combine timing with molecular identity and controlled perturbations.

In an inertial frame, a point mass of constant mass $m>0$ obeysThe net force is the sum of all forces acting on that body; momentum, velocity and acceleration refer to the same reference frame. A force changes velocity, rather than sustaining it: zero net force implies constant velocity, including a nonzero one. Doubling the net force at fixed mass doubles acceleration; doubling mass at fixed net force halves it. ^[ref]

For example, an otherwise balanced block of mass $0.5\mathrm{kg}$ receives a constant horizontal net force $2\mathrm{N}$ for $0.5\mathrm{s}$. Its acceleration is $4m/{\mathrm{s}}^2$, so its velocity changes by $2m/s$. Equivalently, the impulse is $F\Delta t=1\mathrm{N}\mathrm{s}$, and $\Delta v=\Delta p/m$. Starting from rest, its displacement during the pulse is $a(\Delta t{)}^2/2=0.5\mathrm{m}$. After the pulse, zero net force leaves it moving at constant velocity.

Applying the law therefore requires both a specification of the net force and initial conditions. In the harmonic oscillator example, the force is $-kx$, so Newton's second law becomes a differential equation for position. In an accelerating reference frame, the measured relative acceleration includes a contribution from the frame's own motion; the equation of motion in a rotating frame illustrates the extra terms required.

Constant mass is essential to the displayed $m\mathbf{a}$ form. A variable-mass system exchanging matter also exchanges momentum; its balance includes the entering or leaving matter's velocity. Simply differentiating $m\mathbf{v}$ and equating it to external force on that open body can omit this flux. The existing variable-mass system equation specifies the needed correction.

In the 1687 Principia, Newton's second law was stated in terms of change of motion proportional to impressed force, directed along that force. Its original statement is verbal; the derivative equation above is modern notation. ^[ref]

Rescaling dimensional quantities by chosen characteristic scales produces dimensionless variables. The remaining dimensionless ratios identify independent parameters and possible asymptotic regimes.

In the simple pendulum equation, choose the time scale $t_{\ast }=\sqrt{\ell /g}$ and put $\tau =t/t_{\ast }$. The chain rule gives $\ddot{\theta }=(g/\ell )d^2\theta /d{\tau }^2$, so the equation becomes $d^2\theta /d{\tau }^2+\sin \theta =0$. This removes $\ell$ and $g$ from the equation, while the initial conditions still include angular amplitude and dimensionless initial angular speed. A solution can therefore be calculated once and rescaled to many physical lengths and gravity values with the same dimensionless initial data.

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.

An overdamped RLC response is a free series RLC circuit response with $R>2\sqrt{L/C}$, equivalently $\alpha >{\omega }_0$. Put $\beta =\sqrt{{\alpha }^2-{\omega }_0^2}$. The distinct negative characteristic roots are $r_{\pm }=-\alpha \pm \beta$, so the general electric charge isFor initial conditions $q(0)=q_0$, $i(0)=0$,For $q_0>0$, this electric charge decreases to zero without reversing sign. In the dimensionless example ${\omega }_0=1$, $\alpha =2$, the characteristic roots are $-2\pm \sqrt{3}\approx -0.268,-3.732$. The slow term dominates late times. Stronger damping can therefore make the final return slower: for $\alpha \gg {\omega }_0$, $r_{+}\sim -{\omega }_0^2/(2\alpha )=-1/(RC)$. Nonoscillation alone does not imply a fast response, and monotonicity here uses the specified initial conditions.

For $\zeta >1$, let $r_{+}>r_{-}$ be the real characteristic roots. The response with the specified initial conditions isThe generic late response is governed by $r_{+}$, the root nearer zero. For very large damping ratio, $r_{+}\sim -{\omega }_0/(2\zeta )=-k/b$: an increasingly strong brake delays the return to equilibrium. This approximation comes from balancing the spring and viscous damping terms after the fast transient. Choosing $v_0=r_{-}{x}_0$ makes $C_{+}=0$ and leaves only the faster root; that exception requires different initial velocities as $b$ varies.

A partial differential equation relates an unknown function of several independent variables to its partial derivatives. For $u(x,t)$, the notation $u_t=\partial u/\partial t$ means changing time while holding position fixed; $u_{xx}={\partial }^{2}u/\partial x^2$ measures how its spatial slope changes. Unlike an ordinary differential equation, a partial differential equation can describe an entire evolving profile rather than a single trajectory.

The heat equation $u_t=D{u}_{xx}$ smooths spatial variations, the transport equation $u_t+c{u}_x=0$ moves them, and Laplace's equation $\Delta u=0$ describes a spatial equilibrium. Here $D>0$ is a diffusivity, $c$ is a velocity, and $\Delta$ is the Laplacian. The periodic heat and transport comparison makes the difference visible using identical initial conditions.

A partial differential equation usually supplies a local rule rather than a complete answer. For example, every $u(x,t)=A+Bx$ satisfies the heat equation with zero PDE source term; data select which profile is relevant. A complete problem specifies its domain of a function, coefficients, initial conditions, boundary conditions, and meaning of solution of a partial differential equation. A well-posed problem has existence, uniqueness, and continuous dependence on initial data in stated spaces; the equation alone guarantees none of these.^[ref]

Eliminate the juvenile variable from the stage-structured population model. Since $y_{t+1}=2z_t$,The characteristic polynomial of this linear recurrence isso $z_t=A(3/2{)}^t+B(-1{)}^t$. The initial conditions give $z_0=0$ and $z_1=(3/4)10=15/2$, hence $A+B=0$ and $(3/2)A-B=15/2$. ThereforeThe alternating term records the initial imbalance between life stages; the growing $(3/2{)}^t$ term dominates eventually.

The homogeneous characteristic equation is $r^2-r-2=(r-2)(r+1)=0$, so its solutions are $e^{2x}$ and $e^{-x}$. Because the forcing contains the latter exponential, set $y=e^{-x}v$. Direct differentiation reduces the inhomogeneous linear differential equation toA polynomial particular solution is $v=-3x^2-2x$, since $v^{\prime \prime }-3v^{\prime }=-6-3(-6x-2)=18x$. ThereforeThe exponentially growing first term cannot be canceled by the decaying second term, so boundedness at positive infinity requires $C_1=0$. The initial condition $y(0)=1$ sets $C_2=1$. Hence

The mean value theorem states that, if $f$ is continuous on $[a,b]$ and differentiable on $(a,b)$, with $a<b$, there is $c\in (a,b)$ such thatTo deduce it from Rolle's theorem, subtract the secant line: $k(t)=f(t)-f(a)-\frac{f(b)-f(a)}{b-a}(t-a)$. This function has the same continuity and differentiability properties and satisfies $k(a)=k(b)=0$. Rolle's theorem gives an interior point with $k^{\prime }(c)=0$, which is precisely the required equality.

If $h^{\prime }=0$ on $\mathbb{R}$, apply the mean value theorem on any $[x,y]$ with $x<y$. It gives $h(y)-h(x)=(y-x)h^{\prime }(c)=0$. Hence $h$ is constant.

For the differential equation, the product rule givesThe preceding result makes this product constant. Thus all solutions are $f(x)=C{e}^{ax}$, and differentiation verifies every such solution.

For the integral assertion, continuity ensures that $f$ is Riemann integrable on every finite closed interval. Use the usual orientation for an integral with negative upper endpoint. For $h\ne 0$, additivity of the integral givesThe absolute value of this expression is at most the supremum of $|f(t)-f(x)|$ on the segment joining $x$ and $x+h$. Continuity at $x$ makes this bound tend to zero for either sign of $h$. Consequently $F^{\prime }(x)=f(x)$, proving this form of the fundamental theorem of calculus.

Finally, a differentiable solution of the integral equation is continuous, so the fundamental theorem of calculus yields $g^{\prime }=g$; setting $x=0$ gives $g(0)=A$. The differential-equation calculation above now forcesConversely, $A+{\int }_0^{x}A{e}^{t}dt=A+A(e^x-1)=A{e}^x$, so this function satisfies the original integral equation. It is unique because every differentiable solution has already been forced to this form and this initial condition.

The characteristic equation of the homogeneous linear differential equation is $m^2+m-2=(m-1)(m+2)=0$, giving the homogeneous solution $A{e}^t+B{e}^{-2t}$. For the first forcing, a particular solution of the form $C{e}^{-t}$ gives $-2C=1$, soThe initial conditions impose $A+B=1/2$ and $A-2B=-1/2$, whence $A=1/6$, $B=1/3$. Thus the first solution isFor the second forcing, $e^t$ is already a homogeneous solution, so the usual exponential trial must be multiplied by $t$. Substitution of the particular solution $Ct{e}^t$ gives $3C{e}^t=e^t$, hence $C=1/3$. Now $y=A{e}^t+B{e}^{-2t}+t{e}^t/3$, and the initial conditions give $A+B=0$, $A-2B=-1/3$. Therefore the resonant solution isThe additional factor $t$ reflects resonance in a differential equation.

This is a separable differential equation. On the solution with $y<1$, integration from the initial condition givesThe right side lies in $[0,1)$, so taking its sine stays on the required branch. The solution isNear the origin, $1-\sqrt{1-x^2}=x^2/2+O(x^4)$, so $y=x^2/2+O(x^4)$ and the tangent line is horizontal. At the other end, $y\to \sin 1$, whileThus the curve has a vertical limiting tangent line as $x\uparrow 1$, even though $x=1$ is excluded from the domain.

For the direction field, the slope is nonnegative everywhere in the square. The segments are horizontal on $x=0$ and on $y=1$, become steeper as $x$ increases at fixed $y<1$, and become flatter as $y$ increases at fixed $x$. Near $x=1$ they are almost vertical unless $y$ is correspondingly close to one. The excluded corner $(1,1)$ has no single limiting slope: the ratio $(1-y^2)/(1-x^2)$ depends on the path of approach.

An integrating factor for a differential one-form is a nonzero function $\mu (x,y)$ for which $\mu (dy+fdx)$ is an exact differential, say $dF$. Then solutions lie on level sets $F=\text{constant}$. Locally, the compatibility condition is ${\partial }_x\mu ={\partial }_y(\mu f)$.

Here multiplication by $2y{e}^x$, which is nonzero for $y>0$, producesThe initial condition makes this first integral equal to $a^2$. Taking the positive branch yields the solutionon the connected interval containing zero where the radicand is positive. This domain restriction matters because the original differential equation is singular at $y=0$.

Completing the square gives $2x+x^2=(x+1{)}^2-1\ge -1$. As $y>0$, the differential equation consequently impliesAt $x=0$, $y^{\prime }=-a/2<0$. Thus, moving a little to the left from $y(0)=a\ge 1$, the curve rises above $a$. Throughout $y>1$ every segment of the direction field points downwards as $x$ increases, by the inequality above. Tracing backwards therefore makes $y$ rise further, so it cannot return to $y=1$. This graphical barrier argument proves $y>a\ge 1$ and $y^{\prime }<0$ for every $x<0$ on the solution. The explicit first integral also shows this branch extends throughout $x<0$: writing $t=-x>0$, the minimum of $e^t/t^2$ is $e^2/4>1$, so $a^2{e}^t-t^2>0$ for $a\ge 1$. Including the initial point, $y^{\prime }<0$ for all $x\le 0$.

For $a=1$ the curve is $y=\sqrt{e^{-x}-x^2}$ on $(-\infty ,b)$, where the strictly increasing function $x^2{e}^x$ has the unique positive root $x^2{e}^x=1$ at $b\simeq 0.703467$. The curve decreases through $(0,1)$ with slope $-1/2$ and reaches $(b,0)$ only as a limiting endpoint. Its derivative isAs $x\to -\infty$, $y\sim e^{-x/2}$ and $y^{\prime }\sim -\frac{1}{2}{e}^{-x/2}\to -\infty$. As $x\uparrow b$, the numerator tends to $-b^2-2b<0$ and the denominator to zero from above, so again $y^{\prime }\to -\infty$. More precisely, $y\sim \sqrt{b(b+2)(b-x)}$ at the endpoint. Both ends have the requested unbounded negative slope.

For the logistic differential equation, separate variables to obtain $\log |y/(1-ay)|=rt+C$ away from the constant solutions. Applying the positive initial condition givesThe formula also includes the constant solution $y_0=1/a$. Below $1/a$ positive solutions increase; above $1/a$ they decrease. They approach $1/a$ as $t\to \infty$. Equivalently, linearizing $F(y)=ry(1-ay)$ gives $F^{\prime }(0)=r>0$ and $F^{\prime }(1/a)=-r<0$. Thus zero is unstable and $1/a$ is asymptotically stable.

The forward Euler method gives $y_{n+1}=(1+r\delta t)y_n-ra\delta t{y}_n^2$. Set $\lambda =1+r\delta t$ and $u_n=a(\lambda -1)y_n/\lambda$; multiplying the recurrence by this scaling factor gives the logistic mapIts fixed points are $u=0$ and $u_{\ast }=1-1/\lambda$. The fixed-point multipliers are $F^{\prime }(0)=\lambda$ and $F^{\prime }(u_{\ast })=2-\lambda$. Since $\lambda >1$, zero is always unstable. The positive fixed point is stable for $1<\lambda <3$ and unstable for $\lambda >3$. For $\lambda >3$, a small displacement obeys ${ϵ}_{n+1}\simeq (2-\lambda ){ϵ}_n$: its sign alternates and its amplitude grows, giving an oscillatory instability whose period is $2\delta t$.

This is an artifact of Euler discretization of logistic growth. Exact small perturbations decay by the positive factor $e^{-r\delta t}$ per time step; a one-dimensional autonomous differential equation cannot have a nonconstant periodic orbit, since its solutions move monotonically between equilibrium points. The artificial instability occurs when $r\delta t>2$, beyond the stable step range of this discretization.

For $\lambda =3.2$, remove the fixed point factors from $F(F(u))-u$ to obtain ${\lambda }^2{u}^2-\lambda (\lambda +1)u+\lambda +1=0$. The resulting logistic map two-cycle isThese points map to each other. Their periodic-orbit multiplier is $F^{\prime }(u_{-})F^{\prime }(u_{+})=4+2\lambda -{\lambda }^2=0.16$, so the two-cycle attracts nearby iterates. The cobweb shows a small initial departure from $u_{\ast }=0.6875$ growing with alternating sign and then settling into the finite rectangle between $u_{-}$ and $u_{+}$. The time trace shows the corresponding finite-amplitude oscillation.

Write $T=\epsilon t$ and $q(T)=2e^{-T}$. The forcing acts on the fast initial-time scale, so first remove a localized particular solution. Expanding the coefficient there givesA particular choice is $p_1=(t+1/2)e^{-2t}$, henceIts initial value is $1+\epsilon /2$ and its initial derivative is $-2+O({\epsilon }^2)$. The localized forcing determines WKB phase constants; omitting it before imposing the initial conditions would miss the permanent free oscillation.

For the homogeneous part, the slow-time equation is ${\epsilon }^2{u}_{TT}+q^{2}u=0$. Write its WKB approximation as $A(T)[C\cos \Phi +B\sin \Phi ]$, normalizing $A(0)=1$. The leading amplitude transport equation gives $A=e^{T/2}$ and leading phase ${\epsilon }^{-1}{\int }_0^{T}q(s)ds=2(1-e^{-T})/\epsilon$.

To include every $O(\epsilon )$ contribution on bounded slow-time intervals, the first phase correction is needed. If the local phase wavenumber is $q+{\epsilon }^2{q}_2$, the real part of the amplitude-phase equation gives $q_2=A_{TT}/(2qA)=e^T/16$. The WKB phase correction for a slowly varying oscillator, anchored to zero at the initial time, is thereforeThere is no first-order bulk amplitude correction in this normalization. At $t=0$, the homogeneous initial value must be $-\epsilon /2+O({\epsilon }^2)$ and its initial derivative $1+O({\epsilon }^2)$. Since $A_t^{\prime }(0)=\epsilon /2$ and ${\Phi }_t^{\prime }(0)=2+O({\epsilon }^2)$, these give $C=-\epsilon /2+O({\epsilon }^2)$ and $B=1/2+O({\epsilon }^2)$. The requested forced approximation isfor $0\le \epsilon t\le T_0$ with fixed $T_0$. The approximate initial value is exactly one through this order and its derivative is $-1+O({\epsilon }^2)$. The particular-solution residual is uniformly localized and $O({\epsilon }^2)$; the corrected homogeneous residual is $O({\epsilon }^3)$ on the long interval, producing $O({\epsilon }^2)$ solution error after integration over times of order ${\epsilon }^{-1}$.

The WKB approximation requires the relative frequency-change parameterIt therefore fails when $\epsilon e^{\epsilon t}=O(1)$, even though the coefficient has no finite-time zero. This is the Bessel transition for an exponentially decaying oscillator. Introduce the shifted slow time and amplitudeThe exact rescaled equation isFor fixed $\tau$ its right side is smaller than every algebraic power of $\epsilon$. The leading transition problem is consequentlyWith $Z=2e^{-\tau }$ this becomes $Z^2{U}_{0,ZZ}+Z{U}_{0,Z}+Z^2{U}_0=0$, the order-zero Bessel differential equation. This identifies the governing equation without solving it.

The amplitude-phase series also supplies overlap data while $\epsilon e^T\ll 1$, although its bounded-$T$ absolute error estimate is not uniform up to the transition. For matching, combine the earlier sine and cosine terms into $\frac{1}{2}\sin (\Phi -\epsilon )$ through first order. In the overlap $\epsilon \ll e^{\tau }\ll 1$, their phase isto the retained order. Thus amplitude and phase matching as $\tau \to -\infty$ requireThese are oscillatory amplitude-phase conditions, rather than pointwise ratio statements at zeros. The smaller $e^{\tau }/16$ phase term records the first overlap correction if required. The rapidly varying matching constant ${\alpha }_{\epsilon }$ must be retained as data; there is no single epsilon-independent limiting phase. The two matching conditions determine the leading transition solution, which is deliberately left unsolved.