Skip to content

Instantly share code, notes, and snippets.

@gagern
Created June 27, 2013 15:27
Show Gist options
  • Save gagern/5877421 to your computer and use it in GitHub Desktop.
Save gagern/5877421 to your computer and use it in GitHub Desktop.
Formula to compute the locus of parabolas through three given points (ax,ay), (bx,by) and (cx,cy). See http://math.stackexchange.com/a/429858/35416
TO = TermOrder('wdeglex',(100, 10, 10, 1, 1, 1, 1, 1, 1))
PR1.<t, xx, yy, ax, ay, bx, by, cx, cy> = PolynomialRing(QQ, order=TO)
var('mu tau')
A = vector(PR1, [ax, ay, 1])
B = vector(PR1, [bx, by, 1])
C = vector(PR1, [cx, cy, 1])
D = vector(PR1, [1, t, 0])
E1 = A.cross_product(B).column() * C.cross_product(D).row()
E2 = A.cross_product(C).column() * B.cross_product(D).row()
F1 = E1 + E1.transpose()
F2 = E2 + E2.transpose()
mus = solve((F1 + mu*F2).submatrix(0, 0, 2, 2).det(), mu)
assert(len(mus) == 1)
mu1 = mus[0].rhs()
assert((PR1(mu1.denominator())*mu1 - PR1(mu1.numerator())).is_zero())
P = PR1(mu1.denominator())*F1 + PR1(mu1.numerator())*F2
H = vector(PR1, [D[1], -D[0], 0])
L = P*H
M = Matrix(PR1, [
[0, -L[2], L[1]],
[L[2], 0, -L[0]],
[-L[1], L[0], 0]])
N = M.transpose()*P*M
taus = solve((N + tau*M).submatrix(0, 0, 2, 2).det(), tau)
assert((taus[0].rhs() + taus[1].rhs()).is_zero())
assert((PR1(taus[0].rhs()) - taus[0].rhs()).is_zero())
Q = N + PR1(taus[0].rhs())*M
assert(Q[2,2].is_zero())
assert(Q[1,2].is_zero() or Q[2,1].is_zero())
R = Q.row(1) if Q[2,1].is_zero() else Q.column(1)
def degs(v, sym=t):
return [i.degree(sym) for i in v]
assert(degs(R) == [8, 8, 8])
S = vector(PR1, [Ri/gcd(R) for Ri in R])
assert(degs(S) == [7, 7, 7])
# [len(list(Si)) for Si in S]
PR2 = PR1.change_ring(order='deglex')
PR3 = PR2.univariate_ring(t)
S2 = S.change_ring(PR2).change_ring(PR3)
for name, coord in zip(["x", "y", "z"], S2):
for deg, terms in enumerate(coord):
print((name + ' =' if deg == 0 else ') +')
+ " sin(alpha)^{}*cos(alpha)^{}*(".format(deg, coord.degree() - deg))
for coeff, mono in terms:
s = str(coeff*mono)
if s.startswith('-'):
s = ' - ' + s[1:]
else:
s = ' + ' + s
print(s)
print(");")
print("")
x = sin(alpha)^0*cos(alpha)^7*(
- 2*ax^2*by^4
+ 8*ax^2*by^3*cy
- 12*ax^2*by^2*cy^2
+ 8*ax^2*by*cy^3
- 2*ax^2*cy^4
+ 4*ax*ay^2*bx*by^2
- 8*ax*ay^2*bx*by*cy
+ 4*ax*ay^2*bx*cy^2
+ 4*ax*ay^2*by^2*cx
- 8*ax*ay^2*by*cx*cy
+ 4*ax*ay^2*cx*cy^2
- 8*ax*ay*bx*by^2*cy
+ 16*ax*ay*bx*by*cy^2
- 8*ax*ay*bx*cy^3
- 8*ax*ay*by^3*cx
+ 16*ax*ay*by^2*cx*cy
- 8*ax*ay*by*cx*cy^2
+ 4*ax*bx*by^2*cy^2
- 8*ax*bx*by*cy^3
+ 4*ax*bx*cy^4
+ 4*ax*by^4*cx
- 8*ax*by^3*cx*cy
+ 4*ax*by^2*cx*cy^2
- 2*ay^4*bx^2
+ 4*ay^4*bx*cx
- 2*ay^4*cx^2
+ 8*ay^3*bx^2*cy
- 8*ay^3*bx*by*cx
- 8*ay^3*bx*cx*cy
+ 8*ay^3*by*cx^2
- 12*ay^2*bx^2*cy^2
+ 4*ay^2*bx*by^2*cx
+ 16*ay^2*bx*by*cx*cy
+ 4*ay^2*bx*cx*cy^2
- 12*ay^2*by^2*cx^2
+ 8*ay*bx^2*cy^3
- 8*ay*bx*by^2*cx*cy
- 8*ay*bx*by*cx*cy^2
+ 8*ay*by^3*cx^2
- 2*bx^2*cy^4
+ 4*bx*by^2*cx*cy^2
- 2*by^4*cx^2
) + sin(alpha)^1*cos(alpha)^6*(
- 8*ax^2*ay*bx*by^2
+ 16*ax^2*ay*bx*by*cy
- 8*ax^2*ay*bx*cy^2
- 8*ax^2*ay*by^2*cx
+ 16*ax^2*ay*by*cx*cy
- 8*ax^2*ay*cx*cy^2
+ 8*ax^2*bx*by^3
- 16*ax^2*bx*by^2*cy
+ 8*ax^2*bx*by*cy^2
+ 8*ax^2*by^2*cx*cy
- 16*ax^2*by*cx*cy^2
+ 8*ax^2*cx*cy^3
+ 8*ax*ay^3*bx^2
- 16*ax*ay^3*bx*cx
+ 8*ax*ay^3*cx^2
- 8*ax*ay^2*bx^2*by
- 16*ax*ay^2*bx^2*cy
+ 24*ax*ay^2*bx*by*cx
+ 24*ax*ay^2*bx*cx*cy
- 16*ax*ay^2*by*cx^2
- 8*ax*ay^2*cx^2*cy
+ 16*ax*ay*bx^2*by*cy
+ 8*ax*ay*bx^2*cy^2
+ 24*ax*ay*bx*by^2*cx
- 96*ax*ay*bx*by*cx*cy
+ 24*ax*ay*bx*cx*cy^2
+ 8*ax*ay*by^2*cx^2
+ 16*ax*ay*by*cx^2*cy
- 8*ax*bx^2*by*cy^2
- 16*ax*bx*by^3*cx
+ 24*ax*bx*by^2*cx*cy
+ 24*ax*bx*by*cx*cy^2
- 16*ax*bx*cx*cy^3
- 8*ax*by^2*cx^2*cy
- 8*ay^2*bx^2*by*cx
+ 8*ay^2*bx^2*cx*cy
+ 8*ay^2*bx*by*cx^2
- 8*ay^2*bx*cx^2*cy
+ 16*ay*bx^2*by*cx*cy
- 16*ay*bx^2*cx*cy^2
- 16*ay*bx*by^2*cx^2
+ 16*ay*bx*by*cx^2*cy
- 8*bx^2*by*cx*cy^2
+ 8*bx^2*cx*cy^3
+ 8*bx*by^3*cx^2
- 8*bx*by^2*cx^2*cy
) + sin(alpha)^2*cos(alpha)^5*(
+ 4*ax^3*bx*by^2
- 8*ax^3*bx*by*cy
+ 4*ax^3*bx*cy^2
+ 4*ax^3*by^2*cx
- 8*ax^3*by*cx*cy
+ 4*ax^3*cx*cy^2
- 12*ax^2*ay^2*bx^2
+ 24*ax^2*ay^2*bx*cx
- 12*ax^2*ay^2*cx^2
+ 16*ax^2*ay*bx^2*by
+ 8*ax^2*ay*bx^2*cy
- 24*ax^2*ay*bx*by*cx
- 24*ax^2*ay*bx*cx*cy
+ 8*ax^2*ay*by*cx^2
+ 16*ax^2*ay*cx^2*cy
- 12*ax^2*bx^2*by^2
+ 8*ax^2*bx^2*by*cy
- 8*ax^2*bx^2*cy^2
- 4*ax^2*bx*by^2*cx
+ 32*ax^2*bx*by*cx*cy
- 4*ax^2*bx*cx*cy^2
- 4*ax^2*by^4
+ 16*ax^2*by^3*cy
- 8*ax^2*by^2*cx^2
- 24*ax^2*by^2*cy^2
+ 8*ax^2*by*cx^2*cy
+ 16*ax^2*by*cy^3
- 12*ax^2*cx^2*cy^2
- 4*ax^2*cy^4
+ 4*ax*ay^2*bx^3
- 4*ax*ay^2*bx^2*cx
+ 8*ax*ay^2*bx*by^2
- 16*ax*ay^2*bx*by*cy
- 4*ax*ay^2*bx*cx^2
+ 8*ax*ay^2*bx*cy^2
+ 8*ax*ay^2*by^2*cx
- 16*ax*ay^2*by*cx*cy
+ 4*ax*ay^2*cx^3
+ 8*ax*ay^2*cx*cy^2
- 8*ax*ay*bx^3*cy
- 24*ax*ay*bx^2*by*cx
+ 32*ax*ay*bx^2*cx*cy
- 16*ax*ay*bx*by^2*cy
+ 32*ax*ay*bx*by*cx^2
+ 32*ax*ay*bx*by*cy^2
- 24*ax*ay*bx*cx^2*cy
- 16*ax*ay*bx*cy^3
- 16*ax*ay*by^3*cx
+ 32*ax*ay*by^2*cx*cy
- 8*ax*ay*by*cx^3
- 16*ax*ay*by*cx*cy^2
+ 4*ax*bx^3*cy^2
+ 24*ax*bx^2*by^2*cx
- 24*ax*bx^2*by*cx*cy
- 4*ax*bx^2*cx*cy^2
- 4*ax*bx*by^2*cx^2
+ 8*ax*bx*by^2*cy^2
- 24*ax*bx*by*cx^2*cy
- 16*ax*bx*by*cy^3
+ 24*ax*bx*cx^2*cy^2
+ 8*ax*bx*cy^4
+ 8*ax*by^4*cx
- 16*ax*by^3*cx*cy
+ 4*ax*by^2*cx^3
+ 8*ax*by^2*cx*cy^2
- 4*ay^4*bx^2
+ 8*ay^4*bx*cx
+ 2*ay^4*by^2
- 4*ay^4*by*cy
- 4*ay^4*cx^2
+ 2*ay^4*cy^2
+ 16*ay^3*bx^2*cy
- 16*ay^3*bx*by*cx
- 16*ay^3*bx*cx*cy
- 4*ay^3*by^3
+ 4*ay^3*by^2*cy
+ 16*ay^3*by*cx^2
+ 4*ay^3*by*cy^2
- 4*ay^3*cy^3
+ 4*ay^2*bx^3*cx
- 8*ay^2*bx^2*cx^2
- 24*ay^2*bx^2*cy^2
+ 8*ay^2*bx*by^2*cx
+ 32*ay^2*bx*by*cx*cy
+ 4*ay^2*bx*cx^3
+ 8*ay^2*bx*cx*cy^2
+ 2*ay^2*by^4
+ 4*ay^2*by^3*cy
- 24*ay^2*by^2*cx^2
- 12*ay^2*by^2*cy^2
+ 4*ay^2*by*cy^3
+ 2*ay^2*cy^4
- 8*ay*bx^3*cx*cy
+ 8*ay*bx^2*by*cx^2
+ 8*ay*bx^2*cx^2*cy
+ 16*ay*bx^2*cy^3
- 16*ay*bx*by^2*cx*cy
- 8*ay*bx*by*cx^3
- 16*ay*bx*by*cx*cy^2
- 4*ay*by^4*cy
+ 16*ay*by^3*cx^2
+ 4*ay*by^3*cy^2
+ 4*ay*by^2*cy^3
- 4*ay*by*cy^4
+ 4*bx^3*cx*cy^2
- 12*bx^2*by^2*cx^2
+ 16*bx^2*by*cx^2*cy
- 12*bx^2*cx^2*cy^2
- 4*bx^2*cy^4
+ 4*bx*by^2*cx^3
+ 8*bx*by^2*cx*cy^2
- 4*by^4*cx^2
+ 2*by^4*cy^2
- 4*by^3*cy^3
+ 2*by^2*cy^4
) + sin(alpha)^3*cos(alpha)^4*(
+ 8*ax^3*ay*bx^2
- 16*ax^3*ay*bx*cx
+ 8*ax^3*ay*cx^2
- 8*ax^3*bx^2*by
+ 8*ax^3*bx*by*cx
+ 8*ax^3*bx*cx*cy
- 8*ax^3*cx^2*cy
- 8*ax^2*ay*bx^3
+ 8*ax^2*ay*bx^2*cx
- 16*ax^2*ay*bx*by^2
+ 32*ax^2*ay*bx*by*cy
+ 8*ax^2*ay*bx*cx^2
- 16*ax^2*ay*bx*cy^2
- 16*ax^2*ay*by^2*cx
+ 32*ax^2*ay*by*cx*cy
- 8*ax^2*ay*cx^3
- 16*ax^2*ay*cx*cy^2
+ 8*ax^2*bx^3*by
+ 8*ax^2*bx^2*by*cx
- 16*ax^2*bx^2*cx*cy
+ 16*ax^2*bx*by^3
- 32*ax^2*bx*by^2*cy
- 16*ax^2*bx*by*cx^2
+ 16*ax^2*bx*by*cy^2
+ 8*ax^2*bx*cx^2*cy
+ 16*ax^2*by^2*cx*cy
- 32*ax^2*by*cx*cy^2
+ 8*ax^2*cx^3*cy
+ 16*ax^2*cx*cy^3
+ 16*ax*ay^3*bx^2
- 32*ax*ay^3*bx*cx
- 8*ax*ay^3*by^2
+ 16*ax*ay^3*by*cy
+ 16*ax*ay^3*cx^2
- 8*ax*ay^3*cy^2
- 16*ax*ay^2*bx^2*by
- 32*ax*ay^2*bx^2*cy
+ 48*ax*ay^2*bx*by*cx
+ 48*ax*ay^2*bx*cx*cy
+ 12*ax*ay^2*by^3
- 12*ax*ay^2*by^2*cy
- 32*ax*ay^2*by*cx^2
- 12*ax*ay^2*by*cy^2
- 16*ax*ay^2*cx^2*cy
+ 12*ax*ay^2*cy^3
+ 8*ax*ay*bx^3*cx
+ 32*ax*ay*bx^2*by*cy
- 16*ax*ay*bx^2*cx^2
+ 16*ax*ay*bx^2*cy^2
+ 48*ax*ay*bx*by^2*cx
- 192*ax*ay*bx*by*cx*cy
+ 8*ax*ay*bx*cx^3
+ 48*ax*ay*bx*cx*cy^2
- 4*ax*ay*by^4
- 8*ax*ay*by^3*cy
+ 16*ax*ay*by^2*cx^2
+ 24*ax*ay*by^2*cy^2
+ 32*ax*ay*by*cx^2*cy
- 8*ax*ay*by*cy^3
- 4*ax*ay*cy^4
- 16*ax*bx^3*by*cx
+ 8*ax*bx^3*cx*cy
+ 8*ax*bx^2*by*cx^2
- 16*ax*bx^2*by*cy^2
+ 8*ax*bx^2*cx^2*cy
- 32*ax*bx*by^3*cx
+ 48*ax*bx*by^2*cx*cy
+ 8*ax*bx*by*cx^3
+ 48*ax*bx*by*cx*cy^2
- 16*ax*bx*cx^3*cy
- 32*ax*bx*cx*cy^3
+ 4*ax*by^4*cy
- 4*ax*by^3*cy^2
- 16*ax*by^2*cx^2*cy
- 4*ax*by^2*cy^3
+ 4*ax*by*cy^4
- 4*ay^4*bx*by
+ 4*ay^4*bx*cy
+ 4*ay^4*by*cx
- 4*ay^4*cx*cy
+ 12*ay^3*bx*by^2
- 8*ay^3*bx*by*cy
- 4*ay^3*bx*cy^2
- 4*ay^3*by^2*cx
- 8*ay^3*by*cx*cy
+ 12*ay^3*cx*cy^2
- 16*ay^2*bx^2*by*cx
+ 16*ay^2*bx^2*cx*cy
- 8*ay^2*bx*by^3
- 12*ay^2*bx*by^2*cy
+ 16*ay^2*bx*by*cx^2
+ 24*ay^2*bx*by*cy^2
- 16*ay^2*bx*cx^2*cy
- 4*ay^2*bx*cy^3
- 4*ay^2*by^3*cx
+ 24*ay^2*by^2*cx*cy
- 12*ay^2*by*cx*cy^2
- 8*ay^2*cx*cy^3
+ 32*ay*bx^2*by*cx*cy
- 32*ay*bx^2*cx*cy^2
+ 16*ay*bx*by^3*cy
- 32*ay*bx*by^2*cx^2
- 12*ay*bx*by^2*cy^2
+ 32*ay*bx*by*cx^2*cy
- 8*ay*bx*by*cy^3
+ 4*ay*bx*cy^4
+ 4*ay*by^4*cx
- 8*ay*by^3*cx*cy
- 12*ay*by^2*cx*cy^2
+ 16*ay*by*cx*cy^3
+ 8*bx^3*by*cx^2
- 8*bx^3*cx^2*cy
- 8*bx^2*by*cx^3
- 16*bx^2*by*cx*cy^2
+ 8*bx^2*cx^3*cy
+ 16*bx^2*cx*cy^3
+ 16*bx*by^3*cx^2
- 8*bx*by^3*cy^2
- 16*bx*by^2*cx^2*cy
+ 12*bx*by^2*cy^3
- 4*bx*by*cy^4
- 4*by^4*cx*cy
+ 12*by^3*cx*cy^2
- 8*by^2*cx*cy^3
) + sin(alpha)^4*cos(alpha)^3*(
- 2*ax^4*bx^2
+ 4*ax^4*bx*cx
- 2*ax^4*cx^2
+ 4*ax^3*bx^3
- 4*ax^3*bx^2*cx
+ 8*ax^3*bx*by^2
- 16*ax^3*bx*by*cy
- 4*ax^3*bx*cx^2
+ 8*ax^3*bx*cy^2
+ 8*ax^3*by^2*cx
- 16*ax^3*by*cx*cy
+ 4*ax^3*cx^3
+ 8*ax^3*cx*cy^2
- 24*ax^2*ay^2*bx^2
+ 48*ax^2*ay^2*bx*cx
+ 12*ax^2*ay^2*by^2
- 24*ax^2*ay^2*by*cy
- 24*ax^2*ay^2*cx^2
+ 12*ax^2*ay^2*cy^2
+ 32*ax^2*ay*bx^2*by
+ 16*ax^2*ay*bx^2*cy
- 48*ax^2*ay*bx*by*cx
- 48*ax^2*ay*bx*cx*cy
- 12*ax^2*ay*by^3
+ 12*ax^2*ay*by^2*cy
+ 16*ax^2*ay*by*cx^2
+ 12*ax^2*ay*by*cy^2
+ 32*ax^2*ay*cx^2*cy
- 12*ax^2*ay*cy^3
- 2*ax^2*bx^4
- 4*ax^2*bx^3*cx
- 24*ax^2*bx^2*by^2
+ 16*ax^2*bx^2*by*cy
+ 12*ax^2*bx^2*cx^2
- 16*ax^2*bx^2*cy^2
- 8*ax^2*bx*by^2*cx
+ 64*ax^2*bx*by*cx*cy
- 4*ax^2*bx*cx^3
- 8*ax^2*bx*cx*cy^2
+ 12*ax^2*by^3*cy
- 16*ax^2*by^2*cx^2
- 24*ax^2*by^2*cy^2
+ 16*ax^2*by*cx^2*cy
+ 12*ax^2*by*cy^3
- 2*ax^2*cx^4
- 24*ax^2*cx^2*cy^2
+ 16*ax*ay^3*bx*by
- 16*ax*ay^3*bx*cy
- 16*ax*ay^3*by*cx
+ 16*ax*ay^3*cx*cy
+ 8*ax*ay^2*bx^3
- 8*ax*ay^2*bx^2*cx
- 32*ax*ay^2*bx*by^2
+ 16*ax*ay^2*bx*by*cy
- 8*ax*ay^2*bx*cx^2
+ 16*ax*ay^2*bx*cy^2
+ 16*ax*ay^2*by^2*cx
+ 16*ax*ay^2*by*cx*cy
+ 8*ax*ay^2*cx^3
- 32*ax*ay^2*cx*cy^2
- 16*ax*ay*bx^3*cy
- 48*ax*ay*bx^2*by*cx
+ 64*ax*ay*bx^2*cx*cy
+ 16*ax*ay*bx*by^3
+ 16*ax*ay*bx*by^2*cy
+ 64*ax*ay*bx*by*cx^2
- 32*ax*ay*bx*by*cy^2
- 48*ax*ay*bx*cx^2*cy
- 32*ax*ay*by^2*cx*cy
- 16*ax*ay*by*cx^3
+ 16*ax*ay*by*cx*cy^2
+ 16*ax*ay*cx*cy^3
+ 4*ax*bx^4*cx
- 4*ax*bx^3*cx^2
+ 8*ax*bx^3*cy^2
+ 48*ax*bx^2*by^2*cx
- 48*ax*bx^2*by*cx*cy
- 4*ax*bx^2*cx^3
- 8*ax*bx^2*cx*cy^2
- 16*ax*bx*by^3*cy
- 8*ax*bx*by^2*cx^2
+ 16*ax*bx*by^2*cy^2
- 48*ax*bx*by*cx^2*cy
+ 4*ax*bx*cx^4
+ 48*ax*bx*cx^2*cy^2
+ 8*ax*by^2*cx^3
+ 16*ax*by^2*cx*cy^2
- 16*ax*by*cx*cy^3
- 12*ay^3*bx^2*by
+ 12*ay^3*bx^2*cy
+ 12*ay^3*by*cx^2
- 12*ay^3*cx^2*cy
+ 8*ay^2*bx^3*cx
+ 12*ay^2*bx^2*by^2
+ 12*ay^2*bx^2*by*cy
- 16*ay^2*bx^2*cx^2
- 24*ay^2*bx^2*cy^2
+ 16*ay^2*bx*by^2*cx
- 32*ay^2*bx*by*cx*cy
+ 8*ay^2*bx*cx^3
+ 16*ay^2*bx*cx*cy^2
- 24*ay^2*by^2*cx^2
+ 12*ay^2*by*cx^2*cy
+ 12*ay^2*cx^2*cy^2
- 16*ay*bx^3*cx*cy
- 24*ay*bx^2*by^2*cy
+ 16*ay*bx^2*by*cx^2
+ 12*ay*bx^2*by*cy^2
+ 16*ay*bx^2*cx^2*cy
+ 12*ay*bx^2*cy^3
- 16*ay*bx*by^3*cx
+ 16*ay*bx*by^2*cx*cy
- 16*ay*bx*by*cx^3
+ 16*ay*bx*by*cx*cy^2
- 16*ay*bx*cx*cy^3
+ 12*ay*by^3*cx^2
+ 12*ay*by^2*cx^2*cy
- 24*ay*by*cx^2*cy^2
- 2*bx^4*cx^2
+ 4*bx^3*cx^3
+ 8*bx^3*cx*cy^2
- 24*bx^2*by^2*cx^2
+ 12*bx^2*by^2*cy^2
+ 32*bx^2*by*cx^2*cy
- 12*bx^2*by*cy^3
- 2*bx^2*cx^4
- 24*bx^2*cx^2*cy^2
+ 16*bx*by^3*cx*cy
+ 8*bx*by^2*cx^3
- 32*bx*by^2*cx*cy^2
+ 16*bx*by*cx*cy^3
- 12*by^3*cx^2*cy
+ 12*by^2*cx^2*cy^2
) + sin(alpha)^5*cos(alpha)^2*(
+ 16*ax^3*ay*bx^2
- 32*ax^3*ay*bx*cx
- 8*ax^3*ay*by^2
+ 16*ax^3*ay*by*cy
+ 16*ax^3*ay*cx^2
- 8*ax^3*ay*cy^2
- 16*ax^3*bx^2*by
+ 16*ax^3*bx*by*cx
+ 16*ax^3*bx*cx*cy
+ 4*ax^3*by^3
- 4*ax^3*by^2*cy
- 4*ax^3*by*cy^2
- 16*ax^3*cx^2*cy
+ 4*ax^3*cy^3
- 24*ax^2*ay^2*bx*by
+ 24*ax^2*ay^2*bx*cy
+ 24*ax^2*ay^2*by*cx
- 24*ax^2*ay^2*cx*cy
- 16*ax^2*ay*bx^3
+ 16*ax^2*ay*bx^2*cx
+ 28*ax^2*ay*bx*by^2
- 8*ax^2*ay*bx*by*cy
+ 16*ax^2*ay*bx*cx^2
- 20*ax^2*ay*bx*cy^2
- 20*ax^2*ay*by^2*cx
- 8*ax^2*ay*by*cx*cy
- 16*ax^2*ay*cx^3
+ 28*ax^2*ay*cx*cy^2
+ 16*ax^2*bx^3*by
+ 16*ax^2*bx^2*by*cx
- 32*ax^2*bx^2*cx*cy
- 28*ax^2*bx*by^2*cy
- 32*ax^2*bx*by*cx^2
+ 32*ax^2*bx*by*cy^2
+ 16*ax^2*bx*cx^2*cy
- 4*ax^2*bx*cy^3
- 4*ax^2*by^3*cx
+ 32*ax^2*by^2*cx*cy
- 28*ax^2*by*cx*cy^2
+ 16*ax^2*cx^3*cy
+ 28*ax*ay^2*bx^2*by
- 28*ax*ay^2*bx^2*cy
- 28*ax*ay^2*by*cx^2
+ 28*ax*ay^2*cx^2*cy
+ 16*ax*ay*bx^3*cx
- 24*ax*ay*bx^2*by^2
- 8*ax*ay*bx^2*by*cy
- 32*ax*ay*bx^2*cx^2
+ 32*ax*ay*bx^2*cy^2
+ 16*ax*ay*bx*cx^3
+ 32*ax*ay*by^2*cx^2
- 8*ax*ay*by*cx^2*cy
- 24*ax*ay*cx^2*cy^2
- 32*ax*bx^3*by*cx
+ 16*ax*bx^3*cx*cy
+ 24*ax*bx^2*by^2*cy
+ 16*ax*bx^2*by*cx^2
- 20*ax*bx^2*by*cy^2
+ 16*ax*bx^2*cx^2*cy
- 4*ax*bx^2*cy^3
+ 16*ax*bx*by*cx^3
- 32*ax*bx*cx^3*cy
- 4*ax*by^3*cx^2
- 20*ax*by^2*cx^2*cy
+ 24*ax*by*cx^2*cy^2
+ 4*ay^3*bx^3
- 4*ay^3*bx^2*cx
- 4*ay^3*bx*cx^2
+ 4*ay^3*cx^3
- 8*ay^2*bx^3*by
- 4*ay^2*bx^3*cy
- 20*ay^2*bx^2*by*cx
+ 32*ay^2*bx^2*cx*cy
+ 32*ay^2*bx*by*cx^2
- 20*ay^2*bx*cx^2*cy
- 4*ay^2*by*cx^3
- 8*ay^2*cx^3*cy
+ 16*ay*bx^3*by*cy
- 4*ay*bx^3*cy^2
+ 24*ay*bx^2*by^2*cx
- 8*ay*bx^2*by*cx*cy
- 28*ay*bx^2*cx*cy^2
- 28*ay*bx*by^2*cx^2
- 8*ay*bx*by*cx^2*cy
+ 24*ay*bx*cx^2*cy^2
- 4*ay*by^2*cx^3
+ 16*ay*by*cx^3*cy
+ 16*bx^3*by*cx^2
- 8*bx^3*by*cy^2
- 16*bx^3*cx^2*cy
+ 4*bx^3*cy^3
- 24*bx^2*by^2*cx*cy
- 16*bx^2*by*cx^3
+ 28*bx^2*by*cx*cy^2
+ 16*bx^2*cx^3*cy
+ 28*bx*by^2*cx^2*cy
- 24*bx*by*cx^2*cy^2
+ 4*by^3*cx^3
- 8*by^2*cx^3*cy
) + sin(alpha)^6*cos(alpha)^1*(
- 4*ax^4*bx^2
+ 8*ax^4*bx*cx
+ 2*ax^4*by^2
- 4*ax^4*by*cy
- 4*ax^4*cx^2
+ 2*ax^4*cy^2
+ 16*ax^3*ay*bx*by
- 16*ax^3*ay*bx*cy
- 16*ax^3*ay*by*cx
+ 16*ax^3*ay*cx*cy
+ 8*ax^3*bx^3
- 8*ax^3*bx^2*cx
- 8*ax^3*bx*by^2
- 8*ax^3*bx*cx^2
+ 8*ax^3*bx*cy^2
+ 8*ax^3*by^2*cx
+ 8*ax^3*cx^3
- 8*ax^3*cx*cy^2
- 20*ax^2*ay*bx^2*by
+ 20*ax^2*ay*bx^2*cy
+ 20*ax^2*ay*by*cx^2
- 20*ax^2*ay*cx^2*cy
- 4*ax^2*bx^4
- 8*ax^2*bx^3*cx
+ 20*ax^2*bx^2*by*cy
+ 24*ax^2*bx^2*cx^2
- 20*ax^2*bx^2*cy^2
+ 8*ax^2*bx*by^2*cx
- 16*ax^2*bx*by*cx*cy
- 8*ax^2*bx*cx^3
+ 8*ax^2*bx*cx*cy^2
- 20*ax^2*by^2*cx^2
+ 20*ax^2*by*cx^2*cy
- 4*ax^2*cx^4
- 8*ax*ay^2*bx^3
+ 8*ax*ay^2*bx^2*cx
+ 8*ax*ay^2*bx*cx^2
- 8*ax*ay^2*cx^3
+ 16*ax*ay*bx^3*by
- 16*ax*ay*bx^2*cx*cy
- 16*ax*ay*bx*by*cx^2
+ 16*ax*ay*cx^3*cy
+ 8*ax*bx^4*cx
- 16*ax*bx^3*by*cy
- 8*ax*bx^3*cx^2
+ 8*ax*bx^3*cy^2
- 8*ax*bx^2*cx^3
+ 8*ax*bx^2*cx*cy^2
+ 8*ax*bx*by^2*cx^2
+ 8*ax*bx*cx^4
+ 8*ax*by^2*cx^3
- 16*ax*by*cx^3*cy
+ 2*ay^2*bx^4
+ 8*ay^2*bx^3*cx
- 20*ay^2*bx^2*cx^2
+ 8*ay^2*bx*cx^3
+ 2*ay^2*cx^4
- 4*ay*bx^4*cy
- 16*ay*bx^3*by*cx
+ 20*ay*bx^2*by*cx^2
+ 20*ay*bx^2*cx^2*cy
- 16*ay*bx*cx^3*cy
- 4*ay*by*cx^4
- 4*bx^4*cx^2
+ 2*bx^4*cy^2
+ 16*bx^3*by*cx*cy
+ 8*bx^3*cx^3
- 8*bx^3*cx*cy^2
- 20*bx^2*by*cx^2*cy
- 4*bx^2*cx^4
- 8*bx*by^2*cx^3
+ 16*bx*by*cx^3*cy
+ 2*by^2*cx^4
) + sin(alpha)^7*cos(alpha)^0*(
- 4*ax^4*bx*by
+ 4*ax^4*bx*cy
+ 4*ax^4*by*cx
- 4*ax^4*cx*cy
+ 4*ax^3*bx^2*by
- 4*ax^3*bx^2*cy
- 4*ax^3*by*cx^2
+ 4*ax^3*cx^2*cy
+ 4*ax^2*ay*bx^3
- 4*ax^2*ay*bx^2*cx
- 4*ax^2*ay*bx*cx^2
+ 4*ax^2*ay*cx^3
- 4*ax^2*bx^3*cy
- 4*ax^2*bx^2*by*cx
+ 8*ax^2*bx^2*cx*cy
+ 8*ax^2*bx*by*cx^2
- 4*ax^2*bx*cx^2*cy
- 4*ax^2*by*cx^3
- 4*ax*ay*bx^4
+ 8*ax*ay*bx^2*cx^2
- 4*ax*ay*cx^4
+ 4*ax*bx^4*cy
- 4*ax*bx^2*by*cx^2
- 4*ax*bx^2*cx^2*cy
+ 4*ax*by*cx^4
+ 4*ay*bx^4*cx
- 4*ay*bx^3*cx^2
- 4*ay*bx^2*cx^3
+ 4*ay*bx*cx^4
- 4*bx^4*cx*cy
+ 4*bx^3*cx^2*cy
+ 4*bx^2*by*cx^3
- 4*bx*by*cx^4
);
y = sin(alpha)^0*cos(alpha)^7*(
+ 4*ax*ay^2*by^3
- 4*ax*ay^2*by^2*cy
- 4*ax*ay^2*by*cy^2
+ 4*ax*ay^2*cy^3
- 4*ax*ay*by^4
+ 8*ax*ay*by^2*cy^2
- 4*ax*ay*cy^4
+ 4*ax*by^4*cy
- 4*ax*by^3*cy^2
- 4*ax*by^2*cy^3
+ 4*ax*by*cy^4
- 4*ay^4*bx*by
+ 4*ay^4*bx*cy
+ 4*ay^4*by*cx
- 4*ay^4*cx*cy
+ 4*ay^3*bx*by^2
- 4*ay^3*bx*cy^2
- 4*ay^3*by^2*cx
+ 4*ay^3*cx*cy^2
- 4*ay^2*bx*by^2*cy
+ 8*ay^2*bx*by*cy^2
- 4*ay^2*bx*cy^3
- 4*ay^2*by^3*cx
+ 8*ay^2*by^2*cx*cy
- 4*ay^2*by*cx*cy^2
- 4*ay*bx*by^2*cy^2
+ 4*ay*bx*cy^4
+ 4*ay*by^4*cx
- 4*ay*by^2*cx*cy^2
+ 4*bx*by^2*cy^3
- 4*bx*by*cy^4
- 4*by^4*cx*cy
+ 4*by^3*cx*cy^2
) + sin(alpha)^1*cos(alpha)^6*(
- 8*ax^2*ay*by^3
+ 8*ax^2*ay*by^2*cy
+ 8*ax^2*ay*by*cy^2
- 8*ax^2*ay*cy^3
+ 2*ax^2*by^4
+ 8*ax^2*by^3*cy
- 20*ax^2*by^2*cy^2
+ 8*ax^2*by*cy^3
+ 2*ax^2*cy^4
+ 16*ax*ay^3*bx*by
- 16*ax*ay^3*bx*cy
- 16*ax*ay^3*by*cx
+ 16*ax*ay^3*cx*cy
- 20*ax*ay^2*bx*by^2
+ 20*ax*ay^2*bx*cy^2
+ 20*ax*ay^2*by^2*cx
- 20*ax*ay^2*cx*cy^2
+ 16*ax*ay*bx*by^3
- 16*ax*ay*bx*by*cy^2
- 16*ax*ay*by^2*cx*cy
+ 16*ax*ay*cx*cy^3
- 16*ax*bx*by^3*cy
+ 20*ax*bx*by^2*cy^2
- 4*ax*bx*cy^4
- 4*ax*by^4*cx
+ 20*ax*by^2*cx*cy^2
- 16*ax*by*cx*cy^3
+ 2*ay^4*bx^2
- 4*ay^4*bx*cx
- 4*ay^4*by^2
+ 8*ay^4*by*cy
+ 2*ay^4*cx^2
- 4*ay^4*cy^2
- 8*ay^3*bx^2*by
+ 8*ay^3*bx^2*cy
+ 8*ay^3*by^3
- 8*ay^3*by^2*cy
+ 8*ay^3*by*cx^2
- 8*ay^3*by*cy^2
- 8*ay^3*cx^2*cy
+ 8*ay^3*cy^3
+ 8*ay^2*bx^2*by*cy
- 20*ay^2*bx^2*cy^2
+ 20*ay^2*bx*by^2*cx
- 16*ay^2*bx*by*cx*cy
+ 20*ay^2*bx*cx*cy^2
- 4*ay^2*by^4
- 8*ay^2*by^3*cy
- 20*ay^2*by^2*cx^2
+ 24*ay^2*by^2*cy^2
+ 8*ay^2*by*cx^2*cy
- 8*ay^2*by*cy^3
- 4*ay^2*cy^4
+ 8*ay*bx^2*by*cy^2
+ 8*ay*bx^2*cy^3
- 16*ay*bx*by^3*cx
- 16*ay*bx*cx*cy^3
+ 8*ay*by^4*cy
+ 8*ay*by^3*cx^2
- 8*ay*by^3*cy^2
+ 8*ay*by^2*cx^2*cy
- 8*ay*by^2*cy^3
+ 8*ay*by*cy^4
- 8*bx^2*by*cy^3
+ 2*bx^2*cy^4
+ 16*bx*by^3*cx*cy
- 20*bx*by^2*cx*cy^2
+ 16*bx*by*cx*cy^3
+ 2*by^4*cx^2
- 4*by^4*cy^2
- 8*by^3*cx^2*cy
+ 8*by^3*cy^3
- 4*by^2*cy^4
) + sin(alpha)^2*cos(alpha)^5*(
+ 4*ax^3*by^3
- 4*ax^3*by^2*cy
- 4*ax^3*by*cy^2
+ 4*ax^3*cy^3
- 24*ax^2*ay^2*bx*by
+ 24*ax^2*ay^2*bx*cy
+ 24*ax^2*ay^2*by*cx
- 24*ax^2*ay^2*cx*cy
+ 28*ax^2*ay*bx*by^2
- 28*ax^2*ay*bx*cy^2
- 28*ax^2*ay*by^2*cx
+ 28*ax^2*ay*cx*cy^2
- 8*ax^2*bx*by^3
- 20*ax^2*bx*by^2*cy
+ 32*ax^2*bx*by*cy^2
- 4*ax^2*bx*cy^3
- 4*ax^2*by^3*cx
+ 32*ax^2*by^2*cx*cy
- 20*ax^2*by*cx*cy^2
- 8*ax^2*cx*cy^3
- 8*ax*ay^3*bx^2
+ 16*ax*ay^3*bx*cx
+ 16*ax*ay^3*by^2
- 32*ax*ay^3*by*cy
- 8*ax*ay^3*cx^2
+ 16*ax*ay^3*cy^2
+ 28*ax*ay^2*bx^2*by
- 20*ax*ay^2*bx^2*cy
- 8*ax*ay^2*bx*by*cx
- 8*ax*ay^2*bx*cx*cy
- 16*ax*ay^2*by^3
+ 16*ax*ay^2*by^2*cy
- 20*ax*ay^2*by*cx^2
+ 16*ax*ay^2*by*cy^2
+ 28*ax*ay^2*cx^2*cy
- 16*ax*ay^2*cy^3
- 24*ax*ay*bx^2*by^2
+ 32*ax*ay*bx^2*cy^2
- 8*ax*ay*bx*by^2*cx
- 8*ax*ay*bx*cx*cy^2
+ 16*ax*ay*by^3*cy
+ 32*ax*ay*by^2*cx^2
- 32*ax*ay*by^2*cy^2
+ 16*ax*ay*by*cy^3
- 24*ax*ay*cx^2*cy^2
+ 24*ax*bx^2*by^2*cy
- 28*ax*bx^2*by*cy^2
- 4*ax*bx^2*cy^3
+ 16*ax*bx*by^3*cx
- 8*ax*bx*by^2*cx*cy
- 8*ax*bx*by*cx*cy^2
+ 16*ax*bx*cx*cy^3
- 4*ax*by^3*cx^2
- 28*ax*by^2*cx^2*cy
+ 24*ax*by*cx^2*cy^2
+ 4*ay^3*bx^3
- 4*ay^3*bx^2*cx
- 16*ay^3*bx*by^2
+ 16*ay^3*bx*by*cy
- 4*ay^3*bx*cx^2
+ 16*ay^3*by*cx*cy
+ 4*ay^3*cx^3
- 16*ay^3*cx*cy^2
- 4*ay^2*bx^3*cy
- 28*ay^2*bx^2*by*cx
+ 32*ay^2*bx^2*cx*cy
+ 16*ay^2*bx*by^3
+ 16*ay^2*bx*by^2*cy
+ 32*ay^2*bx*by*cx^2
- 32*ay^2*bx*by*cy^2
- 28*ay^2*bx*cx^2*cy
- 32*ay^2*by^2*cx*cy
- 4*ay^2*by*cx^3
+ 16*ay^2*by*cx*cy^2
+ 16*ay^2*cx*cy^3
- 4*ay*bx^3*cy^2
+ 24*ay*bx^2*by^2*cx
- 20*ay*bx^2*cx*cy^2
- 32*ay*bx*by^3*cy
- 20*ay*bx*by^2*cx^2
+ 16*ay*bx*by^2*cy^2
+ 16*ay*bx*by*cy^3
+ 24*ay*bx*cx^2*cy^2
+ 16*ay*by^3*cx*cy
- 4*ay*by^2*cx^3
+ 16*ay*by^2*cx*cy^2
- 32*ay*by*cx*cy^3
+ 4*bx^3*cy^3
- 24*bx^2*by^2*cx*cy
+ 28*bx^2*by*cx*cy^2
- 8*bx^2*cx*cy^3
- 8*bx*by^3*cx^2
+ 16*bx*by^3*cy^2
+ 28*bx*by^2*cx^2*cy
- 16*bx*by^2*cy^3
- 24*bx*by*cx^2*cy^2
+ 4*by^3*cx^3
- 16*by^3*cx*cy^2
+ 16*by^2*cx*cy^3
) + sin(alpha)^3*cos(alpha)^4*(
+ 16*ax^3*ay*bx*by
- 16*ax^3*ay*bx*cy
- 16*ax^3*ay*by*cx
+ 16*ax^3*ay*cx*cy
- 12*ax^3*bx*by^2
+ 12*ax^3*bx*cy^2
+ 12*ax^3*by^2*cx
- 12*ax^3*cx*cy^2
+ 12*ax^2*ay^2*bx^2
- 24*ax^2*ay^2*bx*cx
- 24*ax^2*ay^2*by^2
+ 48*ax^2*ay^2*by*cy
+ 12*ax^2*ay^2*cx^2
- 24*ax^2*ay^2*cy^2
- 32*ax^2*ay*bx^2*by
+ 16*ax^2*ay*bx^2*cy
+ 16*ax^2*ay*bx*by*cx
+ 16*ax^2*ay*bx*cx*cy
+ 8*ax^2*ay*by^3
- 8*ax^2*ay*by^2*cy
+ 16*ax^2*ay*by*cx^2
- 8*ax^2*ay*by*cy^2
- 32*ax^2*ay*cx^2*cy
+ 8*ax^2*ay*cy^3
+ 12*ax^2*bx^2*by^2
+ 16*ax^2*bx^2*by*cy
- 24*ax^2*bx^2*cy^2
+ 12*ax^2*bx*by^2*cx
- 32*ax^2*bx*by*cx*cy
+ 12*ax^2*bx*cx*cy^2
+ 8*ax^2*by^3*cy
- 24*ax^2*by^2*cx^2
- 16*ax^2*by^2*cy^2
+ 16*ax^2*by*cx^2*cy
+ 8*ax^2*by*cy^3
+ 12*ax^2*cx^2*cy^2
- 12*ax*ay^2*bx^3
+ 12*ax*ay^2*bx^2*cx
+ 32*ax*ay^2*bx*by^2
- 48*ax*ay^2*bx*by*cy
+ 12*ax*ay^2*bx*cx^2
+ 16*ax*ay^2*bx*cy^2
+ 16*ax*ay^2*by^2*cx
- 48*ax*ay^2*by*cx*cy
- 12*ax*ay^2*cx^3
+ 32*ax*ay^2*cx*cy^2
+ 16*ax*ay*bx^3*by
+ 16*ax*ay*bx^2*by*cx
- 32*ax*ay*bx^2*cx*cy
- 48*ax*ay*bx*by^2*cy
- 32*ax*ay*bx*by*cx^2
+ 64*ax*ay*bx*by*cy^2
+ 16*ax*ay*bx*cx^2*cy
- 16*ax*ay*bx*cy^3
- 16*ax*ay*by^3*cx
+ 64*ax*ay*by^2*cx*cy
- 48*ax*ay*by*cx*cy^2
+ 16*ax*ay*cx^3*cy
- 16*ax*bx^3*by*cy
+ 12*ax*bx^3*cy^2
- 24*ax*bx^2*by^2*cx
+ 16*ax*bx^2*by*cx*cy
+ 12*ax*bx^2*cx*cy^2
+ 12*ax*bx*by^2*cx^2
+ 16*ax*bx*by^2*cy^2
+ 16*ax*bx*by*cx^2*cy
- 16*ax*bx*by*cy^3
- 24*ax*bx*cx^2*cy^2
- 16*ax*by^3*cx*cy
+ 12*ax*by^2*cx^3
+ 16*ax*by^2*cx*cy^2
- 16*ax*by*cx^3*cy
- 2*ay^4*by^2
+ 4*ay^4*by*cy
- 2*ay^4*cy^2
+ 8*ay^3*bx^2*by
+ 8*ay^3*bx^2*cy
- 16*ay^3*bx*by*cx
- 16*ay^3*bx*cx*cy
+ 4*ay^3*by^3
- 4*ay^3*by^2*cy
+ 8*ay^3*by*cx^2
- 4*ay^3*by*cy^2
+ 8*ay^3*cx^2*cy
+ 4*ay^3*cy^3
+ 12*ay^2*bx^3*cx
- 24*ay^2*bx^2*by^2
- 8*ay^2*bx^2*by*cy
- 24*ay^2*bx^2*cx^2
- 16*ay^2*bx^2*cy^2
+ 16*ay^2*bx*by^2*cx
+ 64*ay^2*bx*by*cx*cy
+ 12*ay^2*bx*cx^3
+ 16*ay^2*bx*cx*cy^2
- 2*ay^2*by^4
- 4*ay^2*by^3*cy
- 16*ay^2*by^2*cx^2
+ 12*ay^2*by^2*cy^2
- 8*ay^2*by*cx^2*cy
- 4*ay^2*by*cy^3
- 24*ay^2*cx^2*cy^2
- 2*ay^2*cy^4
- 16*ay*bx^3*by*cx
+ 48*ay*bx^2*by^2*cy
+ 16*ay*bx^2*by*cx^2
- 8*ay*bx^2*by*cy^2
+ 16*ay*bx^2*cx^2*cy
+ 8*ay*bx^2*cy^3
- 48*ay*bx*by^2*cx*cy
- 48*ay*bx*by*cx*cy^2
- 16*ay*bx*cx^3*cy
+ 4*ay*by^4*cy
+ 8*ay*by^3*cx^2
- 4*ay*by^3*cy^2
- 8*ay*by^2*cx^2*cy
- 4*ay*by^2*cy^3
+ 48*ay*by*cx^2*cy^2
+ 4*ay*by*cy^4
+ 16*bx^3*by*cx*cy
- 12*bx^3*cx*cy^2
+ 12*bx^2*by^2*cx^2
- 24*bx^2*by^2*cy^2
- 32*bx^2*by*cx^2*cy
+ 8*bx^2*by*cy^3
+ 12*bx^2*cx^2*cy^2
- 12*bx*by^2*cx^3
+ 32*bx*by^2*cx*cy^2
+ 16*bx*by*cx^3*cy
- 2*by^4*cy^2
+ 8*by^3*cx^2*cy
+ 4*by^3*cy^3
- 24*by^2*cx^2*cy^2
- 2*by^2*cy^4
) + sin(alpha)^4*cos(alpha)^3*(
- 4*ax^4*bx*by
+ 4*ax^4*bx*cy
+ 4*ax^4*by*cx
- 4*ax^4*cx*cy
- 8*ax^3*ay*bx^2
+ 16*ax^3*ay*bx*cx
+ 16*ax^3*ay*by^2
- 32*ax^3*ay*by*cy
- 8*ax^3*ay*cx^2
+ 16*ax^3*ay*cy^2
+ 12*ax^3*bx^2*by
- 4*ax^3*bx^2*cy
- 8*ax^3*bx*by*cx
- 8*ax^3*bx*cx*cy
- 4*ax^3*by*cx^2
+ 12*ax^3*cx^2*cy
+ 12*ax^2*ay*bx^3
- 12*ax^2*ay*bx^2*cx
- 16*ax^2*ay*bx*by^2
+ 48*ax^2*ay*bx*by*cy
- 12*ax^2*ay*bx*cx^2
- 32*ax^2*ay*bx*cy^2
- 32*ax^2*ay*by^2*cx
+ 48*ax^2*ay*by*cx*cy
+ 12*ax^2*ay*cx^3
- 16*ax^2*ay*cx*cy^2
- 8*ax^2*bx^3*by
- 4*ax^2*bx^3*cy
- 12*ax^2*bx^2*by*cx
+ 24*ax^2*bx^2*cx*cy
- 16*ax^2*bx*by^2*cy
+ 24*ax^2*bx*by*cx^2
+ 16*ax^2*bx*by*cy^2
- 12*ax^2*bx*cx^2*cy
+ 16*ax^2*by^2*cx*cy
- 4*ax^2*by*cx^3
- 16*ax^2*by*cx*cy^2
- 8*ax^2*cx^3*cy
+ 8*ax*ay^3*by^2
- 16*ax*ay^3*by*cy
+ 8*ax*ay^3*cy^2
- 16*ax*ay^2*bx^2*by
- 16*ax*ay^2*bx^2*cy
+ 32*ax*ay^2*bx*by*cx
+ 32*ax*ay^2*bx*cx*cy
- 8*ax*ay^2*by^3
+ 8*ax*ay^2*by^2*cy
- 16*ax*ay^2*by*cx^2
+ 8*ax*ay^2*by*cy^2
- 16*ax*ay^2*cx^2*cy
- 8*ax*ay^2*cy^3
- 4*ax*ay*bx^4
- 8*ax*ay*bx^3*cx
+ 48*ax*ay*bx^2*by*cy
+ 24*ax*ay*bx^2*cx^2
+ 16*ax*ay*bx^2*cy^2
+ 32*ax*ay*bx*by^2*cx
- 192*ax*ay*bx*by*cx*cy
- 8*ax*ay*bx*cx^3
+ 32*ax*ay*bx*cx*cy^2
+ 8*ax*ay*by^3*cy
+ 16*ax*ay*by^2*cx^2
- 16*ax*ay*by^2*cy^2
+ 48*ax*ay*by*cx^2*cy
+ 8*ax*ay*by*cy^3
- 4*ax*ay*cx^4
+ 4*ax*bx^4*cy
+ 16*ax*bx^3*by*cx
- 8*ax*bx^3*cx*cy
- 12*ax*bx^2*by*cx^2
- 32*ax*bx^2*by*cy^2
- 12*ax*bx^2*cx^2*cy
+ 32*ax*bx*by^2*cx*cy
- 8*ax*bx*by*cx^3
+ 32*ax*bx*by*cx*cy^2
+ 16*ax*bx*cx^3*cy
- 32*ax*by^2*cx^2*cy
+ 4*ax*by*cx^4
- 8*ay^3*bx*by^2
+ 8*ay^3*bx*by*cy
+ 8*ay^3*by*cx*cy
- 8*ay^3*cx*cy^2
+ 16*ay^2*bx^3*by
- 32*ay^2*bx^2*by*cx
+ 16*ay^2*bx^2*cx*cy
+ 8*ay^2*bx*by^3
+ 8*ay^2*bx*by^2*cy
+ 16*ay^2*bx*by*cx^2
- 16*ay^2*bx*by*cy^2
- 32*ay^2*bx*cx^2*cy
- 16*ay^2*by^2*cx*cy
+ 8*ay^2*by*cx*cy^2
+ 16*ay^2*cx^3*cy
+ 8*ay^2*cx*cy^3
+ 4*ay*bx^4*cx
- 32*ay*bx^3*by*cy
- 4*ay*bx^3*cx^2
+ 48*ay*bx^2*by*cx*cy
- 4*ay*bx^2*cx^3
- 16*ay*bx^2*cx*cy^2
- 16*ay*bx*by^3*cy
- 16*ay*bx*by^2*cx^2
+ 8*ay*bx*by^2*cy^2
+ 48*ay*bx*by*cx^2*cy
+ 8*ay*bx*by*cy^3
+ 4*ay*bx*cx^4
+ 8*ay*by^3*cx*cy
+ 8*ay*by^2*cx*cy^2
- 32*ay*by*cx^3*cy
- 16*ay*by*cx*cy^3
- 4*bx^4*cx*cy
- 8*bx^3*by*cx^2
+ 16*bx^3*by*cy^2
+ 12*bx^3*cx^2*cy
+ 12*bx^2*by*cx^3
- 16*bx^2*by*cx*cy^2
- 8*bx^2*cx^3*cy
+ 8*bx*by^3*cy^2
- 16*bx*by^2*cx^2*cy
- 8*bx*by^2*cy^3
- 4*bx*by*cx^4
- 8*by^3*cx*cy^2
+ 16*by^2*cx^3*cy
+ 8*by^2*cx*cy^3
) + sin(alpha)^5*cos(alpha)^2*(
+ 2*ax^4*bx^2
- 4*ax^4*bx*cx
- 4*ax^4*by^2
+ 8*ax^4*by*cy
+ 2*ax^4*cx^2
- 4*ax^4*cy^2
- 4*ax^3*bx^3
+ 4*ax^3*bx^2*cx
- 16*ax^3*bx*by*cy
+ 4*ax^3*bx*cx^2
+ 16*ax^3*bx*cy^2
+ 16*ax^3*by^2*cx
- 16*ax^3*by*cx*cy
- 4*ax^3*cx^3
- 12*ax^2*ay^2*by^2
+ 24*ax^2*ay^2*by*cy
- 12*ax^2*ay^2*cy^2
+ 8*ax^2*ay*bx^2*by
+ 8*ax^2*ay*bx^2*cy
- 16*ax^2*ay*bx*by*cx
- 16*ax^2*ay*bx*cx*cy
+ 4*ax^2*ay*by^3
- 4*ax^2*ay*by^2*cy
+ 8*ax^2*ay*by*cx^2
- 4*ax^2*ay*by*cy^2
+ 8*ax^2*ay*cx^2*cy
+ 4*ax^2*ay*cy^3
+ 2*ax^2*bx^4
+ 4*ax^2*bx^3*cx
+ 8*ax^2*bx^2*by*cy
- 12*ax^2*bx^2*cx^2
- 24*ax^2*bx^2*cy^2
+ 32*ax^2*bx*by*cx*cy
+ 4*ax^2*bx*cx^3
+ 4*ax^2*by^3*cy
- 24*ax^2*by^2*cx^2
- 8*ax^2*by^2*cy^2
+ 8*ax^2*by*cx^2*cy
+ 4*ax^2*by*cy^3
+ 2*ax^2*cx^4
+ 16*ax*ay^2*bx*by^2
- 24*ax*ay^2*bx*by*cy
+ 8*ax*ay^2*bx*cy^2
+ 8*ax*ay^2*by^2*cx
- 24*ax*ay^2*by*cx*cy
+ 16*ax*ay^2*cx*cy^2
- 16*ax*ay*bx^3*cy
- 16*ax*ay*bx^2*by*cx
+ 32*ax*ay*bx^2*cx*cy
- 24*ax*ay*bx*by^2*cy
+ 32*ax*ay*bx*by*cx^2
+ 32*ax*ay*bx*by*cy^2
- 16*ax*ay*bx*cx^2*cy
- 8*ax*ay*bx*cy^3
- 8*ax*ay*by^3*cx
+ 32*ax*ay*by^2*cx*cy
- 16*ax*ay*by*cx^3
- 24*ax*ay*by*cx*cy^2
- 4*ax*bx^4*cx
+ 4*ax*bx^3*cx^2
+ 16*ax*bx^3*cy^2
- 16*ax*bx^2*by*cx*cy
+ 4*ax*bx^2*cx^3
+ 8*ax*bx*by^2*cy^2
- 16*ax*bx*by*cx^2*cy
- 8*ax*bx*by*cy^3
- 4*ax*bx*cx^4
- 8*ax*by^3*cx*cy
+ 16*ax*by^2*cx^3
+ 8*ax*by^2*cx*cy^2
+ 4*ay^3*bx^2*by
+ 4*ay^3*bx^2*cy
- 8*ay^3*bx*by*cx
- 8*ay^3*bx*cx*cy
+ 4*ay^3*by*cx^2
+ 4*ay^3*cx^2*cy
- 4*ay^2*bx^4
+ 16*ay^2*bx^3*cx
- 12*ay^2*bx^2*by^2
- 4*ay^2*bx^2*by*cy
- 24*ay^2*bx^2*cx^2
- 8*ay^2*bx^2*cy^2
+ 8*ay^2*bx*by^2*cx
+ 32*ay^2*bx*by*cx*cy
+ 16*ay^2*bx*cx^3
+ 8*ay^2*bx*cx*cy^2
- 8*ay^2*by^2*cx^2
- 4*ay^2*by*cx^2*cy
- 4*ay^2*cx^4
- 12*ay^2*cx^2*cy^2
+ 8*ay*bx^4*cy
- 16*ay*bx^3*cx*cy
+ 24*ay*bx^2*by^2*cy
+ 8*ay*bx^2*by*cx^2
- 4*ay*bx^2*by*cy^2
+ 8*ay*bx^2*cx^2*cy
+ 4*ay*bx^2*cy^3
- 24*ay*bx*by^2*cx*cy
- 16*ay*bx*by*cx^3
- 24*ay*bx*by*cx*cy^2
+ 4*ay*by^3*cx^2
- 4*ay*by^2*cx^2*cy
+ 8*ay*by*cx^4
+ 24*ay*by*cx^2*cy^2
+ 2*bx^4*cx^2
- 4*bx^4*cy^2
- 4*bx^3*cx^3
- 12*bx^2*by^2*cy^2
+ 8*bx^2*by*cx^2*cy
+ 4*bx^2*by*cy^3
+ 2*bx^2*cx^4
+ 16*bx*by^2*cx*cy^2
+ 4*by^3*cx^2*cy
- 4*by^2*cx^4
- 12*by^2*cx^2*cy^2
) + sin(alpha)^6*cos(alpha)^1*(
+ 8*ax^3*ay*by^2
- 16*ax^3*ay*by*cy
+ 8*ax^3*ay*cy^2
- 8*ax^2*ay*bx*by^2
+ 24*ax^2*ay*bx*by*cy
- 16*ax^2*ay*bx*cy^2
- 16*ax^2*ay*by^2*cx
+ 24*ax^2*ay*by*cx*cy
- 8*ax^2*ay*cx*cy^2
- 8*ax^2*bx*by^2*cy
+ 8*ax^2*bx*by*cy^2
+ 8*ax^2*by^2*cx*cy
- 8*ax^2*by*cx*cy^2
- 8*ax*ay^2*bx^2*by
- 8*ax*ay^2*bx^2*cy
+ 16*ax*ay^2*bx*by*cx
+ 16*ax*ay^2*bx*cx*cy
- 8*ax*ay^2*by*cx^2
- 8*ax*ay^2*cx^2*cy
+ 24*ax*ay*bx^2*by*cy
+ 8*ax*ay*bx^2*cy^2
+ 16*ax*ay*bx*by^2*cx
- 96*ax*ay*bx*by*cx*cy
+ 16*ax*ay*bx*cx*cy^2
+ 8*ax*ay*by^2*cx^2
+ 24*ax*ay*by*cx^2*cy
- 16*ax*bx^2*by*cy^2
+ 16*ax*bx*by^2*cx*cy
+ 16*ax*bx*by*cx*cy^2
- 16*ax*by^2*cx^2*cy
+ 8*ay^2*bx^3*by
- 16*ay^2*bx^2*by*cx
+ 8*ay^2*bx^2*cx*cy
+ 8*ay^2*bx*by*cx^2
- 16*ay^2*bx*cx^2*cy
+ 8*ay^2*cx^3*cy
- 16*ay*bx^3*by*cy
+ 24*ay*bx^2*by*cx*cy
- 8*ay*bx^2*cx*cy^2
- 8*ay*bx*by^2*cx^2
+ 24*ay*bx*by*cx^2*cy
- 16*ay*by*cx^3*cy
+ 8*bx^3*by*cy^2
- 8*bx^2*by*cx*cy^2
- 8*bx*by^2*cx^2*cy
+ 8*by^2*cx^3*cy
) + sin(alpha)^7*cos(alpha)^0*(
- 2*ax^4*by^2
+ 4*ax^4*by*cy
- 2*ax^4*cy^2
- 8*ax^3*bx*by*cy
+ 8*ax^3*bx*cy^2
+ 8*ax^3*by^2*cx
- 8*ax^3*by*cx*cy
+ 4*ax^2*ay*bx^2*by
+ 4*ax^2*ay*bx^2*cy
- 8*ax^2*ay*bx*by*cx
- 8*ax^2*ay*bx*cx*cy
+ 4*ax^2*ay*by*cx^2
+ 4*ax^2*ay*cx^2*cy
+ 4*ax^2*bx^2*by*cy
- 12*ax^2*bx^2*cy^2
+ 16*ax^2*bx*by*cx*cy
- 12*ax^2*by^2*cx^2
+ 4*ax^2*by*cx^2*cy
- 8*ax*ay*bx^3*cy
- 8*ax*ay*bx^2*by*cx
+ 16*ax*ay*bx^2*cx*cy
+ 16*ax*ay*bx*by*cx^2
- 8*ax*ay*bx*cx^2*cy
- 8*ax*ay*by*cx^3
+ 8*ax*bx^3*cy^2
- 8*ax*bx^2*by*cx*cy
- 8*ax*bx*by*cx^2*cy
+ 8*ax*by^2*cx^3
- 2*ay^2*bx^4
+ 8*ay^2*bx^3*cx
- 12*ay^2*bx^2*cx^2
+ 8*ay^2*bx*cx^3
- 2*ay^2*cx^4
+ 4*ay*bx^4*cy
- 8*ay*bx^3*cx*cy
+ 4*ay*bx^2*by*cx^2
+ 4*ay*bx^2*cx^2*cy
- 8*ay*bx*by*cx^3
+ 4*ay*by*cx^4
- 2*bx^4*cy^2
+ 4*bx^2*by*cx^2*cy
- 2*by^2*cx^4
);
z = sin(alpha)^0*cos(alpha)^7*(
+ 8*ax*ay^2*by^2
- 16*ax*ay^2*by*cy
+ 8*ax*ay^2*cy^2
- 8*ax*ay*by^3
+ 8*ax*ay*by^2*cy
+ 8*ax*ay*by*cy^2
- 8*ax*ay*cy^3
+ 8*ax*by^3*cy
- 16*ax*by^2*cy^2
+ 8*ax*by*cy^3
- 8*ay^3*bx*by
+ 8*ay^3*bx*cy
+ 8*ay^3*by*cx
- 8*ay^3*cx*cy
+ 8*ay^2*bx*by^2
+ 8*ay^2*bx*by*cy
- 16*ay^2*bx*cy^2
- 16*ay^2*by^2*cx
+ 8*ay^2*by*cx*cy
+ 8*ay^2*cx*cy^2
- 16*ay*bx*by^2*cy
+ 8*ay*bx*by*cy^2
+ 8*ay*bx*cy^3
+ 8*ay*by^3*cx
+ 8*ay*by^2*cx*cy
- 16*ay*by*cx*cy^2
+ 8*bx*by^2*cy^2
- 8*bx*by*cy^3
- 8*by^3*cx*cy
+ 8*by^2*cx*cy^2
) + sin(alpha)^1*cos(alpha)^6*(
- 16*ax^2*ay*by^2
+ 32*ax^2*ay*by*cy
- 16*ax^2*ay*cy^2
+ 8*ax^2*by^3
- 8*ax^2*by^2*cy
- 8*ax^2*by*cy^2
+ 8*ax^2*cy^3
+ 8*ax*ay^2*bx*by
- 8*ax*ay^2*bx*cy
- 8*ax*ay^2*by*cx
+ 8*ax*ay^2*cx*cy
+ 8*ax*ay*bx*by^2
- 32*ax*ay*bx*by*cy
+ 24*ax*ay*bx*cy^2
+ 24*ax*ay*by^2*cx
- 32*ax*ay*by*cx*cy
+ 8*ax*ay*cx*cy^2
- 8*ax*bx*by^2*cy
+ 24*ax*bx*by*cy^2
- 16*ax*bx*cy^3
- 16*ax*by^3*cx
+ 24*ax*by^2*cx*cy
- 8*ax*by*cx*cy^2
+ 8*ay^3*bx^2
- 16*ay^3*bx*cx
+ 8*ay^3*cx^2
- 16*ay^2*bx^2*by
- 8*ay^2*bx^2*cy
+ 24*ay^2*bx*by*cx
+ 24*ay^2*bx*cx*cy
- 8*ay^2*by*cx^2
- 16*ay^2*cx^2*cy
+ 32*ay*bx^2*by*cy
- 8*ay*bx^2*cy^2
- 8*ay*bx*by^2*cx
- 32*ay*bx*by*cx*cy
- 8*ay*bx*cx*cy^2
- 8*ay*by^2*cx^2
+ 32*ay*by*cx^2*cy
- 16*bx^2*by*cy^2
+ 8*bx^2*cy^3
+ 8*bx*by^2*cx*cy
+ 8*bx*by*cx*cy^2
+ 8*by^3*cx^2
- 16*by^2*cx^2*cy
) + sin(alpha)^2*cos(alpha)^5*(
+ 8*ax^3*by^2
- 16*ax^3*by*cy
+ 8*ax^3*cy^2
+ 8*ax^2*ay*bx*by
- 8*ax^2*ay*bx*cy
- 8*ax^2*ay*by*cx
+ 8*ax^2*ay*cx*cy
- 16*ax^2*bx*by^2
+ 24*ax^2*bx*by*cy
- 8*ax^2*bx*cy^2
- 8*ax^2*by^2*cx
+ 24*ax^2*by*cx*cy
- 16*ax^2*cx*cy^2
- 16*ax*ay^2*bx^2
+ 32*ax*ay^2*bx*cx
+ 16*ax*ay^2*by^2
- 32*ax*ay^2*by*cy
- 16*ax*ay^2*cx^2
+ 16*ax*ay^2*cy^2
+ 8*ax*ay*bx^2*by
+ 24*ax*ay*bx^2*cy
- 32*ax*ay*bx*by*cx
- 32*ax*ay*bx*cx*cy
- 16*ax*ay*by^3
+ 16*ax*ay*by^2*cy
+ 24*ax*ay*by*cx^2
+ 16*ax*ay*by*cy^2
+ 8*ax*ay*cx^2*cy
- 16*ax*ay*cy^3
- 8*ax*bx^2*by*cy
- 8*ax*bx^2*cy^2
+ 32*ax*bx*by^2*cx
- 32*ax*bx*by*cx*cy
+ 32*ax*bx*cx*cy^2
+ 16*ax*by^3*cy
- 8*ax*by^2*cx^2
- 32*ax*by^2*cy^2
- 8*ax*by*cx^2*cy
+ 16*ax*by*cy^3
- 16*ay^3*bx*by
+ 16*ay^3*bx*cy
+ 16*ay^3*by*cx
- 16*ay^3*cx*cy
+ 8*ay^2*bx^3
- 8*ay^2*bx^2*cx
+ 16*ay^2*bx*by^2
+ 16*ay^2*bx*by*cy
- 8*ay^2*bx*cx^2
- 32*ay^2*bx*cy^2
- 32*ay^2*by^2*cx
+ 16*ay^2*by*cx*cy
+ 8*ay^2*cx^3
+ 16*ay^2*cx*cy^2
- 16*ay*bx^3*cy
- 8*ay*bx^2*by*cx
+ 24*ay*bx^2*cx*cy
- 32*ay*bx*by^2*cy
+ 24*ay*bx*by*cx^2
+ 16*ay*bx*by*cy^2
- 8*ay*bx*cx^2*cy
+ 16*ay*bx*cy^3
+ 16*ay*by^3*cx
+ 16*ay*by^2*cx*cy
- 16*ay*by*cx^3
- 32*ay*by*cx*cy^2
+ 8*bx^3*cy^2
+ 8*bx^2*by*cx*cy
- 16*bx^2*cx*cy^2
- 16*bx*by^2*cx^2
+ 16*bx*by^2*cy^2
+ 8*bx*by*cx^2*cy
- 16*bx*by*cy^3
- 16*by^3*cx*cy
+ 8*by^2*cx^3
+ 16*by^2*cx*cy^2
) + sin(alpha)^3*cos(alpha)^4*(
- 8*ax^3*bx*by
+ 8*ax^3*bx*cy
+ 8*ax^3*by*cx
- 8*ax^3*cx*cy
+ 8*ax^2*ay*bx^2
- 16*ax^2*ay*bx*cx
- 32*ax^2*ay*by^2
+ 64*ax^2*ay*by*cy
+ 8*ax^2*ay*cx^2
- 32*ax^2*ay*cy^2
+ 8*ax^2*bx^2*by
- 16*ax^2*bx^2*cy
+ 8*ax^2*bx*by*cx
+ 8*ax^2*bx*cx*cy
+ 16*ax^2*by^3
- 16*ax^2*by^2*cy
- 16*ax^2*by*cx^2
- 16*ax^2*by*cy^2
+ 8*ax^2*cx^2*cy
+ 16*ax^2*cy^3
+ 16*ax*ay^2*bx*by
- 16*ax*ay^2*bx*cy
- 16*ax*ay^2*by*cx
+ 16*ax*ay^2*cx*cy
- 8*ax*ay*bx^3
+ 8*ax*ay*bx^2*cx
+ 16*ax*ay*bx*by^2
- 64*ax*ay*bx*by*cy
+ 8*ax*ay*bx*cx^2
+ 48*ax*ay*bx*cy^2
+ 48*ax*ay*by^2*cx
- 64*ax*ay*by*cx*cy
- 8*ax*ay*cx^3
+ 16*ax*ay*cx*cy^2
+ 8*ax*bx^3*cy
- 16*ax*bx^2*by*cx
+ 8*ax*bx^2*cx*cy
- 16*ax*bx*by^2*cy
+ 8*ax*bx*by*cx^2
+ 48*ax*bx*by*cy^2
- 16*ax*bx*cx^2*cy
- 32*ax*bx*cy^3
- 32*ax*by^3*cx
+ 48*ax*by^2*cx*cy
+ 8*ax*by*cx^3
- 16*ax*by*cx*cy^2
+ 16*ay^3*bx^2
- 32*ay^3*bx*cx
+ 16*ay^3*cx^2
- 32*ay^2*bx^2*by
- 16*ay^2*bx^2*cy
+ 48*ay^2*bx*by*cx
+ 48*ay^2*bx*cx*cy
- 16*ay^2*by*cx^2
- 32*ay^2*cx^2*cy
+ 8*ay*bx^3*cx
+ 64*ay*bx^2*by*cy
- 16*ay*bx^2*cx^2
- 16*ay*bx^2*cy^2
- 16*ay*bx*by^2*cx
- 64*ay*bx*by*cx*cy
+ 8*ay*bx*cx^3
- 16*ay*bx*cx*cy^2
- 16*ay*by^2*cx^2
+ 64*ay*by*cx^2*cy
- 8*bx^3*cx*cy
+ 8*bx^2*by*cx^2
- 32*bx^2*by*cy^2
+ 8*bx^2*cx^2*cy
+ 16*bx^2*cy^3
+ 16*bx*by^2*cx*cy
- 8*bx*by*cx^3
+ 16*bx*by*cx*cy^2
+ 16*by^3*cx^2
- 32*by^2*cx^2*cy
) + sin(alpha)^4*cos(alpha)^3*(
+ 16*ax^3*by^2
- 32*ax^3*by*cy
+ 16*ax^3*cy^2
+ 16*ax^2*ay*bx*by
- 16*ax^2*ay*bx*cy
- 16*ax^2*ay*by*cx
+ 16*ax^2*ay*cx*cy
- 32*ax^2*bx*by^2
+ 48*ax^2*bx*by*cy
- 16*ax^2*bx*cy^2
- 16*ax^2*by^2*cx
+ 48*ax^2*by*cx*cy
- 32*ax^2*cx*cy^2
- 32*ax*ay^2*bx^2
+ 64*ax*ay^2*bx*cx
+ 8*ax*ay^2*by^2
- 16*ax*ay^2*by*cy
- 32*ax*ay^2*cx^2
+ 8*ax*ay^2*cy^2
+ 16*ax*ay*bx^2*by
+ 48*ax*ay*bx^2*cy
- 64*ax*ay*bx*by*cx
- 64*ax*ay*bx*cx*cy
- 8*ax*ay*by^3
+ 8*ax*ay*by^2*cy
+ 48*ax*ay*by*cx^2
+ 8*ax*ay*by*cy^2
+ 16*ax*ay*cx^2*cy
- 8*ax*ay*cy^3
- 16*ax*bx^2*by*cy
- 16*ax*bx^2*cy^2
+ 64*ax*bx*by^2*cx
- 64*ax*bx*by*cx*cy
+ 64*ax*bx*cx*cy^2
+ 8*ax*by^3*cy
- 16*ax*by^2*cx^2
- 16*ax*by^2*cy^2
- 16*ax*by*cx^2*cy
+ 8*ax*by*cy^3
- 8*ay^3*bx*by
+ 8*ay^3*bx*cy
+ 8*ay^3*by*cx
- 8*ay^3*cx*cy
+ 16*ay^2*bx^3
- 16*ay^2*bx^2*cx
+ 8*ay^2*bx*by^2
+ 8*ay^2*bx*by*cy
- 16*ay^2*bx*cx^2
- 16*ay^2*bx*cy^2
- 16*ay^2*by^2*cx
+ 8*ay^2*by*cx*cy
+ 16*ay^2*cx^3
+ 8*ay^2*cx*cy^2
- 32*ay*bx^3*cy
- 16*ay*bx^2*by*cx
+ 48*ay*bx^2*cx*cy
- 16*ay*bx*by^2*cy
+ 48*ay*bx*by*cx^2
+ 8*ay*bx*by*cy^2
- 16*ay*bx*cx^2*cy
+ 8*ay*bx*cy^3
+ 8*ay*by^3*cx
+ 8*ay*by^2*cx*cy
- 32*ay*by*cx^3
- 16*ay*by*cx*cy^2
+ 16*bx^3*cy^2
+ 16*bx^2*by*cx*cy
- 32*bx^2*cx*cy^2
- 32*bx*by^2*cx^2
+ 8*bx*by^2*cy^2
+ 16*bx*by*cx^2*cy
- 8*bx*by*cy^3
- 8*by^3*cx*cy
+ 16*by^2*cx^3
+ 8*by^2*cx*cy^2
) + sin(alpha)^5*cos(alpha)^2*(
- 16*ax^3*bx*by
+ 16*ax^3*bx*cy
+ 16*ax^3*by*cx
- 16*ax^3*cx*cy
+ 16*ax^2*ay*bx^2
- 32*ax^2*ay*bx*cx
- 16*ax^2*ay*by^2
+ 32*ax^2*ay*by*cy
+ 16*ax^2*ay*cx^2
- 16*ax^2*ay*cy^2
+ 16*ax^2*bx^2*by
- 32*ax^2*bx^2*cy
+ 16*ax^2*bx*by*cx
+ 16*ax^2*bx*cx*cy
+ 8*ax^2*by^3
- 8*ax^2*by^2*cy
- 32*ax^2*by*cx^2
- 8*ax^2*by*cy^2
+ 16*ax^2*cx^2*cy
+ 8*ax^2*cy^3
+ 8*ax*ay^2*bx*by
- 8*ax*ay^2*bx*cy
- 8*ax*ay^2*by*cx
+ 8*ax*ay^2*cx*cy
- 16*ax*ay*bx^3
+ 16*ax*ay*bx^2*cx
+ 8*ax*ay*bx*by^2
- 32*ax*ay*bx*by*cy
+ 16*ax*ay*bx*cx^2
+ 24*ax*ay*bx*cy^2
+ 24*ax*ay*by^2*cx
- 32*ax*ay*by*cx*cy
- 16*ax*ay*cx^3
+ 8*ax*ay*cx*cy^2
+ 16*ax*bx^3*cy
- 32*ax*bx^2*by*cx
+ 16*ax*bx^2*cx*cy
- 8*ax*bx*by^2*cy
+ 16*ax*bx*by*cx^2
+ 24*ax*bx*by*cy^2
- 32*ax*bx*cx^2*cy
- 16*ax*bx*cy^3
- 16*ax*by^3*cx
+ 24*ax*by^2*cx*cy
+ 16*ax*by*cx^3
- 8*ax*by*cx*cy^2
+ 8*ay^3*bx^2
- 16*ay^3*bx*cx
+ 8*ay^3*cx^2
- 16*ay^2*bx^2*by
- 8*ay^2*bx^2*cy
+ 24*ay^2*bx*by*cx
+ 24*ay^2*bx*cx*cy
- 8*ay^2*by*cx^2
- 16*ay^2*cx^2*cy
+ 16*ay*bx^3*cx
+ 32*ay*bx^2*by*cy
- 32*ay*bx^2*cx^2
- 8*ay*bx^2*cy^2
- 8*ay*bx*by^2*cx
- 32*ay*bx*by*cx*cy
+ 16*ay*bx*cx^3
- 8*ay*bx*cx*cy^2
- 8*ay*by^2*cx^2
+ 32*ay*by*cx^2*cy
- 16*bx^3*cx*cy
+ 16*bx^2*by*cx^2
- 16*bx^2*by*cy^2
+ 16*bx^2*cx^2*cy
+ 8*bx^2*cy^3
+ 8*bx*by^2*cx*cy
- 16*bx*by*cx^3
+ 8*bx*by*cx*cy^2
+ 8*by^3*cx^2
- 16*by^2*cx^2*cy
) + sin(alpha)^6*cos(alpha)^1*(
+ 8*ax^3*by^2
- 16*ax^3*by*cy
+ 8*ax^3*cy^2
+ 8*ax^2*ay*bx*by
- 8*ax^2*ay*bx*cy
- 8*ax^2*ay*by*cx
+ 8*ax^2*ay*cx*cy
- 16*ax^2*bx*by^2
+ 24*ax^2*bx*by*cy
- 8*ax^2*bx*cy^2
- 8*ax^2*by^2*cx
+ 24*ax^2*by*cx*cy
- 16*ax^2*cx*cy^2
- 16*ax*ay^2*bx^2
+ 32*ax*ay^2*bx*cx
- 16*ax*ay^2*cx^2
+ 8*ax*ay*bx^2*by
+ 24*ax*ay*bx^2*cy
- 32*ax*ay*bx*by*cx
- 32*ax*ay*bx*cx*cy
+ 24*ax*ay*by*cx^2
+ 8*ax*ay*cx^2*cy
- 8*ax*bx^2*by*cy
- 8*ax*bx^2*cy^2
+ 32*ax*bx*by^2*cx
- 32*ax*bx*by*cx*cy
+ 32*ax*bx*cx*cy^2
- 8*ax*by^2*cx^2
- 8*ax*by*cx^2*cy
+ 8*ay^2*bx^3
- 8*ay^2*bx^2*cx
- 8*ay^2*bx*cx^2
+ 8*ay^2*cx^3
- 16*ay*bx^3*cy
- 8*ay*bx^2*by*cx
+ 24*ay*bx^2*cx*cy
+ 24*ay*bx*by*cx^2
- 8*ay*bx*cx^2*cy
- 16*ay*by*cx^3
+ 8*bx^3*cy^2
+ 8*bx^2*by*cx*cy
- 16*bx^2*cx*cy^2
- 16*bx*by^2*cx^2
+ 8*bx*by*cx^2*cy
+ 8*by^2*cx^3
) + sin(alpha)^7*cos(alpha)^0*(
- 8*ax^3*bx*by
+ 8*ax^3*bx*cy
+ 8*ax^3*by*cx
- 8*ax^3*cx*cy
+ 8*ax^2*ay*bx^2
- 16*ax^2*ay*bx*cx
+ 8*ax^2*ay*cx^2
+ 8*ax^2*bx^2*by
- 16*ax^2*bx^2*cy
+ 8*ax^2*bx*by*cx
+ 8*ax^2*bx*cx*cy
- 16*ax^2*by*cx^2
+ 8*ax^2*cx^2*cy
- 8*ax*ay*bx^3
+ 8*ax*ay*bx^2*cx
+ 8*ax*ay*bx*cx^2
- 8*ax*ay*cx^3
+ 8*ax*bx^3*cy
- 16*ax*bx^2*by*cx
+ 8*ax*bx^2*cx*cy
+ 8*ax*bx*by*cx^2
- 16*ax*bx*cx^2*cy
+ 8*ax*by*cx^3
+ 8*ay*bx^3*cx
- 16*ay*bx^2*cx^2
+ 8*ay*bx*cx^3
- 8*bx^3*cx*cy
+ 8*bx^2*by*cx^2
+ 8*bx^2*cx^2*cy
- 8*bx*by*cx^3
);
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment