plot ln(erfc(z/sqrt(2))/2), z = 0 .. 30