procedure TUTEST(DATA1:array of real; N1:integer;
                    DATA2:array of real; N2:integer;var T, PROB:real);
var
    AAA,AVE1,AVE2,VAR1,VAR2,DF:real;
begin 
    AVEVAR(DATA1, N1, AVE1, VAR1);
    AVEVAR(DATA2, N2, AVE2, VAR2);
    T:=(AVE1 - AVE2) / Sqrt(VAR1 / N1 + VAR2 / N2);
    AAA:=Sqr(VAR1 / N1 + VAR2 / N2);
    DF:=AAA / (Sqr(VAR1 / N1) / (N1 - 1) + Sqr(VAR2 / N2) / (N2 - 1));
    PROB:=BETAI(0.5 * DF, 0.5, DF / (DF + T*T));
end;