function q = quad_sequential ( a, b, n, f )
%% quad_sequential.m

  q = 0.0;
  for i = 0 : n
    x = ( ( n - i ) * a + i * b ) / n;
    q = q + f ( x );
  end 
  q = q * ( b - a ) / n; 

  return
end
