%
% Kronecker product: X = A*B 
% X = kp(A,B);
%
%
function X = kp(A,B)
   X = [];
   [m,n] = size(A);
   for i=1:m,
      h = [];
      for j=1:n,
         h = [h A(i,j)*B];
      end
      X = [X; h];
   end
end
