\(\int (c+d x^2)^2 \coth ^{-1}(a x) \, dx\) [37]

Optimal result

Integrand size = 14, antiderivative size = 110 \[ \int \left (c+d x^2\right )^2 \coth ^{-1}(a x) \, dx=\frac {d \left (10 a^2 c+3 d\right ) x^2}{30 a^3}+\frac {d^2 x^4}{20 a}+c^2 x \coth ^{-1}(a x)+\frac {2}{3} c d x^3 \coth ^{-1}(a x)+\frac {1}{5} d^2 x^5 \coth ^{-1}(a x)+\frac {\left (15 a^4 c^2+10 a^2 c d+3 d^2\right ) \log \left (1-a^2 x^2\right )}{30 a^5} \]

[Out]

1/30*d*(10*a^2*c+3*d)*x^2/a^3+1/20*d^2*x^4/a+c^2*x*arccoth(a*x)+2/3*c*d*x^3*arccoth(a*x)+1/5*d^2*x^5*arccoth(a
*x)+1/30*(15*a^4*c^2+10*a^2*c*d+3*d^2)*ln(-a^2*x^2+1)/a^5

Rubi [A] (verified)

Time = 0.10 (sec) , antiderivative size = 110, normalized size of antiderivative = 1.00, number of steps used = 5, number of rules used = 5, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.357, Rules used = {200, 6124, 1608, 1261, 712} \[ \int \left (c+d x^2\right )^2 \coth ^{-1}(a x) \, dx=\frac {d x^2 \left (10 a^2 c+3 d\right )}{30 a^3}+\frac {\left (15 a^4 c^2+10 a^2 c d+3 d^2\right ) \log \left (1-a^2 x^2\right )}{30 a^5}+c^2 x \coth ^{-1}(a x)+\frac {2}{3} c d x^3 \coth ^{-1}(a x)+\frac {1}{5} d^2 x^5 \coth ^{-1}(a x)+\frac {d^2 x^4}{20 a} \]

[In]

Int[(c + d*x^2)^2*ArcCoth[a*x],x]

[Out]

(d*(10*a^2*c + 3*d)*x^2)/(30*a^3) + (d^2*x^4)/(20*a) + c^2*x*ArcCoth[a*x] + (2*c*d*x^3*ArcCoth[a*x])/3 + (d^2*
x^5*ArcCoth[a*x])/5 + ((15*a^4*c^2 + 10*a^2*c*d + 3*d^2)*Log[1 - a^2*x^2])/(30*a^5)

Rule 200

Int[((a_) + (b_.)*(x_)^(n_))^(p_), x_Symbol] :> Int[ExpandIntegrand[(a + b*x^n)^p, x], x] /; FreeQ[{a, b}, x]
&& IGtQ[n, 0] && IGtQ[p, 0]

Rule 712

Int[((d_.) + (e_.)*(x_))^(m_)*((a_.) + (b_.)*(x_) + (c_.)*(x_)^2)^(p_.), x_Symbol] :> Int[ExpandIntegrand[(d +
e*x)^m*(a + b*x + c*x^2)^p, x], x] /; FreeQ[{a, b, c, d, e, m}, x] && NeQ[b^2 - 4*a*c, 0] && NeQ[c*d^2 - b*d*
e + a*e^2, 0] && NeQ[2*c*d - b*e, 0] && IntegerQ[p] && (GtQ[p, 0] || (EqQ[a, 0] && IntegerQ[m]))

Rule 1261

Int[(x_)*((d_) + (e_.)*(x_)^2)^(q_.)*((a_) + (b_.)*(x_)^2 + (c_.)*(x_)^4)^(p_.), x_Symbol] :> Dist[1/2, Subst[
Int[(d + e*x)^q*(a + b*x + c*x^2)^p, x], x, x^2], x] /; FreeQ[{a, b, c, d, e, p, q}, x]

Rule 1608

Int[(u_.)*((a_.)*(x_)^(p_.) + (b_.)*(x_)^(q_.) + (c_.)*(x_)^(r_.))^(n_.), x_Symbol] :> Int[u*x^(n*p)*(a + b*x^
(q - p) + c*x^(r - p))^n, x] /; FreeQ[{a, b, c, p, q, r}, x] && IntegerQ[n] && PosQ[q - p] && PosQ[r - p]

