% Berechnet Knoten und Gewichte der Gauss'schen Quadratur im
% Intervall [-1,1].
% @param n Anzahl Stuetzstellen -1
% @return Stuetzstellen x, Gewichte a
function [x,a] = gaussQuad(n)
	b = 1:n; b = b./sqrt(4*b.*b-1);
	J = diag(b,-1)+diag(b,1); [ev,ew]=eig(J);
	for i=1:n+1
		ev(:,i) = ev(:,i)./norm(ev(:,i));
	end
	x=diag(ew); a=(2*(ev(1,:).*ev(1,:)));
end