\(\int \frac {(a+b x)^{5/2} (A+B x)}{x^{5/2}} \, dx\) [310]

Optimal result

Integrand size = 20, antiderivative size = 135 \[ \int \frac {(a+b x)^{5/2} (A+B x)}{x^{5/2}} \, dx=\frac {5}{4} b (4 A b+3 a B) \sqrt {x} \sqrt {a+b x}-\frac {2 (5 A b+3 a B) (a+b x)^{3/2}}{3 \sqrt {x}}+\frac {1}{2} b B \sqrt {x} (a+b x)^{3/2}-\frac {2 A (a+b x)^{5/2}}{3 x^{3/2}}+\frac {5}{4} a \sqrt {b} (4 A b+3 a B) \text {arctanh}\left (\frac {\sqrt {b} \sqrt {x}}{\sqrt {a+b x}}\right ) \] Output:

5/4*b*(4*A*b+3*B*a)*x^(1/2)*(b*x+a)^(1/2)-2/3*(5*A*b+3*B*a)*(b*x+a)^(3/2)/ 
x^(1/2)+1/2*b*B*x^(1/2)*(b*x+a)^(3/2)-2/3*A*(b*x+a)^(5/2)/x^(3/2)+5/4*a*b^ 
(1/2)*(4*A*b+3*B*a)*arctanh(b^(1/2)*x^(1/2)/(b*x+a)^(1/2))
 

Mathematica [A] (verified)

Time = 0.34 (sec) , antiderivative size = 108, normalized size of antiderivative = 0.80 \[ \int \frac {(a+b x)^{5/2} (A+B x)}{x^{5/2}} \, dx=\frac {\sqrt {a+b x} \left (6 b^2 x^2 (2 A+B x)-8 a^2 (A+3 B x)+a b x (-56 A+27 B x)\right )}{12 x^{3/2}}+\frac {5}{2} a \sqrt {b} (4 A b+3 a B) \text {arctanh}\left (\frac {\sqrt {b} \sqrt {x}}{-\sqrt {a}+\sqrt {a+b x}}\right ) \] Input:

Integrate[((a + b*x)^(5/2)*(A + B*x))/x^(5/2),x]
 

Output:

(Sqrt[a + b*x]*(6*b^2*x^2*(2*A + B*x) - 8*a^2*(A + 3*B*x) + a*b*x*(-56*A + 
 27*B*x)))/(12*x^(3/2)) + (5*a*Sqrt[b]*(4*A*b + 3*a*B)*ArcTanh[(Sqrt[b]*Sq 
rt[x])/(-Sqrt[a] + Sqrt[a + b*x])])/2
 

Rubi [A] (verified)

Time = 0.21 (sec) , antiderivative size = 127, normalized size of antiderivative = 0.94, number of steps used = 7, number of rules used = 6, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.300, Rules used = {87, 57, 60, 60, 65, 219}

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 \frac {(a+b x)^{5/2} (A+B x)}{x^{5/2}} \, dx\)

\(\Big \downarrow \) 87

\(\displaystyle \frac {(3 a B+4 A b) \int \frac {(a+b x)^{5/2}}{x^{3/2}}dx}{3 a}-\frac {2 A (a+b x)^{7/2}}{3 a x^{3/2}}\)

\(\Big \downarrow \) 57

\(\displaystyle \frac {(3 a B+4 A b) \left (5 b \int \frac {(a+b x)^{3/2}}{\sqrt {x}}dx-\frac {2 (a+b x)^{5/2}}{\sqrt {x}}\right )}{3 a}-\frac {2 A (a+b x)^{7/2}}{3 a x^{3/2}}\)

\(\Big \downarrow \) 60

\(\displaystyle \frac {(3 a B+4 A b) \left (5 b \left (\frac {3}{4} a \int \frac {\sqrt {a+b x}}{\sqrt {x}}dx+\frac {1}{2} \sqrt {x} (a+b x)^{3/2}\right )-\frac {2 (a+b x)^{5/2}}{\sqrt {x}}\right )}{3 a}-\frac {2 A (a+b x)^{7/2}}{3 a x^{3/2}}\)

\(\Big \downarrow \) 60

