x=c(3,4,5,1,6)
y=3:7
xy= cbind(x,y) #creamos una matriz con "x" en la primera columna e "y" en la segunda
xy 
xy + xy #suma de matrices 
xy * xy #multiplicación de matrices 
#similar para otros operadores (-,/,`^)
rowSums (xy) #sumas por filas
colSums (xy) #sumas por columnas
rowMeans(xy) #medias por filas
colMeans (xy)#medias por columnas