function mcoef = gegenbauer_to_monomial ( n, alpha, gcoef )

%*****************************************************************************80
%
%% gegenbauer_to_monomial() converts from Gegenbauer to monomial form.
%
%  Licensing:
%
%    This code is distributed under the MIT license.
%
%  Modified:
%
%    26 February 2024
%
%  Author:
%
%    John Burkardt.
%
%  Input:
%
%    integer n: the order of the polynomial.
%
%    real alpha: the Gegenbauer parameter.
%
%    real gcoef(1:n+1): the Gegenbauer coefficients of the polynomial.
% 
%  Output:
%
%    real mcoef(1:n+1): the monomial coefficients of the polynomial.
%
  if ( n < 0 )
    return
  end
%
%  Force gcoef to be a column vector!
%
  gcoef = gcoef(:);
%
%  Get matrix.
%
  A = gegenbauer_to_monomial_matrix ( n + 1, alpha );

  mcoef = A * gcoef;

  return
end
