Introduce optFlag so that one can choose not to build the optimization-related cell of matrices called B. It is too computationally costly and should probably be done in a different way.

function A = kr(varargin)
    n = nargin;

    A = 1;

    for i = 1:n
        A = kron(A,varargin{i});
    end
end