general/geometry/intervalquadrature.m

Go to the documentation of this file.

function res = intervalquadrature(poldeg,func,varargin)
%function res = intervalquadrature(poldeg,func,varargin)
% integration of function func over reference interval == unit interval.
% by Gaussian quadrature exactly integrating polynoms of poldeg.
%
% arguments:
%  poldeg:  degree of polynomial the quadrature rule exactly approximates
%           (0-23)
%  func:    is a function getting a local coordinate vector (1d) and giving a
%           (vectorial or scalar) result
%  varargin: optional further arguments for function
%
% return values:
%  res:     result of quadrature
%

% Bernard Haasdonk 27.8.2009

[ points, weights ] = get_quadrature_weights_1d(poldeg);

res = quadrature(weights,points,func,varargin{:});

return