function value = poly_norm ( c, a, b )

%*****************************************************************************80
%
%% poly_norm evaluates the L2 norm of a polynomial over [a,b].
%
%  Input:
%
%    real c(d+1): the polynomial coefficients.
%
%    real a, b: the endpoints of the interval.
%
%  Output:
%
%    real value: the L2 norm of the polynomial over [a,b].
%
  c2 = poly_square ( c );
  c2a = poly_anti_value ( c2, a );
  c2b = poly_anti_value ( c2, b );
  value = sqrt ( c2b - c2a );

  return
end

