\(\int (d+e x) (a+b \arctan (c x^2))^2 \, dx\) [25]

Optimal result
Mathematica [B] (warning: unable to verify)
Rubi [A] (verified)
Maple [F]
Fricas [F]
Sympy [F]
Maxima [F]
Giac [F]
Mupad [F(-1)]
Reduce [F]

Optimal result

Integrand size = 18, antiderivative size = 1325 \[ \int (d+e x) \left (a+b \arctan \left (c x^2\right )\right )^2 \, dx =\text {Too large to display} \] Output:

(-1)^(3/4)*b^2*d*polylog(2,1-2/(1+(-1)^(1/4)*c^(1/2)*x))/c^(1/2)+(-1)^(3/4 
)*b^2*d*polylog(2,1-2/(1-(-1)^(1/4)*c^(1/2)*x))/c^(1/2)+(-1)^(1/4)*b^2*d*p 
olylog(2,1-2/(1+(-1)^(3/4)*c^(1/2)*x))/c^(1/2)-(-1)^(1/4)*b^2*d*arctanh((- 
1)^(3/4)*c^(1/2)*x)^2/c^(1/2)+(-1)^(3/4)*b^2*d*arctan((-1)^(3/4)*c^(1/2)*x 
)^2/c^(1/2)+(-1)^(1/4)*b^2*d*polylog(2,1-2/(1-(-1)^(3/4)*c^(1/2)*x))/c^(1/ 
2)-1/2*(-1)^(1/4)*b^2*d*polylog(2,1+2^(1/2)*((-1)^(3/4)+c^(1/2)*x)/(1+(-1) 
^(3/4)*c^(1/2)*x))/c^(1/2)-1/2*(-1)^(3/4)*b^2*d*polylog(2,1-2^(1/2)*((-1)^ 
(1/4)+c^(1/2)*x)/(1+(-1)^(1/4)*c^(1/2)*x))/c^(1/2)-1/2*(-1)^(3/4)*b^2*d*po 
lylog(2,1+(-1+I)*(1+(-1)^(3/4)*c^(1/2)*x)/(1+(-1)^(1/4)*c^(1/2)*x))/c^(1/2 
)-1/2*(-1)^(1/4)*b^2*d*polylog(2,1-(1+I)*(1+(-1)^(1/4)*c^(1/2)*x)/(1+(-1)^ 
(3/4)*c^(1/2)*x))/c^(1/2)-I*a*b*d*x*ln(1+I*c*x^2)+2*(-1)^(3/4)*a*b*d*arcta 
nh((-1)^(3/4)*c^(1/2)*x)/c^(1/2)-2*(-1)^(1/4)*b^2*d*arctanh((-1)^(3/4)*c^( 
1/2)*x)*ln(2/(1+(-1)^(3/4)*c^(1/2)*x))/c^(1/2)+b*e*(a+b*arctan(c*x^2))*ln( 
2/(1+I*c*x^2))/c+1/2*I*e*(a+b*arctan(c*x^2))^2/c+I*a*b*d*x*ln(1-I*c*x^2)+2 
*(-1)^(1/4)*b^2*d*arctanh((-1)^(3/4)*c^(1/2)*x)*ln(2/(1-(-1)^(3/4)*c^(1/2) 
*x))/c^(1/2)-2*(-1)^(3/4)*a*b*d*arctan((-1)^(3/4)*c^(1/2)*x)/c^(1/2)-2*(-1 
)^(1/4)*b^2*d*arctan((-1)^(3/4)*c^(1/2)*x)*ln(2/(1+(-1)^(1/4)*c^(1/2)*x))/ 
c^(1/2)+2*(-1)^(1/4)*b^2*d*arctan((-1)^(3/4)*c^(1/2)*x)*ln(2/(1-(-1)^(1/4) 
*c^(1/2)*x))/c^(1/2)+(-1)^(1/4)*b^2*d*arctanh((-1)^(3/4)*c^(1/2)*x)*ln(-2^ 
(1/2)*((-1)^(3/4)+c^(1/2)*x)/(1+(-1)^(3/4)*c^(1/2)*x))/c^(1/2)+(-1)^(1/...
 

Mathematica [B] (warning: unable to verify)

Both result and optimal contain complex but leaf count is larger than twice the leaf count of optimal. \(4850\) vs. \(2(1325)=2650\).

Time = 26.49 (sec) , antiderivative size = 4850, normalized size of antiderivative = 3.66 \[ \int (d+e x) \left (a+b \arctan \left (c x^2\right )\right )^2 \, dx=\text {Result too large to show} \] Input:

Integrate[(d + e*x)*(a + b*ArcTan[c*x^2])^2,x]
 

Output:

a^2*d*x + (a^2*e*x^2)/2 + (a*b*d*Sqrt[c*x^2]*(2*Sqrt[c*x^2]*ArcTan[c*x^2] 
- Sqrt[2]*(ArcTan[(-1 + c*x^2)/(Sqrt[2]*Sqrt[c*x^2])] - ArcTanh[(Sqrt[2]*S 
qrt[c*x^2])/(1 + c*x^2)])))/(c*x) + (a*b*e*(c*x^2*ArcTan[c*x^2] + Log[1/Sq 
rt[1 + c^2*x^4]]))/c + (b^2*e*((-I)*ArcTan[c*x^2]^2 + c*x^2*ArcTan[c*x^2]^ 
2 + 2*ArcTan[c*x^2]*Log[1 + E^((2*I)*ArcTan[c*x^2])] - I*PolyLog[2, -E^((2 
*I)*ArcTan[c*x^2])]))/(2*c) + (b^2*d*Sqrt[c*x^2]*(2*Sqrt[c*x^2]*ArcTan[c*x 
^2]^2 + Sqrt[2]*ArcTan[c*x^2]*(2*ArcTan[1 - Sqrt[2]*Sqrt[c*x^2]] - 2*ArcTa 
n[1 + Sqrt[2]*Sqrt[c*x^2]] - Log[1 + c*x^2 - Sqrt[2]*Sqrt[c*x^2]] + Log[1 
+ c*x^2 + Sqrt[2]*Sqrt[c*x^2]]) - Sqrt[2]*((ArcTan[1 - Sqrt[2]*Sqrt[c*x^2] 
] + ArcTan[1 + Sqrt[2]*Sqrt[c*x^2]])*Log[1 + c*x^2 - Sqrt[2]*Sqrt[c*x^2]] 
- (ArcTan[1 - Sqrt[2]*Sqrt[c*x^2]] + ArcTan[1 + Sqrt[2]*Sqrt[c*x^2]])*Log[ 
1 + c*x^2 + Sqrt[2]*Sqrt[c*x^2]] + (Sqrt[c*x^2]*(1 + (1 - Sqrt[2]*Sqrt[c*x 
^2])^2)^(3/2)*(2*(-5*ArcTan[2 + I]*ArcTan[1 - Sqrt[2]*Sqrt[c*x^2]] + 4*Arc 
Tan[1 - Sqrt[2]*Sqrt[c*x^2]]^2 + ((1 + 2*I)*Sqrt[1 + I]*ArcTan[1 - Sqrt[2] 
*Sqrt[c*x^2]]^2)/E^(I*ArcTan[2 + I]) + ((1 - 2*I)*Sqrt[1 - I]*ArcTan[1 - S 
qrt[2]*Sqrt[c*x^2]]^2)/E^ArcTanh[1 + 2*I] - (5*I)*ArcTan[1 - Sqrt[2]*Sqrt[ 
c*x^2]]*ArcTanh[1 + 2*I] + (5*I)*(-ArcTan[2 + I] + ArcTan[1 - Sqrt[2]*Sqrt 
[c*x^2]])*Log[1 - E^((2*I)*(-ArcTan[2 + I] + ArcTan[1 - Sqrt[2]*Sqrt[c*x^2 
]]))] + 5*((-I)*ArcTan[1 - Sqrt[2]*Sqrt[c*x^2]] + ArcTanh[1 + 2*I])*Log[1 
- E^((2*I)*ArcTan[1 - Sqrt[2]*Sqrt[c*x^2]] - 2*ArcTanh[1 + 2*I])] + (5*...
 

Rubi [A] (verified)

Time = 2.51 (sec) , antiderivative size = 1325, normalized size of antiderivative = 1.00, number of steps used = 2, number of rules used = 2, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.111, Rules used = {5397, 2009}

Below are the steps used by Rubi to obtain the solution. The rule number used for the transformation is given above next to the arrow. The rules definitions used are listed below.

\(\displaystyle \int (d+e x) \left (a+b \arctan \left (c x^2\right )\right )^2 \, dx\)

\(\Big \downarrow \) 5397

\(\displaystyle \int \left (d \left (a+b \arctan \left (c x^2\right )\right )^2+e x \left (a+b \arctan \left (c x^2\right )\right )^2\right )dx\)

\(\Big \downarrow \) 2009

\(\displaystyle d x a^2-\frac {2 (-1)^{3/4} b d \arctan \left ((-1)^{3/4} \sqrt {c} x\right ) a}{\sqrt {c}}+\frac {2 (-1)^{3/4} b d \text {arctanh}\left ((-1)^{3/4} \sqrt {c} x\right ) a}{\sqrt {c}}+i b d x \log \left (1-i c x^2\right ) a-i b d x \log \left (i c x^2+1\right ) a+\frac {(-1)^{3/4} b^2 d \arctan \left ((-1)^{3/4} \sqrt {c} x\right )^2}{\sqrt {c}}+\frac {1}{2} e x^2 \left (a+b \arctan \left (c x^2\right )\right )^2+\frac {i e \left (a+b \arctan \left (c x^2\right )\right )^2}{2 c}-\frac {\sqrt [4]{-1} b^2 d \text {arctanh}\left ((-1)^{3/4} \sqrt {c} x\right )^2}{\sqrt {c}}-\frac {1}{4} b^2 d x \log ^2\left (1-i c x^2\right )-\frac {1}{4} b^2 d x \log ^2\left (i c x^2+1\right )+\frac {2 \sqrt [4]{-1} b^2 d \arctan \left ((-1)^{3/4} \sqrt {c} x\right ) \log \left (\frac {2}{1-\sqrt [4]{-1} \sqrt {c} x}\right )}{\sqrt {c}}-\frac {2 \sqrt [4]{-1} b^2 d \arctan \left ((-1)^{3/4} \sqrt {c} x\right ) \log \left (\frac {2}{\sqrt [4]{-1} \sqrt {c} x+1}\right )}{\sqrt {c}}+\frac {\sqrt [4]{-1} b^2 d \arctan \left ((-1)^{3/4} \sqrt {c} x\right ) \log \left (\frac {\sqrt {2} \left (\sqrt {c} x+\sqrt [4]{-1}\right )}{\sqrt [4]{-1} \sqrt {c} x+1}\right )}{\sqrt {c}}+\frac {2 \sqrt [4]{-1} b^2 d \text {arctanh}\left ((-1)^{3/4} \sqrt {c} x\right ) \log \left (\frac {2}{1-(-1)^{3/4} \sqrt {c} x}\right )}{\sqrt {c}}-\frac {2 \sqrt [4]{-1} b^2 d \text {arctanh}\left ((-1)^{3/4} \sqrt {c} x\right ) \log \left (\frac {2}{(-1)^{3/4} \sqrt {c} x+1}\right )}{\sqrt {c}}+\frac {\sqrt [4]{-1} b^2 d \text {arctanh}\left ((-1)^{3/4} \sqrt {c} x\right ) \log \left (-\frac {\sqrt {2} \left (\sqrt {c} x+(-1)^{3/4}\right )}{(-1)^{3/4} \sqrt {c} x+1}\right )}{\sqrt {c}}+\frac {\sqrt [4]{-1} b^2 d \text {arctanh}\left ((-1)^{3/4} \sqrt {c} x\right ) \log \left (\frac {(1+i) \left (\sqrt [4]{-1} \sqrt {c} x+1\right )}{(-1)^{3/4} \sqrt {c} x+1}\right )}{\sqrt {c}}+\frac {\sqrt [4]{-1} b^2 d \arctan \left ((-1)^{3/4} \sqrt {c} x\right ) \log \left (\frac {(1-i) \left ((-1)^{3/4} \sqrt {c} x+1\right )}{\sqrt [4]{-1} \sqrt {c} x+1}\right )}{\sqrt {c}}+\frac {\sqrt [4]{-1} b^2 d \arctan \left ((-1)^{3/4} \sqrt {c} x\right ) \log \left (1-i c x^2\right )}{\sqrt {c}}-\frac {\sqrt [4]{-1} b^2 d \text {arctanh}\left ((-1)^{3/4} \sqrt {c} x\right ) \log \left (1-i c x^2\right )}{\sqrt {c}}+\frac {b e \left (a+b \arctan \left (c x^2\right )\right ) \log \left (\frac {2}{i c x^2+1}\right )}{c}-\frac {\sqrt [4]{-1} b^2 d \arctan \left ((-1)^{3/4} \sqrt {c} x\right ) \log \left (i c x^2+1\right )}{\sqrt {c}}+\frac {\sqrt [4]{-1} b^2 d \text {arctanh}\left ((-1)^{3/4} \sqrt {c} x\right ) \log \left (i c x^2+1\right )}{\sqrt {c}}+\frac {1}{2} b^2 d x \log \left (1-i c x^2\right ) \log \left (i c x^2+1\right )+\frac {(-1)^{3/4} b^2 d \operatorname {PolyLog}\left (2,1-\frac {2}{1-\sqrt [4]{-1} \sqrt {c} x}\right )}{\sqrt {c}}+\frac {(-1)^{3/4} b^2 d \operatorname {PolyLog}\left (2,1-\frac {2}{\sqrt [4]{-1} \sqrt {c} x+1}\right )}{\sqrt {c}}-\frac {(-1)^{3/4} b^2 d \operatorname {PolyLog}\left (2,1-\frac {\sqrt {2} \left (\sqrt {c} x+\sqrt [4]{-1}\right )}{\sqrt [4]{-1} \sqrt {c} x+1}\right )}{2 \sqrt {c}}+\frac {\sqrt [4]{-1} b^2 d \operatorname {PolyLog}\left (2,1-\frac {2}{1-(-1)^{3/4} \sqrt {c} x}\right )}{\sqrt {c}}+\frac {\sqrt [4]{-1} b^2 d \operatorname {PolyLog}\left (2,1-\frac {2}{(-1)^{3/4} \sqrt {c} x+1}\right )}{\sqrt {c}}-\frac {\sqrt [4]{-1} b^2 d \operatorname {PolyLog}\left (2,\frac {\sqrt {2} \left (\sqrt {c} x+(-1)^{3/4}\right )}{(-1)^{3/4} \sqrt {c} x+1}+1\right )}{2 \sqrt {c}}-\frac {\sqrt [4]{-1} b^2 d \operatorname {PolyLog}\left (2,1-\frac {(1+i) \left (\sqrt [4]{-1} \sqrt {c} x+1\right )}{(-1)^{3/4} \sqrt {c} x+1}\right )}{2 \sqrt {c}}-\frac {(-1)^{3/4} b^2 d \operatorname {PolyLog}\left (2,1-\frac {(1-i) \left ((-1)^{3/4} \sqrt {c} x+1\right )}{\sqrt [4]{-1} \sqrt {c} x+1}\right )}{2 \sqrt {c}}+\frac {i b^2 e \operatorname {PolyLog}\left (2,1-\frac {2}{i c x^2+1}\right )}{2 c}\)

Input:

Int[(d + e*x)*(a + b*ArcTan[c*x^2])^2,x]
 

Output:

a^2*d*x - (2*(-1)^(3/4)*a*b*d*ArcTan[(-1)^(3/4)*Sqrt[c]*x])/Sqrt[c] + ((-1 
)^(3/4)*b^2*d*ArcTan[(-1)^(3/4)*Sqrt[c]*x]^2)/Sqrt[c] + ((I/2)*e*(a + b*Ar 
cTan[c*x^2])^2)/c + (e*x^2*(a + b*ArcTan[c*x^2])^2)/2 + (2*(-1)^(3/4)*a*b* 
d*ArcTanh[(-1)^(3/4)*Sqrt[c]*x])/Sqrt[c] - ((-1)^(1/4)*b^2*d*ArcTanh[(-1)^ 
(3/4)*Sqrt[c]*x]^2)/Sqrt[c] + (2*(-1)^(1/4)*b^2*d*ArcTan[(-1)^(3/4)*Sqrt[c 
]*x]*Log[2/(1 - (-1)^(1/4)*Sqrt[c]*x)])/Sqrt[c] - (2*(-1)^(1/4)*b^2*d*ArcT 
an[(-1)^(3/4)*Sqrt[c]*x]*Log[2/(1 + (-1)^(1/4)*Sqrt[c]*x)])/Sqrt[c] + ((-1 
)^(1/4)*b^2*d*ArcTan[(-1)^(3/4)*Sqrt[c]*x]*Log[(Sqrt[2]*((-1)^(1/4) + Sqrt 
[c]*x))/(1 + (-1)^(1/4)*Sqrt[c]*x)])/Sqrt[c] + (2*(-1)^(1/4)*b^2*d*ArcTanh 
[(-1)^(3/4)*Sqrt[c]*x]*Log[2/(1 - (-1)^(3/4)*Sqrt[c]*x)])/Sqrt[c] - (2*(-1 
)^(1/4)*b^2*d*ArcTanh[(-1)^(3/4)*Sqrt[c]*x]*Log[2/(1 + (-1)^(3/4)*Sqrt[c]* 
x)])/Sqrt[c] + ((-1)^(1/4)*b^2*d*ArcTanh[(-1)^(3/4)*Sqrt[c]*x]*Log[-((Sqrt 
[2]*((-1)^(3/4) + Sqrt[c]*x))/(1 + (-1)^(3/4)*Sqrt[c]*x))])/Sqrt[c] + ((-1 
)^(1/4)*b^2*d*ArcTanh[(-1)^(3/4)*Sqrt[c]*x]*Log[((1 + I)*(1 + (-1)^(1/4)*S 
qrt[c]*x))/(1 + (-1)^(3/4)*Sqrt[c]*x)])/Sqrt[c] + ((-1)^(1/4)*b^2*d*ArcTan 
[(-1)^(3/4)*Sqrt[c]*x]*Log[((1 - I)*(1 + (-1)^(3/4)*Sqrt[c]*x))/(1 + (-1)^ 
(1/4)*Sqrt[c]*x)])/Sqrt[c] + I*a*b*d*x*Log[1 - I*c*x^2] + ((-1)^(1/4)*b^2* 
d*ArcTan[(-1)^(3/4)*Sqrt[c]*x]*Log[1 - I*c*x^2])/Sqrt[c] - ((-1)^(1/4)*b^2 
*d*ArcTanh[(-1)^(3/4)*Sqrt[c]*x]*Log[1 - I*c*x^2])/Sqrt[c] - (b^2*d*x*Log[ 
1 - I*c*x^2]^2)/4 + (b*e*(a + b*ArcTan[c*x^2])*Log[2/(1 + I*c*x^2)])/c ...
 

Defintions of rubi rules used

rule 2009
Int[u_, x_Symbol] :> Simp[IntSum[u, x], x] /; SumQ[u]
 

rule 5397
Int[((a_.) + ArcTan[(c_.)*(x_)^(n_)]*(b_.))^(p_)*((d_) + (e_.)*(x_))^(m_.), 
 x_Symbol] :> Int[ExpandIntegrand[(a + b*ArcTan[c*x^n])^p, (d + e*x)^m, x], 
 x] /; FreeQ[{a, b, c, d, e, n}, x] && IGtQ[p, 1] && IGtQ[m, 0]
 
Maple [F]

\[\int \left (e x +d \right ) {\left (a +b \arctan \left (c \,x^{2}\right )\right )}^{2}d x\]

Input:

int((e*x+d)*(a+b*arctan(c*x^2))^2,x)
 

Output:

int((e*x+d)*(a+b*arctan(c*x^2))^2,x)
 

Fricas [F]

\[ \int (d+e x) \left (a+b \arctan \left (c x^2\right )\right )^2 \, dx=\int { {\left (e x + d\right )} {\left (b \arctan \left (c x^{2}\right ) + a\right )}^{2} \,d x } \] Input:

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

Output:

integral(a^2*e*x + a^2*d + (b^2*e*x + b^2*d)*arctan(c*x^2)^2 + 2*(a*b*e*x 
+ a*b*d)*arctan(c*x^2), x)
 

Sympy [F]

\[ \int (d+e x) \left (a+b \arctan \left (c x^2\right )\right )^2 \, dx=\int \left (a + b \operatorname {atan}{\left (c x^{2} \right )}\right )^{2} \left (d + e x\right )\, dx \] Input:

integrate((e*x+d)*(a+b*atan(c*x**2))**2,x)
 

Output:

Integral((a + b*atan(c*x**2))**2*(d + e*x), x)
 

Maxima [F]

\[ \int (d+e x) \left (a+b \arctan \left (c x^2\right )\right )^2 \, dx=\int { {\left (e x + d\right )} {\left (b \arctan \left (c x^{2}\right ) + a\right )}^{2} \,d x } \] Input:

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

Output:

12*b^2*c^2*e*integrate(1/16*x^5*arctan(c*x^2)^2/(c^2*x^4 + 1), x) + b^2*c^ 
2*e*integrate(1/16*x^5*log(c^2*x^4 + 1)^2/(c^2*x^4 + 1), x) + 12*b^2*c^2*d 
*integrate(1/16*x^4*arctan(c*x^2)^2/(c^2*x^4 + 1), x) + 4*b^2*c^2*e*integr 
ate(1/16*x^5*log(c^2*x^4 + 1)/(c^2*x^4 + 1), x) + b^2*c^2*d*integrate(1/16 
*x^4*log(c^2*x^4 + 1)^2/(c^2*x^4 + 1), x) + 8*b^2*c^2*d*integrate(1/16*x^4 
*log(c^2*x^4 + 1)/(c^2*x^4 + 1), x) + 1/2*a^2*e*x^2 + 1/8*b^2*e*arctan(c*x 
^2)^3/c - 8*b^2*c*e*integrate(1/16*x^3*arctan(c*x^2)/(c^2*x^4 + 1), x) - 1 
6*b^2*c*d*integrate(1/16*x^2*arctan(c*x^2)/(c^2*x^4 + 1), x) - 1/2*(c*(2*s 
qrt(2)*arctan(1/2*sqrt(2)*(2*c*x + sqrt(2)*sqrt(c))/sqrt(c))/c^(3/2) + 2*s 
qrt(2)*arctan(1/2*sqrt(2)*(2*c*x - sqrt(2)*sqrt(c))/sqrt(c))/c^(3/2) - sqr 
t(2)*log(c*x^2 + sqrt(2)*sqrt(c)*x + 1)/c^(3/2) + sqrt(2)*log(c*x^2 - sqrt 
(2)*sqrt(c)*x + 1)/c^(3/2)) - 4*x*arctan(c*x^2))*a*b*d + a^2*d*x + b^2*e*i 
ntegrate(1/16*x*log(c^2*x^4 + 1)^2/(c^2*x^4 + 1), x) + 12*b^2*d*integrate( 
1/16*arctan(c*x^2)^2/(c^2*x^4 + 1), x) + b^2*d*integrate(1/16*log(c^2*x^4 
+ 1)^2/(c^2*x^4 + 1), x) + 1/2*(2*c*x^2*arctan(c*x^2) - log(c^2*x^4 + 1))* 
a*b*e/c + 1/8*(b^2*e*x^2 + 2*b^2*d*x)*arctan(c*x^2)^2 - 1/32*(b^2*e*x^2 + 
2*b^2*d*x)*log(c^2*x^4 + 1)^2
 

Giac [F]

\[ \int (d+e x) \left (a+b \arctan \left (c x^2\right )\right )^2 \, dx=\int { {\left (e x + d\right )} {\left (b \arctan \left (c x^{2}\right ) + a\right )}^{2} \,d x } \] Input:

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

Output:

integrate((e*x + d)*(b*arctan(c*x^2) + a)^2, x)
 

Mupad [F(-1)]

Timed out. \[ \int (d+e x) \left (a+b \arctan \left (c x^2\right )\right )^2 \, dx=\int {\left (a+b\,\mathrm {atan}\left (c\,x^2\right )\right )}^2\,\left (d+e\,x\right ) \,d x \] Input:

int((a + b*atan(c*x^2))^2*(d + e*x),x)
 

Output:

int((a + b*atan(c*x^2))^2*(d + e*x), x)
 

Reduce [F]

\[ \int (d+e x) \left (a+b \arctan \left (c x^2\right )\right )^2 \, dx=\frac {4 \sqrt {c}\, \sqrt {2}\, \mathit {atan} \left (\frac {\sqrt {c}\, \sqrt {2}-2 c x}{\sqrt {c}\, \sqrt {2}}\right ) a b d +2 \sqrt {c}\, \sqrt {2}\, \mathit {atan} \left (c \,x^{2}\right ) a b d +4 \mathit {atan} \left (c \,x^{2}\right ) a b c d x +2 \mathit {atan} \left (c \,x^{2}\right ) a b c e \,x^{2}-\sqrt {c}\, \sqrt {2}\, \mathrm {log}\left (-\sqrt {c}\, \sqrt {2}\, x +c \,x^{2}+1\right ) a b d +\sqrt {c}\, \sqrt {2}\, \mathrm {log}\left (\sqrt {c}\, \sqrt {2}\, x +c \,x^{2}+1\right ) a b d +2 \left (\int \mathit {atan} \left (c \,x^{2}\right )^{2}d x \right ) b^{2} c d +2 \left (\int \mathit {atan} \left (c \,x^{2}\right )^{2} x d x \right ) b^{2} c e -\mathrm {log}\left (-\sqrt {c}\, \sqrt {2}\, x +c \,x^{2}+1\right ) a b e -\mathrm {log}\left (\sqrt {c}\, \sqrt {2}\, x +c \,x^{2}+1\right ) a b e +2 a^{2} c d x +a^{2} c e \,x^{2}}{2 c} \] Input:

int((e*x+d)*(a+b*atan(c*x^2))^2,x)
 

Output:

(4*sqrt(c)*sqrt(2)*atan((sqrt(c)*sqrt(2) - 2*c*x)/(sqrt(c)*sqrt(2)))*a*b*d 
 + 2*sqrt(c)*sqrt(2)*atan(c*x**2)*a*b*d + 4*atan(c*x**2)*a*b*c*d*x + 2*ata 
n(c*x**2)*a*b*c*e*x**2 - sqrt(c)*sqrt(2)*log( - sqrt(c)*sqrt(2)*x + c*x**2 
 + 1)*a*b*d + sqrt(c)*sqrt(2)*log(sqrt(c)*sqrt(2)*x + c*x**2 + 1)*a*b*d + 
2*int(atan(c*x**2)**2,x)*b**2*c*d + 2*int(atan(c*x**2)**2*x,x)*b**2*c*e - 
log( - sqrt(c)*sqrt(2)*x + c*x**2 + 1)*a*b*e - log(sqrt(c)*sqrt(2)*x + c*x 
**2 + 1)*a*b*e + 2*a**2*c*d*x + a**2*c*e*x**2)/(2*c)