Rule 6124

Int[((a_.) + ArcCoth[(c_.)*(x_)]*(b_.))*((d_.) + (e_.)*(x_)^2)^(q_.), x_Symbol] :> With[{u = IntHide[(d + e*x^
2)^q, x]}, Dist[a + b*ArcCoth[c*x], u, x] - Dist[b*c, Int[u/(1 - c^2*x^2), x], x]] /; FreeQ[{a, b, c, d, e}, x
] && (IntegerQ[q] || ILtQ[q + 1/2, 0])

Rubi steps \begin{align*} \text {integral}& = c^2 x \coth ^{-1}(a x)+\frac {2}{3} c d x^3 \coth ^{-1}(a x)+\frac {1}{5} d^2 x^5 \coth ^{-1}(a x)-a \int \frac {c^2 x+\frac {2}{3} c d x^3+\frac {d^2 x^5}{5}}{1-a^2 x^2} \, dx \\ & = c^2 x \coth ^{-1}(a x)+\frac {2}{3} c d x^3 \coth ^{-1}(a x)+\frac {1}{5} d^2 x^5 \coth ^{-1}(a x)-a \int \frac {x \left (c^2+\frac {2}{3} c d x^2+\frac {d^2 x^4}{5}\right )}{1-a^2 x^2} \, dx \\ & = c^2 x \coth ^{-1}(a x)+\frac {2}{3} c d x^3 \coth ^{-1}(a x)+\frac {1}{5} d^2 x^5 \coth ^{-1}(a x)-\frac {1}{2} a \text {Subst}\left (\int \frac {c^2+\frac {2 c d x}{3}+\frac {d^2 x^2}{5}}{1-a^2 x} \, dx,x,x^2\right ) \\ & = c^2 x \coth ^{-1}(a x)+\frac {2}{3} c d x^3 \coth ^{-1}(a x)+\frac {1}{5} d^2 x^5 \coth ^{-1}(a x)-\frac {1}{2} a \text {Subst}\left (\int \left (-\frac {d \left (10 a^2 c+3 d\right )}{15 a^4}-\frac {d^2 x}{5 a^2}+\frac {-15 a^4 c^2-10 a^2 c d-3 d^2}{15 a^4 \left (-1+a^2 x\right )}\right ) \, dx,x,x^2\right ) \\ & = \frac {d \left (10 a^2 c+3 d\right ) x^2}{30 a^3}+\frac {d^2 x^4}{20 a}+c^2 x \coth ^{-1}(a x)+\frac {2}{3} c d x^3 \coth ^{-1}(a x)+\frac {1}{5} d^2 x^5 \coth ^{-1}(a x)+\frac {\left (15 a^4 c^2+10 a^2 c d+3 d^2\right ) \log \left (1-a^2 x^2\right )}{30 a^5} \\ \end{align*}

Mathematica [A] (verified)

Time = 0.03 (sec) , antiderivative size = 98, normalized size of antiderivative = 0.89 \[ \int \left (c+d x^2\right )^2 \coth ^{-1}(a x) \, dx=\frac {a^2 d x^2 \left (6 d+a^2 \left (20 c+3 d x^2\right )\right )+4 a^5 x \left (15 c^2+10 c d x^2+3 d^2 x^4\right ) \coth ^{-1}(a x)+\left (30 a^4 c^2+20 a^2 c d+6 d^2\right ) \log \left (1-a^2 x^2\right )}{60 a^5} \]

[In]

Integrate[(c + d*x^2)^2*ArcCoth[a*x],x]

[Out]

(a^2*d*x^2*(6*d + a^2*(20*c + 3*d*x^2)) + 4*a^5*x*(15*c^2 + 10*c*d*x^2 + 3*d^2*x^4)*ArcCoth[a*x] + (30*a^4*c^2
+ 20*a^2*c*d + 6*d^2)*Log[1 - a^2*x^2])/(60*a^5)

Maple [A] (verified)

Time = 0.48 (sec) , antiderivative size = 105, normalized size of antiderivative = 0.95

