x <- c(2,3,7,10) y <- c(7,1,2,30) x-mean(x) (x-mean(x))*(y-mean(y)) sum((x-mean(x))*(y-mean(y)))