I am trying to upper bound the following function and find its growth rate:
\begin{equation} \psi(y) \stackrel{\triangle}{=} \int_{0}^{\infty}\exp\left(-\frac{(y-x)^2}{2(\sigma_0^2 + \sigma_1^2 x)}\right)f(x)\,dx,~y>0, \end{equation} where $f(x)>0$ satisfies $\int_{0}^{+\infty}f(x)\,dx = 1$ (it is a pdf of the random variable $X$). I would like to find the growth rate of $\psi(y)$ for large values of $y$, i.e., would like to know if $\psi(y)$ grows like $O(y^{-k})$, for some $k>1$.
Here are my attempts so far: It is straightforward to show that the function $\exp\left(-\frac{(y-x)^2}{2(\sigma_0^2 + \sigma_1^2 x)}\right)\in[0,1]$ is increasing in $x$ for $y>x>0$ and is decreasing in $x$ for $0<y<x$.
Then one can break the bounds of the integral to $[0,y/2], \, [y/2,3y/2], \, [3y/2, y^{3/2}],\,[y^{3/2},\infty)$ and then evaluate each integral. Doing so will give the following:
\begin{align} \int_{0}^{y/2}\exp\left(-\frac{(y-x)^2}{2(\sigma_0^2 + \sigma_1^2 x)}\right)f(x)\,dx &= O(y^{-k}), \,k>1\\ \int_{y/2}^{3y/2}\exp\left(-\frac{(y-x)^2}{2(\sigma_0^2 + \sigma_1^2 x)}\right)f(x)\,dx &\stackrel{?}{=}O(y^{-k}),\, k>1\\ \int_{3y/2}^{y^{3/2}}\exp\left(-\frac{(y-x)^2}{2(\sigma_0^2 + \sigma_1^2 x)}\right)f(x)\,dx &= O(y^{-k}), \,k>1\\ \int_{y^{3/2}}^{\infty}\exp\left(-\frac{(y-x)^2}{2(\sigma_0^2 + \sigma_1^2 x)}\right)f(x)\,dx &= O(y^{-k}), \,k>1 \end{align}
Therefore, if we can prove that the second integral also grows like $O(y^{-k}),\,k>1$, then we are done. However, I have not been able to show this and am stuck. Any help would be highly appreciated!