A new, numerical framework for the approximation of solutions to matrix valued Riemann–Hilbert problems is developed, based on a recent method for the homogeneous Painlevé II Riemann–Hilbert problem. We demonstrate its effectiveness by computing solutions to other Painlevé transcendents. An implementation in Mathematica is made available online.