In this paper, we prove that the Gauss-Picard modular group PU(3, 1; Z[i]) in three complex dimensions can be generated by five given transformations: two Heisenberg translations, two Heisenberg rotations and an involution. Indeed, our method may work for the other higher dimensional Euclidean-Picard modular groups.