Square-root solution for the discrete-time lyapunov equation