function p=testeqmean(x,y); 
s1=sum((x-mean(x)).^2);
s2=sum((y-mean(y)).^2);
s2=(s1+s2)/(length(x)+length(y)-2);
t=abs(mean(x)-mean(y))/(s2/length(x)+s2/length(y));
p=1-tp(t,length(x)+length(y)-2);
