Reconstruction of the density matrix as a constrained optimization problem

We present a numerical algorithm for the maximum-likelihood estimation of t he density matrix, and apply it to the homodyne tomography of a single-mode radiation field. The algorithm is based on a specific form of the Gauss de composition for positive definite Hermitian matrices. Results from Monte Ca rlo simulated experiments are presented.