It is well-known that for bounded functions of one variable, Darboux integrability (via upper and lower sums) is equivalent to Riemann integrability (via Riemann sums.) I tried to generalize this to functions of two variables, and I encountered some technical issues; consider, for example, a bounded function in two variables defined on a unit disk; the lower sum is built from the little squares inside the disk and infimum of function values, while the upper sum is built from the little squares intersecting the disk and supremum of function values; (Notice that there are more little squares for the upper sum than for the lower sum.) then one defines the Darboux integrability as the lower sum and upper sum being equal in some limit as the little squares become smaller; to make this sandwiching to work via monotone convergence theorem, one needs to assume that the function is nonnegative. A general bounded function is Darboux integrable if it is the difference between two nonnegative Darboux integrable functions.
As for the Riemann integrability, the Riemann sum is built from the little squares inside the disk and function values at an arbitrary point in each little square. It is straightforward to show that Darboux integrability implies Riemann integrability. What I'm having trouble with is the converse: Does Riemann integrability imply Darboux integrability? I can prove this only if the function is nonnegative; since upper and lower sums are used for nonnegative functions only, the usual technique cannot be applied for signed functions. Is there a way around it? Or, could it be that Riemann integrability does not imply Darboux integrability for functions of more than one variable?