function ccii_asymptotic_test ( )

%*****************************************************************************80
%
%% ccii_asymptotic_test(): asymptotic error for Boyd's Hermite quadrature.
%
%  Discussion:
%
%    The integral being approximated is
%      integral ( 0 <= x < +oo ) f(x) exp ( -x^2 ) dx
%    or
%     integral ( 0 <= x < +oo ) f(x) exp ( -x^2/2 ) dx
%    or
%      integral ( 0 <= x < +oo ) f(x) dx
%
%  Licensing:
%
%    This code is distributed under the MIT license.
%
%  Modified:
%
%    25 May 2014
%
%  Author:
%
%    John Burkardt
%
  fprintf ( 1, '\n' );
  fprintf ( 1, 'CCII_ASYMPTOTIC_TEST\n' );
  fprintf ( 1, '  Examine asymptotic quadrature error for Boyd''s\n' );
  fprintf ( 1, '  quadrature rule for\n' );
  fprintf ( 1, '    I(f) = integral ( -oo < x < +oo ) f(x) dx\n' );
  fprintf ( 1, '  with\n' );
  fprintf ( 1, '    f(x) = exp ( -px ) / (1+exp(-qx))\n' );

  p = 1.0;
  q = 3.0;

  n_min = 5;
  n_inc = 5;
  n_max = 50;
  option = 1;
  f_integral = pi / q / sin ( pi * p / q );
  ell = 1.0;

  ccii_asymptotic ( n_min, n_inc, n_max, option, ...
    @(x) exp(-x)./(1+exp(-3*x)), f_integral, ell );

  return
end
