mas01mc@285: function L = lshL(w,k,delta)
mas01mc@285: if nargin<3, delta=0.01;end
mas01mc@285: if nargin<2, k=10;end
mas01mc@285: if nargin<1, w=4;end
mas01mc@285: 
mas01mc@287: P1=lshP(w,1);
mas01mc@287: L = ceil(log(1./delta)/-log(1-P1^k));
mas01mc@285: 
mas01mc@285: endfunction