matthiasm@8: function s = qnorm(inmat, varargin)
matthiasm@8: dim = 1;
matthiasm@8: if nargin > 1
matthiasm@8:     q = varargin{1};
matthiasm@8:     if nargin > 2
matthiasm@8:         dim = varargin{2};
matthiasm@8:     end
matthiasm@8: else
matthiasm@8:     q = 2;
matthiasm@8: end
matthiasm@8: if ~(q == Inf)
matthiasm@8:     s = sum(inmat.^q,dim).^(1/q);
matthiasm@8: else
matthiasm@8:     s = max(inmat,[],dim);
matthiasm@8: end