function res = rmscant2( t,x,s,th )

som = 0;
l = length(s);
xx = length(x);

for k1 = 2:xx
    for k = 1:l
        for kk = 1:l
        som = som +x(k1)^2*x(k1-1)*sinc((s(k)-s(kk))*t)*exp(-x(k1)*abs(th(k)+th(kk)));
    end
    end
end;


res = som/(l^2*xx^2);


end
