Random matrices with external source and multiple orthogonal polynomials