#include<stdio.h>
#include <math.h>
#include <time.h>
#include<stdlib.h>

#include "basestats.h"
#include "ttest.h"

float tTest(int n, int arr1[n], 
            int m, int arr2[m]) 
{ 
    float mean1 = calculateMean(n, arr1); 
    float mean2 = calculateMean(m, arr2); 
    float sd1 = calculateSD(n, arr1); 
    float sd2 = calculateSD(m, arr2); 
  
    // Formula to find t-test 
    // of two set of data. 
    float t_test = (mean1 - mean2) / sqrt((sd1 * sd1) 
                              / n + (sd2 * sd2) / m); 
    return t_test; 
} 

int main(void)
{
    srand(time(NULL));
    int r1[1000];
    int r2[1000];
    int lr1, mr1;
    int lr2, mr2;
    printf("Range 1: low, high: ");
    scanf("%d %d", &lr1, &mr1);
    printf("Range 2: low, high: ");
    scanf("%d %d", &lr2, &mr2);
    for (int i=0; i<1000; i++)
    {
	r1[i]=randrange(lr1, mr1);
	r2[i]=randrange(lr2, mr2);
    }

    printf("T-test: %6.3f\n", tTest(1000, r1, 1000, r2));
}