Compute the error function.
x
2 /
erf (x) = ----- |exp (-t²) dt
√π /
0
Accuracy: The result is a tight enclosure.
erf (infsup (1)) ⇒ ans ⊂ [0.8427, 0.84271]
See also: @infsup/erfc.
Package: interval