function testlambdap
nn=5:50;
for n=nn
    for m=1:1000
        a{m}=100*randn(n,1);
        L=tril(randn(n)); for i=1:n, L(i,i)=1; end
        D=diag(rand(n,1));
        Q{m}=L'*D*L;
    end
    tic;
    for m=1:1000
        [N,f]=lambda(a{m},Q{m},2);
    end
    mt(n)=toc;
    disp(sprintf('n=%2d,t=%10.3fms',n,mt(n)))
end
gut('newfigg','','',[600,400]);
ggt('newaxes','',[0.08,0.1,0.84,0.8]);
plot(nn,mt(nn),'.-')
xlabel('N : Number of Ambiguities')
ylabel('Execution Time (msec)')
title('LAMBDA/MLAMBDA PERFORMANCE')
