Merged in feature/dirac_discr (pull request #17) Add multi-d dirac discretization with tests Approved-by: Vidar Stiernström <vidar.stiernstrom@it.uu.se> Approved-by: Martin Almquist <malmquist@stanford.edu> Approved-by: Jonatan Werpers <jonatan.werpers@it.uu.se>

function A = kr(varargin)
    n = nargin;

    A = 1;

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