function res = thouless9( x )

pp = P(9);
som = 0;

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

res = -2^(-9)*som;


end
