       



fromair=[2.31035 2.31026 2.31024 2.31012 2.31027 ...
      2.31017 2.30986 2.31010 2.31001 ...
      2.31024 2.31010 2.31028];
 fromchem =[2.30143 2.29890 2.29816 2.30182 ...
      2.29869 2.29940 ...
      2.29849 2.29889];
 na = length(fromair); nc = length(fromchem);
 ma= mean(fromair)
%ans =2.3102

 mc=  mean(fromchem)
%ans =2.2995

 s2a = var(fromair)
%ans =1.8752e-008

 s2c = var(fromchem)
%ans =1.9022e-006

 f=s2a/s2c  %0.0099

pv1 = fcdf(f, na-1, nc-1)  %4.1626e-009

 t0 = (ma - mc)/sqrt(s2a/na + s2c/nc) %21.8598
 
 t100 = (ma - mc - 1/100)/sqrt(s2a/na + s2c/nc)  %1.4189
 
 
 
figure(1)
bihist( fromair, fromchem, 100)
%print -depsc 'C:\Springer\Two\Twoeps\Argon.eps'