#include<SAMPLER.H>
#include <stdio.h>
#include <math.h>
#include <stdlib.h>
double x[100];
int i,n;
double mean;
double std;
void meanstd( double *x, int length, double mean, double std_dev)
{
int i;
double sum_x, sum_sq;
sum_x=0.0;
sum_sq=0.0;
for( i=0; i<length; i++)
{
sum_x=sum_x+x[i];
sum_sq=sum_sq+x[i]*x[i];
}
mean=sum_x/length;
std_dev=sqrt((sum_sq-pow(sum_x, 2)/length)/(length-1));
}
int main (void)
{
SAMPLE;
mean=0;
std=0;
n=100;
for(i=0; i<n; i++)
x[i]=rand()%100;
meanstd(x,n, mean, std);
SAMPLE;
return 0;
}