function result = cartprod(varargin)
    %% result = cartprod(varargin)
    %
    % Inputs
    % varargin -
    %
    % Outputs
    % result -

    c = cell(1, numel(varargin));
    [c{:}] = ndgrid( varargin{:} );
    result = cell2mat( cellfun(@(v)v(:), c, 'UniformOutput',false) );