Polyhedron/external/eigen-3.4.0/doc/snippets/MatrixBase_col.cpp

Matrix3d m = Matrix3d::Identity();
m.col(1) = Vector3d(4,5,6);
cout << m << endl;