Help with int(x), getting erf in answer

I am trying to take the definite integral of e^(-x^2)*sin(x) from 0 to inf and keep getting a weird answer. Here is my code and the result:

>> syms x real

>> int(exp(-x^2)*sin(x),x,0,inf)

ans =

-(pi^(1/2)*erf(1i/2)*exp(-1/4)*1i)/2

Any help on what I am doing wrong would be much appreciated. Thanks.

[John D'Errico]

And, why do you think there is a problem? erf is well defined.

vpa(erf(sym('1i/2')))
ans =
0.61495209469651098083968118562364i

Answer 1

[Star Strider]

You’re not doing anything wrong. See the definition of erf and the reason should be straightforward.