method	result	size
parts	\(\frac {d^{2} x^{5} \operatorname {arccoth}\left (a x \right )}{5}+\frac {2 c d \,x^{3} \operatorname {arccoth}\left (a x \right )}{3}+c^{2} x \,\operatorname {arccoth}\left (a x \right )+\frac {a \left (\frac {d \left (\frac {3}{2} a^{2} d \,x^{4}+10 a^{2} c \,x^{2}+3 d \,x^{2}\right )}{2 a^{4}}+\frac {\left (15 a^{4} c^{2}+10 a^{2} c d +3 d^{2}\right ) \ln \left (a^{2} x^{2}-1\right )}{2 a^{6}}\right )}{15}\)	\(105\)
derivativedivides	\(\frac {\operatorname {arccoth}\left (a x \right ) c^{2} a x +\frac {2 a \,\operatorname {arccoth}\left (a x \right ) c d \,x^{3}}{3}+\frac {a \,\operatorname {arccoth}\left (a x \right ) d^{2} x^{5}}{5}+\frac {5 c \,a^{4} d \,x^{2}+\frac {3 d^{2} a^{4} x^{4}}{4}+\frac {3 a^{2} d^{2} x^{2}}{2}+\frac {\left (15 a^{4} c^{2}+10 a^{2} c d +3 d^{2}\right ) \ln \left (a x -1\right )}{2}-\frac {\left (-15 a^{4} c^{2}-10 a^{2} c d -3 d^{2}\right ) \ln \left (a x +1\right )}{2}}{15 a^{4}}}{a}\)	\(137\)
default	\(\frac {\operatorname {arccoth}\left (a x \right ) c^{2} a x +\frac {2 a \,\operatorname {arccoth}\left (a x \right ) c d \,x^{3}}{3}+\frac {a \,\operatorname {arccoth}\left (a x \right ) d^{2} x^{5}}{5}+\frac {5 c \,a^{4} d \,x^{2}+\frac {3 d^{2} a^{4} x^{4}}{4}+\frac {3 a^{2} d^{2} x^{2}}{2}+\frac {\left (15 a^{4} c^{2}+10 a^{2} c d +3 d^{2}\right ) \ln \left (a x -1\right )}{2}-\frac {\left (-15 a^{4} c^{2}-10 a^{2} c d -3 d^{2}\right ) \ln \left (a x +1\right )}{2}}{15 a^{4}}}{a}\)	\(137\)
parallelrisch	\(-\frac {-12 x^{5} \operatorname {arccoth}\left (a x \right ) a^{5} d^{2}-40 x^{3} \operatorname {arccoth}\left (a x \right ) a^{5} c d -3 d^{2} a^{4} x^{4}-60 c^{2} x \,\operatorname {arccoth}\left (a x \right ) a^{5}-20 c \,a^{4} d \,x^{2}-60 \ln \left (a x -1\right ) a^{4} c^{2}-60 \,\operatorname {arccoth}\left (a x \right ) a^{4} c^{2}-6 a^{2} d^{2} x^{2}-40 \ln \left (a x -1\right ) a^{2} c d -40 \,\operatorname {arccoth}\left (a x \right ) a^{2} c d -20 a^{2} c d -12 \ln \left (a x -1\right ) d^{2}-12 \,\operatorname {arccoth}\left (a x \right ) d^{2}-6 d^{2}}{60 a^{5}}\)	\(163\)
risch	\(\left (\frac {1}{10} d^{2} x^{5}+\frac {1}{3} c d \,x^{3}+\frac {1}{2} c^{2} x \right ) \ln \left (a x +1\right )-\frac {d^{2} x^{5} \ln \left (a x -1\right )}{10}-\frac {c d \,x^{3} \ln \left (a x -1\right )}{3}+\frac {d^{2} x^{4}}{20 a}-\frac {c^{2} x \ln \left (a x -1\right )}{2}+\frac {c d \,x^{2}}{3 a}+\frac {\ln \left (a^{2} x^{2}-1\right ) c^{2}}{2 a}+\frac {5 c^{2}}{9 a}+\frac {d^{2} x^{2}}{10 a^{3}}+\frac {\ln \left (a^{2} x^{2}-1\right ) c d}{3 a^{3}}+\frac {c d}{3 a^{3}}+\frac {\ln \left (a^{2} x^{2}-1\right ) d^{2}}{10 a^{5}}+\frac {d^{2}}{20 a^{5}}\)	\(178\)

[In]

int((d*x^2+c)^2*arccoth(a*x),x,method=_RETURNVERBOSE)

[Out]

