cartprod.m

function result = cartprod( varargin )
    %% result = cartprod( varargin )
    %
    % Inputs
    % varargin - Each additional parameter is the possible range of each column.
    %
    % Outputs
    % result - The cartesian product of all parameters

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