title: Using linspace but including endpoints in Runge-Kutta method
body: <p>I've been given a stepsize h that to be used in the implementation of a Runge-Kutta method for higher order ODEs. My problem comes when dividing up the interval from the starting time t_0 to the final time t_f. I thought of using <code>N = ceil((t_f-t_i)/h)</code>, and then using</p>
  
  <pre><code>t = linspace(t_0, t_f, N)
  </code></pre>
  
  <p>But I want to keep the points used by the Runge-Kutta algorithm spaced by h for most of the process, is there a way to include the endpoint t_f while keeping the stepsize at <code>h</code> for the first n-1 steps? I tried using</p>
  
  <pre><code>t = t_0:h:t_f
  </code></pre>
  
  <p>But this does not always include the endpoint t_f. </p>