erf( <var> )
erf(<var>) = 2/sqrt(pi) * integral from 0 to <var> of exp(-x*x) dx
set-double d = erf(0.5) + 1