I have searched the site for posts regarding Darboux integrability $\implies$ Riemann integrability, but haven't found any that specifically adress this question.
My definition of Darboux integrability: Let $f$ be defined and bounded on $[a,b]$, then $f$ is Darboux integrable if for all $\epsilon >0$ there exists a partition $P$ of $[a,b]$ such that $U(f,P)-L(f,P)<\epsilon$ (where $U$ and $L$ are the upper and lower Riemann sums respectively).
My definition of Riemann integrability: Let $f$ be defined and bounded on $[a,b]$, then $f$ is Riemann integrable if $$\lim_{N\to\infty} \sum\limits_{k=1}^{N} f(c_k)(x_{k}-x_{k-1})$$ has the same limit for all sequences of partitions $P_N$ and all choices of $c_k\in[x_{k-1},x_{k}]$.
If my definitions are correct, it seems that Darboux integrability only requires one partition to fulfil the epsilon-inequality, whereas Riemann integrability requires all sequences of partitions to be fulfilled. How can this lead to an implication nevertheless?
If $U(f,P)-L(f,P)<\epsilon$, for a given partition $P$, then the inequality holds for $any$ refinement of $P$, since if $P\subseteq P'$, we have
$L(f,P)\le L(f,P')\le U(f,P')\le U(f,P).$
And since for any partition $P=\{a,x_1,\cdots, x_{n-2},b\}$,
and any $\{x^*_i:x_i\le x^*_i\le x_{i+1};\ 1\le i\le n\},$
we have $L(f,P)\le \sum^{n-1}_{i=1}f(x^*_i)(x_{i+1}-x_i)\le U(f,P),$
it follows that Darboux integrability implies Riemann integrability.
edit: Assuming that the OP wants to show that if bounded function $f$ on $[a, b]$ is Darboux integrable then for each given $\epsilon> 0$, there is a $\delta > 0$ such that
mesh$P < \delta \Rightarrow U(f, P) − L(f, P) < \epsilon,$ here is a sketch:
Suppose we have a partition $P_0=\{a = x_0 < \cdots < x_m = b\}$ such that $U(f,P_0)-L(f,P_0)<\epsilon.$ Let $P=\{a = t_0 < t_1 < \cdots < t_n = b\}$ be any other partition, fine enough so that at most one member of $P_0$ lies between any two members of $P$.
Now, form $Q=P\cup P_0$ and note that in the difference $L(f,Q)-L(f,P)$, the only terms that remain correspond to elements of $Q$ of the form $[t_{n_k}\le x_{n_k}\le t_{n_{k+1}}]$.
Then, with $m_k,m'_k,m''_k$ the mimima of $f$ on $[t_{n_k},x_{n_k}], [x_{n_k}, t_{n_{k+1}}]$ and $[t_{n_k}, t_{n_{k+1}}],\ $ respectively and $M$ an upper bound on $|f|$,
$L(f, Q) − L(f, P)=\sum _k[(m_k(x_{n_k}-t_{n_k})+m'_k(t_{n_{k+1}}-x_{n_k}))-m''_k(t_{n_{k+1}}-t_{n_k})]\le \sum_k[M(t_{n_{k+1}}-t_{n_k})-m_k''((t_{n_{k+1}}-t_{n_k}))]\le MK(\text{mesh P})$ where $K$ is the upper limit of the sum.
Similarly, $U(f,P)-U(f,Q)\le MK(\text{mesh P})$.
The result now follows from the string of inequalities:
$U(f, P) − L(f, P) ≤ U(f, P) − U(f, Q) + U(f, Q) − L(f, Q) + L(f, Q) − L(f, P)$.