#######################################
# MM Matematica - matrici.

# Creazione di una matrice con
# elementi elencati per colonne
# se pc=T.
Mm = function (dati,col,righe,pc=F)
matrix(dati,ncol=col,
  nrow=righe,byrow=!pc)

Mm.autovalori = function (A,simm=F)
{eigen(A,only.values=T,
  symmetric=simm)$values}

Mm.autovettori = function (A,simm=F)
{eigen(A,symmetric=simm)$vectors}

# Lista delle colonne di A.
Mm.col = function (A)
lapply(1:ncol(A),
  function (j) A[,j])

Mm.gershgorin = function (A,i)
{centro=A[i,i]
raggio=sum(abs(A[i,]))-abs(centro)
c(centro,raggio)}

# Lista delle righe di A.
Mm.righe = function (A)
lapply(1:nrow(A),
  function (i) A[i,])

Mm.traccia = function (A)
{if (length(A)>1) sum(diag(A))
else c(A)}