I have this code to define a function that works to do a kronecker product but after the is not workinig
kron <- function(A,B){ nfila = nrow(A)*nrow(B) ncolum