Numerical Solution of Riemann–Hilbert Problems: Random Matrix Theory and Orthogonal Polynomials