function mischeck(m) while 1 A = pnormdim(rand(m),1); B = pnormdim(rand(m),1); lam = rand; AB = lam*A + (1-lam)*B; f1 = lambda2(A); f2 = lambda2(B); f3 = lambda2(AB); good = f3 + 1e-9 < lam*f1 + (1-lam)*f2 if ~good keyboard end end return for t=1:100 m = 2; %a = randn(m,1); %a(m+1) = 1-sum(a); %b = randn(m,1); %b(m+1) = 1-sum(b); a = rand(m+1,1); a=a/sum(a); b = rand(m+1,1); b=b/sum(b); s1 = sum(abs(a-b)); s2 = 2*sum(abs(a(1:m)-b(1:m))); d = abs(s1-s2); if d<=1e-10 'good' else keyboard end end return