  In this paper we calculate, in the large N limit, the eigenvalue density of
an infinite product of random unitary matrices, each of them generated by a
random hermitian matrix. This is equivalent to solving unitary diffusion
generated by a hamiltonian random in time. We find that the result is universal
and depends only on the second moment of the generator of the stochastic
evolution. We find indications of critical behavior (eigenvalue spacing scaling
like $1/N^{3/4}$) close to $\theta=\pi$ for a specific critical evolution time
$t_c$.
