Add two vectors

std::vector<double> add(const std::vector<double>& a, const std::vector<double>& b) {
  std::vector<double> c(a.size());
  for (int i = 0; i < a.size(); i++) {
    c[i] = a[i] + b[i];
  }
  return c;
}