function value = hump_int ( a, b )

%*****************************************************************************80
%
%% hump_int evaluates the definite integral of the hump() function.
%
%  Licensing:
%
%    This code is distributed under the GNU LGPL license.
%
%  Modified:
%
%    06 October 2019
%
%  Author:
%
%    John Burkardt
%
%  Input:
%
%    real a, b: the endpoints of the interval.
%
%  Output:
%
%    real value(): the value of the definite integral of the humps
%    function from a to b.
%
  va = ( 1.0 / 0.1 ) * atan ( ( a - 0.3 ) / 0.1 ) ...
     + ( 1.0 / 0.2 ) * atan ( ( a - 0.9 ) / 0.2 ) ...
     - 6.0 * a;

  vb = ( 1.0 / 0.1 ) * atan ( ( b - 0.3 ) / 0.1 ) ...
     + ( 1.0 / 0.2 ) * atan ( ( b - 0.9 ) / 0.2 ) ...
     - 6.0 * b;

  value = vb - va;

  return
end

