%% quad_parfor_call.m
  w = 4;

  a = 0.0;
  b = 2.0;
  n = 21;
  f = @(x) (4*x+1)*(11-4*x)/8;

  matlabpool ( 'open', 'local', w );
  tic
  q = quad_parfor ( a, b, n, f );
  toc
  matlabpool close

  fprintf ( 1, '  The integral estimate is %f\n', q );