\(\displaystyle \frac {(3 a B+4 A b) \left (5 b \left (\frac {3}{4} a \left (\frac {1}{2} a \int \frac {1}{\sqrt {x} \sqrt {a+b x}}dx+\sqrt {x} \sqrt {a+b x}\right )+\frac {1}{2} \sqrt {x} (a+b x)^{3/2}\right )-\frac {2 (a+b x)^{5/2}}{\sqrt {x}}\right )}{3 a}-\frac {2 A (a+b x)^{7/2}}{3 a x^{3/2}}\)

\(\Big \downarrow \) 65

\(\displaystyle \frac {(3 a B+4 A b) \left (5 b \left (\frac {3}{4} a \left (a \int \frac {1}{1-\frac {b x}{a+b x}}d\frac {\sqrt {x}}{\sqrt {a+b x}}+\sqrt {x} \sqrt {a+b x}\right )+\frac {1}{2} \sqrt {x} (a+b x)^{3/2}\right )-\frac {2 (a+b x)^{5/2}}{\sqrt {x}}\right )}{3 a}-\frac {2 A (a+b x)^{7/2}}{3 a x^{3/2}}\)

\(\Big \downarrow \) 219

\(\displaystyle \frac {(3 a B+4 A b) \left (5 b \left (\frac {3}{4} a \left (\frac {a \text {arctanh}\left (\frac {\sqrt {b} \sqrt {x}}{\sqrt {a+b x}}\right )}{\sqrt {b}}+\sqrt {x} \sqrt {a+b x}\right )+\frac {1}{2} \sqrt {x} (a+b x)^{3/2}\right )-\frac {2 (a+b x)^{5/2}}{\sqrt {x}}\right )}{3 a}-\frac {2 A (a+b x)^{7/2}}{3 a x^{3/2}}\)

Input:

Int[((a + b*x)^(5/2)*(A + B*x))/x^(5/2),x]
 

Output:

(-2*A*(a + b*x)^(7/2))/(3*a*x^(3/2)) + ((4*A*b + 3*a*B)*((-2*(a + b*x)^(5/ 
2))/Sqrt[x] + 5*b*((Sqrt[x]*(a + b*x)^(3/2))/2 + (3*a*(Sqrt[x]*Sqrt[a + b* 
x] + (a*ArcTanh[(Sqrt[b]*Sqrt[x])/Sqrt[a + b*x]])/Sqrt[b]))/4)))/(3*a)
 

Defintions of rubi rules used

rule 57

Int[((a_.) + (b_.)*(x_))^(m_)*((c_.) + (d_.)*(x_))^(n_), x_Symbol] :> Simp[ 
(a + b*x)^(m + 1)*((c + d*x)^n/(b*(m + 1))), x] - Simp[d*(n/(b*(m + 1))) 
Int[(a + b*x)^(m + 1)*(c + d*x)^(n - 1), x], x] /; FreeQ[{a, b, c, d}, x] & 
& GtQ[n, 0] && LtQ[m, -1] &&  !(IntegerQ[n] &&  !IntegerQ[m]) &&  !(ILeQ[m 
+ n + 2, 0] && (FractionQ[m] || GeQ[2*n + m + 1, 0])) && IntLinearQ[a, b, c 
, d, m, n, x]
 

rule 60

Int[((a_.) + (b_.)*(x_))^(m_)*((c_.) + (d_.)*(x_))^(n_), x_Symbol] :> Simp[ 
(a + b*x)^(m + 1)*((c + d*x)^n/(b*(m + n + 1))), x] + Simp[n*((b*c - a*d)/( 
b*(m + n + 1)))   Int[(a + b*x)^m*(c + d*x)^(n - 1), x], x] /; FreeQ[{a, b, 
 c, d}, x] && GtQ[n, 0] && NeQ[m + n + 1, 0] &&  !(IGtQ[m, 0] && ( !Integer 
Q[n] || (GtQ[m, 0] && LtQ[m - n, 0]))) &&  !ILtQ[m + n + 2, 0] && IntLinear 
Q[a, b, c, d, m, n, x]
 

rule 65

Int[1/(Sqrt[(b_.)*(x_)]*Sqrt[(c_) + (d_.)*(x_)]), x_Symbol] :> Simp[2   Sub 
st[Int[1/(b - d*x^2), x], x, Sqrt[b*x]/Sqrt[c + d*x]], x] /; FreeQ[{b, c, d 
}, x] &&  !GtQ[c, 0]
 

