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