function c2 = poly_power ( c1, k )

%*****************************************************************************80
%
%% poly_power returns the coefficients of the k-th power of a polynomial.
%
%  Input:
%
%    real c1(d+1): the coefficients of the polynomial.
%
%    integer k: the power to which the polynomial is to be raised.
%
%  Output:
%
%    real c2(d2+1): the k-th power of the polynomial.
%
  c2 = 1.0;
  for i = 1 : k
    c2 = poly_multiply ( c1, c2 );
  end
  return
end

