procedure TPTEST(DATA1, DATA2:array of real; N:integer;var T,PROB:real);
var
    COV,DF,AVE1,VAR1,AVE2,VAR2,SD:real;  J:integer;
begin
    AVEVAR(DATA1, N, AVE1, VAR1);
    AVEVAR(DATA2, N, AVE2, VAR2);
    COV:=0;
    For J:=1 To N do
        COV:=COV + (DATA1[J] - AVE1) * (DATA2[J] - AVE2);
    DF:=N - 1;
    COV:=COV / DF;
    SD:=Sqrt((VAR1 + VAR2 - 2  * COV) / N);
    T:=(AVE1 - AVE2) / SD;
    PROB:=BETAI(0.5 * DF, 0.5, DF / (DF + T*T));
end;