I've come across the following generalization of Cauchy integral formula in higher dimensions.
Let $h \in \mathcal{O}(U)$, with $U \subset_{\text{open}}$ Let $(p_1, p_2) \in U$. There exists a bidisk $D(p_1, r_1) \times D(p_2, r_2) \subset U$, and we have\begin{align}{ f(p_1,p_2) = \frac {1}{(2\pi i)^2}}\oint\limits_{\partial D(p_1, r_1) \times \partial D(p_2, r_2)} {\frac{f(z,\xi)}{(z - p_1)(\xi - p_2)}} d\xi\,dz \end{align}
I understand that by applying Cauchy formula in one dimension, one has \begin{align*} f(p_1,p_2)={\frac {1}{(2\pi i)^2}}\oint \limits_{\partial D(p_1, r_1)_{\circlearrowleft}}{\frac {1}{z-p_1}} \Big(\oint_{\partial D(p_2, r_2)_{\circlearrowleft}} \frac{f(z,\xi)}{\xi - p_2} d\xi \Big)\,dz\\ \end{align*}
Nevertheless, I don't know how should one understand a double complex integral (i.e complex surface integral), as presented in the theorem. Most books I've encountered don't seem to go much over the details other that invoking compactedness and continuity with Fubini's theorem to go from the iterated line-integral to the formula. While I would sense this indeed is the right thing to do, it seems irrelevant as I don't know how to define such higher complex integral.