//  C = Kronecker product of A and B
<m, n> = size(A);
for i = 1:m, ...
   ci = a(i,1)*B; ...
   for j = 2:n, ci = <ci a(i,j)*B>; end ...
   if i = 1, C = ci; else, C = <C; ci>;
