m = function r = corco(X,Y)
sumxy=0;
sumy=0;
sumx=0;
sumxx=0;
sumyy=0;
n = length(X);
for ii= 1:n
    sumx = sumx + X(ii);
    sumy =sumy + Y(ii);
    sumxy = sumxy + X(ii)*Y(ii);
    sumxx = sumxx+X(ii)*X(ii);
    sumyy = sumyy + Y(ii)*Y(ii);
end
p=n*sumxx - sumx*sumx;
q=n*sumyy - sumy*sumy;
r = (n*sumxy - sumx*sumy)/sqrt(p*q);