%% Example 3.5
% Gauss - Hermite quadrature

n   = 7;
[x,w] = gauss_her(n);
I   = dot(w,cos(x))