function xi = gnodes(n)

% this generates the appropriate
% gaussian nodes on the interval [-1,1]

switch n
	case 1	% this means n=1
		xi = [0];
	case 2	% this means n=2
		xi = [-0.577350269189626,0.577350269189626];
	case 3	% this means n=3
		xi = [-0.774596669241483,0,0.774596669241483];
	case 4	% this means n=4
		xi = [-0.861136311594053, -0.339981043584856, ...
			0.339981043584856, 0.861136311594053];
	case 5	% this means n=5
		xi = [-0.906179845938664, -0.538469310105683, 0,...
			0.538469310105683, 0.906179845938664];
	case 6	% this means n=6
		xi = [-0.932469514203152, -0.661209386466265, ...
			-0.238619186083197, 0.238619186083197, ...
			0.661209386466265, 0.932469514203152];
end