Autovalori = function (A)
{eigen(A,only.values=T)$values}

Autovettori = function (A)
{eigen(A)$vectors}

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

Matrice = function
  (dati, righe=NULL, col=NULL)
{if (!is.null(righe))
matrix(dati,nrow=righe,byrow=T)
else if (!is.null(col))
matrix(dati,ncol=col)
else NULL}

Traccia = function (A)
{if (length(A)>1) sum(diag(A))
else c(A)}