#include <iostream>
#include <boost/numeric/ublas/matrix.hpp>
#include <boost/numeric/ublas/io.hpp>
namespace ublas = boost::numeric::ublas;
using namespace std;
int main()
{
ublas::matrix<double> v(3,3);
ublas::matrix<double> u(3,3);
v(0,0) = 3.0; v(0,1) = 0.0; v(0,2) = 4.0;
v(1,0) = 3.0; v(1,1) = 0.0; v(1,2) = 4.0;
v(2,0) = 3.0; v(2,1) = 0.0; v(2,2) = 4.0;
u(0,0) = 3.0; u(0,1) = 0.0; u(0,2) = 4.0;
u(1,0) = 3.0; u(1,1) = 0.0; u(1,2) = 4.0;
u(2,0) = 3.0; u(2,1) = 0.0; u(2,2) = 4.0;
ublas::matrix<double> v1 = v + u;
ublas::matrix<double> v2 = v - u;
ublas::matrix<double> v3 = v * 2.0;
ublas::matrix<double> v4 = v / 3.0;
ublas::matrix<double> v5 = prod(v, u) ;
cout << "v1 : " << v1 << endl;
}