Solution of Problem 22 in section 4 in chapter 12 of Royden (real analysis third edition (2)) Here is the problem:
Let $(X, \mathcal{S}, \mu)$ and $(Y, \mathcal{T}, \lambda)$ be $\sigma$-finite measure spaces. Suppose that $g: X \rightarrow \mathbb{R}$ is a $\mu-$integrable function and that $h: Y \rightarrow \mathbb{R}$ is a $\lambda-$integrable function. Define $f: X \times Y \rightarrow \mathbb{R}$ by $f(x,y) = g(x)h(y).$ Prove that $f$ is $\mu \times \lambda$ integrable and that $$\int_{X\times Y} f d(\mu \times \lambda) = (\int_{X}g d\mu) (\int_{Y} hd\lambda).$$
Here is its solution:



My Question is:
What are the details of calculations of the fourth line "It follows that the result holds for simple functions and thus nonnegative integrable functions", does the case of nonnegative functions follows by Tonelli's theorem (its statement is:  {Let $(X, \mathcal{S}, \mu)$ and $(Y, \mathcal{T}, \lambda)$ be $\sigma$-finite measure spaces. And let $f$ be a $\mathcal{S} \times \mathcal{T}$ measurable function on $X \times Y.$
(a) If $0 \leq f \leq \infty$ and  $\varphi(x) = \int_{Y}f_{x}d\lambda(y), \psi(y) = \int_{X}f^{y}d\mu(x)$ then $\varphi$ is  $\mathcal{S}-$measurable and $\psi$ is  $\mathcal{T}-$measurable
 and $$\int_{X} \varphi d\mu = \int_{X \times Y} f d(\mu \times \lambda) = \int_{Y} \psi d\lambda.\quad \quad (**)$$ )?
Now for the case of simple functions, Here is my trial:
Let $E_{i} \subset X$ and $F_{j} \subset Y$ where both  $E_{i}$ and $F_{j}$ have finit measures for each $i$ and $j$ (I am not sure if this assumption required or it follows directly from that $X$ and $Y$ have $\sigma-$ finite measures) and let $g = \sum^{n} c_{i} \chi_{E_{i} }$ and $h = \sum^{m} d_{j} \chi_{F_{j} }$ then 
$\int_{X\times Y} f d(\mu \times \lambda) = \int_{X\times Y} (\sum^{n} c_{i} \chi_{E_{i}} )(\sum^{m} d_{j} \chi_{F_{j}} ) = \int_{X\times Y} \sum^{n} (c_{i}  \sum^{m} d_{j} \chi_{F_{j}}. \chi_{E_{i}} ) = \int_{X\times Y} \sum^{n} (c_{i}  \sum^{m} d_{j} \chi_{F_{j} \times E_{i}} ) = \int_{X\times Y} (\sum^{n} \sum^{m} c_{i}   d_{j} \chi_{F_{j} \times E_{i}} ).$
But then I am unable to complete. could anyone help me in completing this proof, please?
This question is also problem 10 on pg.423 in Royden "real analysis fourth edition" and I prefer an answer depending on this edition of Royden  
 A: If $g = \sum\limits_k a_k \mathbf{1}_{\mathrm{A}_k}$ and $h = \sum\limits_j b_j \mathbf{1}_{\mathrm{B_j}},$ the indices $k$ and $j$ running on finite sets, then $f = \sum\limits_{(k,j)} a_k b_j \mathbf{1}_{\mathrm{A}_k \times \mathrm{B}_j},$ with the indices $(k, j)$ running on a finite set. Elementary properties of integral (linearity & definition of product measure) give at once
$$
\begin{align*}
\int\limits_{\mathrm{S} \times \mathrm{T}} d(\mu \otimes \lambda)(s,t) f(s,t) &= \sum\limits_{(k,j)}\ \int\limits_{\mathrm{S} \times \mathrm{T}} d(\mu \otimes \lambda)(s,t) \mathbf{1}_{\mathrm{A}_k \times \mathrm{B}_j}(s,t) \\
&= \sum\limits_{(k,j)}\ \mu(\mathrm{A}_k) \lambda(\mathrm{B}_j) = \sum\limits_k \mu(\mathrm{A}_k) \sum\limits_j \lambda(\mathrm{B}_j) \\
&=\int\limits_{\mathrm{S}} d\mu(s) g(s) \int\limits_{\mathrm{T}} d\lambda(t) h(t). \blacksquare
\end{align*}
$$
Now let us assume that $g$ and $h$ are any non-negative measurable functions. There exists increasing sequences $g_n$ and $h_m$  of simple functions converging simply to $g$ and $h,$ respectively. Apply Monotone convergence theorem to get
$$
\begin{align*}
\int\limits_{\mathrm{S} \times \mathrm{T}} d(\mu \otimes \lambda)(s,t) g_n(s) h(t) &= \lim_m \int\limits_{\mathrm{S} \times \mathrm{T}} d(\mu \otimes \lambda)(s,t) g_n(s) h_m(t) \\
&= \int\limits_{\mathrm{S}} d\mu(s) g_n(s) \lim_m \int\limits_{\mathrm{T}} d\lambda(t) h_m(t) \\
&= \int\limits_{\mathrm{S}} d\mu(s) g_n(s) \int\limits_{\mathrm{T}} d\lambda(t) h(t).
\end{align*}$$
Now let $n \to \infty,$ apply MTC again to reach the desired result.
