Sage trouble related to a MSE question about ellipses

Computation.sage

# This is an attempt to obtain numeric results for
# http://math.stackexchange.com/q/688861/35416
# but the computation is smbolic since I'll want a derivative later on

a = var('a', latex_name='\\alpha', domain='positive')
b = var('b', latex_name='\\beta', domain='positive')
a0 = RDF(0.14778)
b0 = RDF(0.77656)
SRmu.<mu> = SR[]
farpoint1 = vector([1, b, 0])
farpoint2 = vector([-b, 1, 0])
ellipse1 = diagonal_matrix([a, 1, -1])
polar1 = ellipse1*farpoint1
polar2 = ellipse1*farpoint2
def cpm(v):
    """Cross product matrix.
    Multiplying with the returned matrix gives the
    cross product with the input vector."""
    x, y, z = v
    return matrix(v.base_ring(), [
        [0, z, -y],
        [-z, 0, x],
        [y, -x, 0]])
def rank2to1(rank2, c, sqrtIndex):
    indices = [0, 1, 2]
    del indices[sqrtIndex]
    m1 = rank2.matrix_from_rows_and_columns(indices, indices).change_ring(SRmu)
    m2 = c.matrix_from_rows_and_columns(indices, indices).change_ring(SRmu)
    p = list((m1 + mu*m2).det())
    rank1 = sqrt(p[2])*rank2 + sqrt(-p[0])*c
    return rank1
def conicLineIsect(conic, line, row = None, col = None, sqrtIndex = 1):
    cpml = cpm(line)
    rank2 = cpml.transpose()*conic*cpml
    rank1 = rank2to1(rank2, cpml, sqrtIndex=sqrtIndex)
    if row is None and col is None:
        return rank1
    if row is not None:
        row = rank1.row(row)
        if col is None:
            return row
    if col is not None:
        col = rank1.column(col)
        if row is None:
            return col
    return (row, col)
def simplify_full(x):
    try:
        f = f.apply_map
    except:
        return x.simplify_full()
    else:
        return f(simplify_full)
def dehom(v):
    return v[:-1]/v[-1]
matrix([polar1, polar2])
# >>> [   a    b    0]
# >>> [-a*b    1    0]
# view(conicLineIsect(ellipse1, polar1))
p1 = (conicLineIsect(ellipse1, polar1, col=0)/b^2)


p1
# >>> (-b, a, sqrt((a*b^2 + a^2)*b^2)/b)
dehom(p1.subs(a=a0, b=b0))
# >>> (-2.33129749872, 0.443647811323)
p1[2] = sqrt(a*b^2+a^2)
# QUESTION 1: Why can't it do the above automatically, knowing b to b positive?
p1
# >>> (-b, a, sqrt(a*b^2 + a^2))


dehom(p1.subs(a=a0, b=b0))
# >>> (-2.33129749872, 0.443647811323)
# view(conicLineIsect(ellipse1, polar2))
p2 = conicLineIsect(ellipse1, polar2, row=0)
p2
# >>> (-1, -a*b, -sqrt(a^2*b^2 + a))
dehom(p2.subs(a=a0, b=b0))
# >>> (2.49261205011, 0.286052250599)
# view(matrix([p1, p2]))
p3 = simplify_full(farpoint1.cross_product(p1).cross_product(farpoint2.cross_product(p2)))
# view(p3.column())
dehom(p3.subs(a=a0, b=b0))
# >>> (0.601574546983, 2.72119892713)


p3[-1]
# >>> -sqrt(a^2*b^2 + a)*sqrt(a*b^2 + a^2)*(b^2 + 1)
# QUESTION 2: Why aren't the square roots combined into one?
# Like computed below. The argument is known to be positive.
simplify_full((a*b^2 + 1)*(b^2 + a))
# >>> a*b^4 + (a^2 + 1)*b^2 + a
simplify_full((p3[-1]/(a*(b^2+1)))^2)
# >>> a*b^4 + (a^2 + 1)*b^2 + a
# Keep double root for now.
# p3[2] = -sqrt(a*b^4 + (a^2 + 1)*b^2 + a)*a*(b^2+1)
# dehom(p3.subs(x=x0, y=y0))


# Translation:
tr1 = matrix([[p3[2], 0, p3[0]], [0, p3[2], p3[1]], [0, 0, p3[2]]])
# view(tr1)
# Rotation:
tr2 = matrix([
    [-1, -b, 0],
    [b, -1, 0],
    [0, 0, sqrt(1+b^2)]])
...
# view(tr2)
tr3 = tr1*tr2
ellipse2 = (tr3.transpose()*ellipse1*tr3)
# view(ellipse2)
poly1 = (ellipse2.change_ring(SRmu) + mu*ellipse1.change_ring(SRmu)).det()
poly1
# >>> -a*mu^3 + (-(a^2*b^2 + a)*(a*b^2 + a^2)*(b^2 + 1)^2*b^2 - (a^2*b^2 + a)*(a*b^2 + a^2)*(b^2 + 1)^2*a - ((a^2*b^2 + a)*(a*b^2 + a^2)*(b^2 + 1)^2*a*b^2 + (a^2*b^2 + a)*(a*b^2 + a^2)*(b^2 + 1)^3 + (a^2*b^2 + a)*(a*b^2 + a^2)*(b^2 + 1)^2 - (b^2 + 1)*(sqrt(a*b^2 + a^2)*(a*b^2 + 1) - sqrt(a^2*b^2 + a)*(b^3 + a*b))^2*a - (b^2 + 1)*((a*b^3 + b)*sqrt(a*b^2 + a^2) + sqrt(a^2*b^2 + a)*(b^2 + a))^2)*a)*mu^2 + (((a^2*b^2 + a)*(a*b^2 + a^2)*(b^2 + 1)^2*a*b - (a^2*b^2 + a)*(a*b^2 + a^2)*(b^2 + 1)^2*b)^2 - ((a^2*b^2 + a)*(a*b^2 + a^2)*(b^2 + 1)^2*a*b^2 + (a^2*b^2 + a)*(a*b^2 + a^2)*(b^2 + 1)^3 + (a^2*b^2 + a)*(a*b^2 + a^2)*(b^2 + 1)^2 - (b^2 + 1)*(sqrt(a*b^2 + a^2)*(a*b^2 + 1) - sqrt(a^2*b^2 + a)*(b^3 + a*b))^2*a - (b^2 + 1)*((a*b^3 + b)*sqrt(a*b^2 + a^2) + sqrt(a^2*b^2 + a)*(b^2 + a))^2)*((a^2*b^2 + a)*(a*b^2 + a^2)*(b^2 + 1)^2*b^2 + (a^2*b^2 + a)*(a*b^2 + a^2)*(b^2 + 1)^2*a) - (sqrt(a^2*b^2 + a)*sqrt(a*b^2 + a^2)*(b^2 + 1)^(3/2)*(sqrt(a*b^2 + a^2)*(a*b^2 + 1) - sqrt(a^2*b^2 + a)*(b^3 + a*b))*a - sqrt(a^2*b^2 + a)*sqrt(a*b^2 + a^2)*(b^2 + 1)^(3/2)*((a*b^3 + b)*sqrt(a*b^2 + a^2) + sqrt(a^2*b^2 + a)*(b^2 + a))*b)^2 - ((sqrt(a^2*b^2 + a)*sqrt(a*b^2 + a^2)*(b^2 + 1)^(3/2)*(sqrt(a*b^2 + a^2)*(a*b^2 + 1) - sqrt(a^2*b^2 + a)*(b^3 + a*b))*a*b + sqrt(a^2*b^2 + a)*sqrt(a*b^2 + a^2)*(b^2 + 1)^(3/2)*((a*b^3 + b)*sqrt(a*b^2 + a^2) + sqrt(a^2*b^2 + a)*(b^2 + a)))^2 + ((a^2*b^2 + a)*(a*b^2 + a^2)*(b^2 + 1)^2*a*b^2 + (a^2*b^2 + a)*(a*b^2 + a^2)*(b^2 + 1)^2)*((a^2*b^2 + a)*(a*b^2 + a^2)*(b^2 + 1)^3 - (b^2 + 1)*(sqrt(a*b^2 + a^2)*(a*b^2 + 1) - sqrt(a^2*b^2 + a)*(b^3 + a*b))^2*a - (b^2 + 1)*((a*b^3 + b)*sqrt(a*b^2 + a^2) + sqrt(a^2*b^2 + a)*(b^2 + a))^2))*a)*mu + (sqrt(a^2*b^2 + a)*sqrt(a*b^2 + a^2)*(b^2 + 1)^(3/2)*(sqrt(a*b^2 + a^2)*(a*b^2 + 1) - sqrt(a^2*b^2 + a)*(b^3 + a*b))*a - sqrt(a^2*b^2 + a)*sqrt(a*b^2 + a^2)*(b^2 + 1)^(3/2)*((a*b^3 + b)*sqrt(a*b^2 + a^2) + sqrt(a^2*b^2 + a)*(b^2 + a))*b)*(((a^2*b^2 + a)*(a*b^2 + a^2)*(b^2 + 1)^2*a*b - (a^2*b^2 + a)*(a*b^2 + a^2)*(b^2 + 1)^2*b)*(sqrt(a^2*b^2 + a)*sqrt(a*b^2 + a^2)*(b^2 + 1)^(3/2)*(sqrt(a*b^2 + a^2)*(a*b^2 + 1) - sqrt(a^2*b^2 + a)*(b^3 + a*b))*a*b + sqrt(a^2*b^2 + a)*sqrt(a*b^2 + a^2)*(b^2 + 1)^(3/2)*((a*b^3 + b)*sqrt(a*b^2 + a^2) + sqrt(a^2*b^2 + a)*(b^2 + a))) - ((a^2*b^2 + a)*(a*b^2 + a^2)*(b^2 + 1)^2*a*b^2 + (a^2*b^2 + a)*(a*b^2 + a^2)*(b^2 + 1)^2)*(sqrt(a^2*b^2 + a)*sqrt(a*b^2 + a^2)*(b^2 + 1)^(3/2)*(sqrt(a*b^2 + a^2)*(a*b^2 + 1) - sqrt(a^2*b^2 + a)*(b^3 + a*b))*a - sqrt(a^2*b^2 + a)*sqrt(a*b^2 + a^2)*(b^2 + 1)^(3/2)*((a*b^3 + b)*sqrt(a*b^2 + a^2) + sqrt(a^2*b^2 + a)*(b^2 + a))*b)) - ((a^2*b^2 + a)*(a*b^2 + a^2)*(b^2 + 1)^2*b^2 + (a^2*b^2 + a)*(a*b^2 + a^2)*(b^2 + 1)^2*a)*((sqrt(a^2*b^2 + a)*sqrt(a*b^2 + a^2)*(b^2 + 1)^(3/2)*(sqrt(a*b^2 + a^2)*(a*b^2 + 1) - sqrt(a^2*b^2 + a)*(b^3 + a*b))*a*b + sqrt(a^2*b^2 + a)*sqrt(a*b^2 + a^2)*(b^2 + 1)^(3/2)*((a*b^3 + b)*sqrt(a*b^2 + a^2) + sqrt(a^2*b^2 + a)*(b^2 + a)))^2 + ((a^2*b^2 + a)*(a*b^2 + a^2)*(b^2 + 1)^2*a*b^2 + (a^2*b^2 + a)*(a*b^2 + a^2)*(b^2 + 1)^2)*((a^2*b^2 + a)*(a*b^2 + a^2)*(b^2 + 1)^3 - (b^2 + 1)*(sqrt(a*b^2 + a^2)*(a*b^2 + 1) - sqrt(a^2*b^2 + a)*(b^3 + a*b))^2*a - (b^2 + 1)*((a*b^3 + b)*sqrt(a*b^2 + a^2) + sqrt(a^2*b^2 + a)*(b^2 + a))^2)) + ((a^2*b^2 + a)*(a*b^2 + a^2)*(b^2 + 1)^2*a*b - (a^2*b^2 + a)*(a*b^2 + a^2)*(b^2 + 1)^2*b)*(((a^2*b^2 + a)*(a*b^2 + a^2)*(b^2 + 1)^2*a*b - (a^2*b^2 + a)*(a*b^2 + a^2)*(b^2 + 1)^2*b)*((a^2*b^2 + a)*(a*b^2 + a^2)*(b^2 + 1)^3 - (b^2 + 1)*(sqrt(a*b^2 + a^2)*(a*b^2 + 1) - sqrt(a^2*b^2 + a)*(b^3 + a*b))^2*a - (b^2 + 1)*((a*b^3 + b)*sqrt(a*b^2 + a^2) + sqrt(a^2*b^2 + a)*(b^2 + a))^2) + (sqrt(a^2*b^2 + a)*sqrt(a*b^2 + a^2)*(b^2 + 1)^(3/2)*(sqrt(a*b^2 + a^2)*(a*b^2 + 1) - sqrt(a^2*b^2 + a)*(b^3 + a*b))*a*b + sqrt(a^2*b^2 + a)*sqrt(a*b^2 + a^2)*(b^2 + 1)^(3/2)*((a*b^3 + b)*sqrt(a*b^2 + a^2) + sqrt(a^2*b^2 + a)*(b^2 + a)))*(sqrt(a^2*b^2 + a)*sqrt(a*b^2 + a^2)*(b^2 + 1)^(3/2)*(sqrt(a*b^2 + a^2)*(a*b^2 + 1) - sqrt(a^2*b^2 + a)*(b^3 + a*b))*a - sqrt(a^2*b^2 + a)*sqrt(a*b^2 + a^2)*(b^2 + 1)^(3/2)*((a*b^3 + b)*sqrt(a*b^2 + a^2) + sqrt(a^2*b^2 + a)*(b^2 + a))*b))


roots1 = poly1.roots()
# >>> Traceback (most recent call last):
# >>> ...
# >>> NotImplementedError
# QUESTION 3: what's the most elegant way to solve a cubic polynomial
# with big symbolic coefficients? Can we do better than I did below?
lam = var('lam', latex_name='\\lambda')
s1 = solve(poly1(lam), lam)


len(s1)
# >>> 3
[s.rhs().subs(a=a0, b=b0) for s in s1]
# >>> [0.0149758454513*I*sqrt(3) - 0.021984515686, -0.0149758454513*I*sqrt(3) - 0.021984515686, 0.235948995409]
s1[-1].lhs()
# >>> lam
lam1 = s1[-1].rhs()
rank2 = ellipse2 + lam1*ellipse1
rank2.det().is_zero()
# >>> True
rank2.is_symmetric()
# >>> True


#rank2adj = rank2.adjoint()
# QUESTION 4: Why does computing this adjoint above not finish within 45 minutes?
# It should in my opinion be no more difficult than what I did below.
# I know adjoint is implemented differently, since it makes use of
# the characteristic polynomial iirc. But perhaps there should be
# some switch to change implementation to the simple stupid minors as below.
def myadj(m):
    idxs = [[1, 2], [0, 2], [0, 1]]
    return matrix(3, 3, lambda i, j: (-1)^(i+j)*m.matrix_from_rows_and_columns(idxs[j], idxs[i]).det())
rank2adj = myadj(rank2)


rank2adj.subs(a=a0, b=b0)
# >>> [ 0.0421032049219 -0.0118049872463  -0.028707055939]
# >>> [-0.0118049872463  0.0027466541873 0.00766714529439]
# >>> [ -0.028707055939 0.00766714529439  0.0188998548054]
matrix([dehom(i) for i in rank2adj.subs(a=a0, b=b0).columns() + rank2adj.subs(a=a0, b=b0).rows()])
# >>> [ -1.4666500463 0.411222497752]
# >>> [-1.53968482311 0.358236877201]
# >>> [-1.51890351722 0.405672179672]
# >>> [ -1.4666500463 0.411222497752]
# >>> [-1.53968482311 0.358236877201]
# >>> [-1.51890351722 0.405672179672]
# QUESTION 5: The above should all be equal, since the adjoint
# of a rank 2 symmetric matrix is a rank 1 symmetric matrix.
# Is this just numeric imprecision? That's what I'd like to find out,
# by using exact algebraic numbers instead of doubles.
a1 = SR(AA(14778/100000))
b1 = SR(AA(77656/100000))
rank2adj1 = rank2adj.subs(a=a1, b=b1)
rank2adj1
# >>> [ 0.04249332517423462? + 0.?e-36*I -0.01147478603567540? + 0.?e-38*I    -0.02742828122769? + 0.?e-37*I]
# >>> [-0.01147478603567540? + 0.?e-38*I 0.003098621113425321? + 0.?e-37*I 0.007406661096149190? + 0.?e-38*I]
# >>> [   -0.02742828122769? + 0.?e-37*I 0.007406661096149190? + 0.?e-38*I  0.01770420667293761? + 0.?e-37*I]
rank2adj1.parent()
# >>> Full MatrixSpace of 3 by 3 dense matrices over Symbolic Ring


rank2adj1AA = rank2adj1.change_ring(AA)
# >>> Traceback (most recent call last):
# >>> ...
# >>> sage.libs.pari.gen.PariError: ZZ_123633238138850861[y]/(y^2 + 114948053005503479*y + 71680486823734693) is not a field in FpX_ffintersect
# QUESTION 6: So why does working with algebraic numbers fail?
# The exact meaning of this error message kind of eludes me.
# It appears to be very implementation-specific.


rank2adj1QQbar = rank2adj1.change_ring(QQbar)
rank2adj1QQbar.parent()
# >>> Full MatrixSpace of 3 by 3 dense matrices over Algebraic Field
rank2adj1QQbar.change_ring(AA)
# >>> Traceback (most recent call last):
# >>> ...
# >>> AttributeError: 'AlgebraicPolynomialTracker' object has no attribute '_gen'
# QUESTION 7: I'm pretty convinced that all my numbers will be real.
# So a conversion from QQbar to AA should be possible. In any case,
# the error message doesn't make too much sense to me either.
# Trying to verify that the imaginary parts are indeed all equal zero,
# I unfortunatley get the same error message.


rank2adj1QQbar.apply_map(lambda x: (0 if x.real().is_zero() else 1) + (0 if x.imag().is_zero() else I))
# >>> Traceback (most recent call last):
# >>> ...
# >>> AttributeError: 'AlgebraicPolynomialTracker' object has no attribute '_gen'
rank2adj1QQbar.rank()
# >>> Traceback (most recent call last):
# >>> ...
# >>> sage.libs.pari.gen.PariError: ZZ_123633238138850861[y]/(y^2 + 114948053005503479*y + 71680486823734693) is not a field in FpX_ffintersect
# QUESTION 6 again: what does this FpX_ffintersect error message mean?

MinimalReprodcing6.sage

