mas01mc@287: function P2 = lshP(w,c,k)
mas01mc@287: if nargin<3, k=1;end
mas01mc@287: if nargin<2, c=1;end
mas01mc@287: if nargin<1, w=4;end
mas01mc@287: 
mas01mc@287: P2 = 1 - 2*normcdf(-w./c) - 2./(sqrt(2*pi)*(w./c)) .* ( 1-exp(-w.^2./(2*c.^2)) );
mas01mc@287: if(k~=1)
mas01mc@287:   P2 = P2.^k;
mas01mc@287: endif
mas01mc@287: endfunction