MATLAB's linspace in C++

linspace.cpp

template <typename T = double>
vector<T> linspace(T a, T b, size_t N) {
  T h = (b - a) / static_cast<T>(N-1);
  vector<T> xs(N);
  typename vector<T>::iterator x;
  T val;
  for (x = xs.begin(), val = a; x != xs.end(); ++x, val += h)
    *x = val;
  return xs;
}

A simple addition would be to set xs.back()=b, if the issue pointed out by @garrison is not a problem.