3 % integration of
function func over reference interval == unit interval.
4 % by Gaussian quadrature exactly integrating polynoms of poldeg.
7 % poldeg: degree of polynomial the quadrature rule exactly approximates
9 % func: is a function getting a local coordinate vector (1d) and giving a
10 % (vectorial or scalar) result
11 % varargin: optional further arguments for function
14 % res: result of quadrature
17 % Bernard Haasdonk 27.8.2009
19 [ points, weights ] = get_quadrature_weights_1d(poldeg);
21 res = quadrature(weights,points,func,varargin{:});