This is part of my code that calculates xi using a class cosmology.xi_DM.
I am very new to C++, so please bear with me! 
double redshift = 0.75;
  string method = "CAMB";
  vector<double> r = {1., 2.};
  double xi = 0.;
  for (int i=0; i<r.size(); i++)
    xi = cosmology.xi_DM(r[i], method, redshift);
    cout << "xi_DM(z=" << redshift << ") = " << xi << endl;
However, when I print it, I only get the value of xi for r = 2.. It does not print xi for r = 1.. Why is this so? 
 
     
    