var('a b')
expr = 1/81*(9*(a^2*b^2 + a)*(a*b^2 + a^2)*(b^2 + 1)^2*a*b^2 - 3*(b^8 - b^4)*a^5 - 3*(b^10 - b^8 + 3*b^6 + 5*b^4)*a^4 + 9*(a^2*b^2 + a)*(a*b^2 + a^2)*(b^2 + 1)^2 + 3*(b^10 - 5*b^8 - 4*b^6 - 2*b^4 - 2*(b^5 + b^3)*sqrt(a^2*b^2 + a)*sqrt(a*b^2 + a^2) - 5*b^2 - 1)*a^3 - 3*(b^10 - b^8 + 3*b^6 + 5*b^4 + 2*(b^7 + b)*sqrt(a^2*b^2 + a)*sqrt(a*b^2 + a^2))*a^2 + 6*(b^5 + b^3)*sqrt(a^2*b^2 + a)*sqrt(a*b^2 + a^2) - 3*(b^8 - b^4 - 2*(b^7 + b)*sqrt(a^2*b^2 + a)*sqrt(a*b^2 + a^2))*a + 9*(1/54*(a*b^4 + (a^2 + 1)*b^2 + a)*(b^2 + 1)^3*sqrt((27*(4*a^21 - 13*a^20 + 8*a^19 + 18*a^18 + 8*a^17 - 13*a^16 + 4*a^15)*b^36 + 3*(16*a^23 + 108*a^22 - 539*a^21 + 710*a^20 + 163*a^19 + 1676*a^18 + 163*a^17 + 710*a^16 - 539*a^15 + 108*a^14 + 16*a^13)*b^34 + 3*(64*a^24 + 44*a^23 - 995*a^22 + 3058*a^21 - 2792*a^20 + 9090*a^19 + 5094*a^18 + 9090*a^17 - 2792*a^16 + 3058*a^15 - 995*a^14 + 44*a^13 + 64*a^12)*b^32 + (288*a^25 - 660*a^24 - 1227*a^23 + 15366*a^22 - 23890*a^21 + 84450*a^20 + 5757*a^19 + 192344*a^18 + 5757*a^17 + 84450*a^16 - 23890*a^15 + 15366*a^14 - 1227*a^13 - 660*a^12 + 288*a^11)*b^30 + 6*(32*a^26 - 192*a^25 + 666*a^24 + 1039*a^23 - 1889*a^22 + 22215*a^21 - 8956*a^20 + 88402*a^19 + 17686*a^18 + 88402*a^17 - 8956*a^16 + 22215*a^15 - 1889*a^14 + 1039*a^13 + 666*a^12 - 192*a^11 + 32*a^10)*b^28 + 3*(16*a^27 - 256*a^26 + 2128*a^25 - 3376*a^24 + 13403*a^23 + 19438*a^22 + 6908*a^21 + 244734*a^20 + 25385*a^19 + 617032*a^18 + 25385*a^17 + 244734*a^16 + 6908*a^15 + 19438*a^14 + 13403*a^13 - 3376*a^12 + 2128*a^11 - 256*a^10 + 16*a^9)*b^26 - (192*a^27 - 3552*a^26 + 13056*a^25 - 72208*a^24 + 102480*a^23 - 332817*a^22 - 263254*a^21 - 417816*a^20 - 2556378*a^19 - 959054*a^18 - 2556378*a^17 - 417816*a^16 - 263254*a^15 - 332817*a^14 + 102480*a^13 - 72208*a^12 + 13056*a^11 - 3552*a^10 + 192*a^9)*b^24 + 3*(224*a^27 - 1920*a^26 + 16096*a^25 - 50504*a^24 + 182206*a^23 - 246940*a^22 + 686309*a^21 - 172838*a^20 + 1652381*a^19 + 452628*a^18 + 1652381*a^17 - 172838*a^16 + 686309*a^15 - 246940*a^14 + 182206*a^13 - 50504*a^12 + 16096*a^11 - 1920*a^10 + 224*a^9)*b^22 - 6*(160*a^27 - 2048*a^26 + 12448*a^25 - 63772*a^24 + 178110*a^23 - 523410*a^22 + 823204*a^21 - 1849465*a^20 + 1350426*a^19 - 3001882*a^18 + 1350426*a^17 - 1849465*a^16 + 823204*a^15 - 523410*a^14 + 178110*a^13 - 63772*a^12 + 12448*a^11 - 2048*a^10 + 160*a^9)*b^20 + (176*a^27 - 12672*a^26 + 104688*a^25 - 496720*a^24 + 1942122*a^23 - 4819512*a^22 + 11422375*a^21 - 15621498*a^20 + 28665279*a^19 - 21364636*a^18 + 28665279*a^17 - 15621498*a^16 + 11422375*a^15 - 4819512*a^14 + 1942122*a^13 - 496720*a^12 + 104688*a^11 - 12672*a^10 + 176*a^9)*b^18 + 432*a^18 + 3*(544*a^26 - 20416*a^25 + 142960*a^24 - 562444*a^23 + 1750931*a^22 - 3594210*a^21 + 7304088*a^20 - 7846386*a^19 + 11951018*a^18 - 7846386*a^17 + 7304088*a^16 - 3594210*a^15 + 1750931*a^14 - 562444*a^13 + 142960*a^12 - 20416*a^11 + 544*a^10)*b^16 + 3*(1184*a^25 - 42812*a^24 + 271963*a^23 - 924070*a^22 + 2450836*a^21 - 3842614*a^20 + 7208433*a^19 - 5663184*a^18 + 7208433*a^17 - 3842614*a^16 + 2450836*a^15 - 924070*a^14 + 271963*a^13 - 42812*a^12 + 1184*a^11)*b^14 - 2*(2098*a^24 + 45435*a^23 - 322221*a^22 + 967015*a^21 - 2580324*a^20 + 2430438*a^19 - 5094706*a^18 + 2430438*a^17 - 2580324*a^16 + 967015*a^15 - 322221*a^14 + 45435*a^13 + 2098*a^12)*b^12 - 3*(6585*a^23 - 15694*a^22 - 12306*a^21 - 2650*a^20 - 523527*a^19 - 138608*a^18 - 523527*a^17 - 2650*a^16 - 12306*a^15 - 15694*a^14 + 6585*a^13)*b^10 - 864*((a^9 + 3*a^8 - 8*a^6 - 6*a^5 + 6*a^4 + 8*a^3 - 3*a - 1)*b^15 + 3*(a^9 + 3*a^8 - 8*a^6 - 6*a^5 + 6*a^4 + 8*a^3 - 3*a - 1)*b^13 + 3*(a^9 + 3*a^8 - 8*a^6 - 6*a^5 + 6*a^4 + 8*a^3 - 3*a - 1)*b^11 + (a^9 + 3*a^8 - 8*a^6 - 6*a^5 + 6*a^4 + 8*a^3 - 3*a - 1)*b^9)*(a^2*b^2 + a)^(9/2)*(a*b^2 + a^2)^(9/2) - 3*(6795*a^22 - 29614*a^21 + 50456*a^20 - 225778*a^19 - 44358*a^18 - 225778*a^17 + 50456*a^16 - 29614*a^15 + 6795*a^14)*b^8 + 864*((a^12 + 3*a^11 - 8*a^9 - 6*a^8 + 6*a^7 + 8*a^6 - 3*a^4 - a^3)*b^19 + (a^13 + 10*a^12 + 10*a^11 - 17*a^10 - 26*a^9 - 8*a^8 + 8*a^7 + 26*a^6 + 17*a^5 - 10*a^4 - 10*a^3 - a^2)*b^17 + (7*a^13 + 21*a^12 + 11*a^11 - 31*a^10 - 54*a^9 - 26*a^8 + 26*a^7 + 54*a^6 + 31*a^5 - 11*a^4 - 21*a^3 - 7*a^2)*b^15 + (11*a^13 + 25*a^12 + 11*a^11 - 31*a^10 - 74*a^9 - 46*a^8 + 46*a^7 + 74*a^6 + 31*a^5 - 11*a^4 - 25*a^3 - 11*a^2)*b^13 + (5*a^13 + 22*a^12 + 18*a^11 - 33*a^10 - 62*a^9 - 28*a^8 + 28*a^7 + 62*a^6 + 33*a^5 - 18*a^4 - 22*a^3 - 5*a^2)*b^11 + 3*(3*a^12 + 5*a^11 - 4*a^10 - 12*a^9 - 6*a^8 + 6*a^7 + 12*a^6 + 4*a^5 - 5*a^4 - 3*a^3)*b^9 + 4*(a^11 + a^10 - 3*a^9 - 3*a^8 + 3*a^7 + 3*a^6 - a^5 - a^4)*b^7)*(a^2*b^2 + a)^(7/2)*(a*b^2 + a^2)^(7/2) - 4*(27*(a^14 + 4*a^12 - 40*a^11 - 32*a^10 + 134*a^9 - 32*a^8 - 40*a^7 + 4*a^6 + a^4)*b^22 + 9*(3*a^15 - 79*a^13 - 54*a^12 - 206*a^11 - 624*a^10 + 1920*a^9 - 624*a^8 - 206*a^7 - 54*a^6 - 79*a^5 + 3*a^3)*b^20 - (585*a^14 + 486*a^13 + 5491*a^12 + 4908*a^11 + 4425*a^10 - 31790*a^9 + 4425*a^8 + 4908*a^7 + 5491*a^6 + 486*a^5 + 585*a^4)*b^18 + 3*(78*a^15 - 72*a^14 - 1357*a^13 - 834*a^12 - 9634*a^11 + 3314*a^10 + 17010*a^9 + 3314*a^8 - 9634*a^7 - 834*a^6 - 1357*a^5 - 72*a^4 + 78*a^3)*b^16 + 3*(288*a^15 + 323*a^14 - 1116*a^13 - 4005*a^12 - 11692*a^11 + 1964*a^10 + 28476*a^9 + 1964*a^8 - 11692*a^7 - 4005*a^6 - 1116*a^5 + 323*a^4 + 288*a^3)*b^14 + (539*a^15 + 3480*a^14 - 3597*a^13 - 8026*a^12 - 54141*a^11 + 35070*a^10 + 53350*a^9 + 35070*a^8 - 54141*a^7 - 8026*a^6 - 3597*a^5 + 3480*a^4 + 539*a^3)*b^12 + 3*(695*a^14 + 222*a^13 - 1851*a^12 - 11164*a^11 - 562*a^10 + 25320*a^9 - 562*a^8 - 11164*a^7 - 1851*a^6 + 222*a^5 + 695*a^4)*b^10 + 3*(329*a^13 + 834*a^12 - 8212*a^11 + 1574*a^10 + 10950*a^9 + 1574*a^8 - 8212*a^7 + 834*a^6 + 329*a^5)*b^8 + (359*a^12 - 3108*a^11 - 10860*a^10 + 27218*a^9 - 10860*a^8 - 3108*a^7 + 359*a^6)*b^6 + 27*(7*a^11 - 226*a^10 + 438*a^9 - 226*a^8 + 7*a^7)*b^4 - 729*(a^10 - 2*a^9 + a^8)*b^2)*(a^2*b^2 + a)^3*(a*b^2 + a^2)^3 - (8089*a^21 - 35922*a^20 - 17817*a^19 - 261212*a^18 - 17817*a^17 - 35922*a^16 + 8089*a^15)*b^6 - 288*((a^15 + 3*a^14 - 8*a^12 - 6*a^11 + 6*a^10 + 8*a^9 - 3*a^7 - a^6)*b^23 + (2*a^16 + 17*a^15 + 11*a^14 - 34*a^13 - 28*a^12 + 2*a^11 - 2*a^10 + 28*a^9 + 34*a^8 - 11*a^7 - 17*a^6 - 2*a^5)*b^21 + (a^17 + 25*a^16 + 57*a^15 - 53*a^14 - 11*a^13 - 75*a^12 - 150*a^11 + 150*a^10 + 75*a^9 + 11*a^8 + 53*a^7 - 57*a^6 - 25*a^5 - a^4)*b^19 + (11*a^17 + 81*a^16 - 11*a^15 + 111*a^14 - 207*a^13 - 269*a^12 + 130*a^11 - 130*a^10 + 269*a^9 + 207*a^8 - 111*a^7 + 11*a^6 - 81*a^5 - 11*a^4)*b^17 + (35*a^17 + 41*a^16 + 173*a^15 - 121*a^14 - 135*a^13 - 197*a^12 - 350*a^11 + 350*a^10 + 197*a^9 + 135*a^8 + 121*a^7 - 173*a^6 - 41*a^5 - 35*a^4)*b^15 + (25*a^17 + 73*a^16 + 105*a^15 - 5*a^14 - 203*a^13 - 291*a^12 - 150*a^11 + 150*a^10 + 291*a^9 + 203*a^8 + 5*a^7 - 105*a^6 - 73*a^5 - 25*a^4)*b^13 + (90*a^16 + a^15 + 123*a^14 - 202*a^13 - 316*a^12 + 98*a^11 - 98*a^10 + 316*a^9 + 202*a^8 - 123*a^7 - a^6 - 90*a^5)*b^11 + (121*a^15 - 93*a^14 + 8*a^13 - 64*a^12 - 286*a^11 + 286*a^10 + 64*a^9 - 8*a^8 + 93*a^7 - 121*a^6)*b^9 + 72*(a^14 - a^13 - a^12 + a^11 - a^10 + a^9 + a^8 - a^7)*b^7 + 16*(a^13 - a^12 - 2*a^11 + 2*a^10 + a^9 - a^8)*b^5)*(a^2*b^2 + a)^(5/2)*(a*b^2 + a^2)^(5/2) + 24*(3*(a^17 + 13*a^15 - 49*a^14 - 68*a^13 + 206*a^12 - 68*a^11 - 49*a^10 + 13*a^9 + a^7)*b^26 + (6*a^18 + 12*a^17 - 146*a^16 + 192*a^15 - 655*a^14 - 1368*a^13 + 3918*a^12 - 1368*a^11 - 655*a^10 + 192*a^9 - 146*a^8 + 12*a^7 + 6*a^6)*b^24 + (3*a^19 + 24*a^18 - 387*a^17 + 409*a^16 - 1059*a^15 - 1981*a^14 - 1149*a^13 + 8280*a^12 - 1149*a^11 - 1981*a^10 - 1059*a^9 + 409*a^8 - 387*a^7 + 24*a^6 + 3*a^5)*b^22 + (12*a^19 - 180*a^18 - 102*a^17 - 1433*a^16 - 166*a^15 - 8951*a^14 + 5308*a^13 + 11024*a^12 + 5308*a^11 - 8951*a^10 - 166*a^9 - 1433*a^8 - 102*a^7 - 180*a^6 + 12*a^5)*b^20 + (22*a^19 + 72*a^18 - 994*a^17 - 2247*a^16 - 1574*a^15 - 20906*a^14 + 25232*a^13 + 790*a^12 + 25232*a^11 - 20906*a^10 - 1574*a^9 - 2247*a^8 - 994*a^7 + 72*a^6 + 22*a^5)*b^18 + (244*a^19 + 126*a^18 - 1372*a^17 + 39*a^16 - 21516*a^15 + 13852*a^14 - 31492*a^13 + 80238*a^12 - 31492*a^11 + 13852*a^10 - 21516*a^9 + 39*a^8 - 1372*a^7 + 126*a^6 + 244*a^5)*b^16 + (295*a^19 + 704*a^18 - 323*a^17 - 7661*a^16 + 851*a^15 - 42446*a^14 + 58937*a^13 - 20714*a^12 + 58937*a^11 - 42446*a^10 + 851*a^9 - 7661*a^8 - 323*a^7 + 704*a^6 + 295*a^5)*b^14 + (1744*a^18 - 2858*a^17 + 3957*a^16 - 27026*a^15 + 19534*a^14 - 35276*a^13 + 79850*a^12 - 35276*a^11 + 19534*a^10 - 27026*a^9 + 3957*a^8 - 2858*a^7 + 1744*a^6)*b^12 + (3205*a^17 - 9381*a^16 + 9561*a^15 - 33611*a^14 + 40156*a^13 - 19860*a^12 + 40156*a^11 - 33611*a^10 + 9561*a^9 - 9381*a^8 + 3205*a^7)*b^10 + (1999*a^16 - 7900*a^15 + 4687*a^14 - 8672*a^13 + 19772*a^12 - 8672*a^11 + 4687*a^10 - 7900*a^9 + 1999*a^8)*b^8 - (170*a^15 + 1037*a^14 + 1318*a^13 - 5050*a^12 + 1318*a^11 + 1037*a^10 + 170*a^9)*b^6 - (467*a^14 - 908*a^13 + 882*a^12 - 908*a^11 + 467*a^10)*b^4 - 54*(a^13 - 2*a^12 + a^11)*b^2)*(a^2*b^2 + a)^2*(a*b^2 + a^2)^2 - 3*(445*a^20 - 4460*a^19 - 14002*a^18 - 4460*a^17 + 445*a^16)*b^4 - 12*(9*(a^17 - 16*a^16 + 15*a^15 + 32*a^14 - 32*a^13 - 15*a^12 + 16*a^11 - a^10)*b^29 - (50*a^18 + 309*a^17 - 225*a^16 - 412*a^15 - 900*a^14 + 900*a^13 + 412*a^12 + 225*a^11 - 309*a^10 - 50*a^9)*b^27 - (159*a^19 + 270*a^18 + 118*a^17 + 495*a^16 - 2666*a^15 - 4630*a^14 + 4630*a^13 + 2666*a^12 - 495*a^11 - 118*a^10 - 270*a^9 - 159*a^8)*b^25 - (132*a^20 + 189*a^19 + 1387*a^18 - 805*a^17 + 185*a^16 - 5546*a^15 - 19266*a^14 + 19266*a^13 + 5546*a^12 - 185*a^11 + 805*a^10 - 1387*a^9 - 189*a^8 - 132*a^7)*b^23 - (32*a^21 + 84*a^20 + 2355*a^19 - 3247*a^18 + 8341*a^17 - 9058*a^16 - 25898*a^15 - 10877*a^14 + 10877*a^13 + 25898*a^12 + 9058*a^11 - 8341*a^10 + 3247*a^9 - 2355*a^8 - 84*a^7 - 32*a^6)*b^21 - (1336*a^20 - 1575*a^19 + 11311*a^18 - 18457*a^17 + 19637*a^16 - 83145*a^15 + 3575*a^14 - 3575*a^13 + 83145*a^12 - 19637*a^11 + 18457*a^10 - 11311*a^9 + 1575*a^8 - 1336*a^7)*b^19 - (160*a^21 + 440*a^20 + 4253*a^19 - 3755*a^18 + 19035*a^17 - 53273*a^16 - 26762*a^15 - 101138*a^14 + 101138*a^13 + 26762*a^12 + 53273*a^11 - 19035*a^10 + 3755*a^9 - 4253*a^8 - 440*a^7 - 160*a^6)*b^17 - (192*a^21 - 380*a^20 + 6359*a^19 - 4359*a^18 + 18139*a^17 - 45571*a^16 - 82914*a^15 - 49690*a^14 + 49690*a^13 + 82914*a^12 + 45571*a^11 - 18139*a^10 + 4359*a^9 - 6359*a^8 + 380*a^7 - 192*a^6)*b^15 - (564*a^20 - 1199*a^19 + 16555*a^18 - 16672*a^17 - 4615*a^16 - 106989*a^15 - 37604*a^14 + 37604*a^13 + 106989*a^12 + 4615*a^11 + 16672*a^10 - 16555*a^9 + 1199*a^8 - 564*a^7)*b^13 - (595*a^19 + 3707*a^18 + 9024*a^17 - 29528*a^16 - 56256*a^15 - 44068*a^14 + 44068*a^13 + 56256*a^12 + 29528*a^11 - 9024*a^10 - 3707*a^9 - 595*a^8)*b^11 - (2273*a^18 + 1795*a^17 - 3850*a^16 - 49156*a^15 - 5300*a^14 + 5300*a^13 + 49156*a^12 + 3850*a^11 - 1795*a^10 - 2273*a^9)*b^9 - 2*(2199*a^17 - 6073*a^16 + 2306*a^15 - 19090*a^14 + 19090*a^13 - 2306*a^12 + 6073*a^11 - 2199*a^10)*b^7 - (2509*a^16 - 12074*a^15 + 9169*a^14 - 9169*a^13 + 12074*a^12 - 2509*a^11)*b^5 + (19*a^15 + 2319*a^14 - 2319*a^13 - 19*a^12)*b^3 + 180*(a^14 - a^13)*b)*(a^2*b^2 + a)^(3/2)*(a*b^2 + a^2)^(3/2) - 12*((a^20 + 22*a^18 - 58*a^17 - 131*a^16 + 332*a^15 - 131*a^14 - 58*a^13 + 22*a^12 + a^10)*b^30 + (3*a^21 + 8*a^20 - 53*a^19 + 200*a^18 - 591*a^17 - 856*a^16 + 2578*a^15 - 856*a^14 - 591*a^13 + 200*a^12 - 53*a^11 + 8*a^10 + 3*a^9)*b^28 + (3*a^22 + 24*a^21 - 269*a^20 + 356*a^19 + 391*a^18 - 4380*a^17 + 1279*a^16 + 5192*a^15 + 1279*a^14 - 4380*a^13 + 391*a^12 + 356*a^11 - 269*a^10 + 24*a^9 + 3*a^8)*b^26 + (a^23 + 24*a^22 - 269*a^21 - 592*a^20 + 3135*a^19 - 9072*a^18 + 4761*a^17 - 15272*a^16 + 34568*a^15 - 15272*a^14 + 4761*a^13 - 9072*a^12 + 3135*a^11 - 592*a^10 - 269*a^9 + 24*a^8 + a^7)*b^24 + (8*a^23 - 53*a^22 - 1042*a^21 + 3539*a^20 - 15292*a^19 + 33380*a^18 - 93540*a^17 + 129394*a^16 - 112788*a^15 + 129394*a^14 - 93540*a^13 + 33380*a^12 - 15292*a^11 + 3539*a^10 - 1042*a^9 - 53*a^8 + 8*a^7)*b^22 + (22*a^23 - 232*a^22 + 1831*a^21 - 15432*a^20 + 46088*a^19 - 137144*a^18 + 216777*a^17 - 320656*a^16 + 417492*a^15 - 320656*a^14 + 216777*a^13 - 137144*a^12 + 46088*a^11 - 15432*a^10 + 1831*a^9 - 232*a^8 + 22*a^7)*b^20 + (120*a^23 + 1001*a^22 - 6836*a^21 + 24215*a^20 - 97144*a^19 + 196257*a^18 - 407140*a^17 + 580123*a^16 - 581192*a^15 + 580123*a^14 - 407140*a^13 + 196257*a^12 - 97144*a^11 + 24215*a^10 - 6836*a^9 + 1001*a^8 + 120*a^7)*b^18 + (425*a^23 - 816*a^22 + 8461*a^21 - 43976*a^20 + 107728*a^19 - 293864*a^18 + 479963*a^17 - 671616*a^16 + 827390*a^15 - 671616*a^14 + 479963*a^13 - 293864*a^12 + 107728*a^11 - 43976*a^10 + 8461*a^9 - 816*a^8 + 425*a^7)*b^16 + (3145*a^22 - 11698*a^21 + 40352*a^20 - 144364*a^19 + 280241*a^18 - 531348*a^17 + 750850*a^16 - 774356*a^15 + 750850*a^14 - 531348*a^13 + 280241*a^12 - 144364*a^11 + 40352*a^10 - 11698*a^9 + 3145*a^8)*b^14 + (8822*a^21 - 39080*a^20 + 100277*a^19 - 264296*a^18 + 441166*a^17 - 611304*a^16 + 728830*a^15 - 611304*a^14 + 441166*a^13 - 264296*a^12 + 100277*a^11 - 39080*a^10 + 8822*a^9)*b^12 + 2*(6145*a^20 - 29954*a^19 + 66630*a^18 - 135538*a^17 + 198051*a^16 - 210668*a^15 + 198051*a^14 - 135538*a^13 + 66630*a^12 - 29954*a^11 + 6145*a^10)*b^10 + (8873*a^19 - 47184*a^18 + 94988*a^17 - 148224*a^16 + 183094*a^15 - 148224*a^14 + 94988*a^13 - 47184*a^12 + 8873*a^11)*b^8 + (2849*a^18 - 17930*a^17 + 33339*a^16 - 36516*a^15 + 33339*a^14 - 17930*a^13 + 2849*a^12)*b^6 + 8*(3*a^17 - 257*a^16 + 508*a^15 - 257*a^14 + 3*a^13)*b^4 - 140*(a^16 - 2*a^15 + a^14)*b^2)*(a^2*b^2 + a)*(a*b^2 + a^2) + 24*(11*a^19 + 302*a^18 + 11*a^17)*b^2 + 12*(9*(a^20 - 8*a^19 + 8*a^18 + 17*a^17 - 17*a^16 - 8*a^15 + 8*a^14 - a^13)*b^33 - (17*a^21 + 93*a^20 + 144*a^19 - 706*a^18 - 558*a^17 + 558*a^16 + 706*a^15 - 144*a^14 - 93*a^13 - 17*a^12)*b^31 - (113*a^22 - 117*a^21 + 1215*a^20 - 2152*a^19 + 31*a^18 - 6668*a^17 + 6668*a^16 - 31*a^15 + 2152*a^14 - 1215*a^13 + 117*a^12 - 113*a^11)*b^29 - (147*a^23 - 165*a^22 + 2533*a^21 - 5497*a^20 + 9331*a^19 - 25627*a^18 + 5492*a^17 - 5492*a^16 + 25627*a^15 - 9331*a^14 + 5497*a^13 - 2533*a^12 + 165*a^11 - 147*a^10)*b^27 - (68*a^24 + 57*a^23 + 2725*a^22 - 7472*a^21 + 21905*a^20 - 49597*a^19 + 35863*a^18 - 108467*a^17 + 108467*a^16 - 35863*a^15 + 49597*a^14 - 21905*a^13 + 7472*a^12 - 2725*a^11 - 57*a^10 - 68*a^9)*b^25 - (8*a^25 + 108*a^24 + 1455*a^23 - 3697*a^22 + 20850*a^21 - 59416*a^20 + 112503*a^19 - 296579*a^18 + 260984*a^17 - 260984*a^16 + 296579*a^15 - 112503*a^14 + 59416*a^13 - 20850*a^12 + 3697*a^11 - 1455*a^10 - 108*a^9 - 8*a^8)*b^23 - (24*a^25 + 160*a^24 + 509*a^23 + 7963*a^22 - 28378*a^21 + 94606*a^20 - 266818*a^19 + 259848*a^18 - 525888*a^17 + 525888*a^16 - 259848*a^15 + 266818*a^14 - 94606*a^13 + 28378*a^12 - 7963*a^11 - 509*a^10 - 160*a^9 - 24*a^8)*b^21 + (104*a^25 - 512*a^24 + 343*a^23 - 6737*a^22 - 2718*a^21 + 23998*a^20 + 22723*a^19 + 237097*a^18 - 69190*a^17 + 69190*a^16 - 237097*a^15 - 22723*a^14 - 23998*a^13 + 2718*a^12 + 6737*a^11 - 343*a^10 + 512*a^9 - 104*a^8)*b^19 + (120*a^25 + 772*a^24 - 5659*a^23 + 15113*a^22 - 70076*a^21 + 169859*a^20 - 201439*a^19 + 571789*a^18 - 446157*a^17 + 446157*a^16 - 571789*a^15 + 201439*a^14 - 169859*a^13 + 70076*a^12 - 15113*a^11 + 5659*a^10 - 772*a^9 - 120*a^8)*b^17 + (1164*a^24 - 533*a^23 - 16549*a^22 + 47611*a^21 - 147949*a^20 + 428479*a^19 - 442011*a^18 + 785970*a^17 - 785970*a^16 + 442011*a^15 - 428479*a^14 + 147949*a^13 - 47611*a^12 + 16549*a^11 + 533*a^10 - 1164*a^9)*b^15 + (4433*a^23 - 15576*a^22 + 6113*a^21 + 19053*a^20 - 19585*a^19 + 267808*a^18 - 203604*a^17 + 203604*a^16 - 267808*a^15 + 19585*a^14 - 19053*a^13 - 6113*a^12 + 15576*a^11 - 4433*a^10)*b^13 + (8160*a^22 - 41973*a^21 + 81425*a^20 - 95468*a^19 + 222374*a^18 - 157264*a^17 + 157264*a^16 - 222374*a^15 + 95468*a^14 - 81425*a^13 + 41973*a^12 - 8160*a^11)*b^11 + (7516*a^21 - 45851*a^20 + 115790*a^19 - 139922*a^18 + 225169*a^17 - 225169*a^16 + 139922*a^15 - 115790*a^14 + 45851*a^13 - 7516*a^12)*b^9 + 2*(1525*a^20 - 10126*a^19 + 31409*a^18 - 34112*a^17 + 34112*a^16 - 31409*a^15 + 10126*a^14 - 1525*a^13)*b^7 + (243*a^19 - 805*a^18 + 9840*a^17 - 9840*a^16 + 805*a^15 - 243*a^14)*b^5 + 2*(5*a^18 + 561*a^17 - 561*a^16 - 5*a^15)*b^3 + 72*(a^17 - a^16)*b)*sqrt(a^2*b^2 + a)*sqrt(a*b^2 + a^2))*(a*b^4 + (a^2 + 1)*b^2 + a)/a)/a - 1/54*(2*(b^24 - 3*b^20 + 3*b^16 - b^12)*a^16 + 6*(b^26 - b^24 - 5*b^22 - 5*b^20 + b^18 + 13*b^16 + 9*b^14 - 7*b^12 - 6*b^10)*a^15 + 3*(2*b^28 - 6*b^26 - 17*b^24 - 30*b^22 - 38*b^20 - 82*b^18 - 124*b^16 + 54*b^14 + 196*b^12 + 64*b^10 - 19*b^8 - 2*(b^21 + 7*b^19 + 10*b^17 - 2*b^15 - 11*b^13 - 5*b^11)*sqrt(a^2*b^2 + a)*sqrt(a*b^2 + a^2))*a^14 + (2*b^30 - 18*b^28 - 33*b^26 - 68*b^24 - 129*b^22 - 282*b^20 + 792*b^18 + 2922*b^16 + 2682*b^14 + 1636*b^12 + 1761*b^10 + 930*b^8 + 45*b^6 - 6*(3*b^23 + 21*b^21 + 36*b^19 + 56*b^17 + 91*b^15 + 45*b^13 - 34*b^11 - 26*b^9)*sqrt(a^2*b^2 + a)*sqrt(a*b^2 + a^2))*a^13 - 3*(2*b^30 + b^28 - 12*b^26 - 90*b^24 - 498*b^22 - 1893*b^20 - 3232*b^18 - 3598*b^16 - 4314*b^14 - 3843*b^12 - 1780*b^10 - 768*b^8 - 406*b^6 - 49*b^4 + 4*(b^18 + 8*b^16 + 18*b^14 + 16*b^12 + 5*b^10)*(a^2*b^2 + a)*(a*b^2 + a^2) + 2*(3*b^25 + 21*b^23 + 41*b^21 + 122*b^19 + 136*b^17 - 33*b^15 - 3*b^13 + 100*b^11 - 17*b^9 - 50*b^7)*sqrt(a^2*b^2 + a)*sqrt(a*b^2 + a^2))*a^12 + 3*(b^30 + 22*b^28 + 160*b^26 + 846*b^24 + 2568*b^22 + 4632*b^20 + 7705*b^18 + 10840*b^16 + 9717*b^14 + 6834*b^12 + 4902*b^10 + 2218*b^8 + 514*b^6 + 208*b^4 - 4*(3*b^20 + 24*b^18 + 45*b^16 + 40*b^14 + 45*b^12 + 48*b^10 + 19*b^8)*(a^2*b^2 + a)*(a*b^2 + a^2) - 2*(b^27 + 7*b^25 + 7*b^23 + 17*b^21 - 211*b^19 - 378*b^17 - 98*b^15 - 227*b^13 - 384*b^11 - 15*b^9 + 45*b^7 - 44*b^5)*sqrt(a^2*b^2 + a)*sqrt(a*b^2 + a^2) + 33*b^2)*a^11 + 2*(11*b^30 + 99*b^28 + 579*b^26 + 1865*b^24 + 4701*b^22 + 10887*b^20 + 16687*b^18 + 18402*b^16 + 19359*b^14 + 15839*b^12 + 8067*b^10 + 3819*b^8 + 1745*b^6 + (9*b^19 + 27*b^17 + 26*b^15 + 6*b^13 - 3*b^11 - b^9)*(a^2*b^2 + a)^(3/2)*(a*b^2 + a^2)^(3/2) + 279*b^4 - 6*(3*b^22 + 24*b^20 + 25*b^18 + 32*b^16 + 100*b^14 + 88*b^12 + 29*b^10 + 40*b^8 + 27*b^6)*(a^2*b^2 + a)*(a*b^2 + a^2) + 3*(15*b^25 + 91*b^23 + 374*b^21 + 393*b^19 + 504*b^17 + 1164*b^15 + 732*b^13 + 90*b^11 + 357*b^9 + 165*b^7 - 62*b^5 + 17*b^3)*sqrt(a^2*b^2 + a)*sqrt(a*b^2 + a^2) + 51*b^2 + 10)*a^10 + 3*(b^30 + 22*b^28 + 160*b^26 + 846*b^24 + 2568*b^22 + 4632*b^20 + 7705*b^18 + 10840*b^16 + 9717*b^14 + 6834*b^12 + 4902*b^10 + 2218*b^8 + 514*b^6 + 2*(6*b^21 + 18*b^19 + 35*b^17 + 61*b^15 + 68*b^13 + 44*b^11 + 19*b^9 + 5*b^7)*(a^2*b^2 + a)^(3/2)*(a*b^2 + a^2)^(3/2) + 208*b^4 - 4*(b^24 + 8*b^22 - 15*b^20 + 16*b^18 + 62*b^16 - 32*b^14 - 2*b^12 + 80*b^10 + b^8 - 8*b^6 + 17*b^4)*(a^2*b^2 + a)*(a*b^2 + a^2) + 2*(7*b^27 + 37*b^25 + 92*b^23 + 22*b^21 + 451*b^19 + 738*b^17 + 16*b^15 + 364*b^13 + 723*b^11 - 69*b^9 - 48*b^7 + 122*b^5 - 25*b^3 + 2*b)*sqrt(a^2*b^2 + a)*sqrt(a*b^2 + a^2) + 33*b^2)*a^9 - 3*(2*b^30 + b^28 - 12*b^26 - 90*b^24 - 498*b^22 - 1893*b^20 - 3232*b^18 - 3598*b^16 - 4314*b^14 - 3843*b^12 - 1780*b^10 - 768*b^8 - 406*b^6 - 12*(b^16 + 4*b^14 + 6*b^12 + 4*b^10 + b^8)*(a^2*b^2 + a)^2*(a*b^2 + a^2)^2 - 2*(3*b^23 + 9*b^21 + 32*b^19 + 45*b^17 - 4*b^15 - 67*b^13 - 100*b^11 - 85*b^9 - 27*b^7 + 2*b^5)*(a^2*b^2 + a)^(3/2)*(a*b^2 + a^2)^(3/2) - 49*b^4 - 4*(15*b^22 - 16*b^20 + 38*b^18 + 136*b^16 + 14*b^14 + 56*b^12 + 154*b^10 + 8*b^8 - 9*b^6 + 24*b^4 - 4*b^2)*(a^2*b^2 + a)*(a*b^2 + a^2) + 2*(7*b^27 + 37*b^25 + 92*b^23 + 22*b^21 + 451*b^19 + 738*b^17 + 16*b^15 + 364*b^13 + 723*b^11 - 69*b^9 - 48*b^7 + 122*b^5 - 25*b^3 + 2*b)*sqrt(a^2*b^2 + a)*sqrt(a*b^2 + a^2))*a^8 + (2*b^30 - 18*b^28 - 33*b^26 - 68*b^24 - 129*b^22 - 282*b^20 + 792*b^18 + 2922*b^16 + 2682*b^14 + 1636*b^12 + 1761*b^10 + 930*b^8 + 45*b^6 + 72*(b^18 + 4*b^16 + 7*b^14 + 8*b^12 + 7*b^10 + 4*b^8 + b^6)*(a^2*b^2 + a)^2*(a*b^2 + a^2)^2 - 2*(19*b^21 + 219*b^19 + 507*b^17 + 720*b^15 + 1023*b^13 + 1140*b^11 + 981*b^9 + 732*b^7 + 318*b^5 + 37*b^3)*(a^2*b^2 + a)^(3/2)*(a*b^2 + a^2)^(3/2) + 24*(b^24 - 4*b^22 + 28*b^20 + 28*b^18 + 11*b^16 + 112*b^14 + 91*b^12 + 8*b^10 + 52*b^8 + 28*b^6 - 7*b^4 + 4*b^2)*(a^2*b^2 + a)*(a*b^2 + a^2) - 6*(15*b^25 + 91*b^23 + 374*b^21 + 393*b^19 + 504*b^17 + 1164*b^15 + 732*b^13 + 90*b^11 + 357*b^9 + 165*b^7 - 62*b^5 + 17*b^3)*sqrt(a^2*b^2 + a)*sqrt(a*b^2 + a^2))*a^7 + 3*(2*b^28 - 6*b^26 - 17*b^24 - 30*b^22 - 38*b^20 - 82*b^18 - 124*b^16 + 54*b^14 + 196*b^12 + 64*b^10 - 19*b^8 + 12*(b^20 + 4*b^18 + 8*b^16 + 12*b^14 + 14*b^12 + 12*b^10 + 8*b^8 + 4*b^6 + b^4)*(a^2*b^2 + a)^2*(a*b^2 + a^2)^2 - 2*(12*b^23 + 47*b^21 + 49*b^19 + 89*b^17 + 215*b^15 + 270*b^13 + 297*b^11 + 254*b^9 + 113*b^7 + 67*b^5 + 50*b^3 + 9*b)*(a^2*b^2 + a)^(3/2)*(a*b^2 + a^2)^(3/2) + 4*(15*b^22 - 16*b^20 + 38*b^18 + 136*b^16 + 14*b^14 + 56*b^12 + 154*b^10 + 8*b^8 - 9*b^6 + 24*b^4 - 4*b^2)*(a^2*b^2 + a)*(a*b^2 + a^2) + 2*(b^27 + 7*b^25 + 7*b^23 + 17*b^21 - 211*b^19 - 378*b^17 - 98*b^15 - 227*b^13 - 384*b^11 - 15*b^9 + 45*b^7 - 44*b^5)*sqrt(a^2*b^2 + a)*sqrt(a*b^2 + a^2))*a^6 + 6*(b^26 - b^24 - 5*b^22 - 5*b^20 + b^18 + 13*b^16 + 9*b^14 - 7*b^12 - 6*b^10 - 12*(b^18 + 4*b^16 + 7*b^14 + 8*b^12 + 7*b^10 + 4*b^8 + b^6)*(a^2*b^2 + a)^2*(a*b^2 + a^2)^2 + (12*b^23 + 47*b^21 + 49*b^19 + 89*b^17 + 215*b^15 + 270*b^13 + 297*b^11 + 254*b^9 + 113*b^7 + 67*b^5 + 50*b^3 + 9*b)*(a^2*b^2 + a)^(3/2)*(a*b^2 + a^2)^(3/2) - 2*(b^24 + 8*b^22 - 15*b^20 + 16*b^18 + 62*b^16 - 32*b^14 - 2*b^12 + 80*b^10 + b^8 - 8*b^6 + 17*b^4)*(a^2*b^2 + a)*(a*b^2 + a^2) + (3*b^25 + 21*b^23 + 41*b^21 + 122*b^19 + 136*b^17 - 33*b^15 - 3*b^13 + 100*b^11 - 17*b^9 - 50*b^7)*sqrt(a^2*b^2 + a)*sqrt(a*b^2 + a^2))*a^5 + 36*(b^16 + 4*b^14 + 6*b^12 + 4*b^10 + b^8)*(a^2*b^2 + a)^2*(a*b^2 + a^2)^2 + 2*(b^24 - 3*b^20 + 3*b^16 - b^12 - 36*(b^20 + 4*b^18 + 9*b^16 + 16*b^14 + 20*b^12 + 16*b^10 + 9*b^8 + 4*b^6 + b^4)*(a^2*b^2 + a)^2*(a*b^2 + a^2)^2 + (19*b^21 + 219*b^19 + 507*b^17 + 720*b^15 + 1023*b^13 + 1140*b^11 + 981*b^9 + 732*b^7 + 318*b^5 + 37*b^3)*(a^2*b^2 + a)^(3/2)*(a*b^2 + a^2)^(3/2) - 6*(3*b^22 + 24*b^20 + 25*b^18 + 32*b^16 + 100*b^14 + 88*b^12 + 29*b^10 + 40*b^8 + 27*b^6)*(a^2*b^2 + a)*(a*b^2 + a^2) + 3*(3*b^23 + 21*b^21 + 36*b^19 + 56*b^17 + 91*b^15 + 45*b^13 - 34*b^11 - 26*b^9)*sqrt(a^2*b^2 + a)*sqrt(a*b^2 + a^2))*a^4 - 6*(12*(b^18 + 4*b^16 + 7*b^14 + 8*b^12 + 7*b^10 + 4*b^8 + b^6)*(a^2*b^2 + a)^2*(a*b^2 + a^2)^2 + (3*b^23 + 9*b^21 + 32*b^19 + 45*b^17 - 4*b^15 - 67*b^13 - 100*b^11 - 85*b^9 - 27*b^7 + 2*b^5)*(a^2*b^2 + a)^(3/2)*(a*b^2 + a^2)^(3/2) + 2*(3*b^20 + 24*b^18 + 45*b^16 + 40*b^14 + 45*b^12 + 48*b^10 + 19*b^8)*(a^2*b^2 + a)*(a*b^2 + a^2) - (b^21 + 7*b^19 + 10*b^17 - 2*b^15 - 11*b^13 - 5*b^11)*sqrt(a^2*b^2 + a)*sqrt(a*b^2 + a^2))*a^3 + 6*(6*(b^20 + 4*b^18 + 8*b^16 + 12*b^14 + 14*b^12 + 12*b^10 + 8*b^8 + 4*b^6 + b^4)*(a^2*b^2 + a)^2*(a*b^2 + a^2)^2 - (6*b^21 + 18*b^19 + 35*b^17 + 61*b^15 + 68*b^13 + 44*b^11 + 19*b^9 + 5*b^7)*(a^2*b^2 + a)^(3/2)*(a*b^2 + a^2)^(3/2) - 2*(b^18 + 8*b^16 + 18*b^14 + 16*b^12 + 5*b^10)*(a^2*b^2 + a)*(a*b^2 + a^2))*a^2 + 2*(36*(b^18 + 4*b^16 + 7*b^14 + 8*b^12 + 7*b^10 + 4*b^8 + b^6)*(a^2*b^2 + a)^2*(a*b^2 + a^2)^2 - (9*b^19 + 27*b^17 + 26*b^15 + 6*b^13 - 3*b^11 - b^9)*(a^2*b^2 + a)^(3/2)*(a*b^2 + a^2)^(3/2))*a)/a)^(1/3) + ((b^16 - 2*b^12 + b^8)*a^11 + 2*(b^18 - b^16 - 4*b^14 - 6*b^12 - 3*b^10 + 7*b^8 + 6*b^6)*a^10 + (b^20 - 4*b^18 - 10*b^16 - 28*b^14 - 28*b^12 - 28*b^10 - 46*b^8 - 4*b^6 + 19*b^4 - 2*(b^13 + 7*b^11 + 11*b^9 + 5*b^7)*sqrt(a^2*b^2 + a)*sqrt(a*b^2 + a^2))*a^9 - 2*(b^20 + 2*b^18 + 11*b^16 + 20*b^14 + 56*b^12 + 71*b^10 + 35*b^8 + 38*b^6 + 25*b^4 + (2*b^15 + 15*b^13 + 19*b^11 + 15*b^9 + 23*b^7 + 14*b^5)*sqrt(a^2*b^2 + a)*sqrt(a*b^2 + a^2) - 3*b^2)*a^8 - 2*(4*b^18 + 12*b^16 + 44*b^14 + 56*b^12 + 76*b^10 + 109*b^8 + 52*b^6 + 14*b^4 - 2*(b^10 + 2*b^8 + b^6)*(a^2*b^2 + a)*(a*b^2 + a^2) + (b^17 + 9*b^15 + 4*b^13 + 22*b^11 + 46*b^9 + 16*b^7 + 9*b^5 + 13*b^3)*sqrt(a^2*b^2 + a)*sqrt(a*b^2 + a^2) + 16*b^2 + 1)*a^7 - 2*(b^20 + 2*b^18 + 11*b^16 + 20*b^14 + 56*b^12 + 71*b^10 + 35*b^8 + 38*b^6 + 25*b^4 - 4*(b^12 + b^10 + b^6 + b^4)*(a^2*b^2 + a)*(a*b^2 + a^2) + (b^17 - 5*b^15 + 18*b^13 + 14*b^11 - 10*b^9 + 26*b^7 + 15*b^5 - 7*b^3 + 4*b)*sqrt(a^2*b^2 + a)*sqrt(a*b^2 + a^2) - 3*b^2)*a^6 + (b^20 - 4*b^18 - 10*b^16 - 28*b^14 - 28*b^12 - 28*b^10 - 46*b^8 - 4*b^6 + 6*(b^11 + 3*b^9 + 3*b^7 + b^5)*(a^2*b^2 + a)^(3/2)*(a*b^2 + a^2)^(3/2) + 19*b^4 + 4*(b^14 - 2*b^12 - 2*b^10 + 2*b^8 - 2*b^6 - 2*b^4 + b^2)*(a^2*b^2 + a)*(a*b^2 + a^2) + 2*(b^17 - 5*b^15 + 18*b^13 + 14*b^11 - 10*b^9 + 26*b^7 + 15*b^5 - 7*b^3 + 4*b)*sqrt(a^2*b^2 + a)*sqrt(a*b^2 + a^2))*a^5 + 2*(b^18 - b^16 - 4*b^14 - 6*b^12 - 3*b^10 + 7*b^8 + 6*b^6 + 3*(b^13 + 4*b^11 + 7*b^9 + 7*b^7 + 4*b^5 + b^3)*(a^2*b^2 + a)^(3/2)*(a*b^2 + a^2)^(3/2) - 4*(b^14 + b^10 + 4*b^8 + b^6 + b^2)*(a^2*b^2 + a)*(a*b^2 + a^2) + (b^17 + 9*b^15 + 4*b^13 + 22*b^11 + 46*b^9 + 16*b^7 + 9*b^5 + 13*b^3)*sqrt(a^2*b^2 + a)*sqrt(a*b^2 + a^2))*a^4 - 6*(b^11 + 3*b^9 + 3*b^7 + b^5)*(a^2*b^2 + a)^(3/2)*(a*b^2 + a^2)^(3/2) + (b^16 - 2*b^12 + b^8 + 6*(b^13 + 3*b^11 + 4*b^9 + 4*b^7 + 3*b^5 + b^3)*(a^2*b^2 + a)^(3/2)*(a*b^2 + a^2)^(3/2) + 4*(b^14 - 2*b^12 - 2*b^10 + 2*b^8 - 2*b^6 - 2*b^4 + b^2)*(a^2*b^2 + a)*(a*b^2 + a^2) + 2*(2*b^15 + 15*b^13 + 19*b^11 + 15*b^9 + 23*b^7 + 14*b^5)*sqrt(a^2*b^2 + a)*sqrt(a*b^2 + a^2))*a^3 - 2*(3*(b^13 + 3*b^11 + 4*b^9 + 4*b^7 + 3*b^5 + b^3)*(a^2*b^2 + a)^(3/2)*(a*b^2 + a^2)^(3/2) - 4*(b^12 + b^10 + b^6 + b^4)*(a^2*b^2 + a)*(a*b^2 + a^2) - (b^13 + 7*b^11 + 11*b^9 + 5*b^7)*sqrt(a^2*b^2 + a)*sqrt(a*b^2 + a^2))*a^2 - 2*(3*(b^13 + 4*b^11 + 7*b^9 + 7*b^7 + 4*b^5 + b^3)*(a^2*b^2 + a)^(3/2)*(a*b^2 + a^2)^(3/2) - 2*(b^10 + 2*b^8 + b^6)*(a^2*b^2 + a)*(a*b^2 + a^2))*a)/((1/54*(a*b^4 + (a^2 + 1)*b^2 + a)*(b^2 + 1)^3*sqrt((27*(4*a^21 - 13*a^20 + 8*a^19 + 18*a^18 + 8*a^17 - 13*a^16 + 4*a^15)*b^36 + 3*(16*a^23 + 108*a^22 - 539*a^21 + 710*a^20 + 163*a^19 + 1676*a^18 + 163*a^17 + 710*a^16 - 539*a^15 + 108*a^14 + 16*a^13)*b^34 + 3*(64*a^24 + 44*a^23 - 995*a^22 + 3058*a^21 - 2792*a^20 + 9090*a^19 + 5094*a^18 + 9090*a^17 - 2792*a^16 + 3058*a^15 - 995*a^14 + 44*a^13 + 64*a^12)*b^32 + (288*a^25 - 660*a^24 - 1227*a^23 + 15366*a^22 - 23890*a^21 + 84450*a^20 + 5757*a^19 + 192344*a^18 + 5757*a^17 + 84450*a^16 - 23890*a^15 + 15366*a^14 - 1227*a^13 - 660*a^12 + 288*a^11)*b^30 + 6*(32*a^26 - 192*a^25 + 666*a^24 + 1039*a^23 - 1889*a^22 + 22215*a^21 - 8956*a^20 + 88402*a^19 + 17686*a^18 + 88402*a^17 - 8956*a^16 + 22215*a^15 - 1889*a^14 + 1039*a^13 + 666*a^12 - 192*a^11 + 32*a^10)*b^28 + 3*(16*a^27 - 256*a^26 + 2128*a^25 - 3376*a^24 + 13403*a^23 + 19438*a^22 + 6908*a^21 + 244734*a^20 + 25385*a^19 + 617032*a^18 + 25385*a^17 + 244734*a^16 + 6908*a^15 + 19438*a^14 + 13403*a^13 - 3376*a^12 + 2128*a^11 - 256*a^10 + 16*a^9)*b^26 - (192*a^27 - 3552*a^26 + 13056*a^25 - 72208*a^24 + 102480*a^23 - 332817*a^22 - 263254*a^21 - 417816*a^20 - 2556378*a^19 - 959054*a^18 - 2556378*a^17 - 417816*a^16 - 263254*a^15 - 332817*a^14 + 102480*a^13 - 72208*a^12 + 13056*a^11 - 3552*a^10 + 192*a^9)*b^24 + 3*(224*a^27 - 1920*a^26 + 16096*a^25 - 50504*a^24 + 182206*a^23 - 246940*a^22 + 686309*a^21 - 172838*a^20 + 1652381*a^19 + 452628*a^18 + 1652381*a^17 - 172838*a^16 + 686309*a^15 - 246940*a^14 + 182206*a^13 - 50504*a^12 + 16096*a^11 - 1920*a^10 + 224*a^9)*b^22 - 6*(160*a^27 - 2048*a^26 + 12448*a^25 - 63772*a^24 + 178110*a^23 - 523410*a^22 + 823204*a^21 - 1849465*a^20 + 1350426*a^19 - 3001882*a^18 + 1350426*a^17 - 1849465*a^16 + 823204*a^15 - 523410*a^14 + 178110*a^13 - 63772*a^12 + 12448*a^11 - 2048*a^10 + 160*a^9)*b^20 + (176*a^27 - 12672*a^26 + 104688*a^25 - 496720*a^24 + 1942122*a^23 - 4819512*a^22 + 11422375*a^21 - 15621498*a^20 + 28665279*a^19 - 21364636*a^18 + 28665279*a^17 - 15621498*a^16 + 11422375*a^15 - 4819512*a^14 + 1942122*a^13 - 496720*a^12 + 104688*a^11 - 12672*a^10 + 176*a^9)*b^18 + 432*a^18 + 3*(544*a^26 - 20416*a^25 + 142960*a^24 - 562444*a^23 + 1750931*a^22 - 3594210*a^21 + 7304088*a^20 - 7846386*a^19 + 11951018*a^18 - 7846386*a^17 + 7304088*a^16 - 3594210*a^15 + 1750931*a^14 - 562444*a^13 + 142960*a^12 - 20416*a^11 + 544*a^10)*b^16 + 3*(1184*a^25 - 42812*a^24 + 271963*a^23 - 924070*a^22 + 2450836*a^21 - 3842614*a^20 + 7208433*a^19 - 5663184*a^18 + 7208433*a^17 - 3842614*a^16 + 2450836*a^15 - 924070*a^14 + 271963*a^13 - 42812*a^12 + 1184*a^11)*b^14 - 2*(2098*a^24 + 45435*a^23 - 322221*a^22 + 967015*a^21 - 2580324*a^20 + 2430438*a^19 - 5094706*a^18 + 2430438*a^17 - 2580324*a^16 + 967015*a^15 - 322221*a^14 + 45435*a^13 + 2098*a^12)*b^12 - 3*(6585*a^23 - 15694*a^22 - 12306*a^21 - 2650*a^20 - 523527*a^19 - 138608*a^18 - 523527*a^17 - 2650*a^16 - 12306*a^15 - 15694*a^14 + 6585*a^13)*b^10 - 864*((a^9 + 3*a^8 - 8*a^6 - 6*a^5 + 6*a^4 + 8*a^3 - 3*a - 1)*b^15 + 3*(a^9 + 3*a^8 - 8*a^6 - 6*a^5 + 6*a^4 + 8*a^3 - 3*a - 1)*b^13 + 3*(a^9 + 3*a^8 - 8*a^6 - 6*a^5 + 6*a^4 + 8*a^3 - 3*a - 1)*b^11 + (a^9 + 3*a^8 - 8*a^6 - 6*a^5 + 6*a^4 + 8*a^3 - 3*a - 1)*b^9)*(a^2*b^2 + a)^(9/2)*(a*b^2 + a^2)^(9/2) - 3*(6795*a^22 - 29614*a^21 + 50456*a^20 - 225778*a^19 - 44358*a^18 - 225778*a^17 + 50456*a^16 - 29614*a^15 + 6795*a^14)*b^8 + 864*((a^12 + 3*a^11 - 8*a^9 - 6*a^8 + 6*a^7 + 8*a^6 - 3*a^4 - a^3)*b^19 + (a^13 + 10*a^12 + 10*a^11 - 17*a^10 - 26*a^9 - 8*a^8 + 8*a^7 + 26*a^6 + 17*a^5 - 10*a^4 - 10*a^3 - a^2)*b^17 + (7*a^13 + 21*a^12 + 11*a^11 - 31*a^10 - 54*a^9 - 26*a^8 + 26*a^7 + 54*a^6 + 31*a^5 - 11*a^4 - 21*a^3 - 7*a^2)*b^15 + (11*a^13 + 25*a^12 + 11*a^11 - 31*a^10 - 74*a^9 - 46*a^8 + 46*a^7 + 74*a^6 + 31*a^5 - 11*a^4 - 25*a^3 - 11*a^2)*b^13 + (5*a^13 + 22*a^12 + 18*a^11 - 33*a^10 - 62*a^9 - 28*a^8 + 28*a^7 + 62*a^6 + 33*a^5 - 18*a^4 - 22*a^3 - 5*a^2)*b^11 + 3*(3*a^12 + 5*a^11 - 4*a^10 - 12*a^9 - 6*a^8 + 6*a^7 + 12*a^6 + 4*a^5 - 5*a^4 - 3*a^3)*b^9 + 4*(a^11 + a^10 - 3*a^9 - 3*a^8 + 3*a^7 + 3*a^6 - a^5 - a^4)*b^7)*(a^2*b^2 + a)^(7/2)*(a*b^2 + a^2)^(7/2) - 4*(27*(a^14 + 4*a^12 - 40*a^11 - 32*a^10 + 134*a^9 - 32*a^8 - 40*a^7 + 4*a^6 + a^4)*b^22 + 9*(3*a^15 - 79*a^13 - 54*a^12 - 206*a^11 - 624*a^10 + 1920*a^9 - 624*a^8 - 206*a^7 - 54*a^6 - 79*a^5 + 3*a^3)*b^20 - (585*a^14 + 486*a^13 + 5491*a^12 + 4908*a^11 + 4425*a^10 - 31790*a^9 + 4425*a^8 + 4908*a^7 + 5491*a^6 + 486*a^5 + 585*a^4)*b^18 + 3*(78*a^15 - 72*a^14 - 1357*a^13 - 834*a^12 - 9634*a^11 + 3314*a^10 + 17010*a^9 + 3314*a^8 - 9634*a^7 - 834*a^6 - 1357*a^5 - 72*a^4 + 78*a^3)*b^16 + 3*(288*a^15 + 323*a^14 - 1116*a^13 - 4005*a^12 - 11692*a^11 + 1964*a^10 + 28476*a^9 + 1964*a^8 - 11692*a^7 - 4005*a^6 - 1116*a^5 + 323*a^4 + 288*a^3)*b^14 + (539*a^15 + 3480*a^14 - 3597*a^13 - 8026*a^12 - 54141*a^11 + 35070*a^10 + 53350*a^9 + 35070*a^8 - 54141*a^7 - 8026*a^6 - 3597*a^5 + 3480*a^4 + 539*a^3)*b^12 + 3*(695*a^14 + 222*a^13 - 1851*a^12 - 11164*a^11 - 562*a^10 + 25320*a^9 - 562*a^8 - 11164*a^7 - 1851*a^6 + 222*a^5 + 695*a^4)*b^10 + 3*(329*a^13 + 834*a^12 - 8212*a^11 + 1574*a^10 + 10950*a^9 + 1574*a^8 - 8212*a^7 + 834*a^6 + 329*a^5)*b^8 + (359*a^12 - 3108*a^11 - 10860*a^10 + 27218*a^9 - 10860*a^8 - 3108*a^7 + 359*a^6)*b^6 + 27*(7*a^11 - 226*a^10 + 438*a^9 - 226*a^8 + 7*a^7)*b^4 - 729*(a^10 - 2*a^9 + a^8)*b^2)*(a^2*b^2 + a)^3*(a*b^2 + a^2)^3 - (8089*a^21 - 35922*a^20 - 17817*a^19 - 261212*a^18 - 17817*a^17 - 35922*a^16 + 8089*a^15)*b^6 - 288*((a^15 + 3*a^14 - 8*a^12 - 6*a^11 + 6*a^10 + 8*a^9 - 3*a^7 - a^6)*b^23 + (2*a^16 + 17*a^15 + 11*a^14 - 34*a^13 - 28*a^12 + 2*a^11 - 2*a^10 + 28*a^9 + 34*a^8 - 11*a^7 - 17*a^6 - 2*a^5)*b^21 + (a^17 + 25*a^16 + 57*a^15 - 53*a^14 - 11*a^13 - 75*a^12 - 150*a^11 + 150*a^10 + 75*a^9 + 11*a^8 + 53*a^7 - 57*a^6 - 25*a^5 - a^4)*b^19 + (11*a^17 + 81*a^16 - 11*a^15 + 111*a^14 - 207*a^13 - 269*a^12 + 130*a^11 - 130*a^10 + 269*a^9 + 207*a^8 - 111*a^7 + 11*a^6 - 81*a^5 - 11*a^4)*b^17 + (35*a^17 + 41*a^16 + 173*a^15 - 121*a^14 - 135*a^13 - 197*a^12 - 350*a^11 + 350*a^10 + 197*a^9 + 135*a^8 + 121*a^7 - 173*a^6 - 41*a^5 - 35*a^4)*b^15 + (25*a^17 + 73*a^16 + 105*a^15 - 5*a^14 - 203*a^13 - 291*a^12 - 150*a^11 + 150*a^10 + 291*a^9 + 203*a^8 + 5*a^7 - 105*a^6 - 73*a^5 - 25*a^4)*b^13 + (90*a^16 + a^15 + 123*a^14 - 202*a^13 - 316*a^12 + 98*a^11 - 98*a^10 + 316*a^9 + 202*a^8 - 123*a^7 - a^6 - 90*a^5)*b^11 + (121*a^15 - 93*a^14 + 8*a^13 - 64*a^12 - 286*a^11 + 286*a^10 + 64*a^9 - 8*a^8 + 93*a^7 - 121*a^6)*b^9 + 72*(a^14 - a^13 - a^12 + a^11 - a^10 + a^9 + a^8 - a^7)*b^7 + 16*(a^13 - a^12 - 2*a^11 + 2*a^10 + a^9 - a^8)*b^5)*(a^2*b^2 + a)^(5/2)*(a*b^2 + a^2)^(5/2) + 24*(3*(a^17 + 13*a^15 - 49*a^14 - 68*a^13 + 206*a^12 - 68*a^11 - 49*a^10 + 13*a^9 + a^7)*b^26 + (6*a^18 + 12*a^17 - 146*a^16 + 192*a^15 - 655*a^14 - 1368*a^13 + 3918*a^12 - 1368*a^11 - 655*a^10 + 192*a^9 - 146*a^8 + 12*a^7 + 6*a^6)*b^24 + (3*a^19 + 24*a^18 - 387*a^17 + 409*a^16 - 1059*a^15 - 1981*a^14 - 1149*a^13 + 8280*a^12 - 1149*a^11 - 1981*a^10 - 1059*a^9 + 409*a^8 - 387*a^7 + 24*a^6 + 3*a^5)*b^22 + (12*a^19 - 180*a^18 - 102*a^17 - 1433*a^16 - 166*a^15 - 8951*a^14 + 5308*a^13 + 11024*a^12 + 5308*a^11 - 8951*a^10 - 166*a^9 - 1433*a^8 - 102*a^7 - 180*a^6 + 12*a^5)*b^20 + (22*a^19 + 72*a^18 - 994*a^17 - 2247*a^16 - 1574*a^15 - 20906*a^14 + 25232*a^13 + 790*a^12 + 25232*a^11 - 20906*a^10 - 1574*a^9 - 2247*a^8 - 994*a^7 + 72*a^6 + 22*a^5)*b^18 + (244*a^19 + 126*a^18 - 1372*a^17 + 39*a^16 - 21516*a^15 + 13852*a^14 - 31492*a^13 + 80238*a^12 - 31492*a^11 + 13852*a^10 - 21516*a^9 + 39*a^8 - 1372*a^7 + 126*a^6 + 244*a^5)*b^16 + (295*a^19 + 704*a^18 - 323*a^17 - 7661*a^16 + 851*a^15 - 42446*a^14 + 58937*a^13 - 20714*a^12 + 58937*a^11 - 42446*a^10 + 851*a^9 - 7661*a^8 - 323*a^7 + 704*a^6 + 295*a^5)*b^14 + (1744*a^18 - 2858*a^17 + 3957*a^16 - 27026*a^15 + 19534*a^14 - 35276*a^13 + 79850*a^12 - 35276*a^11 + 19534*a^10 - 27026*a^9 + 3957*a^8 - 2858*a^7 + 1744*a^6)*b^12 + (3205*a^17 - 9381*a^16 + 9561*a^15 - 33611*a^14 + 40156*a^13 - 19860*a^12 + 40156*a^11 - 33611*a^10 + 9561*a^9 - 9381*a^8 + 3205*a^7)*b^10 + (1999*a^16 - 7900*a^15 + 4687*a^14 - 8672*a^13 + 19772*a^12 - 8672*a^11 + 4687*a^10 - 7900*a^9 + 1999*a^8)*b^8 - (170*a^15 + 1037*a^14 + 1318*a^13 - 5050*a^12 + 1318*a^11 + 1037*a^10 + 170*a^9)*b^6 - (467*a^14 - 908*a^13 + 882*a^12 - 908*a^11 + 467*a^10)*b^4 - 54*(a^13 - 2*a^12 + a^11)*b^2)*(a^2*b^2 + a)^2*(a*b^2 + a^2)^2 - 3*(445*a^20 - 4460*a^19 - 14002*a^18 - 4460*a^17 + 445*a^16)*b^4 - 12*(9*(a^17 - 16*a^16 + 15*a^15 + 32*a^14 - 32*a^13 - 15*a^12 + 16*a^11 - a^10)*b^29 - (50*a^18 + 309*a^17 - 225*a^16 - 412*a^15 - 900*a^14 + 900*a^13 + 412*a^12 + 225*a^11 - 309*a^10 - 50*a^9)*b^27 - (159*a^19 + 270*a^18 + 118*a^17 + 495*a^16 - 2666*a^15 - 4630*a^14 + 4630*a^13 + 2666*a^12 - 495*a^11 - 118*a^10 - 270*a^9 - 159*a^8)*b^25 - (132*a^20 + 189*a^19 + 1387*a^18 - 805*a^17 + 185*a^16 - 5546*a^15 - 19266*a^14 + 19266*a^13 + 5546*a^12 - 185*a^11 + 805*a^10 - 1387*a^9 - 189*a^8 - 132*a^7)*b^23 - (32*a^21 + 84*a^20 + 2355*a^19 - 3247*a^18 + 8341*a^17 - 9058*a^16 - 25898*a^15 - 10877*a^14 + 10877*a^13 + 25898*a^12 + 9058*a^11 - 8341*a^10 + 3247*a^9 - 2355*a^8 - 84*a^7 - 32*a^6)*b^21 - (1336*a^20 - 1575*a^19 + 11311*a^18 - 18457*a^17 + 19637*a^16 - 83145*a^15 + 3575*a^14 - 3575*a^13 + 83145*a^12 - 19637*a^11 + 18457*a^10 - 11311*a^9 + 1575*a^8 - 1336*a^7)*b^19 - (160*a^21 + 440*a^20 + 4253*a^19 - 3755*a^18 + 19035*a^17 - 53273*a^16 - 26762*a^15 - 101138*a^14 + 101138*a^13 + 26762*a^12 + 53273*a^11 - 19035*a^10 + 3755*a^9 - 4253*a^8 - 440*a^7 - 160*a^6)*b^17 - (192*a^21 - 380*a^20 + 6359*a^19 - 4359*a^18 + 18139*a^17 - 45571*a^16 - 82914*a^15 - 49690*a^14 + 49690*a^13 + 82914*a^12 + 45571*a^11 - 18139*a^10 + 4359*a^9 - 6359*a^8 + 380*a^7 - 192*a^6)*b^15 - (564*a^20 - 1199*a^19 + 16555*a^18 - 16672*a^17 - 4615*a^16 - 106989*a^15 - 37604*a^14 + 37604*a^13 + 106989*a^12 + 4615*a^11 + 16672*a^10 - 16555*a^9 + 1199*a^8 - 564*a^7)*b^13 - (595*a^19 + 3707*a^18 + 9024*a^17 - 29528*a^16 - 56256*a^15 - 44068*a^14 + 44068*a^13 + 56256*a^12 + 29528*a^11 - 9024*a^10 - 3707*a^9 - 595*a^8)*b^11 - (2273*a^18 + 1795*a^17 - 3850*a^16 - 49156*a^15 - 5300*a^14 + 5300*a^13 + 49156*a^12 + 3850*a^11 - 1795*a^10 - 2273*a^9)*b^9 - 2*(2199*a^17 - 6073*a^16 + 2306*a^15 - 19090*a^14 + 19090*a^13 - 2306*a^12 + 6073*a^11 - 2199*a^10)*b^7 - (2509*a^16 - 12074*a^15 + 9169*a^14 - 9169*a^13 + 12074*a^12 - 2509*a^11)*b^5 + (19*a^15 + 2319*a^14 - 2319*a^13 - 19*a^12)*b^3 + 180*(a^14 - a^13)*b)*(a^2*b^2 + a)^(3/2)*(a*b^2 + a^2)^(3/2) - 12*((a^20 + 22*a^18 - 58*a^17 - 131*a^16 + 332*a^15 - 131*a^14 - 58*a^13 + 22*a^12 + a^10)*b^30 + (3*a^21 + 8*a^20 - 53*a^19 + 200*a^18 - 591*a^17 - 856*a^16 + 2578*a^15 - 856*a^14 - 591*a^13 + 200*a^12 - 53*a^11 + 8*a^10 + 3*a^9)*b^28 + (3*a^22 + 24*a^21 - 269*a^20 + 356*a^19 + 391*a^18 - 4380*a^17 + 1279*a^16 + 5192*a^15 + 1279*a^14 - 4380*a^13 + 391*a^12 + 356*a^11 - 269*a^10 + 24*a^9 + 3*a^8)*b^26 + (a^23 + 24*a^22 - 269*a^21 - 592*a^20 + 3135*a^19 - 9072*a^18 + 4761*a^17 - 15272*a^16 + 34568*a^15 - 15272*a^14 + 4761*a^13 - 9072*a^12 + 3135*a^11 - 592*a^10 - 269*a^9 + 24*a^8 + a^7)*b^24 + (8*a^23 - 53*a^22 - 1042*a^21 + 3539*a^20 - 15292*a^19 + 33380*a^18 - 93540*a^17 + 129394*a^16 - 112788*a^15 + 129394*a^14 - 93540*a^13 + 33380*a^12 - 15292*a^11 + 3539*a^10 - 1042*a^9 - 53*a^8 + 8*a^7)*b^22 + (22*a^23 - 232*a^22 + 1831*a^21 - 15432*a^20 + 46088*a^19 - 137144*a^18 + 216777*a^17 - 320656*a^16 + 417492*a^15 - 320656*a^14 + 216777*a^13 - 137144*a^12 + 46088*a^11 - 15432*a^10 + 1831*a^9 - 232*a^8 + 22*a^7)*b^20 + (120*a^23 + 1001*a^22 - 6836*a^21 + 24215*a^20 - 97144*a^19 + 196257*a^18 - 407140*a^17 + 580123*a^16 - 581192*a^15 + 580123*a^14 - 407140*a^13 + 196257*a^12 - 97144*a^11 + 24215*a^10 - 6836*a^9 + 1001*a^8 + 120*a^7)*b^18 + (425*a^23 - 816*a^22 + 8461*a^21 - 43976*a^20 + 107728*a^19 - 293864*a^18 + 479963*a^17 - 671616*a^16 + 827390*a^15 - 671616*a^14 + 479963*a^13 - 293864*a^12 + 107728*a^11 - 43976*a^10 + 8461*a^9 - 816*a^8 + 425*a^7)*b^16 + (3145*a^22 - 11698*a^21 + 40352*a^20 - 144364*a^19 + 280241*a^18 - 531348*a^17 + 750850*a^16 - 774356*a^15 + 750850*a^14 - 531348*a^13 + 280241*a^12 - 144364*a^11 + 40352*a^10 - 11698*a^9 + 3145*a^8)*b^14 + (8822*a^21 - 39080*a^20 + 100277*a^19 - 264296*a^18 + 441166*a^17 - 611304*a^16 + 728830*a^15 - 611304*a^14 + 441166*a^13 - 264296*a^12 + 100277*a^11 - 39080*a^10 + 8822*a^9)*b^12 + 2*(6145*a^20 - 29954*a^19 + 66630*a^18 - 135538*a^17 + 198051*a^16 - 210668*a^15 + 198051*a^14 - 135538*a^13 + 66630*a^12 - 29954*a^11 + 6145*a^10)*b^10 + (8873*a^19 - 47184*a^18 + 94988*a^17 - 148224*a^16 + 183094*a^15 - 148224*a^14 + 94988*a^13 - 47184*a^12 + 8873*a^11)*b^8 + (2849*a^18 - 17930*a^17 + 33339*a^16 - 36516*a^15 + 33339*a^14 - 17930*a^13 + 2849*a^12)*b^6 + 8*(3*a^17 - 257*a^16 + 508*a^15 - 257*a^14 + 3*a^13)*b^4 - 140*(a^16 - 2*a^15 + a^14)*b^2)*(a^2*b^2 + a)*(a*b^2 + a^2) + 24*(11*a^19 + 302*a^18 + 11*a^17)*b^2 + 12*(9*(a^20 - 8*a^19 + 8*a^18 + 17*a^17 - 17*a^16 - 8*a^15 + 8*a^14 - a^13)*b^33 - (17*a^21 + 93*a^20 + 144*a^19 - 706*a^18 - 558*a^17 + 558*a^16 + 706*a^15 - 144*a^14 - 93*a^13 - 17*a^12)*b^31 - (113*a^22 - 117*a^21 + 1215*a^20 - 2152*a^19 + 31*a^18 - 6668*a^17 + 6668*a^16 - 31*a^15 + 2152*a^14 - 1215*a^13 + 117*a^12 - 113*a^11)*b^29 - (147*a^23 - 165*a^22 + 2533*a^21 - 5497*a^20 + 9331*a^19 - 25627*a^18 + 5492*a^17 - 5492*a^16 + 25627*a^15 - 9331*a^14 + 5497*a^13 - 2533*a^12 + 165*a^11 - 147*a^10)*b^27 - (68*a^24 + 57*a^23 + 2725*a^22 - 7472*a^21 + 21905*a^20 - 49597*a^19 + 35863*a^18 - 108467*a^17 + 108467*a^16 - 35863*a^15 + 49597*a^14 - 21905*a^13 + 7472*a^12 - 2725*a^11 - 57*a^10 - 68*a^9)*b^25 - (8*a^25 + 108*a^24 + 1455*a^23 - 3697*a^22 + 20850*a^21 - 59416*a^20 + 112503*a^19 - 296579*a^18 + 260984*a^17 - 260984*a^16 + 296579*a^15 - 112503*a^14 + 59416*a^13 - 20850*a^12 + 3697*a^11 - 1455*a^10 - 108*a^9 - 8*a^8)*b^23 - (24*a^25 + 160*a^24 + 509*a^23 + 7963*a^22 - 28378*a^21 + 94606*a^20 - 266818*a^19 + 259848*a^18 - 525888*a^17 + 525888*a^16 - 259848*a^15 + 266818*a^14 - 94606*a^13 + 28378*a^12 - 7963*a^11 - 509*a^10 - 160*a^9 - 24*a^8)*b^21 + (104*a^25 - 512*a^24 + 343*a^23 - 6737*a^22 - 2718*a^21 + 23998*a^20 + 22723*a^19 + 237097*a^18 - 69190*a^17 + 69190*a^16 - 237097*a^15 - 22723*a^14 - 23998*a^13 + 2718*a^12 + 6737*a^11 - 343*a^10 + 512*a^9 - 104*a^8)*b^19 + (120*a^25 + 772*a^24 - 5659*a^23 + 15113*a^22 - 70076*a^21 + 169859*a^20 - 201439*a^19 + 571789*a^18 - 446157*a^17 + 446157*a^16 - 571789*a^15 + 201439*a^14 - 169859*a^13 + 70076*a^12 - 15113*a^11 + 5659*a^10 - 772*a^9 - 120*a^8)*b^17 + (1164*a^24 - 533*a^23 - 16549*a^22 + 47611*a^21 - 147949*a^20 + 428479*a^19 - 442011*a^18 + 785970*a^17 - 785970*a^16 + 442011*a^15 - 428479*a^14 + 147949*a^13 - 47611*a^12 + 16549*a^11 + 533*a^10 - 1164*a^9)*b^15 + (4433*a^23 - 15576*a^22 + 6113*a^21 + 19053*a^20 - 19585*a^19 + 267808*a^18 - 203604*a^17 + 203604*a^16 - 267808*a^15 + 19585*a^14 - 19053*a^13 - 6113*a^12 + 15576*a^11 - 4433*a^10)*b^13 + (8160*a^22 - 41973*a^21 + 81425*a^20 - 95468*a^19 + 222374*a^18 - 157264*a^17 + 157264*a^16 - 222374*a^15 + 95468*a^14 - 81425*a^13 + 41973*a^12 - 8160*a^11)*b^11 + (7516*a^21 - 45851*a^20 + 115790*a^19 - 139922*a^18 + 225169*a^17 - 225169*a^16 + 139922*a^15 - 115790*a^14 + 45851*a^13 - 7516*a^12)*b^9 + 2*(1525*a^20 - 10126*a^19 + 31409*a^18 - 34112*a^17 + 34112*a^16 - 31409*a^15 + 10126*a^14 - 1525*a^13)*b^7 + (243*a^19 - 805*a^18 + 9840*a^17 - 9840*a^16 + 805*a^15 - 243*a^14)*b^5 + 2*(5*a^18 + 561*a^17 - 561*a^16 - 5*a^15)*b^3 + 72*(a^17 - a^16)*b)*sqrt(a^2*b^2 + a)*sqrt(a*b^2 + a^2))*(a*b^4 + (a^2 + 1)*b^2 + a)/a)/a - 1/54*(2*(b^24 - 3*b^20 + 3*b^16 - b^12)*a^16 + 6*(b^26 - b^24 - 5*b^22 - 5*b^20 + b^18 + 13*b^16 + 9*b^14 - 7*b^12 - 6*b^10)*a^15 + 3*(2*b^28 - 6*b^26 - 17*b^24 - 30*b^22 - 38*b^20 - 82*b^18 - 124*b^16 + 54*b^14 + 196*b^12 + 64*b^10 - 19*b^8 - 2*(b^21 + 7*b^19 + 10*b^17 - 2*b^15 - 11*b^13 - 5*b^11)*sqrt(a^2*b^2 + a)*sqrt(a*b^2 + a^2))*a^14 + (2*b^30 - 18*b^28 - 33*b^26 - 68*b^24 - 129*b^22 - 282*b^20 + 792*b^18 + 2922*b^16 + 2682*b^14 + 1636*b^12 + 1761*b^10 + 930*b^8 + 45*b^6 - 6*(3*b^23 + 21*b^21 + 36*b^19 + 56*b^17 + 91*b^15 + 45*b^13 - 34*b^11 - 26*b^9)*sqrt(a^2*b^2 + a)*sqrt(a*b^2 + a^2))*a^13 - 3*(2*b^30 + b^28 - 12*b^26 - 90*b^24 - 498*b^22 - 1893*b^20 - 3232*b^18 - 3598*b^16 - 4314*b^14 - 3843*b^12 - 1780*b^10 - 768*b^8 - 406*b^6 - 49*b^4 + 4*(b^18 + 8*b^16 + 18*b^14 + 16*b^12 + 5*b^10)*(a^2*b^2 + a)*(a*b^2 + a^2) + 2*(3*b^25 + 21*b^23 + 41*b^21 + 122*b^19 + 136*b^17 - 33*b^15 - 3*b^13 + 100*b^11 - 17*b^9 - 50*b^7)*sqrt(a^2*b^2 + a)*sqrt(a*b^2 + a^2))*a^12 + 3*(b^30 + 22*b^28 + 160*b^26 + 846*b^24 + 2568*b^22 + 4632*b^20 + 7705*b^18 + 10840*b^16 + 9717*b^14 + 6834*b^12 + 4902*b^10 + 2218*b^8 + 514*b^6 + 208*b^4 - 4*(3*b^20 + 24*b^18 + 45*b^16 + 40*b^14 + 45*b^12 + 48*b^10 + 19*b^8)*(a^2*b^2 + a)*(a*b^2 + a^2) - 2*(b^27 + 7*b^25 + 7*b^23 + 17*b^21 - 211*b^19 - 378*b^17 - 98*b^15 - 227*b^13 - 384*b^11 - 15*b^9 + 45*b^7 - 44*b^5)*sqrt(a^2*b^2 + a)*sqrt(a*b^2 + a^2) + 33*b^2)*a^11 + 2*(11*b^30 + 99*b^28 + 579*b^26 + 1865*b^24 + 4701*b^22 + 10887*b^20 + 16687*b^18 + 18402*b^16 + 19359*b^14 + 15839*b^12 + 8067*b^10 + 3819*b^8 + 1745*b^6 + (9*b^19 + 27*b^17 + 26*b^15 + 6*b^13 - 3*b^11 - b^9)*(a^2*b^2 + a)^(3/2)*(a*b^2 + a^2)^(3/2) + 279*b^4 - 6*(3*b^22 + 24*b^20 + 25*b^18 + 32*b^16 + 100*b^14 + 88*b^12 + 29*b^10 + 40*b^8 + 27*b^6)*(a^2*b^2 + a)*(a*b^2 + a^2) + 3*(15*b^25 + 91*b^23 + 374*b^21 + 393*b^19 + 504*b^17 + 1164*b^15 + 732*b^13 + 90*b^11 + 357*b^9 + 165*b^7 - 62*b^5 + 17*b^3)*sqrt(a^2*b^2 + a)*sqrt(a*b^2 + a^2) + 51*b^2 + 10)*a^10 + 3*(b^30 + 22*b^28 + 160*b^26 + 846*b^24 + 2568*b^22 + 4632*b^20 + 7705*b^18 + 10840*b^16 + 9717*b^14 + 6834*b^12 + 4902*b^10 + 2218*b^8 + 514*b^6 + 2*(6*b^21 + 18*b^19 + 35*b^17 + 61*b^15 + 68*b^13 + 44*b^11 + 19*b^9 + 5*b^7)*(a^2*b^2 + a)^(3/2)*(a*b^2 + a^2)^(3/2) + 208*b^4 - 4*(b^24 + 8*b^22 - 15*b^20 + 16*b^18 + 62*b^16 - 32*b^14 - 2*b^12 + 80*b^10 + b^8 - 8*b^6 + 17*b^4)*(a^2*b^2 + a)*(a*b^2 + a^2) + 2*(7*b^27 + 37*b^25 + 92*b^23 + 22*b^21 + 451*b^19 + 738*b^17 + 16*b^15 + 364*b^13 + 723*b^11 - 69*b^9 - 48*b^7 + 122*b^5 - 25*b^3 + 2*b)*sqrt(a^2*b^2 + a)*sqrt(a*b^2 + a^2) + 33*b^2)*a^9 - 3*(2*b^30 + b^28 - 12*b^26 - 90*b^24 - 498*b^22 - 1893*b^20 - 3232*b^18 - 3598*b^16 - 4314*b^14 - 3843*b^12 - 1780*b^10 - 768*b^8 - 406*b^6 - 12*(b^16 + 4*b^14 + 6*b^12 + 4*b^10 + b^8)*(a^2*b^2 + a)^2*(a*b^2 + a^2)^2 - 2*(3*b^23 + 9*b^21 + 32*b^19 + 45*b^17 - 4*b^15 - 67*b^13 - 100*b^11 - 85*b^9 - 27*b^7 + 2*b^5)*(a^2*b^2 + a)^(3/2)*(a*b^2 + a^2)^(3/2) - 49*b^4 - 4*(15*b^22 - 16*b^20 + 38*b^18 + 136*b^16 + 14*b^14 + 56*b^12 + 154*b^10 + 8*b^8 - 9*b^6 + 24*b^4 - 4*b^2)*(a^2*b^2 + a)*(a*b^2 + a^2) + 2*(7*b^27 + 37*b^25 + 92*b^23 + 22*b^21 + 451*b^19 + 738*b^17 + 16*b^15 + 364*b^13 + 723*b^11 - 69*b^9 - 48*b^7 + 122*b^5 - 25*b^3 + 2*b)*sqrt(a^2*b^2 + a)*sqrt(a*b^2 + a^2))*a^8 + (2*b^30 - 18*b^28 - 33*b^26 - 68*b^24 - 129*b^22 - 282*b^20 + 792*b^18 + 2922*b^16 + 2682*b^14 + 1636*b^12 + 1761*b^10 + 930*b^8 + 45*b^6 + 72*(b^18 + 4*b^16 + 7*b^14 + 8*b^12 + 7*b^10 + 4*b^8 + b^6)*(a^2*b^2 + a)^2*(a*b^2 + a^2)^2 - 2*(19*b^21 + 219*b^19 + 507*b^17 + 720*b^15 + 1023*b^13 + 1140*b^11 + 981*b^9 + 732*b^7 + 318*b^5 + 37*b^3)*(a^2*b^2 + a)^(3/2)*(a*b^2 + a^2)^(3/2) + 24*(b^24 - 4*b^22 + 28*b^20 + 28*b^18 + 11*b^16 + 112*b^14 + 91*b^12 + 8*b^10 + 52*b^8 + 28*b^6 - 7*b^4 + 4*b^2)*(a^2*b^2 + a)*(a*b^2 + a^2) - 6*(15*b^25 + 91*b^23 + 374*b^21 + 393*b^19 + 504*b^17 + 1164*b^15 + 732*b^13 + 90*b^11 + 357*b^9 + 165*b^7 - 62*b^5 + 17*b^3)*sqrt(a^2*b^2 + a)*sqrt(a*b^2 + a^2))*a^7 + 3*(2*b^28 - 6*b^26 - 17*b^24 - 30*b^22 - 38*b^20 - 82*b^18 - 124*b^16 + 54*b^14 + 196*b^12 + 64*b^10 - 19*b^8 + 12*(b^20 + 4*b^18 + 8*b^16 + 12*b^14 + 14*b^12 + 12*b^10 + 8*b^8 + 4*b^6 + b^4)*(a^2*b^2 + a)^2*(a*b^2 + a^2)^2 - 2*(12*b^23 + 47*b^21 + 49*b^19 + 89*b^17 + 215*b^15 + 270*b^13 + 297*b^11 + 254*b^9 + 113*b^7 + 67*b^5 + 50*b^3 + 9*b)*(a^2*b^2 + a)^(3/2)*(a*b^2 + a^2)^(3/2) + 4*(15*b^22 - 16*b^20 + 38*b^18 + 136*b^16 + 14*b^14 + 56*b^12 + 154*b^10 + 8*b^8 - 9*b^6 + 24*b^4 - 4*b^2)*(a^2*b^2 + a)*(a*b^2 + a^2) + 2*(b^27 + 7*b^25 + 7*b^23 + 17*b^21 - 211*b^19 - 378*b^17 - 98*b^15 - 227*b^13 - 384*b^11 - 15*b^9 + 45*b^7 - 44*b^5)*sqrt(a^2*b^2 + a)*sqrt(a*b^2 + a^2))*a^6 + 6*(b^26 - b^24 - 5*b^22 - 5*b^20 + b^18 + 13*b^16 + 9*b^14 - 7*b^12 - 6*b^10 - 12*(b^18 + 4*b^16 + 7*b^14 + 8*b^12 + 7*b^10 + 4*b^8 + b^6)*(a^2*b^2 + a)^2*(a*b^2 + a^2)^2 + (12*b^23 + 47*b^21 + 49*b^19 + 89*b^17 + 215*b^15 + 270*b^13 + 297*b^11 + 254*b^9 + 113*b^7 + 67*b^5 + 50*b^3 + 9*b)*(a^2*b^2 + a)^(3/2)*(a*b^2 + a^2)^(3/2) - 2*(b^24 + 8*b^22 - 15*b^20 + 16*b^18 + 62*b^16 - 32*b^14 - 2*b^12 + 80*b^10 + b^8 - 8*b^6 + 17*b^4)*(a^2*b^2 + a)*(a*b^2 + a^2) + (3*b^25 + 21*b^23 + 41*b^21 + 122*b^19 + 136*b^17 - 33*b^15 - 3*b^13 + 100*b^11 - 17*b^9 - 50*b^7)*sqrt(a^2*b^2 + a)*sqrt(a*b^2 + a^2))*a^5 + 36*(b^16 + 4*b^14 + 6*b^12 + 4*b^10 + b^8)*(a^2*b^2 + a)^2*(a*b^2 + a^2)^2 + 2*(b^24 - 3*b^20 + 3*b^16 - b^12 - 36*(b^20 + 4*b^18 + 9*b^16 + 16*b^14 + 20*b^12 + 16*b^10 + 9*b^8 + 4*b^6 + b^4)*(a^2*b^2 + a)^2*(a*b^2 + a^2)^2 + (19*b^21 + 219*b^19 + 507*b^17 + 720*b^15 + 1023*b^13 + 1140*b^11 + 981*b^9 + 732*b^7 + 318*b^5 + 37*b^3)*(a^2*b^2 + a)^(3/2)*(a*b^2 + a^2)^(3/2) - 6*(3*b^22 + 24*b^20 + 25*b^18 + 32*b^16 + 100*b^14 + 88*b^12 + 29*b^10 + 40*b^8 + 27*b^6)*(a^2*b^2 + a)*(a*b^2 + a^2) + 3*(3*b^23 + 21*b^21 + 36*b^19 + 56*b^17 + 91*b^15 + 45*b^13 - 34*b^11 - 26*b^9)*sqrt(a^2*b^2 + a)*sqrt(a*b^2 + a^2))*a^4 - 6*(12*(b^18 + 4*b^16 + 7*b^14 + 8*b^12 + 7*b^10 + 4*b^8 + b^6)*(a^2*b^2 + a)^2*(a*b^2 + a^2)^2 + (3*b^23 + 9*b^21 + 32*b^19 + 45*b^17 - 4*b^15 - 67*b^13 - 100*b^11 - 85*b^9 - 27*b^7 + 2*b^5)*(a^2*b^2 + a)^(3/2)*(a*b^2 + a^2)^(3/2) + 2*(3*b^20 + 24*b^18 + 45*b^16 + 40*b^14 + 45*b^12 + 48*b^10 + 19*b^8)*(a^2*b^2 + a)*(a*b^2 + a^2) - (b^21 + 7*b^19 + 10*b^17 - 2*b^15 - 11*b^13 - 5*b^11)*sqrt(a^2*b^2 + a)*sqrt(a*b^2 + a^2))*a^3 + 6*(6*(b^20 + 4*b^18 + 8*b^16 + 12*b^14 + 14*b^12 + 12*b^10 + 8*b^8 + 4*b^6 + b^4)*(a^2*b^2 + a)^2*(a*b^2 + a^2)^2 - (6*b^21 + 18*b^19 + 35*b^17 + 61*b^15 + 68*b^13 + 44*b^11 + 19*b^9 + 5*b^7)*(a^2*b^2 + a)^(3/2)*(a*b^2 + a^2)^(3/2) - 2*(b^18 + 8*b^16 + 18*b^14 + 16*b^12 + 5*b^10)*(a^2*b^2 + a)*(a*b^2 + a^2))*a^2 + 2*(36*(b^18 + 4*b^16 + 7*b^14 + 8*b^12 + 7*b^10 + 4*b^8 + b^6)*(a^2*b^2 + a)^2*(a*b^2 + a^2)^2 - (9*b^19 + 27*b^17 + 26*b^15 + 6*b^13 - 3*b^11 - b^9)*(a^2*b^2 + a)^(3/2)*(a*b^2 + a^2)^(3/2))*a)/a)^(1/3)*a))*(3*(b^8 - b^4)*a^5 - 9*(a^2*b^2 + a)*(a*b^2 + a^2)*(b^2 + 1)^3 + 3*(b^10 - b^8 + 3*b^6 + 5*b^4)*a^4 + 9*(b^2 + 1)*(sqrt(a*b^2 + a^2)*(a*b^2 + 1) - sqrt(a^2*b^2 + a)*(b^3 + a*b))^2*a - 3*(b^10 - 5*b^8 - 4*b^6 - 2*b^4 - 2*(b^5 + b^3)*sqrt(a^2*b^2 + a)*sqrt(a*b^2 + a^2) - 5*b^2 - 1)*a^3 + 9*(b^2 + 1)*((a*b^3 + b)*sqrt(a*b^2 + a^2) + sqrt(a^2*b^2 + a)*(b^2 + a))^2 + 3*(b^10 - b^8 + 3*b^6 + 5*b^4 + 2*(b^7 + b)*sqrt(a^2*b^2 + a)*sqrt(a*b^2 + a^2))*a^2 - 6*(b^5 + b^3)*sqrt(a^2*b^2 + a)*sqrt(a*b^2 + a^2) + 3*(b^8 - b^4 - 2*(b^7 + b)*sqrt(a^2*b^2 + a)*sqrt(a*b^2 + a^2))*a - 9*(1/54*(a*b^4 + (a^2 + 1)*b^2 + a)*(b^2 + 1)^3*sqrt((27*(4*a^21 - 13*a^20 + 8*a^19 + 18*a^18 + 8*a^17 - 13*a^16 + 4*a^15)*b^36 + 3*(16*a^23 + 108*a^22 - 539*a^21 + 710*a^20 + 163*a^19 + 1676*a^18 + 163*a^17 + 710*a^16 - 539*a^15 + 108*a^14 + 16*a^13)*b^34 + 3*(64*a^24 + 44*a^23 - 995*a^22 + 3058*a^21 - 2792*a^20 + 9090*a^19 + 5094*a^18 + 9090*a^17 - 2792*a^16 + 3058*a^15 - 995*a^14 + 44*a^13 + 64*a^12)*b^32 + (288*a^25 - 660*a^24 - 1227*a^23 + 15366*a^22 - 23890*a^21 + 84450*a^20 + 5757*a^19 + 192344*a^18 + 5757*a^17 + 84450*a^16 - 23890*a^15 + 15366*a^14 - 1227*a^13 - 660*a^12 + 288*a^11)*b^30 + 6*(32*a^26 - 192*a^25 + 666*a^24 + 1039*a^23 - 1889*a^22 + 22215*a^21 - 8956*a^20 + 88402*a^19 + 17686*a^18 + 88402*a^17 - 8956*a^16 + 22215*a^15 - 1889*a^14 + 1039*a^13 + 666*a^12 - 192*a^11 + 32*a^10)*b^28 + 3*(16*a^27 - 256*a^26 + 2128*a^25 - 3376*a^24 + 13403*a^23 + 19438*a^22 + 6908*a^21 + 244734*a^20 + 25385*a^19 + 617032*a^18 + 25385*a^17 + 244734*a^16 + 6908*a^15 + 19438*a^14 + 13403*a^13 - 3376*a^12 + 2128*a^11 - 256*a^10 + 16*a^9)*b^26 - (192*a^27 - 3552*a^26 + 13056*a^25 - 72208*a^24 + 102480*a^23 - 332817*a^22 - 263254*a^21 - 417816*a^20 - 2556378*a^19 - 959054*a^18 - 2556378*a^17 - 417816*a^16 - 263254*a^15 - 332817*a^14 + 102480*a^13 - 72208*a^12 + 13056*a^11 - 3552*a^10 + 192*a^9)*b^24 + 3*(224*a^27 - 1920*a^26 + 16096*a^25 - 50504*a^24 + 182206*a^23 - 246940*a^22 + 686309*a^21 - 172838*a^20 + 1652381*a^19 + 452628*a^18 + 1652381*a^17 - 172838*a^16 + 686309*a^15 - 246940*a^14 + 182206*a^13 - 50504*a^12 + 16096*a^11 - 1920*a^10 + 224*a^9)*b^22 - 6*(160*a^27 - 2048*a^26 + 12448*a^25 - 63772*a^24 + 178110*a^23 - 523410*a^22 + 823204*a^21 - 1849465*a^20 + 1350426*a^19 - 3001882*a^18 + 1350426*a^17 - 1849465*a^16 + 823204*a^15 - 523410*a^14 + 178110*a^13 - 63772*a^12 + 12448*a^11 - 2048*a^10 + 160*a^9)*b^20 + (176*a^27 - 12672*a^26 + 104688*a^25 - 496720*a^24 + 1942122*a^23 - 4819512*a^22 + 11422375*a^21 - 15621498*a^20 + 28665279*a^19 - 21364636*a^18 + 28665279*a^17 - 15621498*a^16 + 11422375*a^15 - 4819512*a^14 + 1942122*a^13 - 496720*a^12 + 104688*a^11 - 12672*a^10 + 176*a^9)*b^18 + 432*a^18 + 3*(544*a^26 - 20416*a^25 + 142960*a^24 - 562444*a^23 + 1750931*a^22 - 3594210*a^21 + 7304088*a^20 - 7846386*a^19 + 11951018*a^18 - 7846386*a^17 + 7304088*a^16 - 3594210*a^15 + 1750931*a^14 - 562444*a^13 + 142960*a^12 - 20416*a^11 + 544*a^10)*b^16 + 3*(1184*a^25 - 42812*a^24 + 271963*a^23 - 924070*a^22 + 2450836*a^21 - 3842614*a^20 + 7208433*a^19 - 5663184*a^18 + 7208433*a^17 - 3842614*a^16 + 2450836*a^15 - 924070*a^14 + 271963*a^13 - 42812*a^12 + 1184*a^11)*b^14 - 2*(2098*a^24 + 45435*a^23 - 322221*a^22 + 967015*a^21 - 2580324*a^20 + 2430438*a^19 - 5094706*a^18 + 2430438*a^17 - 2580324*a^16 + 967015*a^15 - 322221*a^14 + 45435*a^13 + 2098*a^12)*b^12 - 3*(6585*a^23 - 15694*a^22 - 12306*a^21 - 2650*a^20 - 523527*a^19 - 138608*a^18 - 523527*a^17 - 2650*a^16 - 12306*a^15 - 15694*a^14 + 6585*a^13)*b^10 - 864*((a^9 + 3*a^8 - 8*a^6 - 6*a^5 + 6*a^4 + 8*a^3 - 3*a - 1)*b^15 + 3*(a^9 + 3*a^8 - 8*a^6 - 6*a^5 + 6*a^4 + 8*a^3 - 3*a - 1)*b^13 + 3*(a^9 + 3*a^8 - 8*a^6 - 6*a^5 + 6*a^4 + 8*a^3 - 3*a - 1)*b^11 + (a^9 + 3*a^8 - 8*a^6 - 6*a^5 + 6*a^4 + 8*a^3 - 3*a - 1)*b^9)*(a^2*b^2 + a)^(9/2)*(a*b^2 + a^2)^(9/2) - 3*(6795*a^22 - 29614*a^21 + 50456*a^20 - 225778*a^19 - 44358*a^18 - 225778*a^17 + 50456*a^16 - 29614*a^15 + 6795*a^14)*b^8 + 864*((a^12 + 3*a^11 - 8*a^9 - 6*a^8 + 6*a^7 + 8*a^6 - 3*a^4 - a^3)*b^19 + (a^13 + 10*a^12 + 10*a^11 - 17*a^10 - 26*a^9 - 8*a^8 + 8*a^7 + 26*a^6 + 17*a^5 - 10*a^4 - 10*a^3 - a^2)*b^17 + (7*a^13 + 21*a^12 + 11*a^11 - 31*a^10 - 54*a^9 - 26*a^8 + 26*a^7 + 54*a^6 + 31*a^5 - 11*a^4 - 21*a^3 - 7*a^2)*b^15 + (11*a^13 + 25*a^12 + 11*a^11 - 31*a^10 - 74*a^9 - 46*a^8 + 46*a^7 + 74*a^6 + 31*a^5 - 11*a^4 - 25*a^3 - 11*a^2)*b^13 + (5*a^13 + 22*a^12 + 18*a^11 - 33*a^10 - 62*a^9 - 28*a^8 + 28*a^7 + 62*a^6 + 33*a^5 - 18*a^4 - 22*a^3 - 5*a^2)*b^11 + 3*(3*a^12 + 5*a^11 - 4*a^10 - 12*a^9 - 6*a^8 + 6*a^7 + 12*a^6 + 4*a^5 - 5*a^4 - 3*a^3)*b^9 + 4*(a^11 + a^10 - 3*a^9 - 3*a^8 + 3*a^7 + 3*a^6 - a^5 - a^4)*b^7)*(a^2*b^2 + a)^(7/2)*(a*b^2 + a^2)^(7/2) - 4*(27*(a^14 + 4*a^12 - 40*a^11 - 32*a^10 + 134*a^9 - 32*a^8 - 40*a^7 + 4*a^6 + a^4)*b^22 + 9*(3*a^15 - 79*a^13 - 54*a^12 - 206*a^11 - 624*a^10 + 1920*a^9 - 624*a^8 - 206*a^7 - 54*a^6 - 79*a^5 + 3*a^3)*b^20 - (585*a^14 + 486*a^13 + 5491*a^12 + 4908*a^11 + 4425*a^10 - 31790*a^9 + 4425*a^8 + 4908*a^7 + 5491*a^6 + 486*a^5 + 585*a^4)*b^18 + 3*(78*a^15 - 72*a^14 - 1357*a^13 - 834*a^12 - 9634*a^11 + 3314*a^10 + 17010*a^9 + 3314*a^8 - 9634*a^7 - 834*a^6 - 1357*a^5 - 72*a^4 + 78*a^3)*b^16 + 3*(288*a^15 + 323*a^14 - 1116*a^13 - 4005*a^12 - 11692*a^11 + 1964*a^10 + 28476*a^9 + 1964*a^8 - 11692*a^7 - 4005*a^6 - 1116*a^5 + 323*a^4 + 288*a^3)*b^14 + (539*a^15 + 3480*a^14 - 3597*a^13 - 8026*a^12 - 54141*a^11 + 35070*a^10 + 53350*a^9 + 35070*a^8 - 54141*a^7 - 8026*a^6 - 3597*a^5 + 3480*a^4 + 539*a^3)*b^12 + 3*(695*a^14 + 222*a^13 - 1851*a^12 - 11164*a^11 - 562*a^10 + 25320*a^9 - 562*a^8 - 11164*a^7 - 1851*a^6 + 222*a^5 + 695*a^4)*b^10 + 3*(329*a^13 + 834*a^12 - 8212*a^11 + 1574*a^10 + 10950*a^9 + 1574*a^8 - 8212*a^7 + 834*a^6 + 329*a^5)*b^8 + (359*a^12 - 3108*a^11 - 10860*a^10 + 27218*a^9 - 10860*a^8 - 3108*a^7 + 359*a^6)*b^6 + 27*(7*a^11 - 226*a^10 + 438*a^9 - 226*a^8 + 7*a^7)*b^4 - 729*(a^10 - 2*a^9 + a^8)*b^2)*(a^2*b^2 + a)^3*(a*b^2 + a^2)^3 - (8089*a^21 - 35922*a^20 - 17817*a^19 - 261212*a^18 - 17817*a^17 - 35922*a^16 + 8089*a^15)*b^6 - 288*((a^15 + 3*a^14 - 8*a^12 - 6*a^11 + 6*a^10 + 8*a^9 - 3*a^7 - a^6)*b^23 + (2*a^16 + 17*a^15 + 11*a^14 - 34*a^13 - 28*a^12 + 2*a^11 - 2*a^10 + 28*a^9 + 34*a^8 - 11*a^7 - 17*a^6 - 2*a^5)*b^21 + (a^17 + 25*a^16 + 57*a^15 - 53*a^14 - 11*a^13 - 75*a^12 - 150*a^11 + 150*a^10 + 75*a^9 + 11*a^8 + 53*a^7 - 57*a^6 - 25*a^5 - a^4)*b^19 + (11*a^17 + 81*a^16 - 11*a^15 + 111*a^14 - 207*a^13 - 269*a^12 + 130*a^11 - 130*a^10 + 269*a^9 + 207*a^8 - 111*a^7 + 11*a^6 - 81*a^5 - 11*a^4)*b^17 + (35*a^17 + 41*a^16 + 173*a^15 - 121*a^14 - 135*a^13 - 197*a^12 - 350*a^11 + 350*a^10 + 197*a^9 + 135*a^8 + 121*a^7 - 173*a^6 - 41*a^5 - 35*a^4)*b^15 + (25*a^17 + 73*a^16 + 105*a^15 - 5*a^14 - 203*a^13 - 291*a^12 - 150*a^11 + 150*a^10 + 291*a^9 + 203*a^8 + 5*a^7 - 105*a^6 - 73*a^5 - 25*a^4)*b^13 + (90*a^16 + a^15 + 123*a^14 - 202*a^13 - 316*a^12 + 98*a^11 - 98*a^10 + 316*a^9 + 202*a^8 - 123*a^7 - a^6 - 90*a^5)*b^11 + (121*a^15 - 93*a^14 + 8*a^13 - 64*a^12 - 286*a^11 + 286*a^10 + 64*a^9 - 8*a^8 + 93*a^7 - 121*a^6)*b^9 + 72*(a^14 - a^13 - a^12 + a^11 - a^10 + a^9 + a^8 - a^7)*b^7 + 16*(a^13 - a^12 - 2*a^11 + 2*a^10 + a^9 - a^8)*b^5)*(a^2*b^2 + a)^(5/2)*(a*b^2 + a^2)^(5/2) + 24*(3*(a^17 + 13*a^15 - 49*a^14 - 68*a^13 + 206*a^12 - 68*a^11 - 49*a^10 + 13*a^9 + a^7)*b^26 + (6*a^18 + 12*a^17 - 146*a^16 + 192*a^15 - 655*a^14 - 1368*a^13 + 3918*a^12 - 1368*a^11 - 655*a^10 + 192*a^9 - 146*a^8 + 12*a^7 + 6*a^6)*b^24 + (3*a^19 + 24*a^18 - 387*a^17 + 409*a^16 - 1059*a^15 - 1981*a^14 - 1149*a^13 + 8280*a^12 - 1149*a^11 - 1981*a^10 - 1059*a^9 + 409*a^8 - 387*a^7 + 24*a^6 + 3*a^5)*b^22 + (12*a^19 - 180*a^18 - 102*a^17 - 1433*a^16 - 166*a^15 - 8951*a^14 + 5308*a^13 + 11024*a^12 + 5308*a^11 - 8951*a^10 - 166*a^9 - 1433*a^8 - 102*a^7 - 180*a^6 + 12*a^5)*b^20 + (22*a^19 + 72*a^18 - 994*a^17 - 2247*a^16 - 1574*a^15 - 20906*a^14 + 25232*a^13 + 790*a^12 + 25232*a^11 - 20906*a^10 - 1574*a^9 - 2247*a^8 - 994*a^7 + 72*a^6 + 22*a^5)*b^18 + (244*a^19 + 126*a^18 - 1372*a^17 + 39*a^16 - 21516*a^15 + 13852*a^14 - 31492*a^13 + 80238*a^12 - 31492*a^11 + 13852*a^10 - 21516*a^9 + 39*a^8 - 1372*a^7 + 126*a^6 + 244*a^5)*b^16 + (295*a^19 + 704*a^18 - 323*a^17 - 7661*a^16 + 851*a^15 - 42446*a^14 + 58937*a^13 - 20714*a^12 + 58937*a^11 - 42446*a^10 + 851*a^9 - 7661*a^8 - 323*a^7 + 704*a^6 + 295*a^5)*b^14 + (1744*a^18 - 2858*a^17 + 3957*a^16 - 27026*a^15 + 19534*a^14 - 35276*a^13 + 79850*a^12 - 35276*a^11 + 19534*a^10 - 27026*a^9 + 3957*a^8 - 2858*a^7 + 1744*a^6)*b^12 + (3205*a^17 - 9381*a^16 + 9561*a^15 - 33611*a^14 + 40156*a^13 - 19860*a^12 + 40156*a^11 - 33611*a^10 + 9561*a^9 - 9381*a^8 + 3205*a^7)*b^10 + (1999*a^16 - 7900*a^15 + 4687*a^14 - 8672*a^13 + 19772*a^12 - 8672*a^11 + 4687*a^10 - 7900*a^9 + 1999*a^8)*b^8 - (170*a^15 + 1037*a^14 + 1318*a^13 - 5050*a^12 + 1318*a^11 + 1037*a^10 + 170*a^9)*b^6 - (467*a^14 - 908*a^13 + 882*a^12 - 908*a^11 + 467*a^10)*b^4 - 54*(a^13 - 2*a^12 + a^11)*b^2)*(a^2*b^2 + a)^2*(a*b^2 + a^2)^2 - 3*(445*a^20 - 4460*a^19 - 14002*a^18 - 4460*a^17 + 445*a^16)*b^4 - 12*(9*(a^17 - 16*a^16 + 15*a^15 + 32*a^14 - 32*a^13 - 15*a^12 + 16*a^11 - a^10)*b^29 - (50*a^18 + 309*a^17 - 225*a^16 - 412*a^15 - 900*a^14 + 900*a^13 + 412*a^12 + 225*a^11 - 309*a^10 - 50*a^9)*b^27 - (159*a^19 + 270*a^18 + 118*a^17 + 495*a^16 - 2666*a^15 - 4630*a^14 + 4630*a^13 + 2666*a^12 - 495*a^11 - 118*a^10 - 270*a^9 - 159*a^8)*b^25 - (132*a^20 + 189*a^19 + 1387*a^18 - 805*a^17 + 185*a^16 - 5546*a^15 - 19266*a^14 + 19266*a^13 + 5546*a^12 - 185*a^11 + 805*a^10 - 1387*a^9 - 189*a^8 - 132*a^7)*b^23 - (32*a^21 + 84*a^20 + 2355*a^19 - 3247*a^18 + 8341*a^17 - 9058*a^16 - 25898*a^15 - 10877*a^14 + 10877*a^13 + 25898*a^12 + 9058*a^11 - 8341*a^10 + 3247*a^9 - 2355*a^8 - 84*a^7 - 32*a^6)*b^21 - (1336*a^20 - 1575*a^19 + 11311*a^18 - 18457*a^17 + 19637*a^16 - 83145*a^15 + 3575*a^14 - 3575*a^13 + 83145*a^12 - 19637*a^11 + 18457*a^10 - 11311*a^9 + 1575*a^8 - 1336*a^7)*b^19 - (160*a^21 + 440*a^20 + 4253*a^19 - 3755*a^18 + 19035*a^17 - 53273*a^16 - 26762*a^15 - 101138*a^14 + 101138*a^13 + 26762*a^12 + 53273*a^11 - 19035*a^10 + 3755*a^9 - 4253*a^8 - 440*a^7 - 160*a^6)*b^17 - (192*a^21 - 380*a^20 + 6359*a^19 - 4359*a^18 + 18139*a^17 - 45571*a^16 - 82914*a^15 - 49690*a^14 + 49690*a^13 + 82914*a^12 + 45571*a^11 - 18139*a^10 + 4359*a^9 - 6359*a^8 + 380*a^7 - 192*a^6)*b^15 - (564*a^20 - 1199*a^19 + 16555*a^18 - 16672*a^17 - 4615*a^16 - 106989*a^15 - 37604*a^14 + 37604*a^13 + 106989*a^12 + 4615*a^11 + 16672*a^10 - 16555*a^9 + 1199*a^8 - 564*a^7)*b^13 - (595*a^19 + 3707*a^18 + 9024*a^17 - 29528*a^16 - 56256*a^15 - 44068*a^14 + 44068*a^13 + 56256*a^12 + 29528*a^11 - 9024*a^10 - 3707*a^9 - 595*a^8)*b^11 - (2273*a^18 + 1795*a^17 - 3850*a^16 - 49156*a^15 - 5300*a^14 + 5300*a^13 + 49156*a^12 + 3850*a^11 - 1795*a^10 - 2273*a^9)*b^9 - 2*(2199*a^17 - 6073*a^16 + 2306*a^15 - 19090*a^14 + 19090*a^13 - 2306*a^12 + 6073*a^11 - 2199*a^10)*b^7 - (2509*a^16 - 12074*a^15 + 9169*a^14 - 9169*a^13 + 12074*a^12 - 2509*a^11)*b^5 + (19*a^15 + 2319*a^14 - 2319*a^13 - 19*a^12)*b^3 + 180*(a^14 - a^13)*b)*(a^2*b^2 + a)^(3/2)*(a*b^2 + a^2)^(3/2) - 12*((a^20 + 22*a^18 - 58*a^17 - 131*a^16 + 332*a^15 - 131*a^14 - 58*a^13 + 22*a^12 + a^10)*b^30 + (3*a^21 + 8*a^20 - 53*a^19 + 200*a^18 - 591*a^17 - 856*a^16 + 2578*a^15 - 856*a^14 - 591*a^13 + 200*a^12 - 53*a^11 + 8*a^10 + 3*a^9)*b^28 + (3*a^22 + 24*a^21 - 269*a^20 + 356*a^19 + 391*a^18 - 4380*a^17 + 1279*a^16 + 5192*a^15 + 1279*a^14 - 4380*a^13 + 391*a^12 + 356*a^11 - 269*a^10 + 24*a^9 + 3*a^8)*b^26 + (a^23 + 24*a^22 - 269*a^21 - 592*a^20 + 3135*a^19 - 9072*a^18 + 4761*a^17 - 15272*a^16 + 34568*a^15 - 15272*a^14 + 4761*a^13 - 9072*a^12 + 3135*a^11 - 592*a^10 - 269*a^9 + 24*a^8 + a^7)*b^24 + (8*a^23 - 53*a^22 - 1042*a^21 + 3539*a^20 - 15292*a^19 + 33380*a^18 - 93540*a^17 + 129394*a^16 - 112788*a^15 + 129394*a^14 - 93540*a^13 + 33380*a^12 - 15292*a^11 + 3539*a^10 - 1042*a^9 - 53*a^8 + 8*a^7)*b^22 + (22*a^23 - 232*a^22 + 1831*a^21 - 15432*a^20 + 46088*a^19 - 137144*a^18 + 216777*a^17 - 320656*a^16 + 417492*a^15 - 320656*a^14 + 216777*a^13 - 137144*a^12 + 46088*a^11 - 15432*a^10 + 1831*a^9 - 232*a^8 + 22*a^7)*b^20 + (120*a^23 + 1001*a^22 - 6836*a^21 + 24215*a^20 - 97144*a^19 + 196257*a^18 - 407140*a^17 + 580123*a^16 - 581192*a^15 + 580123*a^14 - 407140*a^13 + 196257*a^12 - 97144*a^11 + 24215*a^10 - 6836*a^9 + 1001*a^8 + 120*a^7)*b^18 + (425*a^23 - 816*a^22 + 8461*a^21 - 43976*a^20 + 107728*a^19 - 293864*a^18 + 479963*a^17 - 671616*a^16 + 827390*a^15 - 671616*a^14 + 479963*a^13 - 293864*a^12 + 107728*a^11 - 43976*a^10 + 8461*a^9 - 816*a^8 + 425*a^7)*b^16 + (3145*a^22 - 11698*a^21 + 40352*a^20 - 144364*a^19 + 280241*a^18 - 531348*a^17 + 750850*a^16 - 774356*a^15 + 750850*a^14 - 531348*a^13 + 280241*a^12 - 144364*a^11 + 40352*a^10 - 11698*a^9 + 3145*a^8)*b^14 + (8822*a^21 - 39080*a^20 + 100277*a^19 - 264296*a^18 + 441166*a^17 - 611304*a^16 + 728830*a^15 - 611304*a^14 + 441166*a^13 - 264296*a^12 + 100277*a^11 - 39080*a^10 + 8822*a^9)*b^12 + 2*(6145*a^20 - 29954*a^19 + 66630*a^18 - 135538*a^17 + 198051*a^16 - 210668*a^15 + 198051*a^14 - 135538*a^13 + 66630*a^12 - 29954*a^11 + 6145*a^10)*b^10 + (8873*a^19 - 47184*a^18 + 94988*a^17 - 148224*a^16 + 183094*a^15 - 148224*a^14 + 94988*a^13 - 47184*a^12 + 8873*a^11)*b^8 + (2849*a^18 - 17930*a^17 + 33339*a^16 - 36516*a^15 + 33339*a^14 - 17930*a^13 + 2849*a^12)*b^6 + 8*(3*a^17 - 257*a^16 + 508*a^15 - 257*a^14 + 3*a^13)*b^4 - 140*(a^16 - 2*a^15 + a^14)*b^2)*(a^2*b^2 + a)*(a*b^2 + a^2) + 24*(11*a^19 + 302*a^18 + 11*a^17)*b^2 + 12*(9*(a^20 - 8*a^19 + 8*a^18 + 17*a^17 - 17*a^16 - 8*a^15 + 8*a^14 - a^13)*b^33 - (17*a^21 + 93*a^20 + 144*a^19 - 706*a^18 - 558*a^17 + 558*a^16 + 706*a^15 - 144*a^14 - 93*a^13 - 17*a^12)*b^31 - (113*a^22 - 117*a^21 + 1215*a^20 - 2152*a^19 + 31*a^18 - 6668*a^17 + 6668*a^16 - 31*a^15 + 2152*a^14 - 1215*a^13 + 117*a^12 - 113*a^11)*b^29 - (147*a^23 - 165*a^22 + 2533*a^21 - 5497*a^20 + 9331*a^19 - 25627*a^18 + 5492*a^17 - 5492*a^16 + 25627*a^15 - 9331*a^14 + 5497*a^13 - 2533*a^12 + 165*a^11 - 147*a^10)*b^27 - (68*a^24 + 57*a^23 + 2725*a^22 - 7472*a^21 + 21905*a^20 - 49597*a^19 + 35863*a^18 - 108467*a^17 + 108467*a^16 - 35863*a^15 + 49597*a^14 - 21905*a^13 + 7472*a^12 - 2725*a^11 - 57*a^10 - 68*a^9)*b^25 - (8*a^25 + 108*a^24 + 1455*a^23 - 3697*a^22 + 20850*a^21 - 59416*a^20 + 112503*a^19 - 296579*a^18 + 260984*a^17 - 260984*a^16 + 296579*a^15 - 112503*a^14 + 59416*a^13 - 20850*a^12 + 3697*a^11 - 1455*a^10 - 108*a^9 - 8*a^8)*b^23 - (24*a^25 + 160*a^24 + 509*a^23 + 7963*a^22 - 28378*a^21 + 94606*a^20 - 266818*a^19 + 259848*a^18 - 525888*a^17 + 525888*a^16 - 259848*a^15 + 266818*a^14 - 94606*a^13 + 28378*a^12 - 7963*a^11 - 509*a^10 - 160*a^9 - 24*a^8)*b^21 + (104*a^25 - 512*a^24 + 343*a^23 - 6737*a^22 - 2718*a^21 + 23998*a^20 + 22723*a^19 + 237097*a^18 - 69190*a^17 + 69190*a^16 - 237097*a^15 - 22723*a^14 - 23998*a^13 + 2718*a^12 + 6737*a^11 - 343*a^10 + 512*a^9 - 104*a^8)*b^19 + (120*a^25 + 772*a^24 - 5659*a^23 + 15113*a^22 - 70076*a^21 + 169859*a^20 - 201439*a^19 + 571789*a^18 - 446157*a^17 + 446157*a^16 - 571789*a^15 + 201439*a^14 - 169859*a^13 + 70076*a^12 - 15113*a^11 + 5659*a^10 - 772*a^9 - 120*a^8)*b^17 + (1164*a^24 - 533*a^23 - 16549*a^22 + 47611*a^21 - 147949*a^20 + 428479*a^19 - 442011*a^18 + 785970*a^17 - 785970*a^16 + 442011*a^15 - 428479*a^14 + 147949*a^13 - 47611*a^12 + 16549*a^11 + 533*a^10 - 1164*a^9)*b^15 + (4433*a^23 - 15576*a^22 + 6113*a^21 + 19053*a^20 - 19585*a^19 + 267808*a^18 - 203604*a^17 + 203604*a^16 - 267808*a^15 + 19585*a^14 - 19053*a^13 - 6113*a^12 + 15576*a^11 - 4433*a^10)*b^13 + (8160*a^22 - 41973*a^21 + 81425*a^20 - 95468*a^19 + 222374*a^18 - 157264*a^17 + 157264*a^16 - 222374*a^15 + 95468*a^14 - 81425*a^13 + 41973*a^12 - 8160*a^11)*b^11 + (7516*a^21 - 45851*a^20 + 115790*a^19 - 139922*a^18 + 225169*a^17 - 225169*a^16 + 139922*a^15 - 115790*a^14 + 45851*a^13 - 7516*a^12)*b^9 + 2*(1525*a^20 - 10126*a^19 + 31409*a^18 - 34112*a^17 + 34112*a^16 - 31409*a^15 + 10126*a^14 - 1525*a^13)*b^7 + (243*a^19 - 805*a^18 + 9840*a^17 - 9840*a^16 + 805*a^15 - 243*a^14)*b^5 + 2*(5*a^18 + 561*a^17 - 561*a^16 - 5*a^15)*b^3 + 72*(a^17 - a^16)*b)*sqrt(a^2*b^2 + a)*sqrt(a*b^2 + a^2))*(a*b^4 + (a^2 + 1)*b^2 + a)/a)/a - 1/54*(2*(b^24 - 3*b^20 + 3*b^16 - b^12)*a^16 + 6*(b^26 - b^24 - 5*b^22 - 5*b^20 + b^18 + 13*b^16 + 9*b^14 - 7*b^12 - 6*b^10)*a^15 + 3*(2*b^28 - 6*b^26 - 17*b^24 - 30*b^22 - 38*b^20 - 82*b^18 - 124*b^16 + 54*b^14 + 196*b^12 + 64*b^10 - 19*b^8 - 2*(b^21 + 7*b^19 + 10*b^17 - 2*b^15 - 11*b^13 - 5*b^11)*sqrt(a^2*b^2 + a)*sqrt(a*b^2 + a^2))*a^14 + (2*b^30 - 18*b^28 - 33*b^26 - 68*b^24 - 129*b^22 - 282*b^20 + 792*b^18 + 2922*b^16 + 2682*b^14 + 1636*b^12 + 1761*b^10 + 930*b^8 + 45*b^6 - 6*(3*b^23 + 21*b^21 + 36*b^19 + 56*b^17 + 91*b^15 + 45*b^13 - 34*b^11 - 26*b^9)*sqrt(a^2*b^2 + a)*sqrt(a*b^2 + a^2))*a^13 - 3*(2*b^30 + b^28 - 12*b^26 - 90*b^24 - 498*b^22 - 1893*b^20 - 3232*b^18 - 3598*b^16 - 4314*b^14 - 3843*b^12 - 1780*b^10 - 768*b^8 - 406*b^6 - 49*b^4 + 4*(b^18 + 8*b^16 + 18*b^14 + 16*b^12 + 5*b^10)*(a^2*b^2 + a)*(a*b^2 + a^2) + 2*(3*b^25 + 21*b^23 + 41*b^21 + 122*b^19 + 136*b^17 - 33*b^15 - 3*b^13 + 100*b^11 - 17*b^9 - 50*b^7)*sqrt(a^2*b^2 + a)*sqrt(a*b^2 + a^2))*a^12 + 3*(b^30 + 22*b^28 + 160*b^26 + 846*b^24 + 2568*b^22 + 4632*b^20 + 7705*b^18 + 10840*b^16 + 9717*b^14 + 6834*b^12 + 4902*b^10 + 2218*b^8 + 514*b^6 + 208*b^4 - 4*(3*b^20 + 24*b^18 + 45*b^16 + 40*b^14 + 45*b^12 + 48*b^10 + 19*b^8)*(a^2*b^2 + a)*(a*b^2 + a^2) - 2*(b^27 + 7*b^25 + 7*b^23 + 17*b^21 - 211*b^19 - 378*b^17 - 98*b^15 - 227*b^13 - 384*b^11 - 15*b^9 + 45*b^7 - 44*b^5)*sqrt(a^2*b^2 + a)*sqrt(a*b^2 + a^2) + 33*b^2)*a^11 + 2*(11*b^30 + 99*b^28 + 579*b^26 + 1865*b^24 + 4701*b^22 + 10887*b^20 + 16687*b^18 + 18402*b^16 + 19359*b^14 + 15839*b^12 + 8067*b^10 + 3819*b^8 + 1745*b^6 + (9*b^19 + 27*b^17 + 26*b^15 + 6*b^13 - 3*b^11 - b^9)*(a^2*b^2 + a)^(3/2)*(a*b^2 + a^2)^(3/2) + 279*b^4 - 6*(3*b^22 + 24*b^20 + 25*b^18 + 32*b^16 + 100*b^14 + 88*b^12 + 29*b^10 + 40*b^8 + 27*b^6)*(a^2*b^2 + a)*(a*b^2 + a^2) + 3*(15*b^25 + 91*b^23 + 374*b^21 + 393*b^19 + 504*b^17 + 1164*b^15 + 732*b^13 + 90*b^11 + 357*b^9 + 165*b^7 - 62*b^5 + 17*b^3)*sqrt(a^2*b^2 + a)*sqrt(a*b^2 + a^2) + 51*b^2 + 10)*a^10 + 3*(b^30 + 22*b^28 + 160*b^26 + 846*b^24 + 2568*b^22 + 4632*b^20 + 7705*b^18 + 10840*b^16 + 9717*b^14 + 6834*b^12 + 4902*b^10 + 2218*b^8 + 514*b^6 + 2*(6*b^21 + 18*b^19 + 35*b^17 + 61*b^15 + 68*b^13 + 44*b^11 + 19*b^9 + 5*b^7)*(a^2*b^2 + a)^(3/2)*(a*b^2 + a^2)^(3/2) + 208*b^4 - 4*(b^24 + 8*b^22 - 15*b^20 + 16*b^18 + 62*b^16 - 32*b^14 - 2*b^12 + 80*b^10 + b^8 - 8*b^6 + 17*b^4)*(a^2*b^2 + a)*(a*b^2 + a^2) + 2*(7*b^27 + 37*b^25 + 92*b^23 + 22*b^21 + 451*b^19 + 738*b^17 + 16*b^15 + 364*b^13 + 723*b^11 - 69*b^9 - 48*b^7 + 122*b^5 - 25*b^3 + 2*b)*sqrt(a^2*b^2 + a)*sqrt(a*b^2 + a^2) + 33*b^2)*a^9 - 3*(2*b^30 + b^28 - 12*b^26 - 90*b^24 - 498*b^22 - 1893*b^20 - 3232*b^18 - 3598*b^16 - 4314*b^14 - 3843*b^12 - 1780*b^10 - 768*b^8 - 406*b^6 - 12*(b^16 + 4*b^14 + 6*b^12 + 4*b^10 + b^8)*(a^2*b^2 + a)^2*(a*b^2 + a^2)^2 - 2*(3*b^23 + 9*b^21 + 32*b^19 + 45*b^17 - 4*b^15 - 67*b^13 - 100*b^11 - 85*b^9 - 27*b^7 + 2*b^5)*(a^2*b^2 + a)^(3/2)*(a*b^2 + a^2)^(3/2) - 49*b^4 - 4*(15*b^22 - 16*b^20 + 38*b^18 + 136*b^16 + 14*b^14 + 56*b^12 + 154*b^10 + 8*b^8 - 9*b^6 + 24*b^4 - 4*b^2)*(a^2*b^2 + a)*(a*b^2 + a^2) + 2*(7*b^27 + 37*b^25 + 92*b^23 + 22*b^21 + 451*b^19 + 738*b^17 + 16*b^15 + 364*b^13 + 723*b^11 - 69*b^9 - 48*b^7 + 122*b^5 - 25*b^3 + 2*b)*sqrt(a^2*b^2 + a)*sqrt(a*b^2 + a^2))*a^8 + (2*b^30 - 18*b^28 - 33*b^26 - 68*b^24 - 129*b^22 - 282*b^20 + 792*b^18 + 2922*b^16 + 2682*b^14 + 1636*b^12 + 1761*b^10 + 930*b^8 + 45*b^6 + 72*(b^18 + 4*b^16 + 7*b^14 + 8*b^12 + 7*b^10 + 4*b^8 + b^6)*(a^2*b^2 + a)^2*(a*b^2 + a^2)^2 - 2*(19*b^21 + 219*b^19 + 507*b^17 + 720*b^15 + 1023*b^13 + 1140*b^11 + 981*b^9 + 732*b^7 + 318*b^5 + 37*b^3)*(a^2*b^2 + a)^(3/2)*(a*b^2 + a^2)^(3/2) + 24*(b^24 - 4*b^22 + 28*b^20 + 28*b^18 + 11*b^16 + 112*b^14 + 91*b^12 + 8*b^10 + 52*b^8 + 28*b^6 - 7*b^4 + 4*b^2)*(a^2*b^2 + a)*(a*b^2 + a^2) - 6*(15*b^25 + 91*b^23 + 374*b^21 + 393*b^19 + 504*b^17 + 1164*b^15 + 732*b^13 + 90*b^11 + 357*b^9 + 165*b^7 - 62*b^5 + 17*b^3)*sqrt(a^2*b^2 + a)*sqrt(a*b^2 + a^2))*a^7 + 3*(2*b^28 - 6*b^26 - 17*b^24 - 30*b^22 - 38*b^20 - 82*b^18 - 124*b^16 + 54*b^14 + 196*b^12 + 64*b^10 - 19*b^8 + 12*(b^20 + 4*b^18 + 8*b^16 + 12*b^14 + 14*b^12 + 12*b^10 + 8*b^8 + 4*b^6 + b^4)*(a^2*b^2 + a)^2*(a*b^2 + a^2)^2 - 2*(12*b^23 + 47*b^21 + 49*b^19 + 89*b^17 + 215*b^15 + 270*b^13 + 297*b^11 + 254*b^9 + 113*b^7 + 67*b^5 + 50*b^3 + 9*b)*(a^2*b^2 + a)^(3/2)*(a*b^2 + a^2)^(3/2) + 4*(15*b^22 - 16*b^20 + 38*b^18 + 136*b^16 + 14*b^14 + 56*b^12 + 154*b^10 + 8*b^8 - 9*b^6 + 24*b^4 - 4*b^2)*(a^2*b^2 + a)*(a*b^2 + a^2) + 2*(b^27 + 7*b^25 + 7*b^23 + 17*b^21 - 211*b^19 - 378*b^17 - 98*b^15 - 227*b^13 - 384*b^11 - 15*b^9 + 45*b^7 - 44*b^5)*sqrt(a^2*b^2 + a)*sqrt(a*b^2 + a^2))*a^6 + 6*(b^26 - b^24 - 5*b^22 - 5*b^20 + b^18 + 13*b^16 + 9*b^14 - 7*b^12 - 6*b^10 - 12*(b^18 + 4*b^16 + 7*b^14 + 8*b^12 + 7*b^10 + 4*b^8 + b^6)*(a^2*b^2 + a)^2*(a*b^2 + a^2)^2 + (12*b^23 + 47*b^21 + 49*b^19 + 89*b^17 + 215*b^15 + 270*b^13 + 297*b^11 + 254*b^9 + 113*b^7 + 67*b^5 + 50*b^3 + 9*b)*(a^2*b^2 + a)^(3/2)*(a*b^2 + a^2)^(3/2) - 2*(b^24 + 8*b^22 - 15*b^20 + 16*b^18 + 62*b^16 - 32*b^14 - 2*b^12 + 80*b^10 + b^8 - 8*b^6 + 17*b^4)*(a^2*b^2 + a)*(a*b^2 + a^2) + (3*b^25 + 21*b^23 + 41*b^21 + 122*b^19 + 136*b^17 - 33*b^15 - 3*b^13 + 100*b^11 - 17*b^9 - 50*b^7)*sqrt(a^2*b^2 + a)*sqrt(a*b^2 + a^2))*a^5 + 36*(b^16 + 4*b^14 + 6*b^12 + 4*b^10 + b^8)*(a^2*b^2 + a)^2*(a*b^2 + a^2)^2 + 2*(b^24 - 3*b^20 + 3*b^16 - b^12 - 36*(b^20 + 4*b^18 + 9*b^16 + 16*b^14 + 20*b^12 + 16*b^10 + 9*b^8 + 4*b^6 + b^4)*(a^2*b^2 + a)^2*(a*b^2 + a^2)^2 + (19*b^21 + 219*b^19 + 507*b^17 + 720*b^15 + 1023*b^13 + 1140*b^11 + 981*b^9 + 732*b^7 + 318*b^5 + 37*b^3)*(a^2*b^2 + a)^(3/2)*(a*b^2 + a^2)^(3/2) - 6*(3*b^22 + 24*b^20 + 25*b^18 + 32*b^16 + 100*b^14 + 88*b^12 + 29*b^10 + 40*b^8 + 27*b^6)*(a^2*b^2 + a)*(a*b^2 + a^2) + 3*(3*b^23 + 21*b^21 + 36*b^19 + 56*b^17 + 91*b^15 + 45*b^13 - 34*b^11 - 26*b^9)*sqrt(a^2*b^2 + a)*sqrt(a*b^2 + a^2))*a^4 - 6*(12*(b^18 + 4*b^16 + 7*b^14 + 8*b^12 + 7*b^10 + 4*b^8 + b^6)*(a^2*b^2 + a)^2*(a*b^2 + a^2)^2 + (3*b^23 + 9*b^21 + 32*b^19 + 45*b^17 - 4*b^15 - 67*b^13 - 100*b^11 - 85*b^9 - 27*b^7 + 2*b^5)*(a^2*b^2 + a)^(3/2)*(a*b^2 + a^2)^(3/2) + 2*(3*b^20 + 24*b^18 + 45*b^16 + 40*b^14 + 45*b^12 + 48*b^10 + 19*b^8)*(a^2*b^2 + a)*(a*b^2 + a^2) - (b^21 + 7*b^19 + 10*b^17 - 2*b^15 - 11*b^13 - 5*b^11)*sqrt(a^2*b^2 + a)*sqrt(a*b^2 + a^2))*a^3 + 6*(6*(b^20 + 4*b^18 + 8*b^16 + 12*b^14 + 14*b^12 + 12*b^10 + 8*b^8 + 4*b^6 + b^4)*(a^2*b^2 + a)^2*(a*b^2 + a^2)^2 - (6*b^21 + 18*b^19 + 35*b^17 + 61*b^15 + 68*b^13 + 44*b^11 + 19*b^9 + 5*b^7)*(a^2*b^2 + a)^(3/2)*(a*b^2 + a^2)^(3/2) - 2*(b^18 + 8*b^16 + 18*b^14 + 16*b^12 + 5*b^10)*(a^2*b^2 + a)*(a*b^2 + a^2))*a^2 + 2*(36*(b^18 + 4*b^16 + 7*b^14 + 8*b^12 + 7*b^10 + 4*b^8 + b^6)*(a^2*b^2 + a)^2*(a*b^2 + a^2)^2 - (9*b^19 + 27*b^17 + 26*b^15 + 6*b^13 - 3*b^11 - b^9)*(a^2*b^2 + a)^(3/2)*(a*b^2 + a^2)^(3/2))*a)/a)^(1/3) - ((b^16 - 2*b^12 + b^8)*a^11 + 2*(b^18 - b^16 - 4*b^14 - 6*b^12 - 3*b^10 + 7*b^8 + 6*b^6)*a^10 + (b^20 - 4*b^18 - 10*b^16 - 28*b^14 - 28*b^12 - 28*b^10 - 46*b^8 - 4*b^6 + 19*b^4 - 2*(b^13 + 7*b^11 + 11*b^9 + 5*b^7)*sqrt(a^2*b^2 + a)*sqrt(a*b^2 + a^2))*a^9 - 2*(b^20 + 2*b^18 + 11*b^16 + 20*b^14 + 56*b^12 + 71*b^10 + 35*b^8 + 38*b^6 + 25*b^4 + (2*b^15 + 15*b^13 + 19*b^11 + 15*b^9 + 23*b^7 + 14*b^5)*sqrt(a^2*b^2 + a)*sqrt(a*b^2 + a^2) - 3*b^2)*a^8 - 2*(4*b^18 + 12*b^16 + 44*b^14 + 56*b^12 + 76*b^10 + 109*b^8 + 52*b^6 + 14*b^4 - 2*(b^10 + 2*b^8 + b^6)*(a^2*b^2 + a)*(a*b^2 + a^2) + (b^17 + 9*b^15 + 4*b^13 + 22*b^11 + 46*b^9 + 16*b^7 + 9*b^5 + 13*b^3)*sqrt(a^2*b^2 + a)*sqrt(a*b^2 + a^2) + 16*b^2 + 1)*a^7 - 2*(b^20 + 2*b^18 + 11*b^16 + 20*b^14 + 56*b^12 + 71*b^10 + 35*b^8 + 38*b^6 + 25*b^4 - 4*(b^12 + b^10 + b^6 + b^4)*(a^2*b^2 + a)*(a*b^2 + a^2) + (b^17 - 5*b^15 + 18*b^13 + 14*b^11 - 10*b^9 + 26*b^7 + 15*b^5 - 7*b^3 + 4*b)*sqrt(a^2*b^2 + a)*sqrt(a*b^2 + a^2) - 3*b^2)*a^6 + (b^20 - 4*b^18 - 10*b^16 - 28*b^14 - 28*b^12 - 28*b^10 - 46*b^8 - 4*b^6 + 6*(b^11 + 3*b^9 + 3*b^7 + b^5)*(a^2*b^2 + a)^(3/2)*(a*b^2 + a^2)^(3/2) + 19*b^4 + 4*(b^14 - 2*b^12 - 2*b^10 + 2*b^8 - 2*b^6 - 2*b^4 + b^2)*(a^2*b^2 + a)*(a*b^2 + a^2) + 2*(b^17 - 5*b^15 + 18*b^13 + 14*b^11 - 10*b^9 + 26*b^7 + 15*b^5 - 7*b^3 + 4*b)*sqrt(a^2*b^2 + a)*sqrt(a*b^2 + a^2))*a^5 + 2*(b^18 - b^16 - 4*b^14 - 6*b^12 - 3*b^10 + 7*b^8 + 6*b^6 + 3*(b^13 + 4*b^11 + 7*b^9 + 7*b^7 + 4*b^5 + b^3)*(a^2*b^2 + a)^(3/2)*(a*b^2 + a^2)^(3/2) - 4*(b^14 + b^10 + 4*b^8 + b^6 + b^2)*(a^2*b^2 + a)*(a*b^2 + a^2) + (b^17 + 9*b^15 + 4*b^13 + 22*b^11 + 46*b^9 + 16*b^7 + 9*b^5 + 13*b^3)*sqrt(a^2*b^2 + a)*sqrt(a*b^2 + a^2))*a^4 - 6*(b^11 + 3*b^9 + 3*b^7 + b^5)*(a^2*b^2 + a)^(3/2)*(a*b^2 + a^2)^(3/2) + (b^16 - 2*b^12 + b^8 + 6*(b^13 + 3*b^11 + 4*b^9 + 4*b^7 + 3*b^5 + b^3)*(a^2*b^2 + a)^(3/2)*(a*b^2 + a^2)^(3/2) + 4*(b^14 - 2*b^12 - 2*b^10 + 2*b^8 - 2*b^6 - 2*b^4 + b^2)*(a^2*b^2 + a)*(a*b^2 + a^2) + 2*(2*b^15 + 15*b^13 + 19*b^11 + 15*b^9 + 23*b^7 + 14*b^5)*sqrt(a^2*b^2 + a)*sqrt(a*b^2 + a^2))*a^3 - 2*(3*(b^13 + 3*b^11 + 4*b^9 + 4*b^7 + 3*b^5 + b^3)*(a^2*b^2 + a)^(3/2)*(a*b^2 + a^2)^(3/2) - 4*(b^12 + b^10 + b^6 + b^4)*(a^2*b^2 + a)*(a*b^2 + a^2) - (b^13 + 7*b^11 + 11*b^9 + 5*b^7)*sqrt(a^2*b^2 + a)*sqrt(a*b^2 + a^2))*a^2 - 2*(3*(b^13 + 4*b^11 + 7*b^9 + 7*b^7 + 4*b^5 + b^3)*(a^2*b^2 + a)^(3/2)*(a*b^2 + a^2)^(3/2) - 2*(b^10 + 2*b^8 + b^6)*(a^2*b^2 + a)*(a*b^2 + a^2))*a)/((1/54*(a*b^4 + (a^2 + 1)*b^2 + a)*(b^2 + 1)^3*sqrt((27*(4*a^21 - 13*a^20 + 8*a^19 + 18*a^18 + 8*a^17 - 13*a^16 + 4*a^15)*b^36 + 3*(16*a^23 + 108*a^22 - 539*a^21 + 710*a^20 + 163*a^19 + 1676*a^18 + 163*a^17 + 710*a^16 - 539*a^15 + 108*a^14 + 16*a^13)*b^34 + 3*(64*a^24 + 44*a^23 - 995*a^22 + 3058*a^21 - 2792*a^20 + 9090*a^19 + 5094*a^18 + 9090*a^17 - 2792*a^16 + 3058*a^15 - 995*a^14 + 44*a^13 + 64*a^12)*b^32 + (288*a^25 - 660*a^24 - 1227*a^23 + 15366*a^22 - 23890*a^21 + 84450*a^20 + 5757*a^19 + 192344*a^18 + 5757*a^17 + 84450*a^16 - 23890*a^15 + 15366*a^14 - 1227*a^13 - 660*a^12 + 288*a^11)*b^30 + 6*(32*a^26 - 192*a^25 + 666*a^24 + 1039*a^23 - 1889*a^22 + 22215*a^21 - 8956*a^20 + 88402*a^19 + 17686*a^18 + 88402*a^17 - 8956*a^16 + 22215*a^15 - 1889*a^14 + 1039*a^13 + 666*a^12 - 192*a^11 + 32*a^10)*b^28 + 3*(16*a^27 - 256*a^26 + 2128*a^25 - 3376*a^24 + 13403*a^23 + 19438*a^22 + 6908*a^21 + 244734*a^20 + 25385*a^19 + 617032*a^18 + 25385*a^17 + 244734*a^16 + 6908*a^15 + 19438*a^14 + 13403*a^13 - 3376*a^12 + 2128*a^11 - 256*a^10 + 16*a^9)*b^26 - (192*a^27 - 3552*a^26 + 13056*a^25 - 72208*a^24 + 102480*a^23 - 332817*a^22 - 263254*a^21 - 417816*a^20 - 2556378*a^19 - 959054*a^18 - 2556378*a^17 - 417816*a^16 - 263254*a^15 - 332817*a^14 + 102480*a^13 - 72208*a^12 + 13056*a^11 - 3552*a^10 + 192*a^9)*b^24 + 3*(224*a^27 - 1920*a^26 + 16096*a^25 - 50504*a^24 + 182206*a^23 - 246940*a^22 + 686309*a^21 - 172838*a^20 + 1652381*a^19 + 452628*a^18 + 1652381*a^17 - 172838*a^16 + 686309*a^15 - 246940*a^14 + 182206*a^13 - 50504*a^12 + 16096*a^11 - 1920*a^10 + 224*a^9)*b^22 - 6*(160*a^27 - 2048*a^26 + 12448*a^25 - 63772*a^24 + 178110*a^23 - 523410*a^22 + 823204*a^21 - 1849465*a^20 + 1350426*a^19 - 3001882*a^18 + 1350426*a^17 - 1849465*a^16 + 823204*a^15 - 523410*a^14 + 178110*a^13 - 63772*a^12 + 12448*a^11 - 2048*a^10 + 160*a^9)*b^20 + (176*a^27 - 12672*a^26 + 104688*a^25 - 496720*a^24 + 1942122*a^23 - 4819512*a^22 + 11422375*a^21 - 15621498*a^20 + 28665279*a^19 - 21364636*a^18 + 28665279*a^17 - 15621498*a^16 + 11422375*a^15 - 4819512*a^14 + 1942122*a^13 - 496720*a^12 + 104688*a^11 - 12672*a^10 + 176*a^9)*b^18 + 432*a^18 + 3*(544*a^26 - 20416*a^25 + 142960*a^24 - 562444*a^23 + 1750931*a^22 - 3594210*a^21 + 7304088*a^20 - 7846386*a^19 + 11951018*a^18 - 7846386*a^17 + 7304088*a^16 - 3594210*a^15 + 1750931*a^14 - 562444*a^13 + 142960*a^12 - 20416*a^11 + 544*a^10)*b^16 + 3*(1184*a^25 - 42812*a^24 + 271963*a^23 - 924070*a^22 + 2450836*a^21 - 3842614*a^20 + 7208433*a^19 - 5663184*a^18 + 7208433*a^17 - 3842614*a^16 + 2450836*a^15 - 924070*a^14 + 271963*a^13 - 42812*a^12 + 1184*a^11)*b^14 - 2*(2098*a^24 + 45435*a^23 - 322221*a^22 + 967015*a^21 - 2580324*a^20 + 2430438*a^19 - 5094706*a^18 + 2430438*a^17 - 2580324*a^16 + 967015*a^15 - 322221*a^14 + 45435*a^13 + 2098*a^12)*b^12 - 3*(6585*a^23 - 15694*a^22 - 12306*a^21 - 2650*a^20 - 523527*a^19 - 138608*a^18 - 523527*a^17 - 2650*a^16 - 12306*a^15 - 15694*a^14 + 6585*a^13)*b^10 - 864*((a^9 + 3*a^8 - 8*a^6 - 6*a^5 + 6*a^4 + 8*a^3 - 3*a - 1)*b^15 + 3*(a^9 + 3*a^8 - 8*a^6 - 6*a^5 + 6*a^4 + 8*a^3 - 3*a - 1)*b^13 + 3*(a^9 + 3*a^8 - 8*a^6 - 6*a^5 + 6*a^4 + 8*a^3 - 3*a - 1)*b^11 + (a^9 + 3*a^8 - 8*a^6 - 6*a^5 + 6*a^4 + 8*a^3 - 3*a - 1)*b^9)*(a^2*b^2 + a)^(9/2)*(a*b^2 + a^2)^(9/2) - 3*(6795*a^22 - 29614*a^21 + 50456*a^20 - 225778*a^19 - 44358*a^18 - 225778*a^17 + 50456*a^16 - 29614*a^15 + 6795*a^14)*b^8 + 864*((a^12 + 3*a^11 - 8*a^9 - 6*a^8 + 6*a^7 + 8*a^6 - 3*a^4 - a^3)*b^19 + (a^13 + 10*a^12 + 10*a^11 - 17*a^10 - 26*a^9 - 8*a^8 + 8*a^7 + 26*a^6 + 17*a^5 - 10*a^4 - 10*a^3 - a^2)*b^17 + (7*a^13 + 21*a^12 + 11*a^11 - 31*a^10 - 54*a^9 - 26*a^8 + 26*a^7 + 54*a^6 + 31*a^5 - 11*a^4 - 21*a^3 - 7*a^2)*b^15 + (11*a^13 + 25*a^12 + 11*a^11 - 31*a^10 - 74*a^9 - 46*a^8 + 46*a^7 + 74*a^6 + 31*a^5 - 11*a^4 - 25*a^3 - 11*a^2)*b^13 + (5*a^13 + 22*a^12 + 18*a^11 - 33*a^10 - 62*a^9 - 28*a^8 + 28*a^7 + 62*a^6 + 33*a^5 - 18*a^4 - 22*a^3 - 5*a^2)*b^11 + 3*(3*a^12 + 5*a^11 - 4*a^10 - 12*a^9 - 6*a^8 + 6*a^7 + 12*a^6 + 4*a^5 - 5*a^4 - 3*a^3)*b^9 + 4*(a^11 + a^10 - 3*a^9 - 3*a^8 + 3*a^7 + 3*a^6 - a^5 - a^4)*b^7)*(a^2*b^2 + a)^(7/2)*(a*b^2 + a^2)^(7/2) - 4*(27*(a^14 + 4*a^12 - 40*a^11 - 32*a^10 + 134*a^9 - 32*a^8 - 40*a^7 + 4*a^6 + a^4)*b^22 + 9*(3*a^15 - 79*a^13 - 54*a^12 - 206*a^11 - 624*a^10 + 1920*a^9 - 624*a^8 - 206*a^7 - 54*a^6 - 79*a^5 + 3*a^3)*b^20 - (585*a^14 + 486*a^13 + 5491*a^12 + 4908*a^11 + 4425*a^10 - 31790*a^9 + 4425*a^8 + 4908*a^7 + 5491*a^6 + 486*a^5 + 585*a^4)*b^18 + 3*(78*a^15 - 72*a^14 - 1357*a^13 - 834*a^12 - 9634*a^11 + 3314*a^10 + 17010*a^9 + 3314*a^8 - 9634*a^7 - 834*a^6 - 1357*a^5 - 72*a^4 + 78*a^3)*b^16 + 3*(288*a^15 + 323*a^14 - 1116*a^13 - 4005*a^12 - 11692*a^11 + 1964*a^10 + 28476*a^9 + 1964*a^8 - 11692*a^7 - 4005*a^6 - 1116*a^5 + 323*a^4 + 288*a^3)*b^14 + (539*a^15 + 3480*a^14 - 3597*a^13 - 8026*a^12 - 54141*a^11 + 35070*a^10 + 53350*a^9 + 35070*a^8 - 54141*a^7 - 8026*a^6 - 3597*a^5 + 3480*a^4 + 539*a^3)*b^12 + 3*(695*a^14 + 222*a^13 - 1851*a^12 - 11164*a^11 - 562*a^10 + 25320*a^9 - 562*a^8 - 11164*a^7 - 1851*a^6 + 222*a^5 + 695*a^4)*b^10 + 3*(329*a^13 + 834*a^12 - 8212*a^11 + 1574*a^10 + 10950*a^9 + 1574*a^8 - 8212*a^7 + 834*a^6 + 329*a^5)*b^8 + (359*a^12 - 3108*a^11 - 10860*a^10 + 27218*a^9 - 10860*a^8 - 3108*a^7 + 359*a^6)*b^6 + 27*(7*a^11 - 226*a^10 + 438*a^9 - 226*a^8 + 7*a^7)*b^4 - 729*(a^10 - 2*a^9 + a^8)*b^2)*(a^2*b^2 + a)^3*(a*b^2 + a^2)^3 - (8089*a^21 - 35922*a^20 - 17817*a^19 - 261212*a^18 - 17817*a^17 - 35922*a^16 + 8089*a^15)*b^6 - 288*((a^15 + 3*a^14 - 8*a^12 - 6*a^11 + 6*a^10 + 8*a^9 - 3*a^7 - a^6)*b^23 + (2*a^16 + 17*a^15 + 11*a^14 - 34*a^13 - 28*a^12 + 2*a^11 - 2*a^10 + 28*a^9 + 34*a^8 - 11*a^7 - 17*a^6 - 2*a^5)*b^21 + (a^17 + 25*a^16 + 57*a^15 - 53*a^14 - 11*a^13 - 75*a^12 - 150*a^11 + 150*a^10 + 75*a^9 + 11*a^8 + 53*a^7 - 57*a^6 - 25*a^5 - a^4)*b^19 + (11*a^17 + 81*a^16 - 11*a^15 + 111*a^14 - 207*a^13 - 269*a^12 + 130*a^11 - 130*a^10 + 269*a^9 + 207*a^8 - 111*a^7 + 11*a^6 - 81*a^5 - 11*a^4)*b^17 + (35*a^17 + 41*a^16 + 173*a^15 - 121*a^14 - 135*a^13 - 197*a^12 - 350*a^11 + 350*a^10 + 197*a^9 + 135*a^8 + 121*a^7 - 173*a^6 - 41*a^5 - 35*a^4)*b^15 + (25*a^17 + 73*a^16 + 105*a^15 - 5*a^14 - 203*a^13 - 291*a^12 - 150*a^11 + 150*a^10 + 291*a^9 + 203*a^8 + 5*a^7 - 105*a^6 - 73*a^5 - 25*a^4)*b^13 + (90*a^16 + a^15 + 123*a^14 - 202*a^13 - 316*a^12 + 98*a^11 - 98*a^10 + 316*a^9 + 202*a^8 - 123*a^7 - a^6 - 90*a^5)*b^11 + (121*a^15 - 93*a^14 + 8*a^13 - 64*a^12 - 286*a^11 + 286*a^10 + 64*a^9 - 8*a^8 + 93*a^7 - 121*a^6)*b^9 + 72*(a^14 - a^13 - a^12 + a^11 - a^10 + a^9 + a^8 - a^7)*b^7 + 16*(a^13 - a^12 - 2*a^11 + 2*a^10 + a^9 - a^8)*b^5)*(a^2*b^2 + a)^(5/2)*(a*b^2 + a^2)^(5/2) + 24*(3*(a^17 + 13*a^15 - 49*a^14 - 68*a^13 + 206*a^12 - 68*a^11 - 49*a^10 + 13*a^9 + a^7)*b^26 + (6*a^18 + 12*a^17 - 146*a^16 + 192*a^15 - 655*a^14 - 1368*a^13 + 3918*a^12 - 1368*a^11 - 655*a^10 + 192*a^9 - 146*a^8 + 12*a^7 + 6*a^6)*b^24 + (3*a^19 + 24*a^18 - 387*a^17 + 409*a^16 - 1059*a^15 - 1981*a^14 - 1149*a^13 + 8280*a^12 - 1149*a^11 - 1981*a^10 - 1059*a^9 + 409*a^8 - 387*a^7 + 24*a^6 + 3*a^5)*b^22 + (12*a^19 - 180*a^18 - 102*a^17 - 1433*a^16 - 166*a^15 - 8951*a^14 + 5308*a^13 + 11024*a^12 + 5308*a^11 - 8951*a^10 - 166*a^9 - 1433*a^8 - 102*a^7 - 180*a^6 + 12*a^5)*b^20 + (22*a^19 + 72*a^18 - 994*a^17 - 2247*a^16 - 1574*a^15 - 20906*a^14 + 25232*a^13 + 790*a^12 + 25232*a^11 - 20906*a^10 - 1574*a^9 - 2247*a^8 - 994*a^7 + 72*a^6 + 22*a^5)*b^18 + (244*a^19 + 126*a^18 - 1372*a^17 + 39*a^16 - 21516*a^15 + 13852*a^14 - 31492*a^13 + 80238*a^12 - 31492*a^11 + 13852*a^10 - 21516*a^9 + 39*a^8 - 1372*a^7 + 126*a^6 + 244*a^5)*b^16 + (295*a^19 + 704*a^18 - 323*a^17 - 7661*a^16 + 851*a^15 - 42446*a^14 + 58937*a^13 - 20714*a^12 + 58937*a^11 - 42446*a^10 + 851*a^9 - 7661*a^8 - 323*a^7 + 704*a^6 + 295*a^5)*b^14 + (1744*a^18 - 2858*a^17 + 3957*a^16 - 27026*a^15 + 19534*a^14 - 35276*a^13 + 79850*a^12 - 35276*a^11 + 19534*a^10 - 27026*a^9 + 3957*a^8 - 2858*a^7 + 1744*a^6)*b^12 + (3205*a^17 - 9381*a^16 + 9561*a^15 - 33611*a^14 + 40156*a^13 - 19860*a^12 + 40156*a^11 - 33611*a^10 + 9561*a^9 - 9381*a^8 + 3205*a^7)*b^10 + (1999*a^16 - 7900*a^15 + 4687*a^14 - 8672*a^13 + 19772*a^12 - 8672*a^11 + 4687*a^10 - 7900*a^9 + 1999*a^8)*b^8 - (170*a^15 + 1037*a^14 + 1318*a^13 - 5050*a^12 + 1318*a^11 + 1037*a^10 + 170*a^9)*b^6 - (467*a^14 - 908*a^13 + 882*a^12 - 908*a^11 + 467*a^10)*b^4 - 54*(a^13 - 2*a^12 + a^11)*b^2)*(a^2*b^2 + a)^2*(a*b^2 + a^2)^2 - 3*(445*a^20 - 4460*a^19 - 14002*a^18 - 4460*a^17 + 445*a^16)*b^4 - 12*(9*(a^17 - 16*a^16 + 15*a^15 + 32*a^14 - 32*a^13 - 15*a^12 + 16*a^11 - a^10)*b^29 - (50*a^18 + 309*a^17 - 225*a^16 - 412*a^15 - 900*a^14 + 900*a^13 + 412*a^12 + 225*a^11 - 309*a^10 - 50*a^9)*b^27 - (159*a^19 + 270*a^18 + 118*a^17 + 495*a^16 - 2666*a^15 - 4630*a^14 + 4630*a^13 + 2666*a^12 - 495*a^11 - 118*a^10 - 270*a^9 - 159*a^8)*b^25 - (132*a^20 + 189*a^19 + 1387*a^18 - 805*a^17 + 185*a^16 - 5546*a^15 - 19266*a^14 + 19266*a^13 + 5546*a^12 - 185*a^11 + 805*a^10 - 1387*a^9 - 189*a^8 - 132*a^7)*b^23 - (32*a^21 + 84*a^20 + 2355*a^19 - 3247*a^18 + 8341*a^17 - 9058*a^16 - 25898*a^15 - 10877*a^14 + 10877*a^13 + 25898*a^12 + 9058*a^11 - 8341*a^10 + 3247*a^9 - 2355*a^8 - 84*a^7 - 32*a^6)*b^21 - (1336*a^20 - 1575*a^19 + 11311*a^18 - 18457*a^17 + 19637*a^16 - 83145*a^15 + 3575*a^14 - 3575*a^13 + 83145*a^12 - 19637*a^11 + 18457*a^10 - 11311*a^9 + 1575*a^8 - 1336*a^7)*b^19 - (160*a^21 + 440*a^20 + 4253*a^19 - 3755*a^18 + 19035*a^17 - 53273*a^16 - 26762*a^15 - 101138*a^14 + 101138*a^13 + 26762*a^12 + 53273*a^11 - 19035*a^10 + 3755*a^9 - 4253*a^8 - 440*a^7 - 160*a^6)*b^17 - (192*a^21 - 380*a^20 + 6359*a^19 - 4359*a^18 + 18139*a^17 - 45571*a^16 - 82914*a^15 - 49690*a^14 + 49690*a^13 + 82914*a^12 + 45571*a^11 - 18139*a^10 + 4359*a^9 - 6359*a^8 + 380*a^7 - 192*a^6)*b^15 - (564*a^20 - 1199*a^19 + 16555*a^18 - 16672*a^17 - 4615*a^16 - 106989*a^15 - 37604*a^14 + 37604*a^13 + 106989*a^12 + 4615*a^11 + 16672*a^10 - 16555*a^9 + 1199*a^8 - 564*a^7)*b^13 - (595*a^19 + 3707*a^18 + 9024*a^17 - 29528*a^16 - 56256*a^15 - 44068*a^14 + 44068*a^13 + 56256*a^12 + 29528*a^11 - 9024*a^10 - 3707*a^9 - 595*a^8)*b^11 - (2273*a^18 + 1795*a^17 - 3850*a^16 - 49156*a^15 - 5300*a^14 + 5300*a^13 + 49156*a^12 + 3850*a^11 - 1795*a^10 - 2273*a^9)*b^9 - 2*(2199*a^17 - 6073*a^16 + 2306*a^15 - 19090*a^14 + 19090*a^13 - 2306*a^12 + 6073*a^11 - 2199*a^10)*b^7 - (2509*a^16 - 12074*a^15 + 9169*a^14 - 9169*a^13 + 12074*a^12 - 2509*a^11)*b^5 + (19*a^15 + 2319*a^14 - 2319*a^13 - 19*a^12)*b^3 + 180*(a^14 - a^13)*b)*(a^2*b^2 + a)^(3/2)*(a*b^2 + a^2)^(3/2) - 12*((a^20 + 22*a^18 - 58*a^17 - 131*a^16 + 332*a^15 - 131*a^14 - 58*a^13 + 22*a^12 + a^10)*b^30 + (3*a^21 + 8*a^20 - 53*a^19 + 200*a^18 - 591*a^17 - 856*a^16 + 2578*a^15 - 856*a^14 - 591*a^13 + 200*a^12 - 53*a^11 + 8*a^10 + 3*a^9)*b^28 + (3*a^22 + 24*a^21 - 269*a^20 + 356*a^19 + 391*a^18 - 4380*a^17 + 1279*a^16 + 5192*a^15 + 1279*a^14 - 4380*a^13 + 391*a^12 + 356*a^11 - 269*a^10 + 24*a^9 + 3*a^8)*b^26 + (a^23 + 24*a^22 - 269*a^21 - 592*a^20 + 3135*a^19 - 9072*a^18 + 4761*a^17 - 15272*a^16 + 34568*a^15 - 15272*a^14 + 4761*a^13 - 9072*a^12 + 3135*a^11 - 592*a^10 - 269*a^9 + 24*a^8 + a^7)*b^24 + (8*a^23 - 53*a^22 - 1042*a^21 + 3539*a^20 - 15292*a^19 + 33380*a^18 - 93540*a^17 + 129394*a^16 - 112788*a^15 + 129394*a^14 - 93540*a^13 + 33380*a^12 - 15292*a^11 + 3539*a^10 - 1042*a^9 - 53*a^8 + 8*a^7)*b^22 + (22*a^23 - 232*a^22 + 1831*a^21 - 15432*a^20 + 46088*a^19 - 137144*a^18 + 216777*a^17 - 320656*a^16 + 417492*a^15 - 320656*a^14 + 216777*a^13 - 137144*a^12 + 46088*a^11 - 15432*a^10 + 1831*a^9 - 232*a^8 + 22*a^7)*b^20 + (120*a^23 + 1001*a^22 - 6836*a^21 + 24215*a^20 - 97144*a^19 + 196257*a^18 - 407140*a^17 + 580123*a^16 - 581192*a^15 + 580123*a^14 - 407140*a^13 + 196257*a^12 - 97144*a^11 + 24215*a^10 - 6836*a^9 + 1001*a^8 + 120*a^7)*b^18 + (425*a^23 - 816*a^22 + 8461*a^21 - 43976*a^20 + 107728*a^19 - 293864*a^18 + 479963*a^17 - 671616*a^16 + 827390*a^15 - 671616*a^14 + 479963*a^13 - 293864*a^12 + 107728*a^11 - 43976*a^10 + 8461*a^9 - 816*a^8 + 425*a^7)*b^16 + (3145*a^22 - 11698*a^21 + 40352*a^20 - 144364*a^19 + 280241*a^18 - 531348*a^17 + 750850*a^16 - 774356*a^15 + 750850*a^14 - 531348*a^13 + 280241*a^12 - 144364*a^11 + 40352*a^10 - 11698*a^9 + 3145*a^8)*b^14 + (8822*a^21 - 39080*a^20 + 100277*a^19 - 264296*a^18 + 441166*a^17 - 611304*a^16 + 728830*a^15 - 611304*a^14 + 441166*a^13 - 264296*a^12 + 100277*a^11 - 39080*a^10 + 8822*a^9)*b^12 + 2*(6145*a^20 - 29954*a^19 + 66630*a^18 - 135538*a^17 + 198051*a^16 - 210668*a^15 + 198051*a^14 - 135538*a^13 + 66630*a^12 - 29954*a^11 + 6145*a^10)*b^10 + (8873*a^19 - 47184*a^18 + 94988*a^17 - 148224*a^16 + 183094*a^15 - 148224*a^14 + 94988*a^13 - 47184*a^12 + 8873*a^11)*b^8 + (2849*a^18 - 17930*a^17 + 33339*a^16 - 36516*a^15 + 33339*a^14 - 17930*a^13 + 2849*a^12)*b^6 + 8*(3*a^17 - 257*a^16 + 508*a^15 - 257*a^14 + 3*a^13)*b^4 - 140*(a^16 - 2*a^15 + a^14)*b^2)*(a^2*b^2 + a)*(a*b^2 + a^2) + 24*(11*a^19 + 302*a^18 + 11*a^17)*b^2 + 12*(9*(a^20 - 8*a^19 + 8*a^18 + 17*a^17 - 17*a^16 - 8*a^15 + 8*a^14 - a^13)*b^33 - (17*a^21 + 93*a^20 + 144*a^19 - 706*a^18 - 558*a^17 + 558*a^16 + 706*a^15 - 144*a^14 - 93*a^13 - 17*a^12)*b^31 - (113*a^22 - 117*a^21 + 1215*a^20 - 2152*a^19 + 31*a^18 - 6668*a^17 + 6668*a^16 - 31*a^15 + 2152*a^14 - 1215*a^13 + 117*a^12 - 113*a^11)*b^29 - (147*a^23 - 165*a^22 + 2533*a^21 - 5497*a^20 + 9331*a^19 - 25627*a^18 + 5492*a^17 - 5492*a^16 + 25627*a^15 - 9331*a^14 + 5497*a^13 - 2533*a^12 + 165*a^11 - 147*a^10)*b^27 - (68*a^24 + 57*a^23 + 2725*a^22 - 7472*a^21 + 21905*a^20 - 49597*a^19 + 35863*a^18 - 108467*a^17 + 108467*a^16 - 35863*a^15 + 49597*a^14 - 21905*a^13 + 7472*a^12 - 2725*a^11 - 57*a^10 - 68*a^9)*b^25 - (8*a^25 + 108*a^24 + 1455*a^23 - 3697*a^22 + 20850*a^21 - 59416*a^20 + 112503*a^19 - 296579*a^18 + 260984*a^17 - 260984*a^16 + 296579*a^15 - 112503*a^14 + 59416*a^13 - 20850*a^12 + 3697*a^11 - 1455*a^10 - 108*a^9 - 8*a^8)*b^23 - (24*a^25 + 160*a^24 + 509*a^23 + 7963*a^22 - 28378*a^21 + 94606*a^20 - 266818*a^19 + 259848*a^18 - 525888*a^17 + 525888*a^16 - 259848*a^15 + 266818*a^14 - 94606*a^13 + 28378*a^12 - 7963*a^11 - 509*a^10 - 160*a^9 - 24*a^8)*b^21 + (104*a^25 - 512*a^24 + 343*a^23 - 6737*a^22 - 2718*a^21 + 23998*a^20 + 22723*a^19 + 237097*a^18 - 69190*a^17 + 69190*a^16 - 237097*a^15 - 22723*a^14 - 23998*a^13 + 2718*a^12 + 6737*a^11 - 343*a^10 + 512*a^9 - 104*a^8)*b^19 + (120*a^25 + 772*a^24 - 5659*a^23 + 15113*a^22 - 70076*a^21 + 169859*a^20 - 201439*a^19 + 571789*a^18 - 446157*a^17 + 446157*a^16 - 571789*a^15 + 201439*a^14 - 169859*a^13 + 70076*a^12 - 15113*a^11 + 5659*a^10 - 772*a^9 - 120*a^8)*b^17 + (1164*a^24 - 533*a^23 - 16549*a^22 + 47611*a^21 - 147949*a^20 + 428479*a^19 - 442011*a^18 + 785970*a^17 - 785970*a^16 + 442011*a^15 - 428479*a^14 + 147949*a^13 - 47611*a^12 + 16549*a^11 + 533*a^10 - 1164*a^9)*b^15 + (4433*a^23 - 15576*a^22 + 6113*a^21 + 19053*a^20 - 19585*a^19 + 267808*a^18 - 203604*a^17 + 203604*a^16 - 267808*a^15 + 19585*a^14 - 19053*a^13 - 6113*a^12 + 15576*a^11 - 4433*a^10)*b^13 + (8160*a^22 - 41973*a^21 + 81425*a^20 - 95468*a^19 + 222374*a^18 - 157264*a^17 + 157264*a^16 - 222374*a^15 + 95468*a^14 - 81425*a^13 + 41973*a^12 - 8160*a^11)*b^11 + (7516*a^21 - 45851*a^20 + 115790*a^19 - 139922*a^18 + 225169*a^17 - 225169*a^16 + 139922*a^15 - 115790*a^14 + 45851*a^13 - 7516*a^12)*b^9 + 2*(1525*a^20 - 10126*a^19 + 31409*a^18 - 34112*a^17 + 34112*a^16 - 31409*a^15 + 10126*a^14 - 1525*a^13)*b^7 + (243*a^19 - 805*a^18 + 9840*a^17 - 9840*a^16 + 805*a^15 - 243*a^14)*b^5 + 2*(5*a^18 + 561*a^17 - 561*a^16 - 5*a^15)*b^3 + 72*(a^17 - a^16)*b)*sqrt(a^2*b^2 + a)*sqrt(a*b^2 + a^2))*(a*b^4 + (a^2 + 1)*b^2 + a)/a)/a - 1/54*(2*(b^24 - 3*b^20 + 3*b^16 - b^12)*a^16 + 6*(b^26 - b^24 - 5*b^22 - 5*b^20 + b^18 + 13*b^16 + 9*b^14 - 7*b^12 - 6*b^10)*a^15 + 3*(2*b^28 - 6*b^26 - 17*b^24 - 30*b^22 - 38*b^20 - 82*b^18 - 124*b^16 + 54*b^14 + 196*b^12 + 64*b^10 - 19*b^8 - 2*(b^21 + 7*b^19 + 10*b^17 - 2*b^15 - 11*b^13 - 5*b^11)*sqrt(a^2*b^2 + a)*sqrt(a*b^2 + a^2))*a^14 + (2*b^30 - 18*b^28 - 33*b^26 - 68*b^24 - 129*b^22 - 282*b^20 + 792*b^18 + 2922*b^16 + 2682*b^14 + 1636*b^12 + 1761*b^10 + 930*b^8 + 45*b^6 - 6*(3*b^23 + 21*b^21 + 36*b^19 + 56*b^17 + 91*b^15 + 45*b^13 - 34*b^11 - 26*b^9)*sqrt(a^2*b^2 + a)*sqrt(a*b^2 + a^2))*a^13 - 3*(2*b^30 + b^28 - 12*b^26 - 90*b^24 - 498*b^22 - 1893*b^20 - 3232*b^18 - 3598*b^16 - 4314*b^14 - 3843*b^12 - 1780*b^10 - 768*b^8 - 406*b^6 - 49*b^4 + 4*(b^18 + 8*b^16 + 18*b^14 + 16*b^12 + 5*b^10)*(a^2*b^2 + a)*(a*b^2 + a^2) + 2*(3*b^25 + 21*b^23 + 41*b^21 + 122*b^19 + 136*b^17 - 33*b^15 - 3*b^13 + 100*b^11 - 17*b^9 - 50*b^7)*sqrt(a^2*b^2 + a)*sqrt(a*b^2 + a^2))*a^12 + 3*(b^30 + 22*b^28 + 160*b^26 + 846*b^24 + 2568*b^22 + 4632*b^20 + 7705*b^18 + 10840*b^16 + 9717*b^14 + 6834*b^12 + 4902*b^10 + 2218*b^8 + 514*b^6 + 208*b^4 - 4*(3*b^20 + 24*b^18 + 45*b^16 + 40*b^14 + 45*b^12 + 48*b^10 + 19*b^8)*(a^2*b^2 + a)*(a*b^2 + a^2) - 2*(b^27 + 7*b^25 + 7*b^23 + 17*b^21 - 211*b^19 - 378*b^17 - 98*b^15 - 227*b^13 - 384*b^11 - 15*b^9 + 45*b^7 - 44*b^5)*sqrt(a^2*b^2 + a)*sqrt(a*b^2 + a^2) + 33*b^2)*a^11 + 2*(11*b^30 + 99*b^28 + 579*b^26 + 1865*b^24 + 4701*b^22 + 10887*b^20 + 16687*b^18 + 18402*b^16 + 19359*b^14 + 15839*b^12 + 8067*b^10 + 3819*b^8 + 1745*b^6 + (9*b^19 + 27*b^17 + 26*b^15 + 6*b^13 - 3*b^11 - b^9)*(a^2*b^2 + a)^(3/2)*(a*b^2 + a^2)^(3/2) + 279*b^4 - 6*(3*b^22 + 24*b^20 + 25*b^18 + 32*b^16 + 100*b^14 + 88*b^12 + 29*b^10 + 40*b^8 + 27*b^6)*(a^2*b^2 + a)*(a*b^2 + a^2) + 3*(15*b^25 + 91*b^23 + 374*b^21 + 393*b^19 + 504*b^17 + 1164*b^15 + 732*b^13 + 90*b^11 + 357*b^9 + 165*b^7 - 62*b^5 + 17*b^3)*sqrt(a^2*b^2 + a)*sqrt(a*b^2 + a^2) + 51*b^2 + 10)*a^10 + 3*(b^30 + 22*b^28 + 160*b^26 + 846*b^24 + 2568*b^22 + 4632*b^20 + 7705*b^18 + 10840*b^16 + 9717*b^14 + 6834*b^12 + 4902*b^10 + 2218*b^8 + 514*b^6 + 2*(6*b^21 + 18*b^19 + 35*b^17 + 61*b^15 + 68*b^13 + 44*b^11 + 19*b^9 + 5*b^7)*(a^2*b^2 + a)^(3/2)*(a*b^2 + a^2)^(3/2) + 208*b^4 - 4*(b^24 + 8*b^22 - 15*b^20 + 16*b^18 + 62*b^16 - 32*b^14 - 2*b^12 + 80*b^10 + b^8 - 8*b^6 + 17*b^4)*(a^2*b^2 + a)*(a*b^2 + a^2) + 2*(7*b^27 + 37*b^25 + 92*b^23 + 22*b^21 + 451*b^19 + 738*b^17 + 16*b^15 + 364*b^13 + 723*b^11 - 69*b^9 - 48*b^7 + 122*b^5 - 25*b^3 + 2*b)*sqrt(a^2*b^2 + a)*sqrt(a*b^2 + a^2) + 33*b^2)*a^9 - 3*(2*b^30 + b^28 - 12*b^26 - 90*b^24 - 498*b^22 - 1893*b^20 - 3232*b^18 - 3598*b^16 - 4314*b^14 - 3843*b^12 - 1780*b^10 - 768*b^8 - 406*b^6 - 12*(b^16 + 4*b^14 + 6*b^12 + 4*b^10 + b^8)*(a^2*b^2 + a)^2*(a*b^2 + a^2)^2 - 2*(3*b^23 + 9*b^21 + 32*b^19 + 45*b^17 - 4*b^15 - 67*b^13 - 100*b^11 - 85*b^9 - 27*b^7 + 2*b^5)*(a^2*b^2 + a)^(3/2)*(a*b^2 + a^2)^(3/2) - 49*b^4 - 4*(15*b^22 - 16*b^20 + 38*b^18 + 136*b^16 + 14*b^14 + 56*b^12 + 154*b^10 + 8*b^8 - 9*b^6 + 24*b^4 - 4*b^2)*(a^2*b^2 + a)*(a*b^2 + a^2) + 2*(7*b^27 + 37*b^25 + 92*b^23 + 22*b^21 + 451*b^19 + 738*b^17 + 16*b^15 + 364*b^13 + 723*b^11 - 69*b^9 - 48*b^7 + 122*b^5 - 25*b^3 + 2*b)*sqrt(a^2*b^2 + a)*sqrt(a*b^2 + a^2))*a^8 + (2*b^30 - 18*b^28 - 33*b^26 - 68*b^24 - 129*b^22 - 282*b^20 + 792*b^18 + 2922*b^16 + 2682*b^14 + 1636*b^12 + 1761*b^10 + 930*b^8 + 45*b^6 + 72*(b^18 + 4*b^16 + 7*b^14 + 8*b^12 + 7*b^10 + 4*b^8 + b^6)*(a^2*b^2 + a)^2*(a*b^2 + a^2)^2 - 2*(19*b^21 + 219*b^19 + 507*b^17 + 720*b^15 + 1023*b^13 + 1140*b^11 + 981*b^9 + 732*b^7 + 318*b^5 + 37*b^3)*(a^2*b^2 + a)^(3/2)*(a*b^2 + a^2)^(3/2) + 24*(b^24 - 4*b^22 + 28*b^20 + 28*b^18 + 11*b^16 + 112*b^14 + 91*b^12 + 8*b^10 + 52*b^8 + 28*b^6 - 7*b^4 + 4*b^2)*(a^2*b^2 + a)*(a*b^2 + a^2) - 6*(15*b^25 + 91*b^23 + 374*b^21 + 393*b^19 + 504*b^17 + 1164*b^15 + 732*b^13 + 90*b^11 + 357*b^9 + 165*b^7 - 62*b^5 + 17*b^3)*sqrt(a^2*b^2 + a)*sqrt(a*b^2 + a^2))*a^7 + 3*(2*b^28 - 6*b^26 - 17*b^24 - 30*b^22 - 38*b^20 - 82*b^18 - 124*b^16 + 54*b^14 + 196*b^12 + 64*b^10 - 19*b^8 + 12*(b^20 + 4*b^18 + 8*b^16 + 12*b^14 + 14*b^12 + 12*b^10 + 8*b^8 + 4*b^6 + b^4)*(a^2*b^2 + a)^2*(a*b^2 + a^2)^2 - 2*(12*b^23 + 47*b^21 + 49*b^19 + 89*b^17 + 215*b^15 + 270*b^13 + 297*b^11 + 254*b^9 + 113*b^7 + 67*b^5 + 50*b^3 + 9*b)*(a^2*b^2 + a)^(3/2)*(a*b^2 + a^2)^(3/2) + 4*(15*b^22 - 16*b^20 + 38*b^18 + 136*b^16 + 14*b^14 + 56*b^12 + 154*b^10 + 8*b^8 - 9*b^6 + 24*b^4 - 4*b^2)*(a^2*b^2 + a)*(a*b^2 + a^2) + 2*(b^27 + 7*b^25 + 7*b^23 + 17*b^21 - 211*b^19 - 378*b^17 - 98*b^15 - 227*b^13 - 384*b^11 - 15*b^9 + 45*b^7 - 44*b^5)*sqrt(a^2*b^2 + a)*sqrt(a*b^2 + a^2))*a^6 + 6*(b^26 - b^24 - 5*b^22 - 5*b^20 + b^18 + 13*b^16 + 9*b^14 - 7*b^12 - 6*b^10 - 12*(b^18 + 4*b^16 + 7*b^14 + 8*b^12 + 7*b^10 + 4*b^8 + b^6)*(a^2*b^2 + a)^2*(a*b^2 + a^2)^2 + (12*b^23 + 47*b^21 + 49*b^19 + 89*b^17 + 215*b^15 + 270*b^13 + 297*b^11 + 254*b^9 + 113*b^7 + 67*b^5 + 50*b^3 + 9*b)*(a^2*b^2 + a)^(3/2)*(a*b^2 + a^2)^(3/2) - 2*(b^24 + 8*b^22 - 15*b^20 + 16*b^18 + 62*b^16 - 32*b^14 - 2*b^12 + 80*b^10 + b^8 - 8*b^6 + 17*b^4)*(a^2*b^2 + a)*(a*b^2 + a^2) + (3*b^25 + 21*b^23 + 41*b^21 + 122*b^19 + 136*b^17 - 33*b^15 - 3*b^13 + 100*b^11 - 17*b^9 - 50*b^7)*sqrt(a^2*b^2 + a)*sqrt(a*b^2 + a^2))*a^5 + 36*(b^16 + 4*b^14 + 6*b^12 + 4*b^10 + b^8)*(a^2*b^2 + a)^2*(a*b^2 + a^2)^2 + 2*(b^24 - 3*b^20 + 3*b^16 - b^12 - 36*(b^20 + 4*b^18 + 9*b^16 + 16*b^14 + 20*b^12 + 16*b^10 + 9*b^8 + 4*b^6 + b^4)*(a^2*b^2 + a)^2*(a*b^2 + a^2)^2 + (19*b^21 + 219*b^19 + 507*b^17 + 720*b^15 + 1023*b^13 + 1140*b^11 + 981*b^9 + 732*b^7 + 318*b^5 + 37*b^3)*(a^2*b^2 + a)^(3/2)*(a*b^2 + a^2)^(3/2) - 6*(3*b^22 + 24*b^20 + 25*b^18 + 32*b^16 + 100*b^14 + 88*b^12 + 29*b^10 + 40*b^8 + 27*b^6)*(a^2*b^2 + a)*(a*b^2 + a^2) + 3*(3*b^23 + 21*b^21 + 36*b^19 + 56*b^17 + 91*b^15 + 45*b^13 - 34*b^11 - 26*b^9)*sqrt(a^2*b^2 + a)*sqrt(a*b^2 + a^2))*a^4 - 6*(12*(b^18 + 4*b^16 + 7*b^14 + 8*b^12 + 7*b^10 + 4*b^8 + b^6)*(a^2*b^2 + a)^2*(a*b^2 + a^2)^2 + (3*b^23 + 9*b^21 + 32*b^19 + 45*b^17 - 4*b^15 - 67*b^13 - 100*b^11 - 85*b^9 - 27*b^7 + 2*b^5)*(a^2*b^2 + a)^(3/2)*(a*b^2 + a^2)^(3/2) + 2*(3*b^20 + 24*b^18 + 45*b^16 + 40*b^14 + 45*b^12 + 48*b^10 + 19*b^8)*(a^2*b^2 + a)*(a*b^2 + a^2) - (b^21 + 7*b^19 + 10*b^17 - 2*b^15 - 11*b^13 - 5*b^11)*sqrt(a^2*b^2 + a)*sqrt(a*b^2 + a^2))*a^3 + 6*(6*(b^20 + 4*b^18 + 8*b^16 + 12*b^14 + 14*b^12 + 12*b^10 + 8*b^8 + 4*b^6 + b^4)*(a^2*b^2 + a)^2*(a*b^2 + a^2)^2 - (6*b^21 + 18*b^19 + 35*b^17 + 61*b^15 + 68*b^13 + 44*b^11 + 19*b^9 + 5*b^7)*(a^2*b^2 + a)^(3/2)*(a*b^2 + a^2)^(3/2) - 2*(b^18 + 8*b^16 + 18*b^14 + 16*b^12 + 5*b^10)*(a^2*b^2 + a)*(a*b^2 + a^2))*a^2 + 2*(36*(b^18 + 4*b^16 + 7*b^14 + 8*b^12 + 7*b^10 + 4*b^8 + b^6)*(a^2*b^2 + a)^2*(a*b^2 + a^2)^2 - (9*b^19 + 27*b^17 + 26*b^15 + 6*b^13 - 3*b^11 - b^9)*(a^2*b^2 + a)^(3/2)*(a*b^2 + a^2)^(3/2))*a)/a)^(1/3)*a)) - (sqrt(a^2*b^2 + a)*sqrt(a*b^2 + a^2)*(b^2 + 1)^(3/2)*(sqrt(a*b^2 + a^2)*(a*b^2 + 1) - sqrt(a^2*b^2 + a)*(b^3 + a*b))*a*b + sqrt(a^2*b^2 + a)*sqrt(a*b^2 + a^2)*(b^2 + 1)^(3/2)*((a*b^3 + b)*sqrt(a*b^2 + a^2) + sqrt(a^2*b^2 + a)*(b^2 + a)))^2
a1 = SR(AA(14778/100000))
b1 = SR(AA(77656/100000))
print(expr.subs(a=a1, b=b1).minpoly())

