s=mmstate(varargin)

% mmstate   Density matrix for the maximally mixed state 
%    mmstate(n,d) creates an n-qudit maximally mixed state where d is the 
%    dimension of the qudits. If d is omitted then it is taken to be 2.
%    If argument n is omitted than the default is taken to be
%    the value of global variable N.

function s=mmstate(varargin)

if isempty(varargin)
    global N;
    d=2;
elseif length(varargin)==1
    N=varargin{1};
    d=2;
elseif length(varargin)==2
    N=varargin{1};
    d=varargin{2};
else
    error('Wrong number of input arguments.')
end %if

s=eye(d^N)/d^N;