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