\[ \int \frac {e^{c+b^2 x^2} \text {Erfc}(b x)}{x^3} \, dx \]
Optimal antiderivative \[ \frac {b^{2} {\mathrm e}^{c} \expIntegral \left (b^{2} x^{2}\right )}{2}-\frac {{\mathrm e}^{b^{2} x^{2}+c} \mathrm {erfc}\left (b x \right )}{2 x^{2}}+\frac {b \,{\mathrm e}^{c}}{x \sqrt {\pi }}-\frac {2 b^{3} {\mathrm e}^{c} x \hypergeom \left (\left [\frac {1}{2}, 1\right ], \left [\frac {3}{2}, \frac {3}{2}\right ], b^{2} x^{2}\right )}{\sqrt {\pi }} \]
command
integrate(exp(b**2*x**2+c)*erfc(b*x)/x**3,x)
Sympy 1.10.1 under Python 3.10.4 output
\[ \frac {b^{2} e^{c} \operatorname {Ei}{\left (b^{2} x^{2} \right )}}{2} + \frac {2 b e^{c} {{}_{2}F_{2}\left (\begin {matrix} - \frac {1}{2}, 1 \\ \frac {1}{2}, \frac {3}{2} \end {matrix}\middle | {b^{2} x^{2}} \right )}}{\sqrt {\pi } x} - \frac {e^{c} e^{b^{2} x^{2}}}{2 x^{2}} \]
Sympy 1.8 under Python 3.8.8 output
\[ \text {Exception raised: AttributeError} \]________________________________________________________________________________________