#include <math.h>
void moment(float data[], int n, float *ave, float *var)

{
    void nrerror(char error_text[]);
    int j;
    float ep=0.0,s,p;
    if (n <= 1) nrerror("n must be at least 2 in moment");
                                              
    s=0.0;
    for (j=1;j<=n;j++) s += data[j];
      *ave=s/n;
    variance = sqrt(pow((s - avg), 2)/(n - 1));