n = 1000;

av1 = 0;
av2 = 1;
n1 = 30;
n2 = 30;
sigma1 = 1;
sigma2 = 2;

ps = length(n);
for I = 1:n
    
    % rng(123)
    x1 = randn(1, n1)*sigma1 + av1;
    x2 = randn(1, n1)*sigma2 + av2;

    [~, p] = ttest(x1, x2);
    ps(I) = log(p);

    % pairs(table(x1', x2'))
end

histogram(ps)
xline(log(0.05));