samer@4: function A=pnormalise(p,A,i)
samer@4: % pnormalise - normalise A relative to p-norm in given array dimension
samer@4: %
samer@4: % pnormalise :: nonneg, [D:[[1 E]]], [E] -> [D].
samer@4: if nargin<3, i=1; end
samer@4: A=divnorm(@(x)pnorm(p,x,i),A);
samer@4: 	
samer@4: