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