Gaussian quadrature for multiple orthogonal polynomials