  a = 0.0;
  b = 2.0;
  n = 21;
  x = linspace ( a, b, n );
  y = 2 * ( x + 0.25 ) .* ( 2.75 - x );
  xc(1:n-1) = 0.5 * ( x(1:n-1) + x(2:n) );
  yc = 2 * ( xc + 0.25 ) .* ( 2.75 - xc );
  clf
  hold on
  bar ( xc(1:5), yc(1:5), 0.75, 'r' )
  bar ( xc(6:10), yc(6:10), 0.75, 'g' )
  bar ( xc(11:15), yc(11:15), 0.75, 'b' )
  bar ( xc(16:20), yc(16:20), 0.75, 'm' )
  grid on
  axis ( [0, 2, 0, 5.0 ] )
  plot ( x, y, 'k', 'linewidth', 3 )
  hold off

  print ( '-dpng', 'quad_picture_parallel.png' )

