%% quad_sequential_call.m
  a = 0.0;
  b = 2.0;
  n = 21;
  f = @(x) (4*x+1)*(11-4*x)/8;
  tic
  q = quad_sequential ( a, b, n, f );
  toc
  fprintf ( 1, '  The integral estimate is %f\n', q );
