1. change %*% to tcrossprod