Thanks. However, it is not easy as what you thought. The problem is: x \in (0,\infinity). If we take z=x-i b and we have dz=dx, but the Hermite polynomial becomes
H[2n+1, z+ib] and integral interval becomes (-i b,\infinity). Using the mathematica, it still cannot find the solutions.