如何编写下面这个R程序

参考下如下代码：

est<-function(x,y){
  mux=mean(x)
  muy=mean(y)
  sx=0
  sy=0
  sxy=0
  for (i in 1:length(x)){
  sx=sx+(x[i]-mux)**2
  sy=sy+(y[i]-muy)**2
  sxy=sxy+(x[i]-mux)*(y[i]-muy) }
  res<-list(mu=c(mux,muy),s=matrix(data=list(sx,sxy,sxy,sy),dimnames=list(c('sx','sxy','sxy','sy'))))
  return(res)
}
x=c(150,149,131,143,133,140)
y=c(40,42,35,38,30,43)
est(x,y)