function res = coefc( n )

N = 1000;

som = 0;

for k = 1:N
    som = som + (1/N)*exp(-2*pi*n*((k-1)/N)*complex(0,1))*corrmu(3^((k-1)/N-1))*2^(1-(k-1)/N);
end;

res = som;


end