1/5*d^2*x^5*arccoth(a*x)+2/3*c*d*x^3*arccoth(a*x)+c^2*x*arccoth(a*x)+1/15*a*(1/2*d/a^4*(3/2*a^2*d*x^4+10*a^2*c
*x^2+3*d*x^2)+1/2*(15*a^4*c^2+10*a^2*c*d+3*d^2)/a^6*ln(a^2*x^2-1))

Fricas [A] (verification not implemented)

none

Time = 0.25 (sec) , antiderivative size = 118, normalized size of antiderivative = 1.07 \[ \int \left (c+d x^2\right )^2 \coth ^{-1}(a x) \, dx=\frac {3 \, a^{4} d^{2} x^{4} + 2 \, {\left (10 \, a^{4} c d + 3 \, a^{2} d^{2}\right )} x^{2} + 2 \, {\left (15 \, a^{4} c^{2} + 10 \, a^{2} c d + 3 \, d^{2}\right )} \log \left (a^{2} x^{2} - 1\right ) + 2 \, {\left (3 \, a^{5} d^{2} x^{5} + 10 \, a^{5} c d x^{3} + 15 \, a^{5} c^{2} x\right )} \log \left (\frac {a x + 1}{a x - 1}\right )}{60 \, a^{5}} \]

[In]

integrate((d*x^2+c)^2*arccoth(a*x),x, algorithm="fricas")

[Out]

1/60*(3*a^4*d^2*x^4 + 2*(10*a^4*c*d + 3*a^2*d^2)*x^2 + 2*(15*a^4*c^2 + 10*a^2*c*d + 3*d^2)*log(a^2*x^2 - 1) +
2*(3*a^5*d^2*x^5 + 10*a^5*c*d*x^3 + 15*a^5*c^2*x)*log((a*x + 1)/(a*x - 1)))/a^5

Sympy [C] (verification not implemented)

Result contains complex when optimal does not.

Time = 0.36 (sec) , antiderivative size = 182, normalized size of antiderivative = 1.65 \[ \int \left (c+d x^2\right )^2 \coth ^{-1}(a x) \, dx=\begin {cases} c^{2} x \operatorname {acoth}{\left (a x \right )} + \frac {2 c d x^{3} \operatorname {acoth}{\left (a x \right )}}{3} + \frac {d^{2} x^{5} \operatorname {acoth}{\left (a x \right )}}{5} + \frac {c^{2} \log {\left (x - \frac {1}{a} \right )}}{a} + \frac {c^{2} \operatorname {acoth}{\left (a x \right )}}{a} + \frac {c d x^{2}}{3 a} + \frac {d^{2} x^{4}}{20 a} + \frac {2 c d \log {\left (x - \frac {1}{a} \right )}}{3 a^{3}} + \frac {2 c d \operatorname {acoth}{\left (a x \right )}}{3 a^{3}} + \frac {d^{2} x^{2}}{10 a^{3}} + \frac {d^{2} \log {\left (x - \frac {1}{a} \right )}}{5 a^{5}} + \frac {d^{2} \operatorname {acoth}{\left (a x \right )}}{5 a^{5}} & \text {for}\: a \neq 0 \\\frac {i \pi \left (c^{2} x + \frac {2 c d x^{3}}{3} + \frac {d^{2} x^{5}}{5}\right )}{2} & \text {otherwise} \end {cases} \]

[In]

integrate((d*x**2+c)**2*acoth(a*x),x)

[Out]

Piecewise((c**2*x*acoth(a*x) + 2*c*d*x**3*acoth(a*x)/3 + d**2*x**5*acoth(a*x)/5 + c**2*log(x - 1/a)/a + c**2*a
coth(a*x)/a + c*d*x**2/(3*a) + d**2*x**4/(20*a) + 2*c*d*log(x - 1/a)/(3*a**3) + 2*c*d*acoth(a*x)/(3*a**3) + d*
*2*x**2/(10*a**3) + d**2*log(x - 1/a)/(5*a**5) + d**2*acoth(a*x)/(5*a**5), Ne(a, 0)), (I*pi*(c**2*x + 2*c*d*x*
*3/3 + d**2*x**5/5)/2, True))

Maxima [A] (verification not implemented)

none

