We propose a numerical method to obtain complex amplitude distribution of a three-dimensional (3D) object from a digital hologram. The method consists of two processes. The first process is to measure simultaneously a hologram of the 3D object and an object intensity distribution by two image sensors. These intensity distributions give us the amplitude and absolute value of phase of the 3D object at the image sensor plane. The second process is the determination of phase distribution by a proposed iterative process based on the criterion that the reconstructed 3D object is in focus and its conjugate reconstruction is out of focus. Numerical and experimental results show the effectiveness of the proposed method.