HP Prime bug for approximate(erfc(z)), z as (0+b*i), or erf(b*i)

[Albert Chan]

Cas> a := 3.16845945394*i
Cas> int(exp(-x*x), x=0..a)       → 0.886226925453-51940029439.9*i
Cas> Ans * 2/sqrt(pi)                → 1.-58608047158.3*i
Cas> erf(a)                               → 1.-58608047158.3*i 

Try numerical integration instead (double exponential method, quad6)

Cas> quad(x -> exp(-x*x), 0, a)       → [3833.13077361*i, 1.16978024776e−12]
Cas> Ans[1] * 2/sqrt(pi)                   → 4325.2249097*i == erf(3.16845945394*i)