% SFEJER  Symbolic version of fejer.m
% 
function uv=sfejer(dig,N)
digits(dig);
n=N:-1:1; m=1:floor(N/2);
th=vpa((2*n-1)*pi./(2*N)); 
uv(:,1)=vpa(cos(th'));
for k=N:-1:1
  s=sum(cos(2*m*th(k))./(4*(m.^2)-1));
  uv(k,2)=2*(1-2*s)/N;
end
