Suppose $f(x,y) : [0,1] \times [0,1] \to \mathbb{R}_{\ge 0}$ is a function, not necessarily continuous, nor necessarily Riemann integrable on $[0,1] \times[0,1]$. By existence of an integral here we will also admit also the case in which it is unbounded.
Moreover assume that:
- for each $x \in [0,1]$ the Riemann integral $\int_0^1 f(x,y) dy$ exists;
- for each $y \in [0,1]$ the Riemann integral $\int_0^1 f(x,y) dx$ exists.
Does this implies that the Riemann integral of $F(x) = \int_0^1 f(x,y) dy$ and $G(y) = \int_0^1 f(x,y) dx$ exist?
The latter are nothing else than the iterated integrals of $f(x,y)$, i.e. $\int_0^1 F(x)dx = \int_0^1 \int_0^1 f(x,y) dy dx$, $\int_0^1 G(y)dy = \int_0^1 \int_0^1 f(x,y) dxdy$.
Observe that I cannot apply Fubini's theorem since I do not know whether $f(x,y)$ is integrable in $[0,1] \times [0,1]$. As far as I know Tonelli's theorem only apply to Lebesgue integrals.
Attempt 1: proving that the answer is yes
I have been trying to prove this by contraddiction. First suppose that $f$ is bounded and that, for example, $F(x)$ is non-Riemann integrable. Then its set of discontinuity points $X^*$ must have non-zero Lebesgue measure. I would like to show that this implies that at least for one $y^*$ $f(x, y^*)$ must have a non-zero measure set of discontinuity points, therefore $\int_0^1 f(x,y^*) dx$ does not exist. However, my attempts have failed so far because I do not know how to link the discontinuities of $F(x)$ and those of $f(x,y)$ for $x$ or $y$ fixed.
Attempt 2: finding a counter-example
I know that there are sets $S$ (such as the Sierpinsky set) that have non-zero measure on the plane but such that every projection on the axis has zero measure. However this is not enough, because, for example, taking the characteristic function of $S$ as a counter-example I always end up with iterated integrals that exist and are zero. I also experimented a bit with two-variables version of Thomae's function, but also in that case I always end up with iterated integrals that exist.