doc/snippets/Cwise_erf.cpp

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