function weights = gauss_weights(Nq,h)

  if Nq==2
    weights = [h/2 h/2 h/2 h/2];
  else
    weights(1) = 49*h/(12*(18+sqrt(30)));
    weights(2) = 49*h/(12*(18-sqrt(30)));
    weights(3) = weights(2);
    weights(4) = weights(1);
    weights(5) = weights(1);
    weights(6) = weights(2);
    weights(7) = weights(2);
    weights(8) = weights(1);
  end