N_1=20;
N_2=50;
N_3=100;
N_4=500;
N_5=1000;
N_6=5000;

x=normrnd(100,40,[1,10000]);

N_1_random=datasample(x,N_1);
N_1_mean=mean(N_1_random);
N_1_var=var(N_1_random);

SEM_1 = std(N_1_random)/sqrt(length(N_1_random));
ts_1 = tinv(0.93,length(N_1_random)-1);
CI_1 = mean(N_1_random) + ts_1*SEM_1;

plot(log10(N_1),CI_1);

N_2_random=datasample(x,N_2);
N_2_mean=mean(N_2_random);
N_2_var=var(N_2_random);

SEM_2 = std(N_2_random)/sqrt(length(N_2_random));
ts_2 = tinv(0.93,length(N_2_random)-1);
CI_2 = mean(N_2_random) + ts_2*SEM_2;

N_3_random=datasample(x,N_3);
N_3_mean=mean(N_3_random);
N_3_var=var(N_3_random);

SEM_3 = std(N_3_random)/sqrt(length(N_3_random));
ts_3 = tinv(0.93,length(N_3_random)-1);
CI_3 = mean(N_3_random) + ts_3*SEM_3;

N_4_random=datasample(x,N_4);
N_4_mean=mean(N_4_random);
N_4_var=var(N_4_random);

SEM_4 = std(N_4_random)/sqrt(length(N_4_random));
ts_4 = tinv(0.93,length(N_4_random)-1);
CI_4 = mean(N_4_random) + ts_4*SEM_4;

N_5_random=datasample(x,N_5);
N_5_mean=mean(N_5_random);
N_5_var=var(N_5_random);

SEM_5 = std(N_5_random)/sqrt(length(N_5_random));
ts_5 = tinv(0.93,length(N_5_random)-1);
CI_5 = mean(N_5_random) + ts_5*SEM_5;

N_6_random=datasample(x,N_6);
N_6_mean=mean(N_6_random);
N_6_var=var(N_6_random);

SEM_6 = std(N_6_random)/sqrt(length(N_6_random));
ts_6 = tinv(0.93,length(N_6_random)-1);
CI_6 = mean(N_6_random) + ts_6*SEM_6;