  Several characterizations are given for a square matrix that can be written
as the product of two positive (semidefinite) projections. Based on one of
these characterizations, and the theory of alternating projections, a Matlab
program is written to check the condition and construct the two positive
projections whose product equal to the given matrix, if they exist.