rule 87

Int[((a_.) + (b_.)*(x_))*((c_.) + (d_.)*(x_))^(n_.)*((e_.) + (f_.)*(x_))^(p 
_.), x_] :> Simp[(-(b*e - a*f))*(c + d*x)^(n + 1)*((e + f*x)^(p + 1)/(f*(p 
+ 1)*(c*f - d*e))), x] - Simp[(a*d*f*(n + p + 2) - b*(d*e*(n + 1) + c*f*(p 
+ 1)))/(f*(p + 1)*(c*f - d*e))   Int[(c + d*x)^n*(e + f*x)^(p + 1), x], x] 
/; FreeQ[{a, b, c, d, e, f, n}, x] && LtQ[p, -1] && ( !LtQ[n, -1] || Intege 
rQ[p] ||  !(IntegerQ[n] ||  !(EqQ[e, 0] ||  !(EqQ[c, 0] || LtQ[p, n]))))
 

rule 219

Int[((a_) + (b_.)*(x_)^2)^(-1), x_Symbol] :> Simp[(1/(Rt[a, 2]*Rt[-b, 2]))* 
ArcTanh[Rt[-b, 2]*(x/Rt[a, 2])], x] /; FreeQ[{a, b}, x] && NegQ[a/b] && (Gt 
Q[a, 0] || LtQ[b, 0])
 

Maple [A] (verified)

Time = 0.15 (sec) , antiderivative size = 118, normalized size of antiderivative = 0.87

method	result	size
risch	\(-\frac {\sqrt {b x +a}\, \left (-6 B \,b^{2} x^{3}-12 A \,b^{2} x^{2}-27 B a b \,x^{2}+56 a A b x +24 B \,a^{2} x +8 a^{2} A \right )}{12 x^{\frac {3}{2}}}+\frac {5 a \sqrt {b}\, \left (4 A b +3 B a \right ) \ln \left (\frac {\frac {a}{2}+b x}{\sqrt {b}}+\sqrt {b \,x^{2}+a x}\right ) \sqrt {x \left (b x +a \right )}}{8 \sqrt {x}\, \sqrt {b x +a}}\)	\(118\)
default	\(\frac {\sqrt {b x +a}\, \left (12 B \sqrt {x \left (b x +a \right )}\, b^{\frac {5}{2}} x^{3}+60 A \,b^{2} \ln \left (\frac {2 \sqrt {x \left (b x +a \right )}\, \sqrt {b}+2 b x +a}{2 \sqrt {b}}\right ) a \,x^{2}+24 A \sqrt {x \left (b x +a \right )}\, b^{\frac {5}{2}} x^{2}+45 B b \ln \left (\frac {2 \sqrt {x \left (b x +a \right )}\, \sqrt {b}+2 b x +a}{2 \sqrt {b}}\right ) a^{2} x^{2}+54 B \sqrt {x \left (b x +a \right )}\, b^{\frac {3}{2}} a \,x^{2}-112 A a \,b^{\frac {3}{2}} x \sqrt {x \left (b x +a \right )}-48 B \sqrt {x \left (b x +a \right )}\, \sqrt {b}\, a^{2} x -16 A \,a^{2} \sqrt {x \left (b x +a \right )}\, \sqrt {b}\right )}{24 x^{\frac {3}{2}} \sqrt {x \left (b x +a \right )}\, \sqrt {b}}\)	\(207\)

Input:

int((b*x+a)^(5/2)*(B*x+A)/x^(5/2),x,method=_RETURNVERBOSE)
 

Output:

-1/12*(b*x+a)^(1/2)*(-6*B*b^2*x^3-12*A*b^2*x^2-27*B*a*b*x^2+56*A*a*b*x+24* 
B*a^2*x+8*A*a^2)/x^(3/2)+5/8*a*b^(1/2)*(4*A*b+3*B*a)*ln((1/2*a+b*x)/b^(1/2 
)+(b*x^2+a*x)^(1/2))*(x*(b*x+a))^(1/2)/x^(1/2)/(b*x+a)^(1/2)
 

Fricas [A] (verification not implemented)

