function y=multinom_logistic(x)
%y=multinom_logistic(x) computes probabilities for multinomial logistic
%model (softmax function)
%for nd arrays, sum over dimension nd
%if only one vector, use row vector!


%if size(x,2)==1 error('use row vectors'),end
%if isempty(find(sum(x,1)==0)) error ('one column should consist of zeros'), end
d=ndims(x);
y=exp(x)./repmat(sum(exp(x),d),[ones(1,ndims(x)-1) size(x,d)]);