Hello! Thanks for responding. I'm not using H, a, and z in this loop, no. They are saved and used later in the code. I just needed this function to get their values. The way I understood the gsl to work is that it updated the values of y[] and worked with those, and so I was saving their...