#include <stdio.h>
#include "stat.h"

#define N 1000000

main() {
    int i;
    double x[N];
    double y[N];
    double a, b, c;

    for (i = 0; i < N; i++) {
        x[i] = 2.0*normsim();
        y[i] = 1.0+2.0*x[i]+3.0*normsim();
    }

    /* Should print out approximately  0.0  1.0  2.0  5.0 */
    printf("%lf %lf %lf %lf\n",mean(N,x),mean(N,y),sdev(N,x),sdev(N,y));

    /* Should print out approximately  8.0  0.8 */
    printf("%lf %lf\n",cov(N,x,y),corr(N,x,y));

    /* Should print out approximately  1.0  2.0  3.0 */
    linreg(N,x,y,&a,&b,&c);
    printf("%lf %lf %lf\n",a,b,c);
}