Any sufficiently smooth spatial wave function inside a 2-D box can be expanded in a double sum of the product wave functions, i.e. \begin{equation} \psi(x,y)=\sum_{n=1}^{\infty}\sum_{m=1}^{\infty}\, c_{nm}\; \hbox{eigenfunction}_n(x)\;\hbox{eigenfunction}_m(y) \end{equation} Using your expressions from part (a) above, write out all the terms in this sum out to \(n=3\), \(m=3\). Arrange the terms, conventionally, in terms of increasing energy.
You may find it easier to work in bra/ket notation: \begin{align*} \left|{\psi}\right\rangle &=\sum_{n=1}^{\infty}\sum_{m=1}^{\infty}\, c_{nm}\left|{n}\right\rangle \left|{m}\right\rangle \\ &=\sum_{n=1}^{\infty}\sum_{m=1}^{\infty}\, c_{nm}\left|{nm}\right\rangle \end{align*}