function value = quartic_int ( a, b )

%*****************************************************************************80
%
%% quartic_int returns the exact integral of quartic(x) from a to b.
%
%  Input:
%
%    real a, b: the left and right endpoints of the interval.
%
%  Output:
%
%    real value: the definite integral of quartic(x) from a to b.
%
  quartic_anti = @(x) 0.4 * x.^5 - 4.0 * x.^3 / 3.0 + x.^2 / 2.0 + 20.0 * x;

  value = quartic_anti ( b ) - quartic_anti ( a );

  return
end

