triquad function

Set up Gaussian quadrature points and weights for a triangular domain.