%cdromy
n1=15;
mi_1=6598;
s1=844;
n2=12;
mi_2=6870;
s2=669;
ds=sqrt(1/n1*s1.^2 + 1/n2*s2.^2);
t=(mi_2-mi_1)/ds
df=n1+n2-2
p=1-tcdf(t,25)