clc;
randn('seed',0);
var=10;
N=10;
M=1000;
T=zeros(M,1);
for i=1:M
    x=sqrt(var)*randn(N,1);
    T(i)=mean(x);
end
ngam=100;
gammamin=min(T);
gammamax=max(T);
gamdel=(gammamax-gammamin)/ngam;
gamma=[gammamin:gamdel:gammamax]';
P=zeros(length(gamma),1);Ptrue=P;
for i=1:length(gamma)
    clear Mgam;
    Mgam=find(T>gamma(i));
    P(i)=length(Mgam)/M;
end

Ptrue=Q(gamma/(sqrt(var/N)));
piot(gamma,P,'-',gamma,Ptrue,'--')
xlabel('gamma')
ylabel('P(T>gamma)')
grid