Time = 0.09 (sec) , antiderivative size = 214, normalized size of antiderivative = 1.59 \[ \int \frac {(a+b x)^{5/2} (A+B x)}{x^{5/2}} \, dx=\left [\frac {15 \, {\left (3 \, B a^{2} + 4 \, A a b\right )} \sqrt {b} x^{2} \log \left (2 \, b x + 2 \, \sqrt {b x + a} \sqrt {b} \sqrt {x} + a\right ) + 2 \, {\left (6 \, B b^{2} x^{3} - 8 \, A a^{2} + 3 \, {\left (9 \, B a b + 4 \, A b^{2}\right )} x^{2} - 8 \, {\left (3 \, B a^{2} + 7 \, A a b\right )} x\right )} \sqrt {b x + a} \sqrt {x}}{24 \, x^{2}}, -\frac {15 \, {\left (3 \, B a^{2} + 4 \, A a b\right )} \sqrt {-b} x^{2} \arctan \left (\frac {\sqrt {-b} \sqrt {x}}{\sqrt {b x + a}}\right ) - {\left (6 \, B b^{2} x^{3} - 8 \, A a^{2} + 3 \, {\left (9 \, B a b + 4 \, A b^{2}\right )} x^{2} - 8 \, {\left (3 \, B a^{2} + 7 \, A a b\right )} x\right )} \sqrt {b x + a} \sqrt {x}}{12 \, x^{2}}\right ] \] Input:

integrate((b*x+a)^(5/2)*(B*x+A)/x^(5/2),x, algorithm="fricas")
 

Output:

[1/24*(15*(3*B*a^2 + 4*A*a*b)*sqrt(b)*x^2*log(2*b*x + 2*sqrt(b*x + a)*sqrt 
(b)*sqrt(x) + a) + 2*(6*B*b^2*x^3 - 8*A*a^2 + 3*(9*B*a*b + 4*A*b^2)*x^2 - 
8*(3*B*a^2 + 7*A*a*b)*x)*sqrt(b*x + a)*sqrt(x))/x^2, -1/12*(15*(3*B*a^2 + 
4*A*a*b)*sqrt(-b)*x^2*arctan(sqrt(-b)*sqrt(x)/sqrt(b*x + a)) - (6*B*b^2*x^ 
3 - 8*A*a^2 + 3*(9*B*a*b + 4*A*b^2)*x^2 - 8*(3*B*a^2 + 7*A*a*b)*x)*sqrt(b* 
x + a)*sqrt(x))/x^2]
 

Sympy [A] (verification not implemented)

Time = 6.13 (sec) , antiderivative size = 364, normalized size of antiderivative = 2.70 \[ \int \frac {(a+b x)^{5/2} (A+B x)}{x^{5/2}} \, dx=- \frac {4 A a^{\frac {3}{2}} b}{\sqrt {x} \sqrt {1 + \frac {b x}{a}}} + A \sqrt {a} b^{2} \sqrt {x} \sqrt {1 + \frac {b x}{a}} - \frac {4 A \sqrt {a} b^{2} \sqrt {x}}{\sqrt {1 + \frac {b x}{a}}} - \frac {2 A a^{2} \sqrt {b} \sqrt {\frac {a}{b x} + 1}}{3 x} - \frac {2 A a b^{\frac {3}{2}} \sqrt {\frac {a}{b x} + 1}}{3} + 5 A a b^{\frac {3}{2}} \operatorname {asinh}{\left (\frac {\sqrt {b} \sqrt {x}}{\sqrt {a}} \right )} - \frac {2 B a^{\frac {5}{2}}}{\sqrt {x} \sqrt {1 + \frac {b x}{a}}} + 2 B a^{\frac {3}{2}} b \sqrt {x} \sqrt {1 + \frac {b x}{a}} - \frac {2 B a^{\frac {3}{2}} b \sqrt {x}}{\sqrt {1 + \frac {b x}{a}}} + 4 B a^{2} \sqrt {b} \operatorname {asinh}{\left (\frac {\sqrt {b} \sqrt {x}}{\sqrt {a}} \right )} + 2 B b^{2} \left (\begin {cases} - \frac {a^{2} \left (\begin {cases} \frac {\log {\left (2 \sqrt {b} \sqrt {a + b x} + 2 b \sqrt {x} \right )}}{\sqrt {b}} & \text {for}\: a \neq 0 \\\frac {\sqrt {x} \log {\left (\sqrt {x} \right )}}{\sqrt {b x}} & \text {otherwise} \end {cases}\right )}{8 b} + \frac {a \sqrt {x} \sqrt {a + b x}}{8 b} + \frac {x^{\frac {3}{2}} \sqrt {a + b x}}{4} & \text {for}\: b \neq 0 \\\frac {\sqrt {a} x^{\frac {3}{2}}}{3} & \text {otherwise} \end {cases}\right ) \] Input:

integrate((b*x+a)**(5/2)*(B*x+A)/x**(5/2),x)
 

Output:

-4*A*a**(3/2)*b/(sqrt(x)*sqrt(1 + b*x/a)) + A*sqrt(a)*b**2*sqrt(x)*sqrt(1 
+ b*x/a) - 4*A*sqrt(a)*b**2*sqrt(x)/sqrt(1 + b*x/a) - 2*A*a**2*sqrt(b)*sqr 
t(a/(b*x) + 1)/(3*x) - 2*A*a*b**(3/2)*sqrt(a/(b*x) + 1)/3 + 5*A*a*b**(3/2) 
*asinh(sqrt(b)*sqrt(x)/sqrt(a)) - 2*B*a**(5/2)/(sqrt(x)*sqrt(1 + b*x/a)) + 
 2*B*a**(3/2)*b*sqrt(x)*sqrt(1 + b*x/a) - 2*B*a**(3/2)*b*sqrt(x)/sqrt(1 + 
b*x/a) + 4*B*a**2*sqrt(b)*asinh(sqrt(b)*sqrt(x)/sqrt(a)) + 2*B*b**2*Piecew 
ise((-a**2*Piecewise((log(2*sqrt(b)*sqrt(a + b*x) + 2*b*sqrt(x))/sqrt(b), 
Ne(a, 0)), (sqrt(x)*log(sqrt(x))/sqrt(b*x), True))/(8*b) + a*sqrt(x)*sqrt( 
a + b*x)/(8*b) + x**(3/2)*sqrt(a + b*x)/4, Ne(b, 0)), (sqrt(a)*x**(3/2)/3, 
 True))
 

Maxima [A] (verification not implemented)

Time = 0.05 (sec) , antiderivative size = 191, normalized size of antiderivative = 1.41 \[ \int \frac {(a+b x)^{5/2} (A+B x)}{x^{5/2}} \, dx=\frac {15}{8} \, B a^{2} \sqrt {b} \log \left (2 \, b x + a + 2 \, \sqrt {b x^{2} + a x} \sqrt {b}\right ) + \frac {5}{2} \, A a b^{\frac {3}{2}} \log \left (2 \, b x + a + 2 \, \sqrt {b x^{2} + a x} \sqrt {b}\right ) - \frac {15 \, \sqrt {b x^{2} + a x} B a^{2}}{4 \, x} - \frac {35 \, \sqrt {b x^{2} + a x} A a b}{6 \, x} + \frac {5 \, {\left (b x^{2} + a x\right )}^{\frac {3}{2}} B a}{4 \, x^{2}} - \frac {5 \, \sqrt {b x^{2} + a x} A a^{2}}{6 \, x^{2}} + \frac {{\left (b x^{2} + a x\right )}^{\frac {5}{2}} B}{2 \, x^{3}} - \frac {5 \, {\left (b x^{2} + a x\right )}^{\frac {3}{2}} A a}{6 \, x^{3}} + \frac {{\left (b x^{2} + a x\right )}^{\frac {5}{2}} A}{x^{4}} \] Input:

integrate((b*x+a)^(5/2)*(B*x+A)/x^(5/2),x, algorithm="maxima")
 

Output:

15/8*B*a^2*sqrt(b)*log(2*b*x + a + 2*sqrt(b*x^2 + a*x)*sqrt(b)) + 5/2*A*a* 
b^(3/2)*log(2*b*x + a + 2*sqrt(b*x^2 + a*x)*sqrt(b)) - 15/4*sqrt(b*x^2 + a 
*x)*B*a^2/x - 35/6*sqrt(b*x^2 + a*x)*A*a*b/x + 5/4*(b*x^2 + a*x)^(3/2)*B*a 
/x^2 - 5/6*sqrt(b*x^2 + a*x)*A*a^2/x^2 + 1/2*(b*x^2 + a*x)^(5/2)*B/x^3 - 5 
/6*(b*x^2 + a*x)^(3/2)*A*a/x^3 + (b*x^2 + a*x)^(5/2)*A/x^4
 

