\[ \int (e+f x)^2 \left (a+b \coth ^{-1}(c+d x)\right ) \, dx \]
Optimal antiderivative \[ \frac {b f \left (-c f +d e \right ) x}{d^{2}}+\frac {b \,f^{2} \left (d x +c \right )^{2}}{6 d^{3}}+\frac {\left (f x +e \right )^{3} \left (a +b \,\mathrm {arccoth}\left (d x +c \right )\right )}{3 f}+\frac {b \left (-c f +d e +f \right )^{3} \ln \left (-d x -c +1\right )}{6 d^{3} f}-\frac {b \left (d e -\left (1+c \right ) f \right )^{3} \ln \left (d x +c +1\right )}{6 d^{3} f} \]
command
integrate((f*x+e)^2*(a+b*arccoth(d*x+c)),x, algorithm="giac")
Giac 1.9.0-11 via sagemath 9.6 output
\[ \text {output too large to display} \]
Giac 1.7.0 via sagemath 9.3 output \[ \int {\left (f x + e\right )}^{2} {\left (b \operatorname {arcoth}\left (d x + c\right ) + a\right )}\,{d x} \]_______________________________________________________