void tptest(double data1[],double data2[],int n,double& t,double& prob)
{
	double ave1,ave2,var1,var2;
    avevar(data1, n, ave1, var1);
    avevar(data2, n, ave2, var2);
    double cov = 0.0;
    for (int j = 1; j<=n; j++)
	{
        cov = cov + (data1[j] - ave1) * (data2[j] - ave2);
    }
    int df = n - 1;
    cov = cov / df;
    double sd = sqrt((var1 + var2 - 2.0 * cov) / n);
    t = (ave1 - ave2) / sd;
    prob = betai(0.5 * df, 0.5, df / (df + t * t));
}