I'm trying to get an order bound on the following integral as $\epsilon \to 0$:
$$g(\epsilon) = \int_a^{a+\epsilon} f(\epsilon,v) dv$$ where $f$ is an ugly, but smooth, function when $a\leq |v| <a+\epsilon$. My issue is that when I attempt to expand $$g(\epsilon) = \epsilon f(0,a) + \frac{\epsilon^2}{2}\big(2 f_\epsilon(0,a) + f_v(0,a)\big) + o(\epsilon^2),$$
I run into issues. Namely, some terms go to infinity when evaluated at $\epsilon = 0$. I also notice that, if I look (formally) at the same expansion without evaluating anything in $\epsilon$, the series order decreases with $\epsilon$ as I add more derivatives: first $\epsilon^{-2}$, then $\epsilon^{-4}$, and so on.
Would someone please explain to me what is going on here? I have never done series calculations like this before, so a reference containing theory and/or examples would also be helpful. Thanks,