Kron() calculates the Kroniker
product of two matrices,
C = A⊗B. The (i,j) block of
C is ai, jB. This is a very useful operator in
balanced mixed linear models. See Hocking[#!HK:am!#] for
these uses, and Gram[#!Gm:kp!#] for its applications in
matrix calculus.