Iterative solution of the Lyapunov matrix equation