function res = rmscant( t,n )

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

for k = 1:l
    for kk = 1:l
        som = exp(i(pp(k)-pp(kk))*t)*(thoulessN(pp(k),n)+thoulessN(pp(kk),n))^(-3)
    end
end;

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


end