function cm = centroarea(x,y) cm = sum(x.*y)/sum(y); end