Time = 0.21 (sec) , antiderivative size = 131, normalized size of antiderivative = 1.19 \[ \int \left (c+d x^2\right )^2 \coth ^{-1}(a x) \, dx=\frac {1}{60} \, a {\left (\frac {3 \, a^{2} d^{2} x^{4} + 2 \, {\left (10 \, a^{2} c d + 3 \, d^{2}\right )} x^{2}}{a^{4}} + \frac {2 \, {\left (15 \, a^{4} c^{2} + 10 \, a^{2} c d + 3 \, d^{2}\right )} \log \left (a x + 1\right )}{a^{6}} + \frac {2 \, {\left (15 \, a^{4} c^{2} + 10 \, a^{2} c d + 3 \, d^{2}\right )} \log \left (a x - 1\right )}{a^{6}}\right )} + \frac {1}{15} \, {\left (3 \, d^{2} x^{5} + 10 \, c d x^{3} + 15 \, c^{2} x\right )} \operatorname {arcoth}\left (a x\right ) \]

[In]

integrate((d*x^2+c)^2*arccoth(a*x),x, algorithm="maxima")

[Out]

1/60*a*((3*a^2*d^2*x^4 + 2*(10*a^2*c*d + 3*d^2)*x^2)/a^4 + 2*(15*a^4*c^2 + 10*a^2*c*d + 3*d^2)*log(a*x + 1)/a^
6 + 2*(15*a^4*c^2 + 10*a^2*c*d + 3*d^2)*log(a*x - 1)/a^6) + 1/15*(3*d^2*x^5 + 10*c*d*x^3 + 15*c^2*x)*arccoth(a
*x)

Giac [B] (verification not implemented)

Leaf count of result is larger than twice the leaf count of optimal. 529 vs. \(2 (100) = 200\).

Time = 0.28 (sec) , antiderivative size = 529, normalized size of antiderivative = 4.81 \[ \int \left (c+d x^2\right )^2 \coth ^{-1}(a x) \, dx=\frac {1}{15} \, a {\left (\frac {{\left (15 \, a^{4} c^{2} + 10 \, a^{2} c d + 3 \, d^{2}\right )} \log \left (\frac {{\left | a x + 1 \right |}}{{\left | a x - 1 \right |}}\right )}{a^{6}} - \frac {{\left (15 \, a^{4} c^{2} + 10 \, a^{2} c d + 3 \, d^{2}\right )} \log \left ({\left | \frac {a x + 1}{a x - 1} - 1 \right |}\right )}{a^{6}} + \frac {4 \, {\left (\frac {{\left (5 \, a^{2} c d + 3 \, d^{2}\right )} {\left (a x + 1\right )}^{3}}{{\left (a x - 1\right )}^{3}} - \frac {{\left (10 \, a^{2} c d + 3 \, d^{2}\right )} {\left (a x + 1\right )}^{2}}{{\left (a x - 1\right )}^{2}} + \frac {{\left (5 \, a^{2} c d + 3 \, d^{2}\right )} {\left (a x + 1\right )}}{a x - 1}\right )}}{a^{6} {\left (\frac {a x + 1}{a x - 1} - 1\right )}^{4}} + \frac {{\left (\frac {15 \, {\left (a x + 1\right )}^{4} a^{4} c^{2}}{{\left (a x - 1\right )}^{4}} - \frac {60 \, {\left (a x + 1\right )}^{3} a^{4} c^{2}}{{\left (a x - 1\right )}^{3}} + \frac {90 \, {\left (a x + 1\right )}^{2} a^{4} c^{2}}{{\left (a x - 1\right )}^{2}} - \frac {60 \, {\left (a x + 1\right )} a^{4} c^{2}}{a x - 1} + 15 \, a^{4} c^{2} + \frac {30 \, {\left (a x + 1\right )}^{4} a^{2} c d}{{\left (a x - 1\right )}^{4}} - \frac {60 \, {\left (a x + 1\right )}^{3} a^{2} c d}{{\left (a x - 1\right )}^{3}} + \frac {40 \, {\left (a x + 1\right )}^{2} a^{2} c d}{{\left (a x - 1\right )}^{2}} - \frac {20 \, {\left (a x + 1\right )} a^{2} c d}{a x - 1} + 10 \, a^{2} c d + \frac {15 \, {\left (a x + 1\right )}^{4} d^{2}}{{\left (a x - 1\right )}^{4}} + \frac {30 \, {\left (a x + 1\right )}^{2} d^{2}}{{\left (a x - 1\right )}^{2}} + 3 \, d^{2}\right )} \log \left (-\frac {\frac {\frac {{\left (a x + 1\right )} a}{a x - 1} - a}{a {\left (\frac {a x + 1}{a x - 1} + 1\right )}} + 1}{\frac {\frac {{\left (a x + 1\right )} a}{a x - 1} - a}{a {\left (\frac {a x + 1}{a x - 1} + 1\right )}} - 1}\right )}{a^{6} {\left (\frac {a x + 1}{a x - 1} - 1\right )}^{5}}\right )} \]