Giac [A] (verification not implemented)

Time = 76.46 (sec) , antiderivative size = 175, normalized size of antiderivative = 1.30 \[ \int \frac {(a+b x)^{5/2} (A+B x)}{x^{5/2}} \, dx=-\frac {b^{3} {\left (\frac {15 \, {\left (3 \, B a^{2} + 4 \, A a b\right )} \log \left ({\left | -\sqrt {b x + a} \sqrt {b} + \sqrt {{\left (b x + a\right )} b - a b} \right |}\right )}{b^{\frac {3}{2}}} - \frac {{\left ({\left (3 \, {\left (2 \, {\left (b x + a\right )} B + \frac {3 \, B a^{2} b + 4 \, A a b^{2}}{a b}\right )} {\left (b x + a\right )} - \frac {20 \, {\left (3 \, B a^{3} b + 4 \, A a^{2} b^{2}\right )}}{a b}\right )} {\left (b x + a\right )} + \frac {15 \, {\left (3 \, B a^{4} b + 4 \, A a^{3} b^{2}\right )}}{a b}\right )} \sqrt {b x + a}}{{\left ({\left (b x + a\right )} b - a b\right )}^{\frac {3}{2}}}\right )}}{12 \, {\left | b \right |}} \] Input:

integrate((b*x+a)^(5/2)*(B*x+A)/x^(5/2),x, algorithm="giac")
 

Output:

-1/12*b^3*(15*(3*B*a^2 + 4*A*a*b)*log(abs(-sqrt(b*x + a)*sqrt(b) + sqrt((b 
*x + a)*b - a*b)))/b^(3/2) - ((3*(2*(b*x + a)*B + (3*B*a^2*b + 4*A*a*b^2)/ 
(a*b))*(b*x + a) - 20*(3*B*a^3*b + 4*A*a^2*b^2)/(a*b))*(b*x + a) + 15*(3*B 
*a^4*b + 4*A*a^3*b^2)/(a*b))*sqrt(b*x + a)/((b*x + a)*b - a*b)^(3/2))/abs( 
b)
 

Mupad [F(-1)]

Timed out. \[ \int \frac {(a+b x)^{5/2} (A+B x)}{x^{5/2}} \, dx=\int \frac {\left (A+B\,x\right )\,{\left (a+b\,x\right )}^{5/2}}{x^{5/2}} \,d x \] Input:

int(((A + B*x)*(a + b*x)^(5/2))/x^(5/2),x)
 

Output:

int(((A + B*x)*(a + b*x)^(5/2))/x^(5/2), x)
 

Reduce [B] (verification not implemented)

Time = 0.16 (sec) , antiderivative size = 107, normalized size of antiderivative = 0.79 \[ \int \frac {(a+b x)^{5/2} (A+B x)}{x^{5/2}} \, dx=\frac {-64 \sqrt {x}\, \sqrt {b x +a}\, a^{3}-640 \sqrt {x}\, \sqrt {b x +a}\, a^{2} b x +312 \sqrt {x}\, \sqrt {b x +a}\, a \,b^{2} x^{2}+48 \sqrt {x}\, \sqrt {b x +a}\, b^{3} x^{3}+840 \sqrt {b}\, \mathrm {log}\left (\frac {\sqrt {b x +a}+\sqrt {x}\, \sqrt {b}}{\sqrt {a}}\right ) a^{2} b \,x^{2}+175 \sqrt {b}\, a^{2} b \,x^{2}}{96 x^{2}} \] Input:

int((b*x+a)^(5/2)*(B*x+A)/x^(5/2),x)
 

Output:

( - 64*sqrt(x)*sqrt(a + b*x)*a**3 - 640*sqrt(x)*sqrt(a + b*x)*a**2*b*x + 3 
12*sqrt(x)*sqrt(a + b*x)*a*b**2*x**2 + 48*sqrt(x)*sqrt(a + b*x)*b**3*x**3 
+ 840*sqrt(b)*log((sqrt(a + b*x) + sqrt(x)*sqrt(b))/sqrt(a))*a**2*b*x**2 + 
 175*sqrt(b)*a**2*b*x**2)/(96*x**2)