I am trying to prove the existence of the Riemann integral for continuous functions defined over the closed interval $[a,b].$ I am doing this with Darboux upper and lower integrals. To show that a function $f:[a,b]\to \mathbb{R}$ has a Riemann integral $I=\int_{a}^{b} f(x) dx,$ we must show for every $\varepsilon > 0$ there exists a $\delta >0$ such that any partition $P$ that has partition norm $||\Delta x|| < \delta$ will satisfy $|\sum_{i=1}^{n} f(c_i) \Delta x_i - I| < \varepsilon ,$ where $c_i$ is any element in $[x_{i-1},x_{i}]$ and $n$ is the number of intervals resulting from the partition.
We have a lemma establishing that for any partitions $P$ and $Q$ of $[a,b]$ with $Q \subset P,$ $L(f,Q) \leq L(f,P) \leq I_L \leq I_U \leq U(f,P) \leq U(f,Q)$ where $U$ and $L$ denote the upper and lower sums and $I_L, I_U$ denote the upper and lower integrals, with $I_U = \inf\limits_{P} U(f,P)$ where $P$ runs over all partitions of $[a,b].$
Assuming that $f$ is continuous, I have shown that $I_U$ exists and this is my candidate for what the Riemann integral $I$ of $f$ should be. I know that $f$ is uniformly continuous over its compact domain so given $\varepsilon >0$ I can choose $\delta > 0$ so that $|x-y| < \delta$ implies $|f(x) - f(y)| <\frac{\varepsilon}{3(b-a)}$. Let $P$ be any partition of $[a,b]$ with norm $||\Delta x|| < \delta$. Then by the triangle inequality \begin{align*} | \sum_{i=1}^n f(c_i) \Delta x_i - I_U| & \leq | \sum_{i=1}^n f(c_i) \Delta x_i - U(f,P) | + | U(f,P) - I_U | &\\= |\sum_{i=1}^n (f(c_i) - M_i) \Delta x_i| + | U(f,P) - I_U |, \end{align*} where $M_i$ is the maximum value that $f$ attains over $[x_{i-1},x_i].$ By our choice of $P$ and $\delta$ we are guaranteed that $|f(c_i) - M_i| < \frac{\varepsilon}{3(b-a)}$ so the sum involving this term is less than $\varepsilon /3.$ For the remaining term we have by the infimal nature of $I_U$ that there exists a partition $P_2$ such that $|I_U - U(f,P_2)| < \varepsilon/3,$ and this must be true of $Q= P\cup P_2$ as well by the inequalities in our lemma. Hence we can apply the triangle inequality again to conclude that $$ | U(f,P) - I_U | \leq |U(f,Q) - I_U| + |U(f,P) - U(f,Q)|, $$ where we have said already that $|U(f,Q) - I_U| <\varepsilon/3.$ Now here is where I am stuck. I want to be able to conclude that the remaining term, $|U(f,P) - U(f,Q)|,$ is or can be made arbitrarily small. My issue is that these upper sums are over different partitions so they have differing numbers of terms and differing lengths of each subinterval.
My question is does what I have so far look correct and is it true that this final term I have can be made arbitrarily small?