Integrand size = 26, antiderivative size = 235 \[ \int \frac {x^3}{\sqrt {a^2+2 a b x^3+b^2 x^6}} \, dx=\frac {x \left (a+b x^3\right )}{b \sqrt {a^2+2 a b x^3+b^2 x^6}}+\frac {\sqrt [3]{a} \left (a+b x^3\right ) \arctan \left (\frac {\sqrt [3]{a}-2 \sqrt [3]{b} x}{\sqrt {3} \sqrt [3]{a}}\right )}{\sqrt {3} b^{4/3} \sqrt {a^2+2 a b x^3+b^2 x^6}}-\frac {\sqrt [3]{a} \left (a+b x^3\right ) \log \left (\sqrt [3]{a}+\sqrt [3]{b} x\right )}{3 b^{4/3} \sqrt {a^2+2 a b x^3+b^2 x^6}}+\frac {\sqrt [3]{a} \left (a+b x^3\right ) \log \left (a^{2/3}-\sqrt [3]{a} \sqrt [3]{b} x+b^{2/3} x^2\right )}{6 b^{4/3} \sqrt {a^2+2 a b x^3+b^2 x^6}} \] Output:
x*(b*x^3+a)/b/((b*x^3+a)^2)^(1/2)+1/3*a^(1/3)*(b*x^3+a)*arctan(1/3*(a^(1/3 )-2*b^(1/3)*x)*3^(1/2)/a^(1/3))*3^(1/2)/b^(4/3)/((b*x^3+a)^2)^(1/2)-1/3*a^ (1/3)*(b*x^3+a)*ln(a^(1/3)+b^(1/3)*x)/b^(4/3)/((b*x^3+a)^2)^(1/2)+1/6*a^(1 /3)*(b*x^3+a)*ln(a^(2/3)-a^(1/3)*b^(1/3)*x+b^(2/3)*x^2)/b^(4/3)/((b*x^3+a) ^2)^(1/2)
Time = 1.05 (sec) , antiderivative size = 128, normalized size of antiderivative = 0.54 \[ \int \frac {x^3}{\sqrt {a^2+2 a b x^3+b^2 x^6}} \, dx=\frac {\left (a+b x^3\right ) \left (6 \sqrt [3]{b} x+2 \sqrt {3} \sqrt [3]{a} \arctan \left (\frac {1-\frac {2 \sqrt [3]{b} x}{\sqrt [3]{a}}}{\sqrt {3}}\right )-2 \sqrt [3]{a} \log \left (\sqrt [3]{a}+\sqrt [3]{b} x\right )+\sqrt [3]{a} \log \left (a^{2/3}-\sqrt [3]{a} \sqrt [3]{b} x+b^{2/3} x^2\right )\right )}{6 b^{4/3} \sqrt {\left (a+b x^3\right )^2}} \] Input:
Integrate[x^3/Sqrt[a^2 + 2*a*b*x^3 + b^2*x^6],x]
Output:
((a + b*x^3)*(6*b^(1/3)*x + 2*Sqrt[3]*a^(1/3)*ArcTan[(1 - (2*b^(1/3)*x)/a^ (1/3))/Sqrt[3]] - 2*a^(1/3)*Log[a^(1/3) + b^(1/3)*x] + a^(1/3)*Log[a^(2/3) - a^(1/3)*b^(1/3)*x + b^(2/3)*x^2]))/(6*b^(4/3)*Sqrt[(a + b*x^3)^2])
Time = 0.33 (sec) , antiderivative size = 153, normalized size of antiderivative = 0.65, number of steps used = 12, number of rules used = 11, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.423, Rules used = {1384, 27, 843, 750, 16, 1142, 25, 27, 1082, 217, 1103}
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 {x^3}{\sqrt {a^2+2 a b x^3+b^2 x^6}} \, dx\) |
| \(\Big \downarrow \) 1384 |
| \(\displaystyle \frac {b \left (a+b x^3\right ) \int \frac {x^3}{b \left (b x^3+a\right )}dx}{\sqrt {a^2+2 a b x^3+b^2 x^6}}\) |
| \(\Big \downarrow \) 27 |
| \(\displaystyle \frac {\left (a+b x^3\right ) \int \frac {x^3}{b x^3+a}dx}{\sqrt {a^2+2 a b x^3+b^2 x^6}}\) |
| \(\Big \downarrow \) 843 |
| \(\displaystyle \frac {\left (a+b x^3\right ) \left (\frac {x}{b}-\frac {a \int \frac {1}{b x^3+a}dx}{b}\right )}{\sqrt {a^2+2 a b x^3+b^2 x^6}}\) |
| \(\Big \downarrow \) 750 |
| \(\displaystyle \frac {\left (a+b x^3\right ) \left (\frac {x}{b}-\frac {a \left (\frac {\int \frac {2 \sqrt [3]{a}-\sqrt [3]{b} x}{b^{2/3} x^2-\sqrt [3]{a} \sqrt [3]{b} x+a^{2/3}}dx}{3 a^{2/3}}+\frac {\int \frac {1}{\sqrt [3]{b} x+\sqrt [3]{a}}dx}{3 a^{2/3}}\right )}{b}\right )}{\sqrt {a^2+2 a b x^3+b^2 x^6}}\) |
| \(\Big \downarrow \) 16 |
| \(\displaystyle \frac {\left (a+b x^3\right ) \left (\frac {x}{b}-\frac {a \left (\frac {\int \frac {2 \sqrt [3]{a}-\sqrt [3]{b} x}{b^{2/3} x^2-\sqrt [3]{a} \sqrt [3]{b} x+a^{2/3}}dx}{3 a^{2/3}}+\frac {\log \left (\sqrt [3]{a}+\sqrt [3]{b} x\right )}{3 a^{2/3} \sqrt [3]{b}}\right )}{b}\right )}{\sqrt {a^2+2 a b x^3+b^2 x^6}}\) |
| \(\Big \downarrow \) 1142 |
| \(\displaystyle \frac {\left (a+b x^3\right ) \left (\frac {x}{b}-\frac {a \left (\frac {\frac {3}{2} \sqrt [3]{a} \int \frac {1}{b^{2/3} x^2-\sqrt [3]{a} \sqrt [3]{b} x+a^{2/3}}dx-\frac {\int -\frac {\sqrt [3]{b} \left (\sqrt [3]{a}-2 \sqrt [3]{b} x\right )}{b^{2/3} x^2-\sqrt [3]{a} \sqrt [3]{b} x+a^{2/3}}dx}{2 \sqrt [3]{b}}}{3 a^{2/3}}+\frac {\log \left (\sqrt [3]{a}+\sqrt [3]{b} x\right )}{3 a^{2/3} \sqrt [3]{b}}\right )}{b}\right )}{\sqrt {a^2+2 a b x^3+b^2 x^6}}\) |
| \(\Big \downarrow \) 25 |
| \(\displaystyle \frac {\left (a+b x^3\right ) \left (\frac {x}{b}-\frac {a \left (\frac {\frac {3}{2} \sqrt [3]{a} \int \frac {1}{b^{2/3} x^2-\sqrt [3]{a} \sqrt [3]{b} x+a^{2/3}}dx+\frac {\int \frac {\sqrt [3]{b} \left (\sqrt [3]{a}-2 \sqrt [3]{b} x\right )}{b^{2/3} x^2-\sqrt [3]{a} \sqrt [3]{b} x+a^{2/3}}dx}{2 \sqrt [3]{b}}}{3 a^{2/3}}+\frac {\log \left (\sqrt [3]{a}+\sqrt [3]{b} x\right )}{3 a^{2/3} \sqrt [3]{b}}\right )}{b}\right )}{\sqrt {a^2+2 a b x^3+b^2 x^6}}\) |
| \(\Big \downarrow \) 27 |
| \(\displaystyle \frac {\left (a+b x^3\right ) \left (\frac {x}{b}-\frac {a \left (\frac {\frac {3}{2} \sqrt [3]{a} \int \frac {1}{b^{2/3} x^2-\sqrt [3]{a} \sqrt [3]{b} x+a^{2/3}}dx+\frac {1}{2} \int \frac {\sqrt [3]{a}-2 \sqrt [3]{b} x}{b^{2/3} x^2-\sqrt [3]{a} \sqrt [3]{b} x+a^{2/3}}dx}{3 a^{2/3}}+\frac {\log \left (\sqrt [3]{a}+\sqrt [3]{b} x\right )}{3 a^{2/3} \sqrt [3]{b}}\right )}{b}\right )}{\sqrt {a^2+2 a b x^3+b^2 x^6}}\) |
| \(\Big \downarrow \) 1082 |
| \(\displaystyle \frac {\left (a+b x^3\right ) \left (\frac {x}{b}-\frac {a \left (\frac {\frac {1}{2} \int \frac {\sqrt [3]{a}-2 \sqrt [3]{b} x}{b^{2/3} x^2-\sqrt [3]{a} \sqrt [3]{b} x+a^{2/3}}dx+\frac {3 \int \frac {1}{-\left (1-\frac {2 \sqrt [3]{b} x}{\sqrt [3]{a}}\right )^2-3}d\left (1-\frac {2 \sqrt [3]{b} x}{\sqrt [3]{a}}\right )}{\sqrt [3]{b}}}{3 a^{2/3}}+\frac {\log \left (\sqrt [3]{a}+\sqrt [3]{b} x\right )}{3 a^{2/3} \sqrt [3]{b}}\right )}{b}\right )}{\sqrt {a^2+2 a b x^3+b^2 x^6}}\) |
| \(\Big \downarrow \) 217 |
| \(\displaystyle \frac {\left (a+b x^3\right ) \left (\frac {x}{b}-\frac {a \left (\frac {\frac {1}{2} \int \frac {\sqrt [3]{a}-2 \sqrt [3]{b} x}{b^{2/3} x^2-\sqrt [3]{a} \sqrt [3]{b} x+a^{2/3}}dx-\frac {\sqrt {3} \arctan \left (\frac {1-\frac {2 \sqrt [3]{b} x}{\sqrt [3]{a}}}{\sqrt {3}}\right )}{\sqrt [3]{b}}}{3 a^{2/3}}+\frac {\log \left (\sqrt [3]{a}+\sqrt [3]{b} x\right )}{3 a^{2/3} \sqrt [3]{b}}\right )}{b}\right )}{\sqrt {a^2+2 a b x^3+b^2 x^6}}\) |
| \(\Big \downarrow \) 1103 |
| \(\displaystyle \frac {\left (a+b x^3\right ) \left (\frac {x}{b}-\frac {a \left (\frac {-\frac {\log \left (a^{2/3}-\sqrt [3]{a} \sqrt [3]{b} x+b^{2/3} x^2\right )}{2 \sqrt [3]{b}}-\frac {\sqrt {3} \arctan \left (\frac {1-\frac {2 \sqrt [3]{b} x}{\sqrt [3]{a}}}{\sqrt {3}}\right )}{\sqrt [3]{b}}}{3 a^{2/3}}+\frac {\log \left (\sqrt [3]{a}+\sqrt [3]{b} x\right )}{3 a^{2/3} \sqrt [3]{b}}\right )}{b}\right )}{\sqrt {a^2+2 a b x^3+b^2 x^6}}\) |
Input:
Int[x^3/Sqrt[a^2 + 2*a*b*x^3 + b^2*x^6],x]
Output:
((a + b*x^3)*(x/b - (a*(Log[a^(1/3) + b^(1/3)*x]/(3*a^(2/3)*b^(1/3)) + (-( (Sqrt[3]*ArcTan[(1 - (2*b^(1/3)*x)/a^(1/3))/Sqrt[3]])/b^(1/3)) - Log[a^(2/ 3) - a^(1/3)*b^(1/3)*x + b^(2/3)*x^2]/(2*b^(1/3)))/(3*a^(2/3))))/b))/Sqrt[ a^2 + 2*a*b*x^3 + b^2*x^6]
Int[(c_.)/((a_.) + (b_.)*(x_)), x_Symbol] :> Simp[c*(Log[RemoveContent[a + b*x, x]]/b), x] /; FreeQ[{a, b, c}, x]
Int[(a_)*(Fx_), x_Symbol] :> Simp[a Int[Fx, x], x] /; FreeQ[a, x] && !Ma tchQ[Fx, (b_)*(Gx_) /; FreeQ[b, x]]
Int[((a_) + (b_.)*(x_)^2)^(-1), x_Symbol] :> Simp[(-(Rt[-a, 2]*Rt[-b, 2])^( -1))*ArcTan[Rt[-b, 2]*(x/Rt[-a, 2])], x] /; FreeQ[{a, b}, x] && PosQ[a/b] & & (LtQ[a, 0] || LtQ[b, 0])
Int[((a_) + (b_.)*(x_)^3)^(-1), x_Symbol] :> Simp[1/(3*Rt[a, 3]^2) Int[1/ (Rt[a, 3] + Rt[b, 3]*x), x], x] + Simp[1/(3*Rt[a, 3]^2) Int[(2*Rt[a, 3] - Rt[b, 3]*x)/(Rt[a, 3]^2 - Rt[a, 3]*Rt[b, 3]*x + Rt[b, 3]^2*x^2), x], x] /; FreeQ[{a, b}, x]
Int[((c_.)*(x_))^(m_)*((a_) + (b_.)*(x_)^(n_))^(p_), x_Symbol] :> Simp[c^(n - 1)*(c*x)^(m - n + 1)*((a + b*x^n)^(p + 1)/(b*(m + n*p + 1))), x] - Simp[ a*c^n*((m - n + 1)/(b*(m + n*p + 1))) Int[(c*x)^(m - n)*(a + b*x^n)^p, x] , x] /; FreeQ[{a, b, c, p}, x] && IGtQ[n, 0] && GtQ[m, n - 1] && NeQ[m + n* p + 1, 0] && IntBinomialQ[a, b, c, n, m, p, x]
Int[((a_) + (b_.)*(x_) + (c_.)*(x_)^2)^(-1), x_Symbol] :> With[{q = 1 - 4*S implify[a*(c/b^2)]}, Simp[-2/b Subst[Int[1/(q - x^2), x], x, 1 + 2*c*(x/b )], x] /; RationalQ[q] && (EqQ[q^2, 1] || !RationalQ[b^2 - 4*a*c])] /; Fre eQ[{a, b, c}, x]
Int[((d_) + (e_.)*(x_))/((a_.) + (b_.)*(x_) + (c_.)*(x_)^2), x_Symbol] :> S imp[d*(Log[RemoveContent[a + b*x + c*x^2, x]]/b), x] /; FreeQ[{a, b, c, d, e}, x] && EqQ[2*c*d - b*e, 0]
Int[((d_.) + (e_.)*(x_))/((a_) + (b_.)*(x_) + (c_.)*(x_)^2), x_Symbol] :> S imp[(2*c*d - b*e)/(2*c) Int[1/(a + b*x + c*x^2), x], x] + Simp[e/(2*c) Int[(b + 2*c*x)/(a + b*x + c*x^2), x], x] /; FreeQ[{a, b, c, d, e}, x]
Int[(u_.)*((a_) + (c_.)*(x_)^(n2_.) + (b_.)*(x_)^(n_))^(p_), x_Symbol] :> S imp[(a + b*x^n + c*x^(2*n))^FracPart[p]/(c^IntPart[p]*(b/2 + c*x^n)^(2*Frac Part[p])) Int[u*(b/2 + c*x^n)^(2*p), x], x] /; FreeQ[{a, b, c, n, p}, x] && EqQ[n2, 2*n] && EqQ[b^2 - 4*a*c, 0] && IntegerQ[p - 1/2] && NeQ[u, x^(n - 1)] && NeQ[u, x^(2*n - 1)] && !(EqQ[p, 1/2] && EqQ[u, x^(-2*n - 1)])
Result contains higher order function than in optimal. Order 9 vs. order 3.
Time = 1.81 (sec) , antiderivative size = 74, normalized size of antiderivative = 0.31
| method | result | size |
| risch | \(\frac {x \sqrt {\left (b \,x^{3}+a \right )^{2}}}{\left (b \,x^{3}+a \right ) b}-\frac {\sqrt {\left (b \,x^{3}+a \right )^{2}}\, a \left (\munderset {\textit {\_R} =\operatorname {RootOf}\left (b \,\textit {\_Z}^{3}+a \right )}{\sum }\frac {\ln \left (x -\textit {\_R} \right )}{\textit {\_R}^{2}}\right )}{3 \left (b \,x^{3}+a \right ) b^{2}}\) | \(74\) |
| default | \(\frac {\left (b \,x^{3}+a \right ) \left (6 x b \left (\frac {a}{b}\right )^{\frac {2}{3}}+2 \sqrt {3}\, \arctan \left (\frac {\sqrt {3}\, \left (-2 x +\left (\frac {a}{b}\right )^{\frac {1}{3}}\right )}{3 \left (\frac {a}{b}\right )^{\frac {1}{3}}}\right ) a -2 \ln \left (x +\left (\frac {a}{b}\right )^{\frac {1}{3}}\right ) a +\ln \left (x^{2}-\left (\frac {a}{b}\right )^{\frac {1}{3}} x +\left (\frac {a}{b}\right )^{\frac {2}{3}}\right ) a \right )}{6 \sqrt {\left (b \,x^{3}+a \right )^{2}}\, b^{2} \left (\frac {a}{b}\right )^{\frac {2}{3}}}\) | \(110\) |
Input:
int(x^3/((b*x^3+a)^2)^(1/2),x,method=_RETURNVERBOSE)
Output:
x*((b*x^3+a)^2)^(1/2)/(b*x^3+a)/b-1/3*((b*x^3+a)^2)^(1/2)/(b*x^3+a)/b^2*a* sum(1/_R^2*ln(x-_R),_R=RootOf(_Z^3*b+a))
Time = 0.08 (sec) , antiderivative size = 106, normalized size of antiderivative = 0.45 \[ \int \frac {x^3}{\sqrt {a^2+2 a b x^3+b^2 x^6}} \, dx=\frac {2 \, \sqrt {3} \left (-\frac {a}{b}\right )^{\frac {1}{3}} \arctan \left (\frac {2 \, \sqrt {3} b x \left (-\frac {a}{b}\right )^{\frac {2}{3}} - \sqrt {3} a}{3 \, a}\right ) - \left (-\frac {a}{b}\right )^{\frac {1}{3}} \log \left (x^{2} + x \left (-\frac {a}{b}\right )^{\frac {1}{3}} + \left (-\frac {a}{b}\right )^{\frac {2}{3}}\right ) + 2 \, \left (-\frac {a}{b}\right )^{\frac {1}{3}} \log \left (x - \left (-\frac {a}{b}\right )^{\frac {1}{3}}\right ) + 6 \, x}{6 \, b} \] Input:
integrate(x^3/((b*x^3+a)^2)^(1/2),x, algorithm="fricas")
Output:
1/6*(2*sqrt(3)*(-a/b)^(1/3)*arctan(1/3*(2*sqrt(3)*b*x*(-a/b)^(2/3) - sqrt( 3)*a)/a) - (-a/b)^(1/3)*log(x^2 + x*(-a/b)^(1/3) + (-a/b)^(2/3)) + 2*(-a/b )^(1/3)*log(x - (-a/b)^(1/3)) + 6*x)/b
\[ \int \frac {x^3}{\sqrt {a^2+2 a b x^3+b^2 x^6}} \, dx=\int \frac {x^{3}}{\sqrt {\left (a + b x^{3}\right )^{2}}}\, dx \] Input:
integrate(x**3/((b*x**3+a)**2)**(1/2),x)
Output:
Integral(x**3/sqrt((a + b*x**3)**2), x)
Time = 0.13 (sec) , antiderivative size = 106, normalized size of antiderivative = 0.45 \[ \int \frac {x^3}{\sqrt {a^2+2 a b x^3+b^2 x^6}} \, dx=\frac {x}{b} - \frac {\sqrt {3} a \arctan \left (\frac {\sqrt {3} {\left (2 \, x - \left (\frac {a}{b}\right )^{\frac {1}{3}}\right )}}{3 \, \left (\frac {a}{b}\right )^{\frac {1}{3}}}\right )}{3 \, b^{2} \left (\frac {a}{b}\right )^{\frac {2}{3}}} + \frac {a \log \left (x^{2} - x \left (\frac {a}{b}\right )^{\frac {1}{3}} + \left (\frac {a}{b}\right )^{\frac {2}{3}}\right )}{6 \, b^{2} \left (\frac {a}{b}\right )^{\frac {2}{3}}} - \frac {a \log \left (x + \left (\frac {a}{b}\right )^{\frac {1}{3}}\right )}{3 \, b^{2} \left (\frac {a}{b}\right )^{\frac {2}{3}}} \] Input:
integrate(x^3/((b*x^3+a)^2)^(1/2),x, algorithm="maxima")
Output:
x/b - 1/3*sqrt(3)*a*arctan(1/3*sqrt(3)*(2*x - (a/b)^(1/3))/(a/b)^(1/3))/(b ^2*(a/b)^(2/3)) + 1/6*a*log(x^2 - x*(a/b)^(1/3) + (a/b)^(2/3))/(b^2*(a/b)^ (2/3)) - 1/3*a*log(x + (a/b)^(1/3))/(b^2*(a/b)^(2/3))
Time = 0.13 (sec) , antiderivative size = 143, normalized size of antiderivative = 0.61 \[ \int \frac {x^3}{\sqrt {a^2+2 a b x^3+b^2 x^6}} \, dx=\frac {\left (-\frac {a}{b}\right )^{\frac {1}{3}} \log \left ({\left | x - \left (-\frac {a}{b}\right )^{\frac {1}{3}} \right |}\right ) \mathrm {sgn}\left (b x^{3} + a\right )}{3 \, b} + \frac {x \mathrm {sgn}\left (b x^{3} + a\right )}{b} - \frac {\sqrt {3} \left (-a b^{2}\right )^{\frac {1}{3}} \arctan \left (\frac {\sqrt {3} {\left (2 \, x + \left (-\frac {a}{b}\right )^{\frac {1}{3}}\right )}}{3 \, \left (-\frac {a}{b}\right )^{\frac {1}{3}}}\right ) \mathrm {sgn}\left (b x^{3} + a\right )}{3 \, b^{2}} - \frac {\left (-a b^{2}\right )^{\frac {1}{3}} \log \left (x^{2} + x \left (-\frac {a}{b}\right )^{\frac {1}{3}} + \left (-\frac {a}{b}\right )^{\frac {2}{3}}\right ) \mathrm {sgn}\left (b x^{3} + a\right )}{6 \, b^{2}} \] Input:
integrate(x^3/((b*x^3+a)^2)^(1/2),x, algorithm="giac")
Output:
1/3*(-a/b)^(1/3)*log(abs(x - (-a/b)^(1/3)))*sgn(b*x^3 + a)/b + x*sgn(b*x^3 + a)/b - 1/3*sqrt(3)*(-a*b^2)^(1/3)*arctan(1/3*sqrt(3)*(2*x + (-a/b)^(1/3 ))/(-a/b)^(1/3))*sgn(b*x^3 + a)/b^2 - 1/6*(-a*b^2)^(1/3)*log(x^2 + x*(-a/b )^(1/3) + (-a/b)^(2/3))*sgn(b*x^3 + a)/b^2
Timed out. \[ \int \frac {x^3}{\sqrt {a^2+2 a b x^3+b^2 x^6}} \, dx=\int \frac {x^3}{\sqrt {{\left (b\,x^3+a\right )}^2}} \,d x \] Input:
int(x^3/((a + b*x^3)^2)^(1/2),x)
Output:
int(x^3/((a + b*x^3)^2)^(1/2), x)
Time = 0.19 (sec) , antiderivative size = 78, normalized size of antiderivative = 0.33 \[ \int \frac {x^3}{\sqrt {a^2+2 a b x^3+b^2 x^6}} \, dx=\frac {2 a^{\frac {1}{3}} \sqrt {3}\, \mathit {atan} \left (\frac {a^{\frac {1}{3}}-2 b^{\frac {1}{3}} x}{a^{\frac {1}{3}} \sqrt {3}}\right )+a^{\frac {1}{3}} \mathrm {log}\left (a^{\frac {2}{3}}-b^{\frac {1}{3}} a^{\frac {1}{3}} x +b^{\frac {2}{3}} x^{2}\right )-2 a^{\frac {1}{3}} \mathrm {log}\left (a^{\frac {1}{3}}+b^{\frac {1}{3}} x \right )+6 b^{\frac {1}{3}} x}{6 b^{\frac {4}{3}}} \] Input:
int(x^3/((b*x^3+a)^2)^(1/2),x)
Output:
(2*a**(1/3)*sqrt(3)*atan((a**(1/3) - 2*b**(1/3)*x)/(a**(1/3)*sqrt(3))) + a **(1/3)*log(a**(2/3) - b**(1/3)*a**(1/3)*x + b**(2/3)*x**2) - 2*a**(1/3)*l og(a**(1/3) + b**(1/3)*x) + 6*b**(1/3)*x)/(6*b**(1/3)*b)