This popped up in my research: we have the following mix-effect inequality that $\forall T \geq 1$ \begin{equation}\tag{*} Y(T) - \frac{1}{100T^2}\int_1^T[\alpha^2 + e^{-(T - t)}]Y(t)dt \lesssim \frac{1}{T^2}[\alpha^2 + e^{-T}] Y(1) \end{equation} where $\alpha\in [0,1]$ is a small constant. We denote by $a\lesssim b$ when $a\le Cb$ for some absolute constant $C>0$.
Note we can have a coarsened bound ($\alpha^2 + e^{-T} \lesssim 1$)) which gives $$ Y(T) - \frac{1}{100T^2}\int_1^T Y(t) dt \lesssim \frac{1}{T^2}Y(1) $$ Then the classical Gronwall-type inequality would give (or simply integrating on both sides) $$ Y(T) \lesssim \frac{1}{T^2}Y(1) $$ But this bound is not tight in the dependency on $\alpha$: when $\alpha$ is set as 0 $$ e^T Y(T) - \frac{1}{100T^2} \int_1^T e^t Y(t) dt \lesssim \frac{1}{T^2} Y(1) $$ Integrating on both sides we get $$ e^T Y(T) \lesssim \frac{1}{T^2}Y(1) $$ so there is an exponential decaying prefactor.
Now back to (*), we would conjecture the following fine-grained bound $$ Y(T) \lesssim \frac{\alpha^2 + e^{-T}}{T^2}Y(1) $$ But how do we achieve this (or something similar)?