I am working with Professor Ken McLaughlin to develop numerical techniques for the solution of matrix-valued Riemann-Hilbert problems. We use a modifcation of the trapezoidal method to achieve spectral accuracy and an implementation of the fast multipole method coupled with a multigrid iterative solver to achieve O(n*log(n)) runtime.