function r=corco(x,y)

mx=mean(x);
my=mean(y);
covxy=mean((x-mx).*(y-my));
r=covxy/(std(x)*std(y));