double mean(double* A, int n)
  {
    double sum = 0.0;
    for(int k = 0; k < n; k++) 
    {
      sum = sum + A[k];
    }
    return sum/n;
  }