Array4d v(-0.5,2,0,-7);
cout << v.erfc() << endl;