1 | kron(m1, m2)
|
m1 |
|
m2 |
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 | ##---- Should be DIRECTLY executable !! ----
##-- ==> Define data, use random,
##-- or do help(data=index) for the standard data sets.
## The function is currently defined as
function (m1, m2)
{
m1 <- as.matrix(m1)
m2 <- as.matrix(m2)
kron <- vector(mode = "numeric", length = 0)
for (i in 1:nrow(m1)) {
m3 <- m1[i, 1] * m2
for (j in 2:ncol(m1)) m3 <- cbind(m3, m1[i, j] * m2)
if (i == 1)
kron <- m3
if (i >= 2)
kron <- rbind(kron, m3)
}
kron
}
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.