To prove: Let $f,g:[a,b] \to \mathbb{R}$ be two functions such that $f$ is Riemann integrable and $g$ has an antiderivative $G$ on $[a,b]$. If $$\int_a^b f(x)dx=G(b)-G(a) \text { and } f \le g,$$ then $g$ is also Riemann integrable.
I thought about using Darboux's criterion. I have some trouble proving that $g$ is bounded, but nevertheless, if we consider a sequence of partitions $\Delta_n=(x_0=a,x_1,...,x_n=b)$ such that $||\Delta_n ||\to 0$, we have that $$G(b)-G(a)=\sum_{k=1}^n\left[G(x_k)-G(x_{k-1})\right]\stackrel{\text{MVT}}{=}\sum_{k=1}^n (x_k-x_{k-1})g(c_k)$$ and then I tried to somehow bound the upper and lower Darboux sums, but I wasn't succesful.
If we were lucky enough to have $f$ continuous, we could do this: Define
$$h(x) = G(x)-G(a) -\int_a^x f(t)\,dt.$$
Now $G'=g$ is given. And by the FTC, the derivative of the integral function is $f(x).$ Thus
$$h'(x) =g(x)-f(x)\ge 0.$$
Thus $h$ is nondecreasing on $[a,b].$ But $h(a)=h(b)=0.$ A nondecreasing function that is $0$ at the endpoints must be identically $0.$ Therefore $g\equiv f$ and $g$ is Riemann integrable.