I try for this ...

I want to see the result like this

Plot[(-1.047-1.50133 I)*E^((-0.34381+1.35844 I) t) *(1+ Erf[(0.72714+0.9341
I)* Sqrt(t))+((-1.047+1.50133 I)*E^((-0.34381-1.35844 I) t)*
1+Erf[(0.72714-0.9341 I)*Sqrt(t))]
In seems that your function has an imaginary part of zero for all inputs in range:
expr[t_] := (-1.047 - 1.50133 I)*
E^((-0.34381 + 1.35844 I) t)*(1 + Erf[(0.72714 + 0.9341 I) *Sqrt[t]]) + (-1.047 +
1.50133 I)*E^(( -0.34381 - 1.35844 I) t)*(1 + Erf[(0.72714 - 0.9341 I)*Sqrt[t]])
Table[Im[expr[t]], {t, 0, 15, 0.1}] // Union
{0.}
i, and in Mathematica imaginary unit is a capitalI. – m0nhawk May 23 '15 at 12:49eshould beE. Second,Plotdoes not handle complex numbers. – bbgodfrey May 23 '15 at 12:49