N = 1000000;
results52 = zeros(1, N);
results26 = zeros(1, N);

for i = 1 : N
    cards52 = randperm(52);
    results52(i) = mod(cards52(1), 2) == mod(cards52(2), 2);
    
    cards26 = randperm(52);
    cards26 = cards26(1:26);
    results26(i) = mod(cards26(1), 2) == mod(cards26(2), 2);
end

fprintf('%d trials\n', N);
fprintf('mean 52: %.03f mean 26: %.03f\n', mean(results52), mean(results26));
fprintf('var 52: %.03f var 26: %.03f\n', var(results52), var(results26));