Let U N be a family of N × N independent Haar unitary random matrices and their adjoints, Z N a family of deterministic matrices, P a self-adjoint noncommutative polynomial, i.e. such that for any N , P (U N , Z N ) is self-adjoint, f a smooth function. We prove that for any k, if f is smooth enough, there exist deterministic constantsBesides the constants α P i (f, Z N ) are built explicitly with the help of free probability. As a corollary, we prove that given α < 1/2, for N large enough, every eigenvalue of P (U N , Z N ) is N −α -close to the spectrum of P (u, Z N ) where u is a d-tuple of free Haar unitaries. We also prove the convergence of the norm of any polynomial P (U N ⊗ IM , IN ⊗ Y M ) as long as the family Y M converges strongly and that M ≪ N ln −5/2 (N ).