#include <boost/numeric/ublas/matrix.hpp> using namespace boost::numeric; int main() { ublas::bounded_matrix<double, 3, 3> m, m1, m2; m = m1 * m2; }