[In]

integrate((d*x^2+c)^2*arccoth(a*x),x, algorithm="giac")

[Out]

1/15*a*((15*a^4*c^2 + 10*a^2*c*d + 3*d^2)*log(abs(a*x + 1)/abs(a*x - 1))/a^6 - (15*a^4*c^2 + 10*a^2*c*d + 3*d^
2)*log(abs((a*x + 1)/(a*x - 1) - 1))/a^6 + 4*((5*a^2*c*d + 3*d^2)*(a*x + 1)^3/(a*x - 1)^3 - (10*a^2*c*d + 3*d^
2)*(a*x + 1)^2/(a*x - 1)^2 + (5*a^2*c*d + 3*d^2)*(a*x + 1)/(a*x - 1))/(a^6*((a*x + 1)/(a*x - 1) - 1)^4) + (15*
(a*x + 1)^4*a^4*c^2/(a*x - 1)^4 - 60*(a*x + 1)^3*a^4*c^2/(a*x - 1)^3 + 90*(a*x + 1)^2*a^4*c^2/(a*x - 1)^2 - 60
*(a*x + 1)*a^4*c^2/(a*x - 1) + 15*a^4*c^2 + 30*(a*x + 1)^4*a^2*c*d/(a*x - 1)^4 - 60*(a*x + 1)^3*a^2*c*d/(a*x -
1)^3 + 40*(a*x + 1)^2*a^2*c*d/(a*x - 1)^2 - 20*(a*x + 1)*a^2*c*d/(a*x - 1) + 10*a^2*c*d + 15*(a*x + 1)^4*d^2/
(a*x - 1)^4 + 30*(a*x + 1)^2*d^2/(a*x - 1)^2 + 3*d^2)*log(-(((a*x + 1)*a/(a*x - 1) - a)/(a*((a*x + 1)/(a*x - 1
) + 1)) + 1)/(((a*x + 1)*a/(a*x - 1) - a)/(a*((a*x + 1)/(a*x - 1) + 1)) - 1))/(a^6*((a*x + 1)/(a*x - 1) - 1)^5
))

Mupad [B] (verification not implemented)

Time = 4.40 (sec) , antiderivative size = 115, normalized size of antiderivative = 1.05 \[ \int \left (c+d x^2\right )^2 \coth ^{-1}(a x) \, dx=\frac {a^4\,\left (\frac {c^2\,\ln \left (a^2\,x^2-1\right )}{2}+\frac {d^2\,x^4}{20}+\frac {c\,d\,x^2}{3}\right )+a^2\,\left (\frac {d^2\,x^2}{10}+\frac {c\,d\,\ln \left (a^2\,x^2-1\right )}{3}\right )+\frac {d^2\,\ln \left (a^2\,x^2-1\right )}{10}}{a^5}+c^2\,x\,\mathrm {acoth}\left (a\,x\right )+\frac {d^2\,x^5\,\mathrm {acoth}\left (a\,x\right )}{5}+\frac {2\,c\,d\,x^3\,\mathrm {acoth}\left (a\,x\right )}{3} \]

[In]

int(acoth(a*x)*(c + d*x^2)^2,x)

[Out]

(a^4*((c^2*log(a^2*x^2 - 1))/2 + (d^2*x^4)/20 + (c*d*x^2)/3) + a^2*((d^2*x^2)/10 + (c*d*log(a^2*x^2 - 1))/3) +
(d^2*log(a^2*x^2 - 1))/10)/a^5 + c^2*x*acoth(a*x) + (d^2*x^5*acoth(a*x))/5 + (2*c*d*x^3*acoth(a*x))/3