function res = coefb( n )

N = 5000;

som = 0;

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

res = som;


end

