%gauss nodes
%inputs: a,b : endpoints of interval
%        n : number of nodes
% brad nelson 2/5/12

function [x,w]=gauss_nodes(a,b,n)
[x,w]=gauss(n);
w=w*(b-a)/2;
x=(b-a)*x/2+(b+a)/2;

