function res = rmscant9( t,v,th )


som = 0;
l = length(v);

for k = 1:l
    for kk = 1:l
        som = som +exp(i*(v(k)-v(kk))*t)*exp(-200*(th(k)+th(kk)))*((abs(th(k)+th(kk)))^(-1) + (abs(th(k)+th(kk)))^(-2)+(abs(th(k)+th(kk)))^(-3))*(th(k)*th(kk))^2;
    end
end;

res =som/(l^2);


end
