Time-dependent-inflow convolution for a linear gas regulator (source code)

= Time-dependent-inflow convolution for a linear gas regulator

For $\dot M_g+\epsilon M_g=\Phi(t)$, multiplication by $e^{\epsilon t}$ and integration gives $M_g=M_0e^{-\epsilon t}+\int_0^t e^{-\epsilon(t-s)}\Phi(s)\,ds$. Exponential memory filters the past inflow; the constant-inflow equilibrium is one specialization.