As it stands the left side y′ + p·y is not the derivative of anything you can write down. y′ alone is the derivative of y, and the extra term p·y spoils it. A product μ·y has the derivative μ·y′ + μ′·y, which has the same two-term shape; the only question is whether a μ exists that makes the second term p·μ·y.
Theorem · The integrating factorIf μ′ = p·μ, then (μ·y)′ = μ·(y′ + p·y) for every differentiable y. So y′ + p·y = q implies (μ·y)′ = μ·q.
The product rule, verified exactly on polynomials in Rates of Change and the Derivative, gives (μ·y)′ = μ·y′ + μ′·y.
Replace μ′ by p·μ: μ·y′ + p·μ·y = μ·(y′ + p·y). If y′ + p·y = q the right side is μ·q. Only the product rule and the one condition on μ were used.
□
Two families of p need nothing else. For a constant p = a, μ = e^(a·t), since the rate of e^(a·t) is a times itself, a claim from the earlier courses. For p = a/t with a whole number a, μ = tᵃ, and the check is the power rule alone. It is exact, and the lab prints the result without a logarithm.
p = 2/t, μ = t²
μ′ = 2t
p·μ = (2/t)·t² = 2t
equal, so μ′ = p·μ
Why the exponent is the antiderivative of p
For the general formula μ = e^(∫p dt) the chain rule, verified on polynomials in Rates of Change and the Derivative and applied here to the exponential as a stated claim, gives μ′ = (∫p dt)′·e^(∫p dt) = p·μ, which is the condition. The exponent has to be an antiderivative of p. Using p itself fails at once: for the constant p = 3, e^p = e³ is a constant, its rate is 0, and p·μ = 3·e³ is not.
Which antiderivative does not matter. A different constant added to ∫p dt multiplies μ by a constant factor, such as 5 or e⁷, and that factor multiplies both sides of the equation and cancels. μ = 5t² works as well as μ = t², and the final y is the same.
Example · A case with an exponential factor and a forcingTake y′ + 2y = e^t. Then μ = e^(2t), and (e^(2t)·y)′ = e^(2t)·e^t = e^(3t).
Integrate: e^(2t)·y = (1/3)·e^(3t) + C, so y = (1/3)·e^t + C·e^(−2t). The lab prints exactly this, with C left free because no starting value was given.
For p = a/t the same steps stay in fractions. With y′ + (2/t)·y = t² and y(1) = 1, the factor is t², (t²·y)′ = t⁴, t²·y = t⁵/5 + C, and y = t³/5 + C/t². Substituting to check: y′ = 3t²/5 − (2C)/t³ and (2/t)·y = 2t²/5 + (2C)/t³, which add to t² for every C.
y′ + (2/t)·y = t², y(1) = 1
μ = t²
(t²·y)′ = t⁴
t²·y = t⁵/5 + C
y = t³/5 + C/t²
y(1) = 1/5 + C = 1, so C = 4/5
The lab limits itself to the families above, and says so when it refuses. It takes a constant p, or a/t with a a whole number from −4 to 4; it will not integrate tᵃ·q when that integral has a ln t in it; and for p = a/t it asks for a starting time other than 0, since the equation is not defined at t = 0. These are limits of the lab, and the integrating factor itself has no such limit.