function vegas
nsamples=10000;
count=0;
for i=1:nsamples
    c=randperm(52);
    points=value(c);
    aces=find(points==1);
    points(aces)=points(aces)+10;
    hand=points(1)+points(3);
    if hand==21
        count=count+1;
    end
end
twentyones=count/nsamples
%
function v = value(x)
% Evaluate hand
v = mod(x-1,13)+1;
v = min(v,10);