Sub TPTEST(DATA1(), DATA2(), N, T, PROB)
    Call AVEVAR(DATA1(), N, AVE1, VAR1)
    Call AVEVAR(DATA2(), N, AVE2, VAR2)
    COV = 0#
    For J = 1 To N
        COV = COV + (DATA1(J) - AVE1) * (DATA2(J) - AVE2)
    Next J
    DF = N - 1
    COV = COV / DF
    SD = Sqr((VAR1 + VAR2 - 2# * COV) / N)
    T = (AVE1 - AVE2) / SD
    PROB = BETAI(0.5 * DF, 0.5, DF / (DF + T ^ 2))
End Sub