It was 26 years ago, on 1992, when Fokas, Its and Kitaev, in the context of 2D quantum gravity, discovered that certain Riemann-Hilbert problem was solved in terms of orthogonal polynomials in the real line (OPRL), [1]}. Namely, it was found that the solution of a \(2\times 2\) 
Riemann–Hilbert problem can be expressed in terms of orthogonal polynomials in the real line and its Cauchy transforms. Later, Deift and Zhou combined these ideas with a non-linear steepest descent analysis in a series of papers [2]}, [3]}, [4]}, [5]} which was the seed for a large activity in the field. To mention just a few relevant results let us cite the study of strong asymptotic with applications in random matrix theory, [2]}, [7]}, the analysis of determinantal point processes [8]}, [9]}, [10]}, [11]}, orthogonal Laurent polynomials [12]}, [13]} and Painlevé equations [14]}, [15]}.
