13   float maxNorm = 
mat.colwise().sum().maxCoeff(&maxIndex);
 
   15   std::cout << 
"Maximum sum at position " << maxIndex << std::endl;
 
   17   std::cout << 
"The corresponding vector is: " << std::endl;
 
   18   std::cout << 
mat.col( maxIndex ) << std::endl;
 
   19   std::cout << 
"And its sum is is: " << maxNorm << std::endl;