For continuous a,fa,f, the y′+a(t)y=f(t)y'+a(t)y=f(t) has integrating factor

μ(t)=exp⁡(∫t0ta(s) ds).\mu(t)=\exp\left(\int_{t_0}^ta(s)\,ds\right).

The product rule gives (μy)′=μf(\mu y)'=\mu f, so

y(t)=μ(t)−1(y(t0)+∫t0tμ(s)f(s) ds).y(t)=\mu(t)^{-1}\left(y(t_0)+\int_{t_0}^t\mu(s)f(s)\,ds\right).
Regularity and matrices

For integrable a,fa,f, the same formula gives an absolutely continuous solution and the equation holds almost everywhere. For matrix coefficients at different times, multiplication need not commute; replacing a scalar integral exponential by a matrix exponential requires a commutation condition. The general matrix construction uses a .