#include <boost/numeric/ublas/fwd.hpp>
#include <boost/numeric/ublas/vector.hpp>
#include <boost/numeric/ublas/vector_proxy.hpp>
#include <boost/numeric/ublas/vector_expression.hpp>
#include <boost/numeric/ublas/storage.hpp>
#include <boost/numeric/ublas/io.hpp>
#include <iostream>
using std::cout;
using std::endl;

using boost::numeric::ublas::vector;
using boost::numeric::ublas::vector_reference;
using boost::numeric::ublas::vector_range;
using boost::numeric::ublas::range;

typedef vector<double> dvector;
typedef vector_range<dvector> dvrange;

int main(void)
{
  dvector x(10);
  dvrange xr(x, range(1,5) );  
  vector_reference<dvector> refx(x);
  vector_reference<const dvector> crefx(x);

  vector_reference<dvector> refxr(xr.data());
  vector_reference<const dvector> crefxr(xr.data());

  for(size_t i=0; i<10; ++i)
    x(i) = i;

  cout << x << endl;
  cout << refx << endl;
  cout << crefx << endl;
  cout << refxr << endl;
  cout << crefxr << endl;
}