We numerically and experimentally demonstrate an iterative method to simultaneously reconstruct two unknown interfering wavefronts. A three-dimensional interference pattern is analyzed and then Zernike polynomials and the stochastic parallel gradient descent algorithm are used to expand and calculate wavefronts. (C) 2014 Optical Society of America