"Testing erf"

clear

check(erf(-x) == -erf(x))
check(erf((x,y)) == (erf(x),erf(y)))

check(erf(0) == 0)
check(erf(0.0) == 0)
check(erf(1000.0) == 1)
check(erf(-1000.0) == -1)
