Let $\Omega \subset R^3$ be a bounded open region. It is well known that there exists a smooth surface $\Gamma$ with minimum area and constant mean curvature which is orthogonal to $\partial \Omega$ and divides $\Omega$ into two regions of equal volume. I wonder if there is a numerical algorithm for constructing such $\Gamma$.