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