kprod(diag(2), t(ones(2)))
kprod(t(ones(2)), diag(2))
kprod(diag(2), t(ones(2)), t(ones(2)))
kprod(t(ones(2)), diag(2), t(ones(2)))
kprod(t(ones(2)), t(ones(2)), diag(2))
rbind(
kprod(diag(2), t(ones(2))),
kprod(t(ones(2)), diag(2))
)
Run the code above in your browser using DataLab