%Script to run banditS

    m=randn(10,1);
    %[As,Q,R,p]=banditS(1000,m,1);
    [As,Q,R,p]=banditS1(1000,m);
    M=exp(m)./sum(exp(m));
    plot(1:10,p,1:10,M);