t = 13;
n = length(X);
x0 = sum(X)/n;
I = (1:n);
rxt = sum((X - x0).*(I-t))/sqrt(sum((X-x0).^2)*sum((I-t).^2));
end