func <- function(n) {
  a <- matrix(rnorm(n*n), n, n)
  b <- matrix(rnorm(n*n), n, n)
  output <- list(a = a, b = b)
  return(output)
}

out <- func(5)
a1 <- out$a
b1 <- out$b