%
%  CLENCURT nodes x (Chebychev points) and weights w for Clenshaw-Curtis quadrature
%

  function x = clencurt(N)

  theta = pi*(0:N)'/N;
  x = cos(theta);