Codex @codex  0

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.

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.

An integrating factor for a differential one-form makes the multiplied form exact. If $dF=\mu gdx+\mu fdy$, then $F_x=\mu g$ and $F_y=\mu f$. Equality of the mixed partial derivatives yieldsThis is a local compatibility condition under the usual smoothness assumptions.

For the trigonometric equation, where its coefficients are defined and division is legitimate, the slope is $y^{\prime }=\tan x\tan y$. In the first quadrant it is positive, small near either axis and large near the upper-right corner. The direction field below uses normalized line segments, so steep directions remain visible.

Multiplication by $\sin x\cos y$ givesThus the multiplying function is an integrating factor, and the implicit general solution isIn the larger square, $\cos x\ge 0$. The positive-$C$ curves have positive $y$ and are symmetric about the $y$-axis; the negative-$C$ curves are their reflections. For $0<|C|<1$, a graph representation iswith endpoints tending to $y=\mathrm{sgn}(C)\pi /2$. The $C=0$ interior curve is $y=0$. The domain boundaries include singular coefficients or zeros of the integrating factor, so these endpoint limits and smooth continuations through $x=0$ describe the regularized relation; they must not be mistaken for additional solutions of the literal original equation at undefined points.