title: Using odeint function definition
body: <p>Pretty noob question so please bear with me.</p>
  
  <p>I am following the example given here-->
  <a href="http://www.codeproject.com/Articles/268589/odeint-v2-Solving-ordinary-differential-equations" rel="nofollow">http://www.codeproject.com/Articles/268589/odeint-v2-Solving-ordinary-differential-equations</a></p>
  
  <p>In particular, I am looking at this function:</p>
  
  <pre><code>void lorenz( state_type &amp;x , state_type &amp;dxdt , double t )
  {                                                         
  dxdt[0] = sigma * ( x[1] - x[0] );
  dxdt[1] = R * x[0] - x[1] - x[0] * x[2];
  dxdt[2] = x[0]*x[1] - b * x[2];
   } 
  </code></pre>
  
  <p>In my case, R takes on a series of values (vector with 100 doubles).</p>
  
  <p>odeint is called as:</p>
  
  <pre><code>integrate_const( runge_kutta4&lt; state_type &gt;() , lorenz , x , 0.0 , 10.0 , dt );
  </code></pre>
  
  <p>I would like to do this for each value of R. How can I accomplish this? My knowledge of C++/OOP is limited, but I am willing to learn.</p>
  
  <p>Thank you.</p>