I am trying to prove the following theorem below. I have also found the following theorem here (see Theorem 2.1 that clearly covers this but no proof is given). I don't want to make these evenly spaced intervals for this theorem as they could be any length. I am wanting to prove this formally from start to finish; any help would be greatly appreciated!
Definition: Let $f(x)$ be a continuous function on $[a, b]$. Then, we define $\int_a^b f(x):=\lim_{n\to \infty}\sum_{i=0}^nf(x_i^*)\,\frac{b-a}{n}$.
Theorem: Let $f(x)$ be a continuous function on $[a, b]$. Then, $\int_a^b f(x)=\lim_{\text{mesh} \to 0}\sum_{i=0}^nf(x_i^*)\,\Delta x_i$ as shown in the definition here.
Proof:
\begin{align*} \lim_{\text{mesh} \to 0}\sum_{i=0}^nf(x_i^*)\,\Delta x_i&= \lim_{\text{mesh} \to 0}\sum_{i=0}^nf(x_i^*)\,\Delta x_i\\ &= \lim_{\max\{\Delta x_1, \Delta x_2, \ldots, \Delta x_n\} \to 0}\sum_{i=0}^nf(x_i^*)\,\Delta x_i \text{(CAN'T DO THIS)} \end{align*}
I have made some progress on the following. The limit in the second equation is NOT an Epsilon-Delta. Read the article here and see minimalrho's answer (note order is refinement here). This order makes is needed in regards to partitions to make it be well-defined. Now, look at the second definition here for Riemann Integral. It is equivalent to the first definition of a Riemann integral in that page which in my honest opinion seems easier to use. Now, the first definition kind of looks like the Epsilon-Delta but in NO WAY SHAPE OR FORM is it related. Now, we can use the first definition of Wikipedia in that article. Note the definition where $n\to \infty$ is $F(b)-F(a)$ by the FTOC but that requires the antiderivative to exist.
Since $f$ is uniformly continuous on $[a,b]$, for every $\varepsilon>0$ there exists $\delta>0$ such that $$|f(x)-f(y)|<\frac \varepsilon{4(b-a)}\qquad\text{if}\qquad|x-y|<\delta$$ Then, if $P_1$ and $P_2$ are tagged partitions of $[a,b]$ with $\|P_1\|<\delta$ and $\|P_2\|<\delta$, the corresponding Riemann sums satisfy $$|R(f,P_1)-R(f,P_2)|< \frac\varepsilon{2}$$See $\;\;$ R. Courant, F. John - Introduction to Calculus and Analysis vol. I - 1999, pp. 193-194 .
It follows that, for any sequence $\{Q_n\}$ of tagged partitions whose meshes converge to zero, the sequence $\{R(f,Q_n)\}$ of the corresponding Riemann sums is Cauchy.
If $A$ is its limit, then there exists $\bar N\in\mathbb N$ such that $\quad|R(f,Q_n)-A|<\dfrac\varepsilon{2} \quad$if$\quad n>\bar N$ .
Let also $N\in\mathbb N \;$ be such that $\quad \|Q_n\|<\delta \quad$if$\quad n>N$ .
Consider now any tagged partition $P$ with $\|P\|<\delta$ .
Then, if $n>\max \{N,\bar N\}$, you have $$|R(f,P)-A|\le|R(f,P)-R(f,Q_n)|+|R(f,Q_n)-A|<\varepsilon$$