trasp <- function(A) {
  A = as.matrix(A)       # trasforma un eventuale vettore
                         # in una matrice di una sola colonna
  nr = dim(A)[2]
  nc = dim(A)[1]  
  tA = matrix( rep(NA, nr*nc), c(nr, nc) )
  for(i in 1:nr) {
    for(j in 1:nc) {
      tA[i,j] = A[j,i]
    }
  }
  return(tA)
}