%grupa leczona grupa kontrolna
D=[6.9 6.4
7.6 6.7
7.3 5.4
7.6 8.2
6.8 5.3
7.2 6.6
8.0 5.8
5.5 5.7
5.8 6.2
7.3 7.1
8.2 7.0
6.9 6.9
6.8 5.6
5.7 4.2
8.6 6.8];
%srednia: 
srednie=[7.08 6.26];
wektor=[D(:,1)' D(:,2)'];
A=D(:,1)';
B=D(:,2)';
tosamo=(mean(A)-mean(B))/var([A B]);
liczba=1e4;
sukces=zeros(liczba,1);
for i=1:liczba
    A=randsample(wektor,15,'true');
    B=randsample(wektor,15,'true');
    sukces(i)=(mean(A)-mean(B))/var([A B]);
end
p_1=sum(sukces>=tosamo)/liczba

len=length(wektor);
for i=1:liczba
    permutacja=wektor(randperm(len));
    A=permutacja(1:len/2); 
    B=permutacja(len/2+1:len);
    
    sukces(i)=(mean(A)-mean(B))/var([A B]);
end
p_2=sum(sukces>=tosamo)/liczba