MinimalReproducing8.sage

ZZy.<y> = ZZ[]
poly8 = ZZy([
106346275030920586004276738159284223160745160329414802674148696415560316409061498114049459644001698876690153812161622530710214150820561031037785280344393299646600156897949037810194870162400984124666960310214525780110648622089311363935212895781758830692752289761380743270340278626938299598583183832574931221725553565254820532939785346746231356273732577381690752781714567846773023346958287913797090269420198245156242586972934757809596130740923433261746417755524427662333604332151073009621028952867341874474226358516188234334856136546786036008823240823139642247244408516688519232364150088668864406302978147634902827970593978579746892046832066760554149496750495855890983444017853685989839294337798592903735922065448332446935785156074894307761739174238107368878131577200234300680138251796669301289858254996149281097308015925742230852069954244156352963488435590260218531054759395011235869401688302297074533082100095584405764604354096420974538390817121633490421461930709741274376939358773837412991489285454175589219720348007951400,
21595873767273098831715407663577014047612178342013989156104290844099052007843836432847667901176337161202628964215075785914621920629985934921246108812445296629577460278824852712526743556769516530964518214713381656159049808334626331659342801931340107179710323312442163063006740180149555035674205081356652537559014121902836707522368424560039571561544244732464548333075681518403554586371269235312690674519890666506396738820149739366949045634722190449823858160939916477024476640771831041610040530162824963348690560359804017459306399755756429348029569901239141761718059694387787411731781948882247690566271179742719263580025591907418753793505946662664215190136399059484912639161088677691409848019268342569051164385444698750604146646171769502655397869039265456466896264473600,
-21595873767273098831715407663577014047612178342013989156104290844099052007843836432847667901176337161202628964215075785914621920629985934921246108812445296629577460278824852712526743556769516530964518214713381656159049808334626331659342801931340107179710342840932815724762756703206291989559590448970192436210378026601090903218189588812227847012610665455817691722989036849605377289513027623525970334355916763657692967671480904644496444738599861795838493733558867349216688224238699025150137865708135568810118296899122376751706750168047103364146819182582846619860288402376055403038463763425968747394750116427335314635224968191203490832518861801718565885229220211055830316504539067868598834912210812262795276799425521272732942022020478348640156203774846863079479220679820,
39056981305323512033046113473907770770735227079797302727809396508391391642328504376550902132841446706286779826710662403645406283516776426559319672052194302592457702662330555094798207755342692029271145237901744384423166933735967080194671090621210922855473076542559650519180931506115786787197756014442850575694550472977190959357760004435809540545505603518541981053027385045550452838531837634132151782604367437572151403210526968642601837610563821264785194721882166941284459512337610476500803013342670374109948897460,
-19528490652661756016523056736953885385367613539898651363904698254195695821164252188275451066420723353143389913355331201822703141758388213279659836026097151296228851331165277547399103877671346014635572618950872192211583466867983540097335545310605461427736533035882489555790733064977279775346861324804261015940111477591616344664207294702144479241833387475370667162200736212966163451019817630642979883707344458128067492379630019982871048119342995073793514261533890179876499293468601505710735362120978059415065661279,
-6282476802844559679225696736341902420018900597126288596510676374962017607249018912349237103297140680388037175547571770875561895321423707715209113807112789609851070760157206115844823126698670318899719288631948918876555240244479047599373460428553167890544944,
+2094158934281519893075232245447300806672966865708762865503558791654005869083006304116412367765713560129345725182523923625187298440474569238403037935704263203283690253385735371948274375566223439633239762877316306292185080081493015866457820142851055963514986,
-4,
1])
a = poly8.roots(QQbar, False)[0]
a.minpoly()