function res = thouless( x )

pp = P(4);
som = 0;

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

res = som;
    

end

