double vec_sum(double *xvec, int n) {
double xsum = 0.0;
for (int i = 0; i < n; i++) {xsum += xvec[i];}
return xsum;
}