#include <boost/numeric/ublas/matrix.hpp>
#include <boost/numeric/ublas/io.hpp>
int main () {
  using namespace boost::numeric::ublas;
  using namespace std;  
  typedef double mtype; //real matrices 
  //typedef complex<double> mtype; // complex matr.
  {
    identity_matrix<mtype> m(3);
    cout << "identity matrix  ";
    cout << m << endl;
  }      
  {
    zero_matrix<mtype> m(3,4);
    cout << "    zero matrix  ";
    cout << m << endl;
  }    
  {
    matrix<mtype> m(3,4);
    for (unsigned i=0; i<m.size1();++i)
      for (unsigned j=0; j<m.size2();++j)
	m(i,j) = 3*i+j;
    cout << "      a  matrix  ";
    cout << m << endl;
  }
  return EXIT_SUCCESS; 
}