function ccii_asymptotic ( n_min, n_inc, n_max, option, f, f_integral, ell )

%*****************************************************************************80
%
%% ccii_asymptotic() computes asymptotic Laguerre quadrature errors.
%
%  Discussion:
%
%    The integral being approximated is either:
%      integral ( -oo < x < +oo ) f(x) exp ( - x^2 ) dx
%    or
%      integral ( -oo < x < +oo ) f(x) dx
%
%  Licensing:
%
%    This code is distributed under the MIT license.
%
%  Modified:
%
%    18 May 2014
%
%  Author:
%
%    John Burkardt
%
%  Input:
%
%    integer N_MIN, N_INC, N_MAX, the minimum, increment, and maximum
%    for the number of points in the rule.
%
%    integer OPTION, specifies the form of the integral.
%    0: integral ( -oo < x < +oo ) f(x) exp ( - x^2 ) dx
%    1: integral ( -oo < x < +oo ) f(x) dx
%
%    function pointer F(x), returns the value of the integrand
%    at the N points X.
%
%    real F_INTEGRAL, the exact value of the integral.
%
%    real ELL, the scaling parameter.
%
  fprintf ( 1, '\n' );
  fprintf ( 1, '     N  |Quad error|\n' );
  fprintf ( 1, '\n' );
  for n = n_min : n_inc : n_max
    if ( option == 0 )
      error ( 'CCII_0 not written yet!' );
    elseif ( option == 1 )
      [ x, w ] = ccii_1 ( n, ell );
    end
    fx = f ( x );
    q = w' * fx;
    e = abs ( f_integral - q );
    fprintf ( 1, '  %4d  %8.2e\n', n, e );
  end

  return
end
