clear all

load T4-1.DAT

X=T4_1;

qqplot(X);
n=length(X);

for i=1:n
    q(i)=norminv((i-1/2)/n,0,1);
end 

X_sort=sort(X);
r_Q=corr(q',X_sort)

