function test94
l = 2.0;
sigma = 0.2;
N = 20000;
options = optimset;
m = fminsearch('wast',l*1.3,options,sigma,l,N);

function J = wast(m,sigma,l,N)
len = randn(1,N) * sigma + m;
s = sum(len>=l);
w = sum(len) - s*l;
J = w/s;