I'm reading Lanczos Variational Principles of Mechanics p.124, and following a discussion of how for scleronomic systems we get
$$\sum_{i=1}^{n} p_i\dot q_i - L = const.\tag{53.12}$$
For rheonomic systems it's stated that
$$\delta L=dL-\frac{\partial L}{\partial t}dt = \epsilon\left(\dot L -\frac{\partial L}{\partial t}\right)\tag{53.22}$$ where $\epsilon=dt$, which leads to $$\left[\sum_{i=1}^{n} p_i\dot q_i - L\right]^{t_2}_{t_1} = -\int^{t_2}_{t_1} \frac{\partial L}{\partial t} dt\tag{53.23}$$
However, when I do the variation
$$\delta\int_{t_1}^{t_2} L~dt= \epsilon\int_{t_1}^{t_2} \left(\dot L -\frac{\partial L}{\partial t}\right)dt = \epsilon L|_{t_1}^{t_2} - \int_{t_1}^{t_2}\frac{\partial L}{\partial t}dt$$
I'm getting an extra $\epsilon L|_{t_1}^{t_2}$ term? Any insight on what missing would be greatly appreciated!