With Mathematica, the

** Series** function gives a power series expansion. But I am looking for a different type of expansion. Let me explain :

$Assumptions = Element[x, Reals] && x > 0

f := Exp[x^2]*Erfc[x]

Series[f, {x, 0, 5}]

This gives the result :

SeriesData[x, 0, {

1, (-2) Pi^Rational[-1, 2], 1, Rational[-4, 3] Pi^Rational[-1, 2],

Rational[1, 2], Rational[-8, 15] Pi^Rational[-1, 2]}, 0, 6, 1]

But I am interested in getting an asymptotic expansion of the formĀ given in Abrahamson and Stegun (as shown below). How can I get these type of expansions with Mathematica? I would appreciate any help that I can get.