function res = rmscant10( t,v,th )


som = 0;
l = length(v);

for k = 1:l
    for kk = 1:l
        som = som +sinc((v(k)-v(kk))*t)*exp(-20*(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
