function res = thoulessN( x,n )

pp = P(n);
som = 0;
l = 2^n;

for k = 1:l
    if x == pp(k)
        som = som +0;
    else
    som = som + log(abs(x-pp(k)));
end;

res = 2^(-n)*som;


end
