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;