
% Integrate f(x) on [-1,1] using gll quadrature

close all; clear all; format compact; format longe

for N=1:20;

    [z,w]=zwgll(N);

    f=1./(1.2+cos(5*z)); 

    quad=w'*f 

end;


% Integrate f(x) on [a,b] using gll quadrature:

%%% ...
