This is a theorem in Yeh's Real Analysis :
Let $(X, \mathbf{A}, \mu)$ and $(Y, \mathbf{B}, v)$ be two finite measure spaces. Consider the product measure space $(X \times Y, \sigma(\mathbf{\mathbf { A }} \times \mathbf{B}), \mu \times v) .$ Let $\mathbf{P}$ be the collection of nonnegative functions on $X \times Y$ of the type $$\psi(x, y)=\sum_{j=1}^{k} c_{j} \mathbf{1}_{A_{j}}(x) \mathbf{1}_{B_{j}}(y)$$ where $c_{j}>0, A_{j} \in \mathbf{A},$ and $B_{j} \in \mathbf{B}$ for $j=1, \ldots, k .$ Then forevery nonnegative extended real-valued $\sigma(\mathbf{A} \times \mathbf{B})$ -measurable function $f$ on $X \times Y,$ there exists an increasing sequence $\left(\psi_{n}: n \in \mathbb{N}\right)$ in $\mathbf{P}$ such that $\psi_{n} \uparrow f$ a.e. on $X \times Y .$
The suthor first approximate $f$ using simple functions in the product measure space, and then approximate sets in $\sigma(\mathbf{\mathbf { A }} \times \mathbf{B})$ by sets in $\mathbf{\mathbf { A }} \times \mathbf{B}$.
However, I find that the proof does not completely solve the problem. Instead, it only proves that there exists a sequence of $\psi_{n}$ that converges a.e. to $f$, but they may not be monotone.
Can anyone give a correct proof to this theorem please ?