function [x,w]=gauss_jacobi(N,alpha,beta)
ab=r_jacobi(N, alpha,beta);
xw=gauss(N,ab);
x=xw(:,1); 
w=xw(:,2); 