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