Array3d v(M_PI, M_PI/2, M_PI/3);
cout << v.cos() << endl;