$\bar {\partial } $ -extension of the matrix Riemann–Hilbert method is used to study asymptotics of the polynomials $ P_n(z) $ satisfying orthogonality relations

$$ \begin{align*} \int_{-1}^1 x^lP_n(x)\frac{\rho(x)dx}{\sqrt{1-x^2}}=0, \quad l\in\{0,\ldots,n-1\}, \end{align*} $$

where $ \rho (x) $ is a positive $ m $ times continuously differentiable function on $ [-1,1] $ , $ m\geq 3 $ .

In this work type II Hermite-Padé approximants for a vector of Cauchy transforms of smooth Jacobi-type densities are considered. It is assumed that densities are supported on mutually disjoint intervals (an Angelesco system with complex weights). The formulae of strong asymptotics are derived for any ray sequence of multi-indices.