mas01mc@285: function Rho = lshRho(w,c)
mas01mc@285: if nargin<2, c=1;end
mas01mc@285: if nargin<1, w=4;end
mas01mc@285: 
mas01mc@287: P1 = lshP(w,1);
mas01mc@287: P2 = lshP(w,c);
mas01mc@285: 
mas01mc@287: Rho = log(1./P1) ./ log(1./P2);
mas01mc@285: 
mas01mc@285: endfunction