Orthonormal basis for product $L^2$ space
Solution 1:
The way I would do it is first to show that $\{f_n\times g_m\}$ is orthonormal: indeed, $$ \langle f_n\times g_m, f_s\times g_t\rangle = \langle f_n,f_s\rangle\,\langle g_m,g_t\rangle =\delta_{n,s}\,\delta_{m,t} = \delta_{(n,m),(s,t)}. $$
And then it only remains to show that it is a basis; for that, we can show that its orthogonal complement is zero. So, if $\langle h, f_n\times g_m\rangle = 0 $ for all $n,m$, we have $$ 0=\int_X\left(\int_Y h(x,y)\,g_m(y)\,d\nu(y)\right)\,f_n(x)\,d\mu(x); $$ so, as $n$ is arbitrary, the function $x\mapsto\int_Y h(x,y)\,g_m(y)\,d\nu(y)$ is zero almost everywhere for each $m$. Let $$ E_m=\{x\in X: \int_Y h(x,y)\,g_m(y)\,d\nu(y)\ne0\}. $$ Each $E_m$ is a null-set, and then so is its (countable) union $E$. Outside of $E$, $$ \int_Y h(x,y)\,g_m(y)\,d\nu(y)=0\ \ \ \mbox{for all } m. $$ Thus for each $x\in X\setminus E$, $h(x,y)=0$ almost everywhere. As $|h|^2$ is integrable, its integral agrees with the iterated integrals, so $$ \int_{X\times Y} |h(x,y)|^2\, d(\mu\times\nu)=\int_X\int_Y|h(x,y)|^2\,d\nu(y)\,d\mu(x) =\int_{X\setminus E}\int_Y|h(x,y)|^2\,d\nu(y)\,d\mu(x)=0. $$ So $h=0$ in $L^2(X\times Y)$.