demo_quadratures.m

Go to the documentation of this file.

function qres = demo_quadratures(pdeg,qdeg)
%function qres = demo_quadratures(pdeg,qdeg)
% small script demonstrating interval integration of polynomials
% of degree 'pdeg' with quadratures of degree 'qdeg'
%
% Parameters:
%  pdeg: polynomial degree (default = 25)
%  qdeg: quadrature degree (default = 23)
%
% Return values:
%  qres: a matrix with the integration results

% Bernard Haasdonk 28.8.2009

help demo_quadratures;

if nargin < 2
  qdeg = 23;
end;
if nargin < 1
  pdeg = 25;
end;
qres = zeros(pdeg+1,qdeg+1);
for d=0:pdeg
  f = @(x) x.^d;
  for q = 0:qdeg
    qres(d+1,q+1)=intervalquadrature(q,f);
  end;
end;
disp(qres); % otherwise in the demos command window no output appears.