function p = WelchTest(Mean1, Std1, N1, Mean2, Std2, N2)

%Data = [Mean1 Mean2; Std1 Std2];

a = Std1*Std1/N1;
b = Std2*Std2/N2;

T = (Mean1-Mean2)/sqrt(a+b);
v = (a+b)*(a+b)/((a*a/(N1-1))+(b*b/(N2-1)));

%disp('P two-tailed:')
p = min(1,2*(1-tcdf(abs(T),v)));

%oder
%p = min(1, 2*tcdf((-1)*abs(T),v))

%auf der Internetseite wird v gerundet. Dadurch bekommt man leicht andere
%p-Werte

