function [ we ] = weights( nbig)
%  builds a vector with the quadrature weights
% gaussian quadrature with  4 nodes on each subinterval
%  see formula 4.2 in page 1717 
eta(1)= 0.347854845137454;
eta(2)=  0.652145154862546;
eta(3)=eta(2);
eta(4)=eta(1);
n=nbig/4;
h=1/n;
for i=1:nbig
    ind=mod(i-1,4)+1;
    we(i)=eta(ind)*h;
end
%no=wnodes(4);
%su=0;
%for i=1:nbig
%    for j=1:nbig
  %      su=su+we(i)*we(j)*kernel2d(1,no(4),no(4),no(i),no(j));    
 %   end
%end
end


