 Solution

Let $N_{ij}$ be the observed cell counts, $N_{i+}$ the row totals, $N_{+j}$ the column totals and $N=33009$. The fitted expected frequencies under independence are $E_{ij}=N_{i+}{N}_{+j}/N$. For example the three-child, one-boy cell has expected count $7541\cdot 16472/33009=3763.1$.

The Pearson chi-squared statistic isUnder $H_0$ it has asymptotically a chi-squared distribution with $(2-1)(3-1)=2$ degrees of freedom, assuming the multinomial sampling model and sufficiently large expected counts. All expected counts here are large. The distribution is an asymptotic approximation, rather than an exact finite-sample identity.