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